傅里叶变换

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 ) , − π &lt; x &lt; π , 0 , f(x)=\left\{ \begin{aligned} &amp;sin(3x),-\pi&lt;x&lt;\pi,\\ &amp;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

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

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值