- ベストアンサー
離散フーリエ変換について
離散フーリエ変換によって得られた値についての質問です。 多くのサイトでその値は Σ(k=0～N-1) f(k)exp(-2πkni/N) という式から求められるとあります。 離散フーリエ変換は本来、ある周期関数が、どのくらいの振幅でどのくらいの周波数の波からできているかを調べるために行うものだと思います。 しかし上記の公式から得られるスペクトル（sqrt（Re^2+Im^2)）では振幅の値は得られません。振幅を得るには刻み幅Δ（関数をサンプリングした際の幅）を乗じて Σ(k=0～N-1) f(k)exp(-2πkni/N)*Δ とすれば得られることが分かりました。 最初の公式から得られるスペクトルはなにを表しているのでしょうか？またなぜ刻み幅Δを乗じることで、振幅が求まるのでしょうか？ よろしくお願いします。m(__)m
- taro_kao
- お礼率80% (8/10)
- 数学・算数
- 回答数4
- ありがとう数2
- みんなの回答 （4）
- 専門家の回答
質問者が選んだベストアンサー
パッと見だけど, Δを掛けた式で N → ∞ の極限をとると連続の Fourier 変換にならないかなぁ?
その他の回答 (3)
- guuman
- ベストアンサー率30% (100/331)
逆変換においてNで割るか変換においてNで割るかは流儀の違い なので好きな方を選べばよい DFTが何を洗わしているかが重要 DFTはサンプリングされた周期関数のフーリエ変換の結果であるインパルス列の係数を表す 周期関数をフーリエ変換すると周波数軸上でインパルス列となるが そのインパルス列の大きさを表すのである さらにサンプリングされているからそのインパルス列は周期を持つ
お礼
ありがとうございました。 参考になり助かりました。
- Meowth
- ベストアンサー率35% (130/362)
もうすこし 正規化はどうなっているのですか。 使っているDFTの 離散フーリエ変換と逆離散フーリエ変換 の式をおしえてください。 ｛おそらく フーリエ変換ではNで割らない 逆フーリエ変換ではNで割る｝だとおもいますが。 ΔをかけるとかNで割るとか、 Δをかける前と後のスペクトル（sqrt（Re^2+Im^2)）の 違いがみえてくると思いますが。
補足
書き込みありがとうございます。 使っている離散フーリエ変換の式は 離散フーリエ変換 Σ(k=0～N-1) f(k)exp(-2πkni/N) 離散逆フーリエ変換 （1／Ｎ）Σ(n=0～N-1) f(n)exp(2πkni/N) です。（１／Ｎ）は逆変換の方についています。 （１／Ｎ）を逆変換ではなく、変換側につけますと、逆変換の式は Σ(n=0～N-1) f(n)exp(2πkni/N) になると思います。 スペクトル値として振幅をえるためには、 Σ(k=0～N-1) f(k)exp(-2πkni/N)×Δ としなければなりませんでした。Δなしですと、振幅×100の値がスペクトル値（sqrt（Re^2+Im^2)）としてえられました。 この100倍のスペクトル値はなにを示しているのでしょうか？ また刻み幅Δを乗じることで、なぜ振幅が求まるのでしょうか？ 刻み幅Δは Δ＝サンプリング領域Ｌ／サンプリング数Ｎ です。
- Meowth
- ベストアンサー率35% (130/362)
Σ(k=0～N-1) f(k)exp(-2πkni/N) はなんでしょうか？ Nで割らないのでしょうか 刻み幅ΔとNの関係は？
補足
書き込みありがとうございます。 Σ(k=0～N-1) f(k)exp(-2πkni/N)は離散フーリエ変換の式です。 http://ja.wikipedia.org/wiki/%E9%9B%A2%E6%95%A3%E3%83%95%E3%83%BC%E3%83%AA%E3%82%A8%E5%A4%89%E6%8F%9B （wikiでの離散フーリエ変換の式） Ｎで割るタイプ Σ(k=0～N-1) f(k)exp(-2πkni/N)／Ｎ http://akita-nct.jp/~yamamoto/comp/mathematica/material/fft/fourier_html/node1.html もあるようです。Δ＝刻み幅＝（サンプリングする領域Ｌ）／Ｎ ですから、このＮで割る公式の場合はＬを乗じて Σ(k=0～N-1) f(k)exp(-2πkni/N)／Ｎ×Ｌとし、スペクトルを計算しますと、振幅の値がスペクトルとして表示されました。 Ｎで割らないタイプ（Σ(k=0～N-1) f(k)exp(-2πkni/N)）ですと、振幅をスペクトルとして表示させるには、Δを乗じなければなりませんでした。
関連するQ&A
- 離散フーリエ変換のスペクトルについて
関数f(x)=2sin(πx)をx=0～2まで等間隔1000点でサンプリングし、 離散フーリエ変換 Σ(k=0～N-1) f(k)exp(-2πkni/N) の式から、言語プログラムで計算する式をつくり、1000個の実数Reと虚数Imを得ました。 ピークはもちろん周波数πのときで、スペクトルの値が1000でした。 √（Re^2+Im^2）をスペクトル値、√なしをパワースペクトル値をいうそうですが、元の関数の振幅2とこのスペクトル値とはどのような関係があるのでしょうか？ 異なる正弦波を混ぜれば、スペクトル値を見ることによって振幅の比は分かりますが、スペクトル値と振幅には式的になんらかの関係は存在するのでしょうか？ 波のエネルギーは振幅の2乗になると思い、2^2=4がスペクトル値としてでる事を期待していましたが途方もなく異なる値が出てしまいました。 どうぞよろしくお願いします。
- 締切済み
- 数学・算数
- ある時間関数を離散フーリエ変換して得られるフーリエスペクトルの振幅値に
ある時間関数を離散フーリエ変換して得られるフーリエスペクトルの振幅値について教えて下さい。 今想定している離散フーリエ変換の式は一般的なもので Σ(k=0～N-1) f(k)exp(-2πkni/N) を考えています。 また、離散フーリエ変換して得られるスペクトルは √（Re^2+Im^2) で計算します。 離散フーリエ変換を適用する関数を、 振幅1の直流、及び振幅1で周波数5[Hz]の正弦波とします。 （この2つの信号は別々の信号で合成されていません。） サンプリング周波数を20[Hz]とした場合、サンプリングして得られるデータ列はそれぞれ、 直流： 「1, 1, 1, 1」 正弦波： 「0, 1, 0, -1」 となると想定されます。 （正弦波をサンプリングする場合は位相が関わってきますが、今回は気にしないで下さい。） このデータ列に対して上記の離散フーリエ変換を適用した場合、 得られるフーリエスペクトルの振幅値はそれぞれ、 直流： 「4」（直流のフーリエスペクトルの振幅値値） 正弦波： 「2」（5[Hz]のフーリエスペクトルの振幅値） となります。 （データ点数は上の通り4点） ここで質問なのですが、 離散フーリエ変換して得られるスペクトルの振幅値から元の関数の振幅値を求める場合、 フーリエスペクトルをサンプリングの総データ点数で割ることは数学的に納得できます。 しかしこの例の場合、フーリエスペクトルを総データ点数で割ると、 直流： 「4 -> 1」 正弦波： 「2 -> 0.5」 となってしまい、直流は正しいのですが、正弦波の元の振幅値を正確に求めることは出来ません。 この例の場合、フーリエスペクトルの振幅値から正弦波の振幅値を正しく求めるには、 「フーリエスペクトルの振幅値*2/データ点数」 としてやらなければいけません。 上記のことに関して、なぜこのようになるのかを（2をかける理由を）教えて頂けないでしょうか。 当方、数学についてはあまり詳しくないため、簡単に説明して頂けると幸いです。
- 締切済み
- 数学・算数
- 離散フーリエ変換について
G(n/Nτ)= Σ[k=0->N-1] {τ*f(kτ)*e^(-i2πkn/N)} 上記の式は離散フーリエ変換の式らしいですが、 これを関数化しようとしてつまずいています。 どのように解釈すれば、関数化できるのでしょうか？ 特に、複素数iがよくわかりません。 e^iが点を左に90度回転させるくらいはわかります。 e^(-i2πkn/N)を関数powで表現できなくて困っていますが、多分見当違いだとは思います。 フーリエ級数展開はわかります。 最終的にはFFTを行いたいのですが、 その理解の前に離散フーリエ変換ができないといけないと思っています。 よろしくお願いします。 void GraphClass::ScatteredFourierConvert() { /* N-1 G(n/Nτ)= Σ {τ*f(kτ)*e^(-i2πkn/N)} k=0 k=何番目の値か？ τ=値を読んだ間隔 N=値を読んだ数 */ int Tau=10; int N=RangeSize.cx*2/Tau; double Sum=0; double e=2.7182818; for(int k=0;k<N-1;k++) { Sum+=Tau*Data[k*Tau]*pow(e,-i*2*PI*k*n/N); } }
- ベストアンサー
- C・C++・C#
- フーリエスペクトルの振幅について
ある時間関数を離散フーリエ変換して得られるフーリエスペクトルの 振幅値について教えて下さい。 今想定している離散フーリエ変換の式は一般的なもので Σ(k=0～N-1) f(k)exp(-2πkni/N) を考えています。 また、離散フーリエ変換して得られるスペクトルは √（Re^2+Im^2) で計算します。 離散フーリエ変換を適用する関数を、 振幅1の直流、及び振幅1で周波数5[Hz]の正弦波とします。 （この2つの信号は別々の信号で合成されていません。） サンプリング周波数を20[Hz]とした場合、 サンプリングして得られるデータ列はそれぞれ、 直流： 「1, 1, 1, 1」 正弦波： 「0, 1, 0, -1」 となると想定されます。 （正弦波をサンプリングする場合は位相が関わってきますが、 今回は気にしないで下さい。） このデータ列に対して上記の離散フーリエ変換を適用した場合、 得られるフーリエスペクトルの振幅値はそれぞれ、 直流： 「4」（直流のフーリエスペクトルの振幅値値） 正弦波： 「2」（5[Hz]のフーリエスペクトルの振幅値） となります。 （データ点数は上の通り4点） ここで質問なのですが、 離散フーリエ変換して得られるスペクトルの振幅値から 元の関数の振幅値を求める場合、 フーリエスペクトルをサンプリングの総データ点数で割ることは 数学的に納得できます。 しかしこの例の場合、フーリエスペクトルを総データ点数で割ると、 直流： 「4 -> 1」 正弦波： 「2 -> 0.5」 となってしまい、直流は正しいのですが、 正弦波の元の振幅値を正確に求めることは出来ません。 フーリエスペクトルの振幅値から正弦波の振幅値を正しく求めるには、 「フーリエスペクトルの振幅値*2/データ点数」 としてやらなければいけません。 上記のことに関して、 なぜこのようになるのかを（2をかける理由を）教えて頂けないでしょうか。
- ベストアンサー
- 物理学
- 離散フーリエ変換（ＤＦＴ）の公式について
離散フーリエ変換の公式は、参考書等によりますと色々な記述があります。 今回は２例の違いを教えていただきたいのです。 １）F(u)=1/NΣf(x)・・・・ ２）F(u)=Σf(x)・・・・ との式があります（両式とも詳細は省いて記述しました）。この規格化定数（１/Ｎ）がある公式１）と、ない公式２）があります。 本来の周波数スペクトル（振幅）を表しているのは１）式であると考えていますが いかがでしょうか？ 公式２）を使用して算出した場合には、その値をサンプリング数で除すれば公式１）同じ結果となるのですが、なぜ公式２）が記述してあるのでしょうか？
- ベストアンサー
- 数学・算数
- 離散フーリエ変換の対称定理について
離散フーリエ変換において、 Re[X(k)]=Re[X(N-k)] Im[X(k)]=-Im[X(N-k)] |X(k)|=|X(N-k)| といった対称定理が成り立ちます。 数学的にこれらが成り立つということは、大体理解できたのですが、直感的な意味が分かりません。 僕の理解が正しいとしたら、N=128だとしたら、 N=1の極めて周期の大きな波と、N=127の極めて周期の短い波の、 振幅スペクトル（≒振幅？）が等しいということになると思うのですが、どうしても納得が行きません。 N>64の周波数をフィルターでカットした場合も、対称定理は成り立ち、N>64の周波数の波は出てくる（というか、カットしないと、エイリアシングが生じる）と思うのですが、どういうことなのでしょうか？ N>64の周波数の波は、計算上は出てくるが、実際にはそのような波は存在しないのですか？ あるいは、N>64の周波数同士の波で、上手く相殺されて、逆離散フーリエ変換をしたときには、影響が出ないということなのですか？ また、通常のフーリエ変換では、対称定理は成り立たないと思うのですが、そのことも併せて教えて頂けたらと思います。 よろしくお願いします。
- 締切済み
- 数学・算数
- 離散フーリエ変換について
double DiscreteFourierTransformClass::DiscreteFourierTransformProcess(int n) { double e=2.7182818; double Pai=3.1415926535; int Tau=10; int N=DataLen/Tau; double Sum=0; for(int k=0;k<N;k++) { Sum+=Tau*Data[k*Tau]*pow(e,i*2*Pai*k*n/N); } return Sum; } Tau*Data[k*Tau]=帯 pow(e,i*2*Pai*k*n/N)=ほしい周波数の乗算 ということで、離散フーリエ展開と同じなのでしょうか？ パソコン上で計算するには、離散フーリエ展開で計算して、 FFTというのは、pow(e,i*2*Pai*k*n/N)この回転子が3,4個程度のcos,sinで決まるから その値で計算すれば高速になるよということでしょうか？ 離散フーリエ変換 G(n/N)=Σ[k=0::N-1]{τ*f(kτ)e^-i2πkn/N} k:何番目の値か τ:値を読んだ間隔 N:値を読んだ数 n:基本周波数の何倍か G(n/N)->ほしいnの成分を得る
- ベストアンサー
- C・C++・C#
- 離散フーリエ展開の意味（イメージ）
離散フーリエ展開のよく意味がわかりません。 フーリエ変換は任意の周期関数をsunやcosの関数の和 で書くことですよね？それの複素表示も理解できるのですが…。 DFT（離散フーリエ展開）は関数の値を細分化して、 そしてそれを変えるいうことですか？ それはどのようなことなのでしょうか？ Uj=cos(2π/N)jとVj=sin(2π/N)jの離散フーリエ展開を すると、離散フーリエ変換したUkがk=1、N-1で振動成分をもつとはどういうことですか？ まだ、大学１年で非常に難しくて困っています。 どなたかわかりやすく教えてほしいです。 お願いします。
- ベストアンサー
- 数学・算数
お礼
その一言が私の求めていた答えでした！！ 本当に感謝します。