本文是对斯坦福大学公开课-傅里叶变换及其应用的一些学习总结
对应的课程视频链接为:课程链接
课程用到的讲义下载地址:讲义链接
傅里叶变换的内容对应课程5-8节,讲义第二章。
由傅里叶级数到傅里叶变换
在前面介绍了利用傅里叶级数将一个周期函数表达为由一组正交的函数组线性表出的形式,那么对于那些非周期函数,是否也能将其以这样的形式表示呢?注意到,对于每个非周期函数,其实可以将其看做是周期趋于无限大的一个周期函数,由这一点出发,令
T−>∞
,写出对应的傅里叶级数:
f(t)=∑n=−∞∞cne2πint/T
cn=1T∫T/2−T/2e−2πint/Tf(t)dt
似乎这样就完成了?来考虑一个具体的实例:
现在计算这个非周期函数的傅里叶系数:
cn=1T∫T/2−T/2e−2πint/Tf(t)dt=1T∫1/2−1/2e−2πint/Tdt
很显然,积分式是有界的,当
T−>∞
时,
cn−>0
,这意味着傅里叶级数每一项的系数都趋近于0了,也就是说在这种情况下,傅里叶级数的每一项都消失了!
因此如果采用这样的
cn
用傅里叶级数是没办法表示这样的非周期函数的,那么现在的问题在于
cn−>0
,可以考虑利用一些技巧让
cn
不趋近于0。
对于傅里叶级数,做如下处理:
f(t)=∑n=−∞∞1T⋅Tcne2πint/T=1T∑n=−∞∞Tcne2πint/T
考虑新的傅里叶系数:
Tcn=TT∫T/2−T/2e−2πint/Tf(t)dt=∫T/2−T/2e−2πint/Tf(t)dt
很显然这个新的傅里叶系数的值取决于积分式,当
T−>∞
时,并不是一个趋近于0的式子。这就达到了目的。
原本的傅里叶系数的自变量是
n
,现在把新的傅里叶系数的自变量设为
n/T
,因此在
T−>∞
时,将新的傅里叶系数表达为:
cs=∫∞−∞e−2πinsf(t)dt
更标准的写法为:
f^(s)=∫∞−∞e−2πinsf(t)dt
其中
s=nT
有了新的傅里叶系数,再来看此时的傅里叶级数:
f(t)=1T∑n=−∞∞f^(n/T)e2πint/T
此时,这个式子就是一个积分的定义式!
f(t)=∫∞−∞f^(s)e2πistds
这里将傅里叶级数中先乘一个T再除以一个T实现了两项功能:1.让傅里叶系数不趋近于0 2.让傅里叶级数的求和式变成了一个积分式。因此这个操作是非常巧妙的。
这就得到了傅里叶变换及其逆变换,求解傅里叶系数
f^(s)
的过程叫做傅里叶变换,而由傅里叶系数得到原函数的过程(即最后的这个积分式)叫做傅里叶逆变换。
可以看到,傅里叶变换和傅里叶逆变换都是在傅里叶级数中
T−>∞
时,对其做一定的变形得到的,即傅里叶变换是傅里叶级数的极限形式。
下面做一个傅里叶变换和逆变换的小结:
傅里叶变换:对于一个信号
f(t)
,它的傅里叶变换为:
f^(s)=∫∞−∞f(t)e−2πistdt
这个函数的函数值是复数,当
s=0
时有:
f^(0)=∫∞−∞f(t)dt
f(t)
大多数情况下都是实值函数,因此
f^(0)
大多数情况下都是实数,但是
f^
的其他函数值有可能为复数。
傅里叶变换将信号由时域变换到频域。
傅里叶逆变换:对于一个频域下的信号
g(s)
,它的逆傅里叶变换为:
gˇ(t)=∫∞−∞e2πistg(s)ds
当
t=0
时,有
gˇ(0)=∫∞−∞g(s)ds=∫∞−∞f^(s)ds
注意到等式左侧一般都是实数,而等式右侧是对复数的积分,结果得到的是一个实数。
傅里叶变换将一个信号由频域变换到时域。
因此傅里叶变换和傅里叶逆变换实现了信号在时域和频域的切换。
还需注意帕塞瓦尔定理也不难由傅里叶级数推广到傅里叶变换中:
∫∞−∞|f(t)|2dt=∫∞−∞|f^(s)|2ds
实例:高斯函数的傅里叶变换
最简单的高斯函数形式为:
f(x)=e−x2
,这个函数的形状应该很熟悉了,它的积分为:
∫∞−∞e−x2dx=x−−√
高斯函数是没办法找到它的原函数的,因此这个积分不能按照常规的方法来计算。一般是利用二重积分来证明这个等式。
下面来求标准化的高斯函数的傅里叶变换:
已知高斯函数满足:
∫∞−∞e−πx2dx=1
求
e−πx2
的傅里叶变换:
f^(s)=∫∞−∞e−2πisxe−πx2dx
两边对s求导得:
ddsf^(s)=−2πix∫∞−∞e−2πisxe−πx2dx=i∫∞−∞e−2πisxde−πx2
=i(e−πx2e−2πisx|∞−∞+(2πis)∫∞−∞e−2πisxe−πx2dx)=−2πsf^(s)
解这个常微分方程可得:
f^(s)=f^(0)e−πs2
而
f^(0)=∫∞−∞e−πx2=1
因此
f^(s)=e−πs2
即高斯函数
e−πx2
的傅里叶变换是它本身。
傅里叶变换的性质
首先有两个很直观的性质,对一个时域下的信号先作傅里叶变换再作傅里叶逆变换,得到的应该就是原信号。同理对一个频域下的信号先作逆傅里叶变换再作傅里叶变换,得到的应该也是原信号。
对偶性:
对于一个信号
f(t)
,考虑
f^(−s)
和
fˇ(−t)
具有什么性质:
f^(−s)=∫∞−∞e−2πi(−s)tf(t)dt=∫∞−∞e2πsitf(t)dt=fˇ(s)
fˇ(−t)=∫∞−∞e2πis(−t)f(s)ds=∫∞−∞e−2πistf(s)ds=f^(t)
这两个式子看上去十分错乱,s是频域中的自变量,它怎么可以被作为傅里叶逆变换的自变量呢?同样,t是时域中的自变量,它怎么可以被作为傅里叶变换的自变量呢?
对于这两个等式,比较好的处理方法应该是不要试图从理论层面来解释,而应该就把它们当作是两个数学上的等式,仅此而已,它所表达的就是一种数学上的等量关系,这就足够了。
由这两个等量关系,可以推导出:
f^(s)^=fˇ(−s)^=f(−s)
f^(−s)^=fˇ(s)^=f(s)
即对一个时域下的信号,先作傅里叶变换,再作傅里叶变换,得到的应该是原信号关于y轴对称的反转信号。
时延理论:考虑在时域下的信号在水平方向上平移,那么对应的频域下的信号会怎么变化?即
f^(t+b)
与
f^(t)
有什么关系?
f^(t+b)=∫∞−∞e−2πistf(t+b)dt
设
u=t+b
,则
t=u−b
∫∞−∞e−2πistf(t+b)dt=∫∞−∞e−2πis(u−b)f(u)du=e2πisb∫∞−∞e−2πisuf(u)du=e2πisbf^(s)
注意到频域的振幅是没有发生改变的,
|e2πisb|=1
伸缩理论:考虑在时域下的信号在水平方向上进行压缩或者拉伸,那么对应的频域下的信号会怎么变化?即
f^(at)
与
f^(t)
有什么关系?
当
a>0
时:
f^(s)=∫∞−∞e−2πistf(at)dt
设
u=at
,则
t=u/a
,
dt=du/a
f^(s)=1a∫∞−∞e−2πisu/af(u)du=1af^(sa)
当
a<0
时:
f^(s)=∫∞−∞e−2πistf(at)dt
设
u=at
,则
t=u/a
,
dt=du/a
f^(s)=1a∫−∞∞e−2πisu/af(u)du=−1af^(sa)
因此:若
f(t)
对应
F(s)
,那么
f(at)
就对应
1|a|F(sa)
。
由此不难看出,当信号在时域中进行压缩时,在频域中会进行拉伸并变矮;当信号在时域中进行拉伸时,在频域中会进行压缩并变高。