1. 程式人生 > >第十八講 傅立葉變換

第十八講 傅立葉變換

一,傅立葉級數的複數形式:

  • 一般形式:
  1. f(t)=\frac{a_{0}}{2}+\sum_{n=1}^{\infty }[a_{n}cos(nk_{0}t)+b_{n}sin(nk_{0}t)],週期T=2Lk_{0}=\frac{2\pi }{T}
  2. a_{n}=\frac{1}{L}\int_{0}^{2L }f(t)cos(nk_{0}t)dt,n為任意整數
  3. b_{n}=\frac{1}{L}\int_{0}^{2L }f(t)sin(nk_{0}t)dt,n為任意整數
  • 利用逆向尤拉公式(第十講第二節)展開:
  • a_{n}cos(nk_{0}t)+b_{n}sin(nk_{0}t)=\frac{a_{n}-ib_{n}}{2}e^{ink_{0}t}+\frac{a_{n}+ib_{n}}{2}e^{-ink_{0}t}
  • 原式:f(t)=\frac{a_{0}}{2}+\sum_{n=1}^{\infty }[\frac{a_{n}-ib_{n}}{2}e^{ink_{0}t}+\frac{a_{n}+ib_{n}}{2}e^{-ink_{0}t}]=\frac{a_{0}}{2}+\sum_{n=1}^{\infty }\frac{a_{n}-ib_{n}}{2}e^{ink_{0}t}+\sum_{n=1}^{\infty }\frac{a_{n}+ib_{n}}{2}e^{-ink_{0}t}
  • 第一項:\frac{a_{0}}{2}=\sum_{n=0}^{0}\frac{a_{0}}{2}e^{ink_{0}t}
  • 第三項:當n< 0時,\sum_{n=1}^{\infty }\frac{a_{n}+ib_{n}}{2}e^{-ink_{0}t}=\sum_{n=-\infty }^{-1}\frac{a_{-n}+ib_{-n}}{2}e^{ink_{0}t}
  • 代入原式:f(t)=\sum_{n=0}^{0}\frac{a_{0}}{2}e^{ink_{0}t}+\sum_{n=1}^{\infty }\frac{a_{n}-ib_{n}}{2}e^{ink_{0}t}+\sum_{n=-\infty }^{-1}\frac{a_{-n}+ib_{-n}}{2}e^{ink_{0}t}=\sum_{n=-\infty }^{\infty }c_{n}e^{ink_{0}t}
  • n=0時:c_{n}=\frac{a_{0}}{2}=\frac{1}{2}\cdot \frac{1}{L}\int_{0}^{2L }f(t)cos(0)dt=\frac{1}{2L}\int_{0}^{2L }f(t)dt=\frac{1}{2L}\int_{0}^{2L }f(t)e^{-ink_{0}t}dt
  • n> 0時:c_{n}=\frac{a_{n}-ib_{n}}{2}=\frac{1}{2L}\int_{0}^{2L }f(t)cos(nk_{0}t)dt-i\cdot\frac{1}{2L}\int_{0}^{2L }f(t)sin(nk_{0}t)dt=\frac{1}{2L}\int_{0}^{2L }f(t)[cos(nk_{0}t)-isin(nk_{0}t)]dt=\frac{1}{2L}\int_{0}^{2L }f(t)e^{-ink_{0}t}dt
  • n< 0時:c_{n}=\frac{a_{-n}+ib_{-n}}{2}=\frac{1}{2L}\int_{0}^{2L }f(t)e^{-ink_{0}t}dt
  • 總結:f(t)=\sum_{n=-\infty }^{\infty }c_{n}e^{ink_{0}t}c_{n}=\frac{1}{T}\int_{0}^{T}f(t)e^{-ink_{0}t}dt

二,傅立葉變換:

  • f_{T}(t)表示週期為T的f(t)f_{T}(t)=f(t+T)
  • 對於任意f_{T}(t)=\sum_{n=-\infty }^{\infty }c_{n}e^{ink_{0}t},其中\sum_{n=-\infty }^{\infty }e^{ink_{0}t}都是不變的
  • 唯一改變的是c_{n}=\frac{1}{T}\int_{0}^{T}f_{T}(t)e^{-ink_{0}t}dt=\frac{1}{T}\int_{-L}^{L}f_{T}(t)e^{-ink_{0}t}dt
  • 頻譜如圖(圖中的\omega_{0}k_{0}意思相同):
  • 兩個不同角速度\omega之間的間距\Delta \omega =(n+1)k_{0}-nk_{0}=k_{0},如圖(圖中的\omega_{0}k_{0}意思相同):
  • 當週期T\rightarrow \infty時,周期函式f_{T}(t)拓展到非周期函式f(t)\lim_{T\rightarrow \infty }f_{T}(t)=f(t)
  • 此時間距\Delta \omega =k_{0}=\frac{2\pi }{T}\rightarrow 0,原本離散的角速度\omega變連續了,原本的n從只能取整數變成可以取整個實數
  • \omega =nk_{0}這個角速度公式才能在整個實數軸上成立,如圖(圖中的\omega_{0}k_{0}意思相同):
  • 將傅立葉係數c_{n}代入周期函式f_{T}(t)
  • f_{T}(t)=\sum_{n=-\infty }^{\infty }\frac{1}{T}\int_{-L}^{L}f_{T}(t)e^{-ink_{0}t}dt\cdot e^{ink_{0}t}
  • \frac{1}{T}=\frac{k_{0}}{2\pi }=\frac{\Delta \omega }{2\pi }也代入:
  • f_{T}(t)=\sum_{n=-\infty }^{\infty }\frac{\Delta \omega }{2\pi }\int_{-L}^{L}f_{T}(t)e^{-ink_{0}t}dt\cdot e^{ink_{0}t}
  • T\rightarrow \infty時:\sum_{n=-\infty }^{\infty }\Delta \omega =\int_{-\infty }^{\infty }d\omegaL\rightarrow \infty\omega =nk_{0}
  •  
  • f(t)=\frac{1}{2\pi }\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }f(t)e^{-i\omega t}dt\cdot e^{i\omega t}d\omega,這個式子被稱作逆變換
  • 而其中\int_{-\infty }^{\infty }f(t)e^{-i\omega t}dt=F(\omega ),被稱作傅立葉變換

三,傅立葉變換的應用(濾波):

  • 首先將f(t)進行傅立葉變換,檢視有哪些角速度(角頻率)的波需要去掉,再用濾波器濾掉(濾掉高頻波最簡單的手段就是對傅立葉級數進行積分)如圖(黑線是原函式,藍線是一重積分,綠線是二重積分):