傅里叶变换的原理及matlab实现

1、屁明理2大穹傅里叶变换的原理及matlab实现课程名称:学 院：专 业:年 级:学生姓名: 指导教师: 日 期：数字图像处理信息工程与自动化学院计算机科学与技术09级11111112012-6-10教务处制 TOC o 1-5 h z HYPERLINK l bookmark17 o Current Document 一、傅立叶变化的原理；3 HYPERLINK l bookmark21 o Current Document （1）原理3 HYPERLINK l bookmark25 o Current Document （2）计算方法3 HYPERLINK l bookmark38 o Cu

2、rrent Document 二、傅立叶变换的应用； 3（1）、频谱分析4（2）、数据压缩4（3）、OFDM4 HYPERLINK l bookmark54 o Current Document 三、傅里叶变换的本质；4 HYPERLINK l bookmark64 o Current Document 四、实验内容；8 HYPERLINK l bookmark71 o Current Document 五、傅立叶变换方法；8 HYPERLINK l bookmark78 o Current Document 六、实验结果及分析；8 HYPERLINK l bookmark82 o Curre

3、nt Document 七、傅立叶变换的意义；9 HYPERLINK l bookmark86 o Current Document （1）、傅立叶变换的物理意义9 HYPERLINK l bookmark89 o Current Document （2）、图像傅立叶变换的物理意义10 HYPERLINK l bookmark92 o Current Document 八、总结；11 HYPERLINK l bookmark96 o Current Document 九、附录；11一、傅立叶变化的原理;（1）原理正交级数的展开是其理论基础！将一个在时域收敛的函数展开成一系列不同 频率谐波的叠加

4、，从而达到解决周期函数问题的目的。在此基础上进行推广，从 而可以对一个非周期函数进行时频变换。从分析的角度看，他是用简单的函数去逼近（或代替）复杂函数，从几何的角度 看，它是以一族正交函数为基向量，将函数空间进行正交分解，相应的系数即为 坐标。从变幻的角度的看，他建立了周期函数与序列之间的对应关系；而从物理 意义上看，他将信号分解为一些列的简谐波的复合，从而建立了频谱理论。当然Fourier积分建立在傅氏积分基础上，一个函数除了要满足狄氏条件外，一 般来说还要在积分域上绝对可积，才有古典意义下的傅氏变换。引入衰减因子 eA（-st），从而有了 Laplace变换。（好像走远了）。（2）计算方法

5、连续傅里叶变换将平方可积的函数f（t）表示成复指数函数的积分或级数 形式。89）=月伸）=I 此8这是将频率域的函数F（3）表示为时间域的函数f （t）的积分形式。连续傅里叶变换的逆变换（inverse Fourier transform）为8川）=尸吓&） = / gww8即将时间域的函数f （t）表示为频率域的函数F（3 ）的积分。一般可称函数f （t）为原函数，而称函数F（3）为傅里叶变换的像函数，原函数 和像函数构成一个傅里叶变换对（transform pair）。二、傅立叶变换的应用;DFT在诸多多领域中有着重要应用，下面仅是颉取的几个例子。需要指出的 是，所有DFT的实际应用都依赖

6、于计算离散傅里叶变换及其逆变换的快速算法， 即快速傅里叶变换（快速傅里叶变换（即FFT）是计算离散傅里叶变换及其逆变 换的快速算法。）。、频谱分析DFT是连续傅里叶变换的近似。因此可以对连续信号x(t)均匀采样并截断以得到有限长的离 散序列，对这一序列作离散傅里叶变换，可以分析连续信号x(t)频谱的性质。前面还提到 DFT应用于频谱分析需要注意的两个问题：即采样可能导致信号混叠和截断信号引起的频 谱泄漏。可以通过选择适当的采样频率(见奈奎斯特频率)消减混叠。选择适当的序列长度 并加窗可以抑制频谱泄漏。、数据压缩由于人类感官的分辨能力存在极限，因此很多有损压缩算法利用这一点将语音、 音频、图像、

7、视频等信号的高频部分除去。高频信号对应于信号的细节，滤除高 频信号可以在人类感官可以接受的范围内获得很高的压缩比。这一去除高频分量 的处理就是通过离散傅里叶变换完成的。将时域或空域的信号转换到频域，仅储 存或传输较低频率上的系数，在解压缩端采用逆变换即可重建信号。、 OFDMOFDM (正交频分复用)在宽带无线通信中有重要的应用。这种技术将带宽为N 个等间隔的子载波，可以证明这些子载波相互正交。尤其重要的是，OFDM调制 可以由IDFT实现，而解调可以由DFT实现。OFDM还利用DFT的移位性质，在每 个帧头部加上循环前缀(Cyclic Prefix)，使得只要信道延时小于循环前缀的长 度，就

8、能消除信道延时对传输的影响。三、傅里叶变换的本质；傅里叶变换的公式为F(o) = j f (t)e-网出-s可以把傅里叶变换也成另外一种形式：F(O) = -1 Y f (t), ejo A2兀可以看出，傅里叶变换的本质是内积，三角函数是完备的正交函数集，不同频率的三角 函数的之间的内积为0，只有频率相等的三角函数做内积时，才不为0。Y e jQ1t,e jQ2t A =j ej(r-q2)tdt = 2兀8 (Q Q )下面从公式解释下傅里叶变换的意义因为傅里叶变换的本质是内积，所以f(t)和e加求内积的时候，只有f(t)中频率为的 分量才会有内积的结果，其余分量的内积为0。可以理解为f(t

9、)在e加上的投影，积分值是 时间从负无穷到正无穷的积分，就是把信号每个时间在的分量叠加起来，可以理解为f(t) 在ejt上的投影的叠加，叠加的结果就是频率为的分量，也就形成了频谱。傅里叶逆变换的公式为1专8f (t) =j F(沁仙2兀-8下面从公式分析下傅里叶逆变换的意义傅里叶逆变换就是傅里叶变换的逆过程，在F(o )和ej求内积的时候，F(o )只有t 时刻的分量内积才会有结果，其余时间分量内积结果为0，同样积分值是频率从负无穷到正 无穷的积分，就是把信号在每个频率在t时刻上的分量叠加起来，叠加的结果就是f(t)在t 时刻的值，这就回到了我们观察信号最初的时域。对一个信号做傅里叶变换，然后

10、直接做逆变换，这样做是没有意义的，在傅里叶变换和 傅里叶逆变换之间有一个滤波的过程。将不要的频率分量给滤除掉，然后再做逆变换，就得 到了想要的信号。比如信号中掺杂着噪声信号，可以通过滤波器将噪声信号的频率给去除， 再做傅里叶逆变换，就得到了没有噪声的信号。优点：频率的定位很好，通过对信号的频率分辨率很好，可以清晰的得到信号所包含的 频率成分，也就是频谱。缺点：因为频谱是时间从负无穷到正无穷的叠加，所以，知道某一频率，不能判断，该 频率的时间定位。不能判断某一时间段的频率成分。例子：平稳信号：x(t)=cos(2*pi*5*t)+cos(2*pi*10*t)+cos(2*pi*20*t)+cos

11、(2*pi*50*t)傅里叶变换的结果:由于信号是平稳信号，每处的频率都相等，所以看不到傅里叶变换的缺点。对于非平稳信号：信号是余弦信号，仍然有四个频率分量傅里叶变换的结果:由上图看出知道某一频率，不能判断，该频率的时间定位。不能判断某一时间段的频率成分。四、实验内容;（一）用快速傅立叶变换FFT实现数字图像的傅立叶变换，进一步加深对DFT算法原理 和基本性质的理解（因为FFT只是DF T的一种快速算法，所以FF T的算法结果必然满足DFT 的基本性质）。（二）熟悉FFT算法原理和FFT子程序的应用。（三）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分布 误差及其原因，

12、以便在实际中正确应用FFT。五、傅立叶变换方法；（1）、矩阵形式的傅立叶变换的算法如下:数字图像f的傅立叶正变换：r = A fAt c R数字图像f的傅立叶反变换：f= Atr AC R变换矩阵Ac = Ar = 0 2 3WWW0 12WWWo o oWWW一.WW 2（ N 一 1）W 0Wn-1W2（N-1）:w（ N-1）2一 皿，其中w =e -iN，N为图像的维数。六、实验结果及分析；（一）对原图像进行傅立叶变换，实验结果如图5-1：图5-1分析：图像显示了原图像及其傅立叶频谱。观察傅立叶谱中心对称，在此图像进行傅立叶变换的计算之前被乘以（-1*+，以此增强了灰度级细节。（二）输

13、出彩色图像greens.jpg的傅立叶频谱，实验结果如图5-2：,，己田整叶障立叶地:-图5-2分析：图像显示了原图像和其彩色图像傅立叶频谱。可以看出图像的频率分布是以原点 为圆心，对称分布的。变换之后的图像在原点平移之前四角是低频，最亮，平移之后中间部 分是低频，最亮，亮度大说明低频的能量大（幅角比较大）（三）对彩色图像football.jpg进行二维DCT变换，实验结果如图5-3:图5-3分析：二维DCT变换后的频谱图亮点在左上角。七、傅立叶变换的意义;（1）、傅立叶变换的物理意义傅立叶变换是数字信号处理领域一种很重要的算法。要知道傅立叶变换算 法的意义，首先要了解傅立叶原理的意义。傅立叶

14、原理表明：任何连续测量的时 序或信号，都可以表示为不同频率的正弦波信号的无限叠加。而根据该原理创立 的傅立叶变换算法利用直接测量到的原始信号，以累加方式来计算该信号中不同 正弦波信号的频率、振幅和相位。和傅立叶变换算法对应的是反傅立叶变换算法。 该反变换从本质上说也是一种累加处理，这样就可以将单独改变的正弦波信号转 换成一个信号。因此，可以说，傅立叶变换将原来难以处理的时域信号转换成了 易于分析的频域信号（信号的频谱），可以利用一些工具对这些频域信号进行处 理、加工。最后还可以利用傅立叶反变换将这些频域信号转换成时域信号。从现代数学的眼光来看，傅里叶变换是一种特殊的积分变换。它能将满足一定条

15、件的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域，傅里 叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。在数学 领域，尽管最初傅立叶分析是作为热过程的解析分析的工具，但是其思想方法仍 然具有典型的还原论和分析主义的特征。任意的函数通过一定的分解，都能够 表示为正弦函数的线性组合的形式，而正弦函数在物理上是被充分研究而相对简 单的函数类：1.傅立叶变换是线性算子，若赋予适当的范数，它还是酉算子；2.傅 立叶变换的逆变换容易求出，而且形式与正变换非常类似；3.正弦基函数是微分 运算的本征函数,从而使得线性微分方程的求解可以转化为常系数的代数方程的 求解.在线性时不变

16、杂的卷积运算为简单的乘积运算，从而提供了计算卷积的一种 简单手段；4.离散形式的傅立叶的物理系统内,频率是个不变的性质，从而系统对 于复杂激励的响应可以通过组合其对不同频率正弦信号的响应来获取；5.著名 的卷积定理指出:傅立叶变换可以化复变换可以利用数字计算机快速的算出（其 算法称为快速傅立叶变换算法（FFT）。正是由于上述的良好性质，傅里叶变换在物 理学、数论、组合数学、信号处理、概率、统计、密码学、声学、光学等领域都 有着广泛的应用。（2）、图像傅立叶变换的物理意义图像的频率是表征图像中灰度变化剧烈程度的指标，是灰度在平面空间上的梯 度。如：大面积的沙漠在图像中是一片灰度变化缓慢的区域，对

17、应的频率值很低； 而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域，对应 的频率值较高。傅立叶变换在实际中有非常明显的物理意义，设f是一个能量有 限的模拟信号，则其傅立叶变换就表示f的谱。从纯粹的数学意义上看，傅立叶 变换是将一个函数转换为一系列周期函数来处理的。从物理效果看，傅立叶变换 是将图像从空间域转换到频率域，其逆变换是将图像从频率域转换到空间域。换 句话说，傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布 函数，傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数。傅立叶变换以前，图像（未压缩的位图）是由对在连续空间（现实空间）上的采样得到一系 列点的集

18、合，我们习惯用一个二维矩阵表示空间上各点，则图像可由z=f（x,殊表示。由于 空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就由梯度来表示，这 样我们可以通过观察图像得知物体在三维空间中的对应关系。为什么要提梯度？因为实际上 对图像进行二维傅立叶变换得到频谱图，就是图像梯度的分布图，当然频谱图上的各点与图 像上各点并不存在一一对应的关系，即使在不移频的情况下也是没有。傅立叶频谱图上我们 看到的明暗不一的亮点，实际上图像上某一点与邻域点差异的强弱，即梯度的大小，也即该 点的频率的大小（可以这么理解，图像中的低频部分指低梯度的点，高频部分相反）。一般 来讲，梯度大则该点的亮度强，否

19、则该点亮度弱。这样通过观察傅立叶变换后的频谱图，也 叫功率图，我们首先就可以看出，图像的能量分布，如果频谱图中暗的点数更多，那么实际 图像是比较柔和的（因为各点与邻域差异都不大，梯度相对较小），反之，如果频谱图中亮 的点数多，那么实际图像一定是尖锐的，边界分明且边界两边像素差异较大的。对频谱移频 到原点以后，可以看出图像的频率分布是以原点为圆心，对称分布的。将频谱移频到圆心除 了可以清晰地看出图像频率分布以外，还有一个好处，它可以分离出有周期性规律的干扰信 号，比如正弦干扰，一副带有正弦干扰，移频到原点的频谱图上可以看出除了中心以外还存 在以某一点为中心，对称分布的亮点集合，这个集合就是干扰噪音产生的，这时可以很直观 的通过在该位置放置带阻滤波器消除干扰。八、总结；图像经过二维傅立叶变换后，其变换系数矩阵表明：若变换矩阵Fn原点 设在中心，其频谱能量集中分布在变换系数短阵的中心附近(图中阴影区)。若所 用的二维傅立叶变换矩阵Fn的原点设在左上角，那么图像信号能量将集中