UA OPTI512R 傅立叶光学导论13 傅立叶变换及其性质
UA OPTI512R 傅立葉光學(xué)導(dǎo)論13 傅立葉變換及其性質(zhì)
- 傅立葉變換的定義
- 驗證Fourier變換與其逆變換的定義互為逆變換
- 實值函數(shù)的Fourier變換
- 例:矩形波的Fourier變換
- 傅立葉變換的性質(zhì)
- 常用Fourier變換對
傅立葉變換的定義
上一講介紹了周期函數(shù)的Fourier級數(shù)展開:
fp(x)=∑n=?∞+∞cnej2πnξ0x,fp(x+nT)=fp(x),?n∈Zcn=1T∫t0t0+Tfp(x)e?j2πnξ0xf_p(x)= \sum_{n = -\infty}^{+\infty} c_n e^{j2 \pi n \xi_0 x},f_p(x+nT)=f_p(x),\forall n \in \mathbb{Z} \\ c_n = \frac{1}{T} \int_{t_0}^{t_0+T} f_p(x) e^{-j2 \pi n \xi_0 x}fp?(x)=n=?∞∑+∞?cn?ej2πnξ0?x,fp?(x+nT)=fp?(x),?n∈Zcn?=T1?∫t0?t0?+T?fp?(x)e?j2πnξ0?x
盡管Fourier級數(shù)展開的方法對于周期函數(shù)是非常簡單通用的,但是在光學(xué)中處理波形的時候大部分情況還是非周期的波,這時就需要引入Fourier變換來處理了。對于任意復(fù)值函數(shù)f(x)f(x)f(x),定義它的Fourier變換為
F(ξ)=∫?∞+∞f(x)e?j2πξxdxF(\xi) = \int_{-\infty}^{+\infty} f(x)e^{-j2 \pi \xi x}dxF(ξ)=∫?∞+∞?f(x)e?j2πξxdx
可以發(fā)現(xiàn)它的本質(zhì)與Fourier級數(shù)展開的系數(shù)類似,都是f(x)f(x)f(x)與基函數(shù)e?j2πξxe^{-j 2\pi \xi x}e?j2πξx的內(nèi)積,所以Fourier變換其實就是f(x)f(x)f(x)在基函數(shù)空間{ej2πξx:ξ∈R}\{e^{j 2 \pi \xi x}: \xi \in \mathbb{R}\}{ej2πξx:ξ∈R}上的投影,只是這里的頻率ξ\xiξ可以是任意值,而Fourier級數(shù)展開是把f(x)f(x)f(x)投影到頻率為nξ0n\xi_0nξ0?的基函數(shù)列{ej2πnξ0x:n∈Z}\{e^{j2 \pi n \xi_0 x}: n \in \mathbb{Z}\}{ej2πnξ0?x:n∈Z}上。
已知F(ξ)F(\xi)F(ξ),我們也可以還原出f(x)f(x)f(x),這個操作叫Fourier變換的逆變換(inverse Fourier transform),
f(x)=∫?∞+∞F(ξ)ej2πξxdξf(x) = \int_{-\infty}^{+\infty} F(\xi)e^{j 2 \pi \xi x}d\xif(x)=∫?∞+∞?F(ξ)ej2πξxdξ
這個式子可以類比周期函數(shù)的Fourier級數(shù)展開,{ej2πξx}\{e^{j 2 \pi \xi x}\}{ej2πξx}是基函數(shù),F(ξ)F(\xi)F(ξ)相當(dāng)于基函數(shù)對應(yīng)的系數(shù),稱∣F(ξ)∣|F(\xi)|∣F(ξ)∣為f(x)f(x)f(x)的頻率譜(frequency spectrum)。
驗證Fourier變換與其逆變換的定義互為逆變換
記F\mathcal{F}F為Fourier變換算子,F?1\mathcal{F}^{-1}F?1為Fourier變換的逆變換算子,則
F(ξ)=F[f(x)]f(x)=F?1[F(ξ)]F(\xi) = \mathcal{F}[f(x)] \\ f(x)=\mathcal{F}^{-1}[F(\xi)]F(ξ)=F[f(x)]f(x)=F?1[F(ξ)]
驗證一下這兩個算子在數(shù)學(xué)上確實互為逆算子,即驗證
f(x)=F?1[F(ξ)]=F?1[F[f(x)]]f(x) = \mathcal{F}^{-1}[F(\xi)] = \mathcal{F}^{-1}[\mathcal{F}[f(x)]]f(x)=F?1[F(ξ)]=F?1[F[f(x)]]
好餓啊不想打公式了,直接貼老師的notes好了。。。
第二步到第三步是Fubini定理,第三步到第四步的積分可以使用111的Fourier變換是2πδ(ξ)2 \pi \delta(\xi)2πδ(ξ)的結(jié)論,這個結(jié)論的推導(dǎo)可以參考Masunipur的Field, Force, Energy and Momentum in Classical Electrodynamics第66頁的example 1:
實值函數(shù)的Fourier變換
假設(shè)f(x)=v(x)+jw(x)f(x)=v(x)+jw(x)f(x)=v(x)+jw(x),因為
e?j2πξx=cos?(2πξx)?jsin?(2πξx)e^{-j 2 \pi \xi x} = \cos(2 \pi \xi x)- j \sin (2 \pi \xi x)e?j2πξx=cos(2πξx)?jsin(2πξx)
所以
Re[F(ξ)]=∫?∞+∞[v(x)cos?(2πξx)+w(x)sin?(2πξx)]dxIm[F(ξ)]=∫?∞+∞[w(x)cos?(2πξx)?v(x)sin?(2πξx)]dxRe[F(\xi)] = \int_{-\infty}^{+\infty} \left[ v(x) \cos(2 \pi \xi x)+w(x) \sin(2 \pi \xi x) \right]dx \\ Im[F(\xi)] =\int_{-\infty}^{+\infty} \left[ w(x) \cos(2 \pi \xi x)-v(x) \sin(2 \pi \xi x) \right]dx Re[F(ξ)]=∫?∞+∞?[v(x)cos(2πξx)+w(x)sin(2πξx)]dxIm[F(ξ)]=∫?∞+∞?[w(x)cos(2πξx)?v(x)sin(2πξx)]dx
如果f?(x)=f(x)f^*(x)=f(x)f?(x)=f(x),即f(x)f(x)f(x)為實值函數(shù),則w(x)=0w(x)=0w(x)=0,
Re[F(ξ)]=∫?∞+∞v(x)cos?(2πξx)dxIm[F(ξ)]=∫?∞+∞?v(x)sin?(2πξx)dxRe[F(\xi)] = \int_{-\infty}^{+\infty} v(x) \cos(2 \pi \xi x)dx \\ Im[F(\xi)] =\int_{-\infty}^{+\infty}-v(x) \sin(2 \pi \xi x) dx Re[F(ξ)]=∫?∞+∞?v(x)cos(2πξx)dxIm[F(ξ)]=∫?∞+∞??v(x)sin(2πξx)dx
如果在此基礎(chǔ)上還滿足f(x)=f(?x)f(x)=f(-x)f(x)=f(?x),即f(x)f(x)f(x)為偶函數(shù),則
∫?∞+∞?v(x)sin?(2πξx)dx=0\int_{-\infty}^{+\infty}-v(x) \sin(2 \pi \xi x) dx = 0∫?∞+∞??v(x)sin(2πξx)dx=0
此時
F(ξ)=∫?∞+∞v(x)cos?(2πξx)dxF(\xi) = \int_{-\infty}^{+\infty} v(x) \cos(2 \pi \xi x)dxF(ξ)=∫?∞+∞?v(x)cos(2πξx)dx
這個形式也被稱為余弦變換(cosine transform)。
例:矩形波的Fourier變換
將這個結(jié)果與Fourier級數(shù)展開進行對比:
傅立葉變換的性質(zhì)
性質(zhì)1 如果F(ξ)=F[f(x)]F(\xi)=\mathcal{F}[f(x)]F(ξ)=F[f(x)],則F(?ξ)=F[f(?x)]F(-\xi) = \mathcal{F}[f(-x)]F(?ξ)=F[f(?x)]
性質(zhì)2 如果F(ξ)=F[f(x)]F(\xi)=\mathcal{F}[f(x)]F(ξ)=F[f(x)],則F?(?ξ)=F[f?(x)]F^*(-\xi) = \mathcal{F}[f^*(x)]F?(?ξ)=F[f?(x)]
性質(zhì)3 如果F(ξ)=F[f(x)]F(\xi)=\mathcal{F}[f(x)]F(ξ)=F[f(x)],則f(?ξ)=F[F(x)]f(-\xi) = \mathcal{F}[F(x)]f(?ξ)=F[F(x)]
例 ej2πξ0xe^{j2 \pi \xi_0 x}ej2πξ0?x與δ(ξ?ξ0)\delta(\xi-\xi_0)δ(ξ?ξ0?)是一對Fourier變換對,根據(jù)性質(zhì)三可以直接得到δ(x?x0)\delta(x-x_0)δ(x?x0?)的Fourier變換為e?j2πξ0xe^{-j2 \pi \xi_0 x}e?j2πξ0?x,驗證:
性質(zhì)四 如果F(ξ)=F[f(x)]F(\xi)=\mathcal{F}[f(x)]F(ξ)=F[f(x)],則∣b∣F(bξ)=F[f(x/b)],?b≠0|b|F(b\xi) = \mathcal{F}[f(x/b)],\forall b \ne 0∣b∣F(bξ)=F[f(x/b)],?b?=0
性質(zhì)五 如果F(ξ)=F[f(x)]F(\xi)=\mathcal{F}[f(x)]F(ξ)=F[f(x)],則e?j2πξx0F(ξ)=F[f(x?x0)]e^{-j 2 \pi \xi x_0}F(\xi) = \mathcal{F}[f(x-x_0)]e?j2πξx0?F(ξ)=F[f(x?x0?)]
性質(zhì)六 如果F(ξ)=F[f(x)]F(\xi)=\mathcal{F}[f(x)]F(ξ)=F[f(x)],則(j2πξ)nF(ξ)=F[f(n)(x)](j 2 \pi \xi)^nF(\xi) = \mathcal{F}[f^{(n)}(x)](j2πξ)nF(ξ)=F[f(n)(x)],以n=1n=1n=1為例,
常用Fourier變換對
| 1 | δ(ξ)\delta(\xi)δ(ξ) |
| ej2πξ0xe^{j 2 \pi \xi_0 x}ej2πξ0?x | δ(ξ?ξ0)\delta(\xi-\xi_0)δ(ξ?ξ0?) |
| δ(x)\delta(x)δ(x) | 111 |
| δ(x?x0)\delta(x-x_0)δ(x?x0?) | e?j2πξx0e^{-j 2 \pi \xi x_0}e?j2πξx0? |
| cos?(2πξ0x)\cos(2 \pi \xi_0 x)cos(2πξ0?x) | [δ(ξ?ξ0)+δ(ξ+ξ0)]/2[\delta(\xi-\xi_0)+\delta(\xi+\xi_0)]/2[δ(ξ?ξ0?)+δ(ξ+ξ0?)]/2 |
| sin?(2πξ0x)\sin(2 \pi \xi_0 x)sin(2πξ0?x) | [δ(ξ?ξ0)?δ(ξ+ξ0)]/(2j)[\delta(\xi-\xi_0)-\delta(\xi+\xi_0)]/(2j)[δ(ξ?ξ0?)?δ(ξ+ξ0?)]/(2j) |
| fp(x)=∑n=?∞n=+∞cnej2πnξ0xf_p(x)=\sum_{n=-\infty}^{n=+\infty}c_ne^{j 2 \pi n\xi_0 x}fp?(x)=∑n=?∞n=+∞?cn?ej2πnξ0?x | ∑n=?∞n=+∞cnδ(ξ?nξ0)\sum_{n=-\infty}^{n=+\infty}c_n \delta(\xi-n \xi_0)∑n=?∞n=+∞?cn?δ(ξ?nξ0?) |
| comb(x)=∑n=?∞+∞δ(x?n)comb(x)=\sum_{n=-\infty}^{+\infty}\delta(x-n)comb(x)=∑n=?∞+∞?δ(x?n) | ∑n=?∞+∞e?j2πnξ=comb(ξ)\sum_{n=-\infty}^{+\infty}e^{-j2 \pi n \xi}=comb(\xi)∑n=?∞+∞?e?j2πnξ=comb(ξ) |
| rect(x)rect(x)rect(x) | sinc(ξ)sinc(\xi)sinc(ξ) |
| Gaus(x)Gaus(x)Gaus(x) | Gaus(ξ)Gaus(\xi)Gaus(ξ) |
| e?xstep(x)e^{-x}step(x)e?xstep(x) | 11+j2πξ\frac{1}{1+j 2 \pi \xi}1+j2πξ1? |
第一個的推導(dǎo)很簡單,
第二個就是性質(zhì)三的例,第三個在第一個的基礎(chǔ)上用性質(zhì)三,第四個在第三個的基礎(chǔ)上用性質(zhì)五,第五個和第六個分別寫成
cos?(2πξ0x)=ej2πξ0x+e?j2πξ0x2sin?(2πξ0x)=ej2πξ0x?e?j2πξ0x2j\cos(2\pi \xi_0 x) = \frac{e^{j2 \pi \xi_0 x}+e^{-j2 \pi \xi_0 x}}{2} \\ \sin(2 \pi \xi_0 x) = \frac{e^{j2 \pi \xi_0 x}-e^{-j2 \pi \xi_0 x}}{2j}cos(2πξ0?x)=2ej2πξ0?x+e?j2πξ0?x?sin(2πξ0?x)=2jej2πξ0?x?e?j2πξ0?x?
然后用第二個變換對、性質(zhì)一和線性性即可,第七個說的是周期函數(shù)寫成Fourier級數(shù)展開的形式再做Fourier變換,在第二個的基礎(chǔ)上用線性性得證,第八個涉及comb函數(shù),它有一點特殊,如果直接推導(dǎo)可以得到
但比價巧合的是它也是一個周期為1的函數(shù),所以可以用第七個Fourier變換對,先寫出它的Fourier級數(shù)展開
然后做Fourier變換,
第九個變換對就是上文矩形波那個例子,第十個變換對的推導(dǎo)如下,
第十一個變換對的推導(dǎo)如下:
總結(jié)
以上是生活随笔為你收集整理的UA OPTI512R 傅立叶光学导论13 傅立叶变换及其性质的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: android ListView实现下拉
- 下一篇: 第一次作业--四则运算题目生成程序