深刻理解傅里叶变换资料.docx
《深刻理解傅里叶变换资料.docx》由会员分享,可在线阅读,更多相关《深刻理解傅里叶变换资料.docx(23页珍藏版)》请在冰豆网上搜索。
![深刻理解傅里叶变换资料.docx](https://file1.bdocx.com/fileroot1/2023-1/29/df707f54-80f2-467e-851d-b8aefb2ccd9b/df707f54-80f2-467e-851d-b8aefb2ccd9b1.gif)
深刻理解傅里叶变换资料
理解离散傅立叶变换
(一)------傅立叶变换的由来
一、傅立叶变换的提出
让我们先看看为什么会有傅立叶变换?
傅立叶是一位法国数学家和物理学家的名字,英语原名是JeanBaptisteJosephFourier(1768-1830),Fourier对热传递很感兴趣,于1807年在法国科学学会上发表了一篇论文,运用正弦曲线来描述温度分布,论文里有个在当时具有争议性的决断:
任何连续周期信号可以由一组适当的正弦曲线组合而成。
当时审查这个论文的人,其中有两位是历史上著名的数学家拉格朗日(JosephLouisLagrange,1736-1813)和拉普拉斯(PierreSimondeLaplace,1749-1827),当拉普拉斯和其它审查者投票通过并要发表这个论文时,拉格朗日坚决反对,在近50年的时间里,拉格朗日坚持认为傅立叶的方法无法表示带有棱角的信号,如在方波中出现非连续变化斜率。
法国科学学会屈服于拉格朗日的威望,拒绝了傅立叶的工作,幸运的是,傅立叶还有其它事情可忙,他参加了政治运动,随拿破仑远征埃及,法国大革命后因会被推上断头台而一直在逃避。
直到拉格朗日死后15年这个论文才被发表出来。
拉格朗日是对的:
正弦曲线无法组合成一个带有棱角的信号。
但是,我们可以用正弦曲线来非常逼近地表示它,逼近到两种表示方法不存在能量差别,基于此,傅立叶是对的。
为什么我们要用正弦曲线来代替原来的曲线呢?
如我们也还可以用方波或三角波来代替呀,分解信号的方法是无穷的,但分解信号的目的是为了更加简单地处理原来的信号。
用正余弦来表示原信号会更加简单,因为正余弦拥有原信号所不具有的性质:
正弦曲线保真度。
一个正弦曲线信号输入后,输出的仍是正弦曲线,只有幅度和相位可能发生变化,但是频率和波的形状仍是一样的。
且只有正弦曲线才拥有这样的性质,正因如此我们才不用方波或三角波来表示。
二、傅立叶变换分类
根据原信号的不同类型,我们可以把傅立叶变换分为四种类别:
1
非周期性连续信号
傅立叶变换(FourierTransform)
2
周期性连续信号
傅立叶级数(FourierSeries)
3
非周期性离散信号
离散时域傅立叶变换(DiscreteTimeFourierTransform)
4
周期性离散信号
离散傅立叶变换(DiscreteFourierTransform)
下图是四种原信号图例:
这四种傅立叶变换都是针对正无穷大和负无穷大的信号,即信号的的长度是无穷大的,我们知道这对于计算机处理来说是不可能的,那么有没有针对长度有限的傅立叶变换呢?
没有。
因为正余弦波被定义成从负无穷小到正无穷大,我们无法把一个长度无限的信号组合成长度有限的信号。
面对这种困难,方法是把长度有限的信号表示成长度无限的信号,可以把信号无限地从左右进行延伸,延伸的部分用零来表示,这样,这个信号就可以被看成是非周期性离散信号,我们就可以用到离散时域傅立叶变换的方法。
还有,也可以把信号用复制的方法进行延伸,这样信号就变成了周期性离解信号,这时我们就可以用离散傅立叶变换方法进行变换。
这里我们要学的是离散信号,对于连续信号我们不作讨论,因为计算机只能处理离散的数值信号,我们的最终目的是运用计算机来处理信号的。
但是对于非周期性的信号,我们需要用无穷多不同频率的正弦曲线来表示,这对于计算机来说是不可能实现的。
所以对于离散信号的变换只有离散傅立叶变换(DFT)才能被适用,对于计算机来说只有离散的和有限长度的数据才能被处理,对于其它的变换类型只有在数学演算中才能用到,在计算机面前我们只能用DFT方法,后面我们要理解的也正是DFT方法。
这里要理解的是我们使用周期性的信号目的是为了能够用数学方法来解决问题,至于考虑周期性信号是从哪里得到或怎样得到是无意义的。
每种傅立叶变换都分成实数和复数两种方法,对于实数方法是最好理解的,但是复数方法就相对复杂许多了,需要懂得有关复数的理论知识,不过,如果理解了实数离散傅立叶变换(realDFT),再去理解复数傅立叶就更容易了,所以我们先把复数的傅立叶放到一边去,先来理解实数傅立叶变换,在后面我们会先讲讲关于复数的基本理论,然后在理解了实数傅立叶变换的基础上再来理解复数傅立叶变换。
还有,这里我们所要说的变换(transform)虽然是数学意义上的变换,但跟函数变换是不同的,函数变换是符合一一映射准则的,对于离散数字信号处理(DSP),有许多的变换:
傅立叶变换、拉普拉斯变换、Z变换、希尔伯特变换、离散余弦变换等,这些都扩展了函数变换的定义,允许输入和输出有多种的值,简单地说变换就是把一堆的数据变成另一堆的数据的方法。
三、一个关于实数离散傅立叶变换(RealDFT)的例子
先来看一个变换实例,下图是一个原始信号图像:
这个信号的长度是16,于是可以把这个信号分解9个余弦波和9个正弦波(一个长度为N的信号可以分解成N/2+1个正余弦信号,这是为什么呢?
结合下面的18个正余弦图,我想从计算机处理精度上就不难理解,一个长度为N的信号,最多只能有N/2+1个不同频率,再多的频率就超过了计算机所能所处理的精度范围),如下图:
9个余弦信号:
9个正弦信号:
把以上所有信号相加即可得到原始信号,至于是怎么分别变换出9种不同频率信号的,我们先不急,先看看对于以上的变换结果,在程序中又是该怎么表示的,我们可以看看下面这个示例图:
上图中左边表示时域中的信号,右边是频域信号表示方法,从左向右表示正向转换(ForwardDFT),从右向左表示逆向转换(InverseDFT),用小写x[]表示信号在每个时间点上的幅度值数组,用大写X[]表示每种频率的幅度值数组,因为有N/2+1种频率,所以该数组长度为N/2+1,X[]数组又分两种,一种是表示余弦波的不同频率幅度值:
ReX[],另一种是表示正弦波的不同频率幅度值:
ImX[],Re是实数(Real)的意思,Im是虚数(Imagine)的意思,采用复数的表示方法把正余弦波组合起来进行表示,但这里我们不考虑复数的其它作用,只记住是一种组合方法而已,目的是为了便于表达(在后面我们会知道,复数形式的傅立叶变换长度是N,而不是N/2+1)。
下一节我们将来看一下实数傅立叶变换的具体方法。
理解离散傅立叶变换
(二)------实数形式离散傅立叶变换(RealDFT)
上一节我们看到了一个实数形式离散傅立叶变换的例子,通过这个例子能够让我们先对傅立叶变换有一个较为形象的感性认识,现在就让我们来看看实数形式离散傅立叶变换的正向和逆向是怎么进行变换的。
在此,我们先来看一下频率的多种表示方法。
一、 频域中关于频率的四种表示方法
1、 序号表示方法,根据时域中信号的样本数取0~N/2,用这种方法在程序中使用起来可以更直接地取得每种频率的幅度值,因为频率值跟数组的序号是一一对应的:
X[k],取值范围是0~N/2;
2、 分数表示方法,根据时域中信号的样本数的比例值取0~0.5:
X[ƒ],ƒ=k/N,取值范围是0~N/2;
3、 用弧度值来表示,把ƒ乘以一个2π得到一个弧度值,这种表示方法叫做自然频率(naturalfrequency):
X[ω],ω=2πƒ=2πk/N,取值范围是0~π;
4、 以赫兹(Hz)为单位来表示,这个一般是应用于一些特殊应用,如取样率为10kHz表示每秒有10,000个样本数:
取值范围是0到取样率的一半。
二、 DFT基本函数
ck[i]=cos(2πki/N)
sk[i]=sin(2πki/N)
其中k表示每个正余弦波的频率,如为2表示在0到N长度中存在两个完整的周期,10即有10个周期,如下图:
上图中至于每个波的振幅(amplitude)值(ReX[k],ImX[k])是怎么算出来的,这个是DFT的核心,也是最难理解的部分,我们先来看看如何把分解出来的正余弦波合成原始信号(InverseDFT)。
三、 合成运算方法(RealInverseDFT)
DFT合成等式:
如果有学过傅立叶级数,对这个等式就会有似曾相识的感觉,不错!
这个等式跟傅立叶级数是非常相似的:
当然,差别是肯定是存在的,因为这两个等式是在两个不同条件下运用的,至于怎么证明DFT合成公式,这个我想需要非常强的高等数学理论知识了,这是研究数学的人的工作,对于普通应用者就不需要如此的追根究底了,但是傅立叶级数是好理解的,我们起码可以从傅立叶级数公式中看出DFT合成公式的合理性。
DFT合成等式中的Im
[k]和Re
[k]跟ImX[k]和ReX[k]是不一样的,下面是转换方法:
但k等于0和N/2时,实数部分的计算要用下面的等式:
上面四个式中的N是时域中点的总数,k是从0到N/2的序号。
为什么要这样进行转换呢?
这个可以从频谱密度(spectraldensity)得到理解,如下图就是个频谱图:
这是一个频谱图,横坐标表示频率大小,纵坐标表示振幅大小,原始信号长度为N(这里是32),经DFT转换后得到的17个频率的频谱,频谱密度表示每单位带宽中为多大的振幅,那么带宽是怎么计算出来的呢?
看上图,除了头尾两个,其余点的所占的宽度是2/N,这个宽度便是每个点的带宽,头尾两个点的带宽是1/N,而ImX[k]和ReX[k]表示的是频谱密度,即每一个单位带宽的振幅大小,但Im[k]和Re[k]表示2/N(或1/N)带宽的振幅大小,所以Im[k]和Re[k]分别应当是ImX[k]和ReX[k]的2/N(或1/N)。
频谱密度就象物理中物质密度,原始信号中的每一个点就象是一个混合物,这个混合物是由不同密度的物质组成的,混合物中含有的每种物质的质量是一样的,除了最大和最小两个密度的物质外,这样我们只要把每种物质的密度加起来就可以得到该混合物的密度了,又该混合物的质量是单位质量,所以得到的密度值跟该混合物的质量值是一样的。
至于为什么虚数部分是负数,这是为了跟复数DFT保持一致,这个我们将在后面会知道这是数学计算上的需要(ImX[k]计算时加上了一个负号,Im[k]再加上负号,结果便是正的,等于没有变化)。
如果已经得到了DFT结果,这时要进行逆转换,即合成原始信号,则可按如下步骤进行转换:
1、 先根据上面四个式子计算得出Im[k]和Re[k]的值;
2、 再根据DFT合成等式得到原始信号数据。
下面是用BASIC语言来实现的转换源代码:
100‘DFT逆转换方法
110‘/XX[]数组存储计算结果(时域中的原始信号)
120‘/REX[]数组存储频域中的实数分量,IMX[]为虚分量
130‘
140DIMXX[511]
150DIMREX[256]
160DIMIMX[256]
170‘
180PI=3.14159265
190N%=512
200‘
210GOSUBXXXX‘转到子函数去获取REX[]和IMX[]数据
220‘
230‘
240‘
250FORK%=0TO256
260REX[K%]=REX[K%]/(N%/2)
270IMX[K%]=-IMX[K%]/(N%/2)
280NEXTk%
290‘
300REX[0]=REX[0]/N
310REX[256]=REX[256]/N
320‘
330‘初始化XX[]数组
340FORI%=0TO511
350XX[I%]=0
360NEXTI%
370‘
380‘
390‘
400‘
410‘
420FORK%=0TO256
430FORI%=0TO511
440‘
450XX[I%]=XX[I%]+REX[K%]*COS(2*PI*K%*I%/N%)
460XX[I%]=XX[I%]+IMX[K%]*SIN(2*PI*K%*I%/N%)
470‘
480NEXTI%
490NEXTK%
500‘
510END
上面代码中420至490换成如下形式也许更好理解,但结果都是一样的:
420FORI%=0TO511
430FORK%=0TO256
440‘
450XX[I%]=XX[I%]+REX[K%]*COS(2*PI*K%*I%/N%)
460XX[I%]=XX[I%]+IMX[K%]*SIN(2*PI*K%*I%/N%)
470‘
480NEXTI%
490NEXTK%
四、 分解运算方法(DFT)
有三种完全不同的方法进行DFT:
一种方法是通过联立方程进行求解,从代数的角度看,要从N个已知值求N个未知值,需要N个联立方程,且N个联立方程必须是线性独立的,但这是这种方法计算量非常的大且极其复杂,所以很少被采用;第二种方法是利用信号的相关性(correlation)进行计算,这个是我们后面将要介绍的方法;第三种方法是快速傅立叶变换(FFT),这是一个非常具有创造性和革命性的的方法,因为它大大提高了运算速度,使得傅立叶变换能够在计算机中被广泛应用,但这种算法是根据复数形式的傅立叶变换来实现的,它把N个点的信号分解成长度为N的频域,这个跟我们现在所进行的实域DFT变换不一样,而且这种方法也较难理解,这里我们先不去理解,等先理解了复数DFT后,再来看一下FFT。
有一点很重要,那就是这三种方法所得的变换结果是一样的,经过实践证明,当频域长度为32时,利用相关性方法进行计算效率最好,否则FFT算法效率较高。
现在就让我们来看一下相关性算法。
利用信号的相关性(correlation)可以从噪声背景中检测出已知的信号,我们也可以利用这个方法检测信号波中是否含有某个频率的信号波:
把一个待检测信号波乘以另一个信号波,得到一个新的信号波,再把这个新的信号波所有的点进行相加,从相加的结果就可以判断出这两个信号的相似程度。
如下图:
上面a和b两个图是待检测信号波,图a很明显可以看出是个3个周期的正弦信号波,图b的信号波则看不出是否含有正弦或余弦信号,图c和d都是个3个周期的正弦信号波,图e和f分别是a、b两图跟c、d两图相乘后的结果,图e所有点的平均值是0.5,说明信号a含有振幅为1的正弦信号c,但图f所有点的平均值是0,则说明信号b不含有信号d。
这个就是通过信号相关性来检测是否含有某个信号的方法。
相应地,我也可以通过把输入信号和每一种频率的正余弦信号进行相乘(关联操作),从而得到原始信号与每种频率的关联程度(即总和大小),这个结果便是我们所要的傅立叶变换结果,下面两个等式便是我们所要的计算方法:
第二个式子中加了个负号,是为了保持复数形式的一致,前面我们知道在计算Im[k]时又加了个负号,所以这只是个形式的问题,并没有实际意义,你也可以把负号去掉,并在计算Im[k]时也不加负号。
这里有一点必须明白一个正交的概念:
两个函数相乘,如果结果中的每个点的总和为0,则可认为这两个函数为正交函数。
要确保关联性算法是正确的,则必须使得跟原始信号相乘的信号的函数形式是正交的,我们知道所有的正弦或余弦函数是正交的,这一点我们可以通过简单的高数知识就可以证明它,所以我们可以通过关联的方法把原始信号分离出正余弦信号。
当然,其它的正交函数也是存在的,如:
方波、三角波等形式的脉冲信号,所以原始信号也可被分解成这些信号,但这只是说可以这样做,却是没有用的。
下面是实域傅立叶变换的BASIC语言代码:
到此为止,我们对傅立叶变换便有了感性的认识了吧。
但要记住,这只是在实域上的离散傅立叶变换,其中虽然也用到了复数的形式,但那只是个替代的形式,并无实际意义,现实中一般使用的是复数形式的离散傅立叶变换,且快速傅立叶变换是根据复数离散傅立叶变换来设计算法的,在后面我们先来复习一下有关复数的内容,然后再在理解实域离散傅立叶变换的基础上来理解复数形式的离散傅立叶变换。
理解离散傅立叶变换(三)
------复数
复数扩展了我们一般所能理解的数的概念,复数包含了实数和虚数两部分,利用复数的形式可以把由两个变量表示的表达式变成由一个变量(复变量)来表达,使得处理起来更加自然和方便,我们知道傅立叶变换的结果是由两部分组成的,使用复数形式可以缩短变换表达式,使得我们可以单独处理一个复变量(这个在后一节的描述中我们就可以更加确切地知道),而且快速傅立叶变换正是基于复数形式的,所以几乎所有描述的傅立叶变换形式都是复数的形式。
但是复数的概念超过了我们日常生活中所能理解的概念,要理解复数是较难的,所以我们在理解复数傅立叶变换之前,先来专门复习一下有关复数的知识,这对后面的理解非常重要。
一、 复数的提出
在此,先让我们看一个物理实验:
把一个球从某点向上抛出,然后根据初速度和时间来计算球所在高度,这个方法可以根据下面的式子计算得出:
其中h表示高度,g表示重力加速度(9.8m/s2),v表示初速度,t表示时间。
现在反过来,假如知道了高度,要求计算到这个高度所需要的时间,这时我们又可以通过下式来计算:
经过计算我们可以知道,当高度是3米时,有两个时间点到达该高度:
球向上运动时的时间是0.38秒,球向下运动时的时间是1.62秒。
但是如果高度等于10时,结果又是什么呢?
根据上面的式子可以发现存在对负数进行开平方运算,我们知道这肯定是不现实的。
第一次使用这个不一般的式子的人是意大利数学家GirolamoCardano(1501-1576),两个世纪后,德国伟大数学家CarlFriedrichGause(1777-1855)提出了复数的概念,为后来的应用铺平了道路,他对复数进行这样表示:
复数由实数(real)和虚数(imaginary)两部分组成,虚数中的根号负1用i来表示(在这里我们用j来表示,因为i在电力学中表示电流的意思)。
我们可以把横坐标表示成实数,纵坐标表示成虚数,则坐标中的每个点的向量就可以用复数来表示,如下图:
上图中的ABC三个向量可以表示成如下的式子:
A=2+6j
B=-4–1.5j
C=3–7j
这样子来表达方便之处在于运用一个符号就能把两个原来难以联系起来的数组合起来了,不方便的是我们要分辨哪个是实数和哪个是虚数,我们一般是用Re()和Im()来表示实数和虚数两部分,如:
ReA=2ImA=6
ReB=-4ImB=-1.5
ReC=3ImC=-7
复数之间也可以进行加减乘除运算:
这里有个特殊的地方是j2等于-1,上面第四个式子的计算方法是把分子和分母同时乘以c–dj,这样就可消去分母中的j了。
复数也符合代数运算中的交换律、结合律、分配律:
AB=BA
(A+B)+C=A+(B+C)
A(B+C)=AB+AC
二、复数的极坐标表示形式
前面提到的是运用直角坐标来表示复数,其实更为普遍应用的是极坐标的表示方法,如下图:
上图中的M即是数量积(magnitude),表示从原点到坐标点的距离,θ是相位角(phaseangle),表示从X轴正方向到某个向量的夹角,下面四个式子是计算方法:
我们还可以通过下面的式子进行极坐标到直角坐标的转换:
a+jb=M(cosθ+jsinθ)
上面这个等式中左边是直角坐标表达式,右边是极坐标表达式。
还有一个更为重要的等式——欧拉等式(欧拉是瑞士的著名数学家,LeonhardEuler,1707-1783):
ejx=cosx+jsinx
这个等式可以从下面的级数变换中得到证明:
上面中右边的两个式子分别是cos(x)和sin(x)的泰勒(Taylor)级数。
这样子我们又可以把复数的表达式表示成指数的形式了:
a+jb=Mejθ(这便是复数的两个表达式)
指数形式是数字信号处理中数学方法的支柱,也许是因为用指数形式进行复数的乘除运算极为简单的缘故吧:
三、复数是数学分析中的一个工具
为什么要使用复数呢?
其实它只是个工具而已,就如钉子和锤子的关系,复数就象那锤子,作为一种使用的工具。
我们把要解决的问题表达成复数的形式(因为有些问题用复数的形式进行运算更加方便),然后对复数进行运算,最后再转换回来得到我们所需要的结果。
有两种方法使用复数,一种是用复数进行简单的替换,如前面所说的向量表达式方法和前一节中我们所讨论的实域DFT,另一种是更高级的方法:
数学等价(mathematicalequivalence),复数形式的傅立叶变换用的便是数学等价的方法,但在这里我们先不讨论这种方法,这里我们先来看一下用复数进行替换中的问题。
用复数进行替换的基本思想是:
把所要分析的物理问题转换成复数的形式,其中只是简单地添加一个复数的符号j,当返回到原来的物理问题时,则只是把符号j去掉就可以了。
有一点要明白的是并不是所有问题都可以用复数来表示,必须看用复数进行分析是否适用,有个例子可以看出用复数来替换原来问题的表达方式明显是谬误的:
假设一箱的苹果是5美元,一箱的桔子是10美元,于是我们把它表示成5+10j,有一个星期你买了6箱苹果和2箱桔子,我们又把它表示成6+2j,最后计算总共花的钱是(5+10j)(6+2j)=10+70j,结果是买苹果花了10美元的,买桔子花了70美元,这样的结果明显是错了,所以复数的形式不适合运用于对这种问题的解决。
四、用复数来表示正余弦函数表达式
对于象Mcos(ωt+φ)和Acos(ωt)+Bsin(ωt)表达式,用复数来表示,可以变得非常简洁,对于直角坐标形式可以按如下形式进行转换:
上式中余弦幅值A经变换生成a,正弦幅值B的相反数经变换生成b:
Aa,B-b,但要注意的是,这不是个等式,只是个替换形式而已。
对于极坐标形式可以按如下形式进行转换:
上式中,MM,θ-φ。
这里虚数部分采用负数的形式主要是为了跟复数傅立叶变换表达式保持一致,对于这种替换的方法来表示正余弦,符号的变换没有什么好处,但替换时总会被改变掉符号以跟更高级的等价变换保持形式上的一致。
在离散信号处理中,运用复数形式来表示正余弦波是个常用的技术,这是因为利用复数进行各种运算得到的结果跟原来的正余弦运算结果是一致的,但是,我们要小心使用复数操作,如加、减、乘、除,有些操作是不能用的,如两个正弦信号相加,采用复数形式进行相加,得到的结果跟替换前的直接相加的结果是一样的,但是如果两个正弦信号相乘,则采用复数形式来相乘结果是不一样的。
幸运的是,我们已严格定义了正余弦复数形式的运算操作条件:
1、 参加运算的所有正余弦的频率必须是一样的;
2、 运算操作必须是线性的,如两个正弦信号可以进行相加减,但不能进行乘除,象信号的放大、衰减、高低通滤波等系统都是线性的,象平方、缩短、取限等则不是线性的。
要记住的是卷积和傅立叶分析也只有线性操作才可以进行。