お世話になります。
フーリエ変換の勉強を始めています。
f(t)=cosωt という関数を次式で手計算でフーリエ変換しました。
Cn = 1/T・∫f(t)e^(-jω_nt)dt 積分区間は-T/2~T/2
この場合、結果は
Cn = 1 (ω=ω_n) 振幅(縦軸)=1
0 (ω≠ω_n)
となるのが正しいと思うのですが、計算結果をエクセルでグラフで書いたところ、
Cn = 1 (ω≒ω_n) 振幅(縦軸)=0.5
0 (ω≠ω_n)
特異解(ω=ω_n)
となりました。
手計算は添付に示します。どこがおかしいのかご教示頂けませんでしょうか?
よろしくお願い致します
No.2ベストアンサー
- 回答日時:
ANo.1へのコメントについてです。
> Cn = 1 (ω=ω_n) 振幅(縦軸)=1
> 0 (ω≠ω_n)
>
> となるのが正しいと思う
という思い込みが間違っているんです。そうではなく、
C[n] = C[-n] = 0.5
になるのが正しい。総和Σを取る添字の範囲に注意してANo.1をご再読あれ。
No.1
- 回答日時:
計算なさったのはフーリエ変換
F(ω) = ∫{t=-∞~∞} f(t) exp(-jωt) dt
ではなく、複素フーリエ級数展開ですね。
実数から複素数への関数f(t)が周期Tを持つ滑らかな関数のとき、fの複素フーリエ級数は
f(t) = Σ{n=-∞~∞} C[n] exp(2πj n t/T)
C[n] = (1/T)∫{t=-T/2~T/2} f(t) exp(-2πj n t/T) dt
ここで
○ fが実数から実数への関数である場合、C[n]とC[-n]は共役複素数である。
○ 全てのtについてf(t)がf(-t)の共役複素数である場合、C[n]は実数である。
は容易に証明できます。なので、
○ fが実数から実数への偶関数である場合、C[n]は実数で、C[n] = C[-n] である。
ご質問の場合はこれに該当しますね。従って、
f(t) = Σ{n=-∞~∞} C[n] exp(2πj n t/T)
= C[0] + Σ{n=1~∞} (C[n] exp(2πj n t/T)+ C[-n]exp(-2πj n t/T))
= C[0] + Σ{n=1~∞} C[n] (exp(2πj n t/T)+exp(-2πj n t/T))
= C[0] + 2 Σ{n=1~∞} C[n] cos(2πn t/T)
だから計算なさった結果には何の不思議も無い。
要するに、総和をn=-∞~∞の範囲で取ることをお忘れであり、C[-n]を足していないんです。
あるいは、もしかすると、実数から実数への周期Tを持つ関数f(t)に関する実フーリエ級数
f(t) = c[0] + Σ{n=1~∞} (c[n] cos(2π n t/T) + s[n] sin(2π n t/T) )
と混同なさったのかも知れない。実フーリエ級数と複素フーリエ級数との関係はご自分で検討できるでしょう。
ちなみに、数値計算の場合には、整数から複素数(実数)への周期関数f(n)の離散フーリエ変換 (DFT)が行われます。複素離散フーリエ変換のアルゴリズムのひとつが高速フーリエ変換(FFT)です。なお、フーリエ変換とフーリエ級数展開との関係をキッチリ考えるには超関数の理論を学ぶ必要があります。
ともあれ、「相互に関連のある似たような用語がいろいろあり、それらを区別しなくちゃいけない」というのが回答です。
この回答への補足
遅くなり申し訳ありません。ご回答いただきありがとうございました。
おっしゃるとおり、私がしたのはフーリエ変換ではなく、複素フーリエ級数展開でした。
頂いた回答を読んだのですが、理解しかねるところがあるので、もう少し教えて頂けませんでしょうか?
f(t)=cosωtとしたときの複素フーリエ係数Cnを求めたいのですが、式にそのまま当てはめると
Cn = 1/T・∫f(t)e^(-jω_nt)dt 積分区間は-T/2~T/2
となり、それを解くと
Cn=1/T[1/(ω-ωn)×sin{(ω-ωn)T/2}}+1/(ω+ωn)×sin{(ω+ωn)T/2}]
となりました。ここまでが質問に書いた部分です。
この計算が違うということでしょうか?正しくはどういう式になるのでしょうか?そしてその式の解(複素フーリエ係数)はどうなるのでしょうか?恐れ入りますが、ご教示頂けませんでしょうか?
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 離散フーリエ逆変換が周波数分割数をNにできる理由について 4 2022/09/18 12:56
- 数学 フーリエ変換後の負の周波数成分の扱いについて 4 2022/09/03 10:18
- 物理学 フーリエ変換の振幅について 1 2022/09/04 08:56
- 数学 フーリエ変換、逆変換の「2π」の扱いについて 3 2022/10/07 08:31
- Excel(エクセル) エクセルのSUM関数について 4 2023/04/18 10:37
- 数学 f(x)=e^(-ax+b) のフーリエ変換をフーリエ変換の定義に従って計算せよ。但し、a>0、bは 1 2023/02/06 18:26
- その他(Microsoft Office) Excelで時間計算(負) 8 2023/02/26 05:47
- 数学 f(x)=e^|-ax+b| のフーリエ変換を計算せよ。 但し、a>0、bは定数 この問題についての 2 2023/02/10 16:13
- 数学 数学の質問です。 関数f(t)のフーリエ変換をF(ω)=∫[-∞→∞]f(t)exp(-iωt)dt 1 2023/07/29 01:08
- 数学 フーリエ変換についての質問です。 h(t)=cos(ω0t)×cos(ω1t) のフーリエ変換を教え 1 2022/07/23 17:37
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報