傅里叶级数展开的详细推导和部分证明
正文
傅里葉級數表述為:
f(t)=a0+∑k=1∞{akcos?(2πkT0t)+bksin?(2πkT0t)}f(t) = a_0 + \sum^\infin_{k=1} \left \{ a_k \cos \left( \frac {2\pi k}{T_0}t \right) + b_k \sin \left( \frac{2\pi k}{T_0}t \right)\right \} f(t)=a0?+k=1∑∞?{ak?cos(T0?2πk?t)+bk?sin(T0?2πk?t)}
上式中,T0T_0T0? 是傅里葉展開前原周期函數的一個周期。可以看到三角函數的角頻率 ω\omegaω 為 ω=2πkT0\omega = \frac {2 \pi k}{T_0}ω=T0?2πk?
※三角函數的中 ttt 前面的系數是角頻率 ω\omegaω 它與周期的關系是 ω=2πT\omega = \frac {2\pi}{T}ω=T2π? ,所以可以發現,在原函數的一個周期 T0T_0T0? 中, 當 kkk 變化,角頻率 ω\omegaω 隨之變化,于是周期 TTT , 變成了 T02\frac{T_0}{2}2T0?? 、T03\frac{T_0}{3}3T0?? 、T04\frac{T_0}{4}4T0?? … T0∞\frac{T_0}{\infin}∞T0?? 。
于是,上式就成了無窮個同相位、各種振幅的三角函數之和了。
( ※注意:振幅為 aka_kak? 或 bkb_kbk? )
周期 TTT 隨 kkk 增大不斷減小,即,頻率 fff 和角頻率 ω\omegaω 不斷增大, 最大的周期稱為基本周期 T0T_0T0? ,最小的(角)頻率稱為基本(角)頻率 f0f_0f0? ( ω0\omega_0ω0? ),傅里葉級數的每一份三角函數的頻率都是基本頻率的整數倍。
接下來,級數展開需要求的三角函數的系數,即,傅里葉系數
求常數項 a0a_0a0? (直流成分)
∫?T0/2T0/2f(t)dt=∫?T0/2T0/2{a0+∑k=1∞{akcos?(2πkT0t)+bksin?(2πkT0t)}dt}\int_{-T_0/2}^{T_0/2} f(t)dt = \int_{-T_0/2}^{T_0/2} \left \{ a_0 + \sum^\infin_{k=1} \left \{ a_k \cos \left( \frac {2\pi k}{T_0}t \right) + b_k \sin \left( \frac{2\pi k}{T_0}t \right)\right \}dt \right \} ∫?T0?/2T0?/2?f(t)dt=∫?T0?/2T0?/2?{a0?+k=1∑∞?{ak?cos(T0?2πk?t)+bk?sin(T0?2πk?t)}dt}
由于積分和累加本質上就是求和(積分是連續和,累加是離散和)所以把積分號放進括號里,有:
∫?T0/2T0/2f(t)dt=∫?T0/2T0/2a0dt+∑k=1∞{∫?T0/2T0/2akcos?(2πkT0t)dt+∫?T0/2T0/2bksin?(2πkT0t)dt}\int_{-T_0/2}^{T_0/2} f(t)dt = \int_{-T_0/2}^{T_0/2} a_0dt + \sum^\infin_{k=1} \left \{ \int_{-T_0/2}^{T_0/2} a_k \cos \left( \frac {2\pi k}{T_0}t \right)dt + \int_{-T_0/2}^{T_0/2} b_k \sin \left( \frac{2\pi k}{T_0}t \right)dt \right \} ∫?T0?/2T0?/2?f(t)dt=∫?T0?/2T0?/2?a0?dt+k=1∑∞?{∫?T0?/2T0?/2?ak?cos(T0?2πk?t)dt+∫?T0?/2T0?/2?bk?sin(T0?2πk?t)dt}
在一個周期內,三角函數的定積分等于0,所以右側后項被消掉,得:
∫?T0/2T0/2f(t)dt=∫?T0/2T0/2a0dt=T0a0\int_{-T_0/2}^{T_0/2} f(t)dt = \int_{-T_0/2}^{T_0/2} a_0dt = T_0a_0 ∫?T0?/2T0?/2?f(t)dt=∫?T0?/2T0?/2?a0?dt=T0?a0?
移項,
a0=1T0∫?T0/2T0/2f(t)dta_0 = \frac{1}{T_0} \int_{-T_0/2}^{T_0/2} f(t)dt a0?=T0?1?∫?T0?/2T0?/2?f(t)dt
即直流成分的傅里葉系數 a0a_0a0? 為原函數的1個周期的定積分除以周期,a0a_0a0? 物理意義就是,原函數的時間平均值。
(關于三角函數正交性的一部分推導)
試證:
∫?T0/2T0/2cos?(2πmT0t)cos?(2πnT0t)dt=T02δm,n\int_{-T_0/2}^{T_0/2} \cos( \frac{2 \pi m}{T_0}t) \cos( \frac{2 \pi n}{T_0}t)dt = \frac{T_0}{2} \delta_{m,n} ∫?T0?/2T0?/2?cos(T0?2πm?t)cos(T0?2πn?t)dt=2T0??δm,n?
其中,克羅內克德爾塔為
δm,n={0,if?m≠n1,if?m=n\delta_{m,n}= \begin{cases} 0, & \text{if $m \not = n$} \\ 1, & \text{if $m = n$} \end{cases} δm,n?={0,1,?if?m?=nif?m=n?
證明:
根據積化和差公式: (help)
哥哥 = 哥 + 哥
cos?αcos?β=cos?(α+β)+cos?(α?β)2\cos \alpha \cos \beta = \frac{ \cos( \alpha + \beta) + \cos( \alpha - \beta)}{2} cosαcosβ=2cos(α+β)+cos(α?β)?
有:
∫?T0/2T0/2cos?(2πmT0t)cos?(2πnT0t)dt=∫?T0/2T0/2cos?(2πmT0t+2πnT0t)+cos?(2πmT0t?2πnT0t)2dt=∫?T0/2T0/2cos?(2πmT0t+2πnT0t)2dt+∫?T0/2T0/2cos?(2πmT0t?2πnT0t)2dt\int_{-T_0/2}^{T_0/2} \cos ( \frac{2 \pi m}{T_0}t) \cos ( \frac{2 \pi n}{T_0}t)dt = \int_{-T_0/2}^{T_0/2} \frac{ \cos( \frac{2 \pi m}{T_0}t + \frac{2 \pi n}{T_0}t) + \cos( \frac{2 \pi m}{T_0}t - \frac{2 \pi n}{T_0}t)}{2}dt\\ = \int_{-T_0/2}^{T_0/2} \frac{ \cos( \frac{2 \pi m}{T_0}t + \frac{2 \pi n}{T_0}t)}{2}dt + \int_{-T_0/2}^{T_0/2}\frac{ \cos( \frac{2 \pi m}{T_0}t - \frac{2 \pi n}{T_0}t)}{2}dt ∫?T0?/2T0?/2?cos(T0?2πm?t)cos(T0?2πn?t)dt=∫?T0?/2T0?/2?2cos(T0?2πm?t+T0?2πn?t)+cos(T0?2πm?t?T0?2πn?t)?dt=∫?T0?/2T0?/2?2cos(T0?2πm?t+T0?2πn?t)?dt+∫?T0?/2T0?/2?2cos(T0?2πm?t?T0?2πn?t)?dt
右側第一項分子始終不為0,則一個周期內的積分為0;同理,右側第二項只有當分子cos?(2πmT0t?2πnT0t)\cos( \frac{2 \pi m}{T_0}t - \frac{2 \pi n}{T_0}t)cos(T0?2πm?t?T0?2πn?t)的括號中為0時,該項變為常數項,積分才不為零,即只有 m=nm=nm=n 時,cos?0=1\cos0 = 1cos0=1,等式右側變為:
∫?T0/2T0/212dt\int_{-T_0/2}^{T_0/2} \frac{1}{2}dt ∫?T0?/2T0?/2?21?dt
等于:
T02\frac{T_0}{2} 2T0??
得證。
現在以求系數 a3a_3a3? 為例,在 ∫?T0/2T0/2f(t)dt=∫?T0/2T0/2a0dt+∑k=1∞{∫?T0/2T0/2akcos?(2πkT0t)dt+∫?T0/2T0/2bksin?(2πkT0t)dt}\int_{-T_0/2}^{T_0/2} f(t)dt = \int_{-T_0/2}^{T_0/2} a_0dt + \sum^\infin_{k=1} \left \{ \int_{-T_0/2}^{T_0/2} a_k \cos \left( \frac {2\pi k}{T_0}t \right)dt + \int_{-T_0/2}^{T_0/2} b_k \sin \left( \frac{2\pi k}{T_0}t \right)dt \right \}∫?T0?/2T0?/2?f(t)dt=∫?T0?/2T0?/2?a0?dt+∑k=1∞?{∫?T0?/2T0?/2?ak?cos(T0?2πk?t)dt+∫?T0?/2T0?/2?bk?sin(T0?2πk?t)dt} 的兩邊乘上一個 cos?(2π?3T0)\cos ( \frac{2 \pi \cdot 3}{T_0})cos(T0?2π?3?) ,則:
∫?T0/2T0/2f(t)cos?(2π?3T0)dt=∫?T0/2T0/2a0cos?(2π?3T0)dt+∑k=1∞{∫?T0/2T0/2akcos?(2πkT0t)cos?(2π?3T0)dt+∫?T0/2T0/2bksin?(2πkT0t)cos?(2π?3T0)dt}=0+∫?T0/2T0/2a3?T02dt=a3T02\int_{-T_0/2}^{T_0/2} f(t) \cos( \frac{2 \pi \cdot 3}{T_0}) dt = \int_{-T_0/2}^{T_0/2} a_0 \cos( \frac{2 \pi \cdot 3}{T_0})dt + \sum^\infin_{k=1} \left \{ \int_{-T_0/2}^{T_0/2} a_k \cos \left( \frac {2\pi k}{T_0}t \right) \cos( \frac{2 \pi \cdot 3}{T_0})dt + \int_{-T_0/2}^{T_0/2} b_k \sin \left( \frac{2\pi k}{T_0}t \right) \cos( \frac{2 \pi \cdot 3}{T_0}) dt \right \}\\ = 0 + \int_{-T_0/2}^{T_0/2} a_3 \cdot \frac {T_0}{2}dt = \frac{a_3T_0}{2} ∫?T0?/2T0?/2?f(t)cos(T0?2π?3?)dt=∫?T0?/2T0?/2?a0?cos(T0?2π?3?)dt+k=1∑∞?{∫?T0?/2T0?/2?ak?cos(T0?2πk?t)cos(T0?2π?3?)dt+∫?T0?/2T0?/2?bk?sin(T0?2πk?t)cos(T0?2π?3?)dt}=0+∫?T0?/2T0?/2?a3??2T0??dt=2a3?T0??
于是,移項得:
a3=2T0∫?T0/2T0/2f(t)cos?(2π?3T0)dta_3 = \frac{2}{T_0} \int_{-T_0/2}^{T_0/2} f(t) \cos( \frac{2 \pi \cdot 3}{T_0}) dt a3?=T0?2?∫?T0?/2T0?/2?f(t)cos(T0?2π?3?)dt
以此類推便完成了傅里葉級數展開, 將時域連續的原函數展開成了一大堆頻域離散三角函數的和。
即:
a0=1T0∫?T0/2T0/2f(t)dtak=2T0∫?T0/2T0/2f(t)cos?(2π?kT0)dt,(k=1,2,3...)bk=2T0∫?T0/2T0/2f(t)sin?(2π?kT0)dt,(k=1,2,3...)a_0 = \frac{1}{T_0} \int_{-T_0/2}^{T_0/2} f(t)dt\\ a_k = \frac{2}{T_0} \int_{-T_0/2}^{T_0/2} f(t) \cos( \frac{2 \pi \cdot k}{T_0}) dt, \quad (k=1,2,3...)\\ b_k = \frac{2}{T_0} \int_{-T_0/2}^{T_0/2} f(t) \sin( \frac{2 \pi \cdot k}{T_0}) dt, \quad (k=1,2,3...) a0?=T0?1?∫?T0?/2T0?/2?f(t)dtak?=T0?2?∫?T0?/2T0?/2?f(t)cos(T0?2π?k?)dt,(k=1,2,3...)bk?=T0?2?∫?T0?/2T0?/2?f(t)sin(T0?2π?k?)dt,(k=1,2,3...)
總結
(1)周期為 T0T_0T0? 的周期函數 f(t)f(t)f(t) 能表示為f(t)=a0+∑k=1∞{akcos?(2πkT0t)+bksin?(2πkT0t)}f(t) = a_0 + \sum^\infin_{k=1} \left \{ a_k \cos \left( \frac {2\pi k}{T_0}t \right) + b_k \sin \left( \frac{2\pi k}{T_0}t \right)\right \}f(t)=a0?+∑k=1∞?{ak?cos(T0?2πk?t)+bk?sin(T0?2πk?t)}的形式,這個被稱為傅里葉展開
(2)ak,bka_k, \, b_kak?,bk? 被稱為傅里葉系數
(3)累加起來的這些三角函數周期只存在原函數頻率 f0f_0f0? 的整數倍,即 f=2f0,3f0...f = 2f_0, \, 3f_0 \, ...f=2f0?,3f0?... 不存在其他頻率的正余弦函數。
注意:
1)級數展開右側第一項a0a_0a0?在很多教材中寫成a02\frac{a_0}{2}2a0?? 是為了在傅里葉系數表達式中,把 a0a_0a0? 公式包含進 aka_kak? 的式子中。實際上無論怎么寫這個展開式都一樣,僅僅是為了系數表達式的美觀
2)在推導系數表達式時積分和累加直接互換其實是不嚴謹的,這里的級數展開式的等號實際需要證明其成立,此處并未詳細說明,在DFT中會進一步說明,未完待續…
reference
[1] 如何巧記「積化和差」與「和差化積」公式? - 青話的回答 - 知乎
[2] フーリエ級數
總結
以上是生活随笔為你收集整理的傅里叶级数展开的详细推导和部分证明的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 哪几种情况可以销毁session(一般有
- 下一篇: 编程手册集锦