関数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がスペクトル値としてでる事を期待していましたが途方もなく異なる値が出てしまいました。
どうぞよろしくお願いします。
No.1
- 回答日時:
離散フーリエ変換に関数f(x)=2sin(πx)をかければ、きれいに1つだけ
ピークがでると思います。それにCOSの成分はほとんど0になっているはずです。
指数表示からSIN、COS表示に直して、それぞれの成分をみてください。
スペクトルの値が合わないということですので、ソフトが違っている可能性が
高いです。(データ入力ミスか、ソフト自体のバグか)
直接結果をみることができないので あくまで想像ですが。
ありがとうございます。おっしゃるようにピークがでているのは虚数部分の値で、実数部分(Re)はほとんどゼロです。プログラム自体は下記のサイトから転用したので、間違っている可能性は少ないかもしれません。
http://www.geocities.jp/supermisosan/spectrum.html
サイトでも大きな値がスペクトル値として出ています。
振幅とスペクトルの関係をもう少し考えてみます。
No.2
- 回答日時:
え~, ちょっといろいろ書いてみたんですが,
逆変換を考えてみてください
というアドバイスが最も簡単かつ適切かなぁ.
ちなみに原信号の振幅を A とすると
1000 / N = A / 2
が成り立っています.
回答ありがとうございます。
逆変換するとやはり元の実空間の数値に戻るようです。
考えてみましたら、実空間において関数の振幅の値は、場所によって異なりますから、スペクトル値は振幅の2乗となるからA(最大振幅のみ)^2がスペクトルとなるというのはおかしいかもしれませんね。すべての点における振幅の2乗がスペクトルとなるのかもしれません。
教えていただいた式の1000の部分はスペクトルの大きさ、Nがサンプリング数という理解でよろしいでしょうか??
No.3ベストアンサー
- 回答日時:
#2 です.
はい, 1000 はスペクトルの大きさで, N は点数です.
DFT の基本をきちんとおさえておけばいいんだけど, 自分でも思い出しながら書いてみる:
ω = exp (2πi/N)
とおきます. ω は 1 の原始 N 乗根です.
今考えている関数だと f(k) = 2 sin(2πk/N) = 2(ω^k - ω^(-k)) / (2i) = -i (ω^k - ω^(-k))
です... あ, 離散フーリエ変換の式中の f(k) は 2 sin (2πk/N) ですよね? f(x) = 2 sin πx をそのまま使っていないですよね?
すると, exp (-2πink/N) = ω^(-nk) ですから
Σ (k=0~N-1) f(k) exp (-2πink/N)
= Σ (k=0~N-1) -i (ω^k - ω^(-k)) ω^(-nk)
= Σ (k=0~N-1) -i ω^(1-n)k + Σ (k=0~N-1) i ω^(-(1+n)k))
となりますが, 第1項は 1 - n ≡ 0 mod N のときにのみ 0 でない値となり, 一方第2項は -(1+n) ≡ 0 mod N のときにのみ 0 でない値となります. 順に考えると
1.1 - n ≡ 0 mod N のとき. このとき (0 ≦ n < N より) n = 1 です. そして和の値は
1 + ω^N + ω^(2N) + ... + ω^((N-1)N) = N
となります. 1+n = 2 ですから第2項は 0. つまり n = 1 に対して -iN という値になります.
2.-(1+n) ≡ 0 mod N のとき. このとき n = N-1 で, 和の値は1と同じく N になります. 第1項は 0 になりますから, フーリエ変換の値として iN が得られます.
3.その他のとき. 第1項, 第2項ともに 0 なので 0.
ということで, プログラムで得られた結果と一致します. 振幅の 2 は sin を複素表示したときに分母にある 2 と相殺されていることに注意してください.
ありがとうございました。
お忙しい中、丁寧な解説ありがとうございます。
完璧とまではいきませんが、良く理解することができました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 フーリエ変換後の負の周波数成分の扱いについて 4 2022/09/03 10:18
- 数学 フーリエ変換、逆変換の「2π」の扱いについて 3 2022/10/07 08:31
- 物理学 フーリエ変換の振幅について 1 2022/09/04 08:56
- 数学 離散フーリエ逆変換が周波数分割数をNにできる理由について 4 2022/09/18 12:56
- その他(自然科学) 原子の線スペクトルについて 1 2022/05/07 09:10
- 物理学 打撃力は広範囲の連続周波数スペクトルを有するとあります。 どういうことなのでしょうか。 連続周波数ス 6 2022/04/21 15:50
- 物理学 なぜ、問題文からこの式が出てくるのでしょうか? 5 2023/08/15 15:07
- 工学 周波数fで表現したフーリエ変換の対称性に関する質問です。 1 2022/09/14 12:27
- 物理学 有限の大きさの物質では、周期的境界条件を満たすように格子振動が発生する。もし、満たさない場合、物質の 1 2022/07/05 18:37
- 物理学 風力発電での音 1 2023/04/16 08:55
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
2÷3などの余りについて
-
信頼区間の1.96や1.65ってどこ...
-
マイナスからプラスへ転じた時...
-
ある商品のロス率を5%見込み、...
-
EXCELの分散分析表のP-値が....
-
n^2-20n+91が素数とな...
-
大きな負の値?負の大きな値???
-
中学生で数学です。 √84nが自然...
-
制御工学解ける方!!
-
1/fゆらぎの傾きについて
-
ニューラルネットワークの重み...
-
青チャートI基本例題89
-
エクセルで可視セルにのみ値貼...
-
エックスの値は5cmですか?
-
教えてください。数学Bの二項分...
-
√19-aの値が整数となる正の整数...
-
パーセント点はP1 P5 P10 P25 P...
-
エクセル 3つの値の中からデー...
-
パーセントの出し方を教えて下さい
-
行列です
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
マイナスからプラスへ転じた時...
-
2÷3などの余りについて
-
信頼区間の1.96や1.65ってどこ...
-
Aの値からBの値を除するとは??
-
エクセルで可視セルにのみ値貼...
-
変数とパラメータとは違うもの...
-
ある商品のロス率を5%見込み、...
-
値差の%計算方法について
-
エックスの値は5cmですか?
-
a^2の√=a が成り立たない場合
-
Excelで1つしかない値だけを抽...
-
20'(角度)の計算がわかりま...
-
「an=(n-1)/(n+1)のときlim[n→∞...
-
EXCELの分散分析表のP-値が....
-
「Aに対するBの割合」と「Aに対...
-
10%引いた元の数字を出すには?
-
パーセントの出し方を教えて下さい
-
楕円の外周の計算方法
-
中学生で数学です。 √84nが自然...
-
教えてください。数学Bの二項分...
おすすめ情報