傅里叶变换

    xiaoxiao2022-06-30  162

    UTF8gbsn

    Introduction

    我们回忆一下傅里叶级数。傅里叶级数要求周期性函数。而别展开式实际上是离散频率的 信号。也就是说,展开式中的频率不是连续的。只有离散的信号。而且仅有下面的频率信 号。

    1 2 l , 2 2 l , 3 2 l , . . . , n 2 l , ( n = , 1 , 2 , 3 , . . . . ) \frac{1}{2l},\frac{2}{2l},\frac{3}{2l},...,\frac{n}{2l}, (n=,1,2,3,....) 2l1,2l2,2l3,...,2ln,(n=,1,2,3,....)

    Complex Fourier series

    不过在这之前,我们来看看傅里叶级数的复数表示形式。我们都知道傅里叶级数的一般表 示形式为,傅里叶级数的复数表示形式和我们一般用的常规表示形式。没有任何区别,只 是因为复数表示形式更为简单而已。

    f ( x ) = a 0 2 + ∑ n = 1 ∞ [ a n c o s ( n π x l ) + b n s i n ( n π x l ) ] f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}[a_ncos(\frac{n\pi x}{l})+b_nsin(\frac{n\pi x}{l})] f(x)=2a0+n=1[ancos(lnπx)+bnsin(lnπx)]

    其中 a n , b n a_n,b_n an,bn, 为

    a n = 1 l ∫ − l l f ( x ) c o s ( n π x l ) d x , a_n= \frac{1}{l} \int_{-l}^{l}{f(x)cos(\frac{n\pi x}{l})dx}, an=l1llf(x)cos(lnπx)dx, b n = 1 l ∫ − l l f ( x ) s i n ( n π x l ) d x , b_n= \frac{1}{l} \int_{-l}^{l}{f(x)sin(\frac{n\pi x}{l})dx}, bn=l1llf(x)sin(lnπx)dx,

    利用欧拉公式

    c o s ( t ) = e i t + e − i t 2 , s i n ( t ) = e i t − e − i t 2 i cos(t)=\frac{e^{it}+e^{-it}}{2}, sin(t)=\frac{e^{it}-e^{-it}}{2i} cos(t)=2eit+eit,sin(t)=2ieiteit

    傅里叶级数可以变为。

    a 0 2 + ∑ n = 1 ∞ [ a n − i b n 2 e i n π x l + a n + i b n 2 e − i n π x l ] \frac{a_0}{2}+\sum_{n=1}^{\infty}[\frac{a_n-ib_n}{2}e^{i \frac{n\pi x}{l}}+\frac{a_n+ib_n}{2}e^{-i \frac{n\pi x}{l}}] 2a0+n=1[2anibneilnπx+2an+ibneilnπx]

    c 0 = a 0 2 = 1 2 l ∫ − l l f ( x ) d x c_0=\frac{a_0}{2}=\frac{1}{2l} \int_{-l}^{l}{f(x)dx} c0=2a0=2l1llf(x)dx

    c n = 1 2 l ∫ − l l f ( x ) e − i n π x l d x = a n − i b n 2 c_n=\frac{1}{2l} \int_{-l}^{l}{f(x)e^{-i \frac{n\pi x}{l}}dx}=\frac{a_n-ib_n}{2} cn=2l1llf(x)eilnπxdx=2anibn

    c − n = 1 2 l ∫ − l l f ( x ) e i n π x l d x = a n + i b n 2 c_{-n}=\frac{1}{2l} \int_{-l}^{l}{f(x)e^{i \frac{n\pi x}{l}}dx}=\frac{a_n+ib_n}{2} cn=2l1llf(x)eilnπxdx=2an+ibn

    最后,我们把形式统一写成

    f ( x ) = ∑ n = − ∞ ∞ α n e i n π x l f(x)=\sum_{n=-\infty}^{\infty}\alpha_n e^{\frac{i n \pi x}{l}} f(x)=n=αnelinπx

    α n = 1 2 l ∫ − l l f ( x ) e − i n π x l d x , ( n = 0 , ± 1 , ± 2 , . . . ) \alpha_n=\frac{1}{2l}\int_{-l}^{l}{f(x)e^{-i \frac{n \pi x}{l}}dx},(n=0,\pm 1,\pm 2,...) αn=2l1llf(x)eilnπxdx,(n=0,±1,±2,...)

    周期无穷大

    周期无穷大。也就是函数整个定义在实数轴上。上一小节中,我们用的周期是 l l l. 当周期为无 穷大时。我们表示为 l → ∞ l\rightarrow \infty l

    f ( x ) = lim ⁡ l → ∞ [ ∑ n = − ∞ ∞ ( 1 2 l ∫ − l l f ( t ) e − i n π t / l d t ) e i n π x / l ] = lim ⁡ l → ∞ [ ∑ n = − ∞ ∞ ( 1 2 l ∫ − l l f ( t ) e i n π ( x − t ) / l d t ) ] \left. \begin{aligned} f(x)=& \lim_{l\rightarrow \infty}\left[\sum_{n=-\infty}^{\infty}(\frac{1}{2l}\int_{-l}^{l}{f(t)e^{-in\pi t/l}dt})e^{in\pi x/l}\right]\\ =& \lim_{l\rightarrow \infty}\left[\sum_{n=-\infty}^{\infty}(\frac{1}{2l}\int_{-l}^{l}{f(t)e^{in\pi (x-t)/l}dt})\right] \end{aligned} \right. f(x)==llim[n=(2l1llf(t)einπt/ldt)einπx/l]llim[n=(2l1llf(t)einπ(xt)/ldt)]

    我们先来看看 l → ∞ l\rightarrow \infty l 会产生的一些变化。

    当l为一定数值的时候。原始信号又离散的一些信号合成。 1 2 l , 2 2 l , 3 2 l , 4 2 l , . . . , n 2 l , ( n = , 1 , 2 , 3 , . . . ) \frac{1}{2l}, \frac{2}{2l}, \frac{3}{2l},\frac{4}{2l},...,\frac{n}{2l}, (n=,1,2,3,...) 2l1,2l2,2l3,2l4,...,2ln,(n=,1,2,3,...)

    l → ∞ l\rightarrow \infty l意味着什么?

    我们看到,当 l l l不断增大的时候,同一个区间里面的离散频率变得越来越稠密了。当 l → ∞ l\rightarrow \infty l,我们得到所有的频率成分。也就是说,我们得到了连续的频 率成分。

    Fourier Transform

    接前面的介绍,我们令 λ n = n π l , Δ λ = λ n + 1 − λ n = π l \lambda_n=\frac{n\pi}{l},\Delta \lambda=\lambda_{n+1}-\lambda_n=\frac{\pi}{l} λn=lnπ,Δλ=λn+1λn=lπ我们得到

    f ( x ) = lim ⁡ l → ∞ ∑ n = − ∞ ∞ [ 1 2 π ∫ − l l f ( t ) e λ n i ( x − t ) d t ] Δ λ f(x)=\lim_{l\rightarrow \infty} \sum_{n=-\infty}^{\infty}\left[\frac{1}{2\pi}\int_{-l}^{l}{f(t)e^{\lambda_n i (x-t)}dt} \right]\Delta \lambda f(x)=llimn=[2π1llf(t)eλni(xt)dt]Δλ

    接着我们对当 l → ∞ l\rightarrow \infty l时,我们可以得到下面的等式

    f ( x ) = 1 2 π ∫ − ∞ ∞ ∫ − ∞ ∞ f ( t ) e i λ ( x − t ) d t d λ f(x)=\frac{1}{2\pi}\int_{-\infty}^{\infty}{\int_{-\infty}^{\infty}{f(t)e^{i\lambda(x-t)}dt}d\lambda} f(x)=2π1f(t)eiλ(xt)dtdλ f ( x ) = 1 2 π ∫ − ∞ ∞ ( 1 2 π ∫ − ∞ ∞ f ( t ) e − i λ t d t ) e i λ x d x f(x)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}{\left(\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}{f(t)e^{-i\lambda t}dt}\right)e^{i\lambda x}dx} f(x)=2π 1(2π 1f(t)eiλtdt)eiλxdx 其中我们把

    f ^ ( λ ) = 1 2 π ∫ − ∞ ∞ f ( t ) e i λ t d t \widehat{f}(\lambda)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}{f(t)e^{i\lambda t}dt} f (λ)=2π 1f(t)eiλtdt 叫做傅里叶变换。

    那么现在让我们来看看 f ^ ( λ ) \widehat{f}(\lambda) f (λ) vs α n \alpha_n αn 的比较。 α n = 1 2 l ∫ − l l f ( x ) e − i n π x l d x , ( n = 0 , ± 1 , ± 2 , . . . ) \alpha_n=\frac{1}{2l}\int_{-l}^{l}{f(x)e^{-i \frac{n \pi x}{l}}dx},(n=0,\pm 1,\pm 2,...) αn=2l1llf(x)eilnπxdx,(n=0,±1,±2,...)

    我们可以看出 λ \lambda λ代表的就是频率成分。而 f ^ \widehat{f} f α n \alpha_n αn的取极值的 表示方式。

    Original Form

    让我们来看看如何从傅里叶变换,还原计算出原始的没有引入复数表示形式的傅里叶级数 形式。

    α n = 1 2 l ∫ − l l f ( x ) e − i n π x l d x = a n − i b n 2 \alpha_n=\frac{1}{2l} \int_{-l}^{l}{f(x)e^{-i \frac{n\pi x}{l}}dx}=\frac{a_n-ib_n}{2} αn=2l1llf(x)eilnπxdx=2anibn

    α − n = 1 2 l ∫ − l l f ( x ) e i n π x l d x = a n + i b n 2 \alpha_{-n}=\frac{1}{2l} \int_{-l}^{l}{f(x)e^{i \frac{n\pi x}{l}}dx}=\frac{a_n+ib_n}{2} αn=2l1llf(x)eilnπxdx=2an+ibn

    如何计算 a n , b n a_n,b_n an,bn就十分的简单了。让我们换成傅里叶变换的一般表示形式,就是

    α λ = 1 2 π f ^ ( λ ) = a λ − i b λ 2 \alpha_{\lambda}=\frac{1}{\sqrt{2\pi}}\widehat{f}(\lambda)=\frac{a_{\lambda}-ib_{\lambda}}{2} αλ=2π 1f (λ)=2aλibλ α − λ = 1 2 π f ^ ( − λ ) = a λ + i b λ 2 \alpha_{-\lambda}=\frac{1}{\sqrt{2\pi}}\widehat{f}(-\lambda)=\frac{a_{\lambda}+ib_{\lambda}}{2} αλ=2π 1f (λ)=2aλ+ibλ

    解出 a λ , b λ a_{\lambda},b_{\lambda} aλ,bλ,就可以恢复出 λ \lambda λ成分信号 a λ c o s ( λ x ) + b λ s i n ( λ x ) a_{\lambda}cos(\lambda x)+b_{\lambda}sin(\lambda x) aλcos(λx)+bλsin(λx)

    For Example

    我们具体来看看一个函数的傅里叶变换是个什么样子。

    f ( x ) = { s i n ( 3 x ) , − π < x < π , 0 , f(x)=\left\{ \begin{aligned} &sin(3x),-\pi<x<\pi,\\ &0, \end{aligned} \right. f(x)={sin(3x),π<x<π,0, 他的傅里叶变换函数为

    3 i 2 π s i n ( π λ ) λ 2 − 9 \frac{3i \sqrt{\frac{2}{\pi}}sin(\pi \lambda)}{\lambda^2-9} λ293iπ2 sin(πλ)

    这个函数只有虚部我们把图像画出来,如下。

    The End

    傅里叶变换,和傅里叶级数是紧密联系的。可以说傅里叶变换就是把傅里叶级数取极限。 然后系数的表达式也就是傅里叶变换函数。傅里叶变化,可以让我们观察一个函数


    最新回复(0)