傅里叶变换及应用

《傅里叶变换及应用》由会员分享，可在线阅读，更多相关《傅里叶变换及应用（18页珍藏版）》请在金锄头文库上搜索。

1、傅里叶变换在MATLZB里的应用摘要：在现代数学中，傅里叶变换是一种非常重要的变换，且在数字信号处理中有 着广泛的应用。本文首先介绍了傅里叶变换的基本概念、性质及发展情况；其次，详细 介绍了分离变数法及积分变换法在解数学物理方程中的应用。傅立叶变换将原来难以处 理的时域信号转换成了易于分析的频域信号，再利用傅立叶反变换将这些频域信号转换 成时域信号。应用MATLAB实现信号的谱分析和对信号消噪。关键词：傅里叶变换；MATLAB软件；信号消噪Abstract: In modern mathematics,Fourier transform is a transform is very impor

2、tant ,And has been widely used in digital signal paper first introduces the basic concepts, properties and development situation of Fourier transform ;Secondly, introduces in detail the method of separation of variables and integral transform method in solving equations in Mathematical transformatio

3、n makes the original time domain signal whose analysis is difficult easy, by transforming it into frequency domain signal that can be transformed into time domain signal by inverse transformation of Fourier. Using Mat lab realizes signal spectral analysis and signal denoising.Key word: Fourier trans

4、formation, software of mat lab ,signal denoising1、傅里叶变换的提出及发展 在自然科学和工程技术中为了把较复杂的运算转化为较简单的运算，人们常常采用 所谓变换的方法来达到目的例如在初等数学中，数量的乘积和商可以通过对数变换化 为较简单的加法和减法运算。在工程数学里积分变换能够将分析运算（如微分，积分） 转化为代数运算，正是积分变换这一特性，使得它在微分方程和其它方程的求解中成为 重要方法之一。1804年，法国科学家傅里叶由于当时工业上处理金属的需要，开始从事热流动的 研究他在题为热的解析理论一文中，发展了热流动方程，并且指出如何求解在求 解过程中

5、，他提出了任意周期函数都可以用三角级数来表示的想法。他的这种思想，虽 然缺乏严格的论证，但对近代数学以及物理、工程技术却都产生了深远的影响，成为傅 里叶变换的起源。从现代数学的眼光来看，傅里叶变换是一种特殊的积分变换。它能将满足一定条件 的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域，傅里叶变换具 有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。1傅里叶变换通过对函数的分析来达到对复杂函数的深入理解和研究。最初，傅立叶 分析是作为热过程的解析分析的工具，但是其思想方法仍然具有典型的还原论和分析主 义的特征。“任意”的函数通过一定的分解，都能够表示为正弦函数的线性组合的形

6、式， 而正弦函数在物理上是被充分研究而相对简单的函数类。利用这一点，傅里叶变换可通 过对相对简单的事物的研究来了解复杂事物，而且现代数学发现傅里叶变换具有非常好 的性质：（1）傅里叶变换是线性算子，若赋予适当的范数+它还是酉算子；（2）傅里叶变换的逆变换容易求出，而且形式与正变换非常类似；（3）正弦基函数是微分运算的本征函数，从而使得线性微分方程的求解可以转化 为常系数的代数方程的求解在线性时不变的物理系统内，频率是个不变的性质，从而 系统对于复杂激励的响应可以通过组合其对不同频率正弦信号的响应来获取；（4）著名的卷积定理指出.傅里叶变换可以化复杂的卷积运算为简单的乘积运算， 从而提供了计算卷

7、积的一种简单手段；（5）离散形式的傅里叶变换可以利用数字计算机快速的算出（其算法称为快速傅 里叶变换算法）。（6）正是由于上述的良好性质,傅里叶变换在物理学、数论、组合数学、信号处理、 概率、统计、密码学、声学、光学等领域都有着广泛的应用。2、傅里叶变换的基本概念由傅里叶级数知，一个周期函数可以展开成为傅里叶级数，而一个非周期函数可以 看成某个周期函数其周期趋向于无穷大转化而来。根据这个思路，我们可以得到傅里叶 积分公式及傅里叶积分公式成立的充分条件傅里叶积分定理。傅里叶级数的指数形式/ / 一, 一定理设fT是以T 0 T S为周期的实函数2且在I 2 2丿上满足狄利 f ()克雷条件，即f

8、TV在一个周期上满足：(1)连续或只有有限个第一类间断点；(2)只有有限个极值点.则在连续点处，有f ()= a +(a cos nt + b sin nt)(1)T 2nnn=1a =12 f (t )dt0 T T T 其中2,a = j 2 f (t)cos nw tdt(n = 1,2,) n T T T2 ,b = j 2 f Osin nw tdt(n = 1,2.) n T T T2,f (0+0)+f (0- )在间断点点处，(1)式右端级数收敛于2ez + ei.ez e询2i,.于是cos 0 =sin 0 =又4)=+兰T2、 einwt + einwteinwt ein

9、wta+ bn 2 nn=12iag=-2 + n =1a + ibn einwt + n e inwt2 丿aa 一 ibC = , C = nn令 02 n2a + ibC = nnn2n 1,2,3, , 则4)= T=c + 9 eiwt + c ei2wt + 0 1 2C einwt +ne_ iwt + C e_i 2wt + C einwt +)1 2 n (2)c einwt nn= g(2)式称为傅里叶级数的复指数形式，具有明显的物理意义.容易证明Cn可以合写成一个式子，即c = J 2 f C-indt(n = 0,l,2,) n T - T T-2(3)傅里叶积分任何一

10、个非周期函数f都可看成是由某个周期函数讣、当T+8时转化而来的.即iim fTf (t)由公式(2) 、(3)得Tn=-ge-inOT dTT T2einOt可知limT T+wn=-we-inOT dTeinOt=no,Ao = o -on nn -1o =Ao于是f ( )= lim Jrrin=-wT T +w=limAo t 0 2兀nTnf(t)=艺 J 2 f G)e - ionT dT艺 J： f d-io”T dTn=-weio tneio t AonJ 2 f G )?-io”T dT eiont2兀-T Tlim (o koAonT0 =-w T n n n n = -w(

11、4)注意到当Ao n0,即T 时,e (o ) T卩+w f G2 - io T dT eiontT nn 2兀-w从而按照积分的定义，（4）可以写为：或者f c)=2-卜卜 f G St di eitd(5)2兀 s g傅氏积分公式.公式(5)称为函数定理 若f)在(-8, +8)上满足条件：(1) f l在任一有限区间上满足狄氏条件；(2) f O l在无限区间(-8, +8)上绝对可 卜廿伽f Qf Q t积，即一8收敛，则(5)在的连续点成里；而在的间断点10处f (t + 0)+ f (t - 0)0 0应以2来代替.上述定理称为傅氏积分定理.可以证明，当/)满足傅氏积分定理条件时，

12、公式(5)可以写为三角形式，即f()|fQ 在f(t)连续点处，-J+gJ+gfG)COSSt T)dTd = b(t + 0)+ f(t 0)甘宀兀0 -gI?其匕.I 2周期傅里叶变换描述周期现象的最简单的周期函数是物理学上所说的谐波函数，它由正弦或余弦函数来表示y (t) = A cos( wt + a)而所有函数都可以看做是不同频率的正弦或余弦函数的叠加。下面介绍周期函数的傅里叶变换3。将一个周期为T的函数分解为Fourier级数，其三角形式展开为:f (t) = a +(a cos nt + b sin nt)0nn离散傅里叶变换 但我们在数字资料处理中经常的不是一个函数，而是一个离

13、散的序列。与连续时间 信号的分析类似，对于连续时间信号进行离散Fourier变换，一般可概括为时域采样, 时域截断，频域采样三个步骤，最终导出离散傅立叶变换2对为：X (n)二吕三X （艸n=0n=0,l,2,N-1它通过连续傅立叶变换，将N个时域采样点与N个频域采样点联系起来。3、傅立叶变换的应用傅立叶变换以前，图像（未压缩的位图）是由对在连续空间（现实空间）上的采样 得到一系列点的集合，我们习惯用一个二维矩阵表示空间上各点，则图像可由z=f（x,y） 来表示。由于空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就 由梯度来表示，这样我们可以通过观察图像得知物体在三维空间中的对

14、应关系。为什么 要提梯度因为实际上对图像进行二维傅立叶变换得到频谱图，就是图像梯度的分布图， 当然频谱图上的各点与图像上各点并不存在一一对应的关系，即使在不移频的情况下也 是没有。傅立叶频谱图上我们看到的明暗不一的亮点，实际上图像上某一点与邻域点差 异的强弱，即梯度的大小，也即该点的频率的大小（可以这么理解，图像中的低频部分 指低梯度的点，高频部分相反）。一般来讲，梯度大则该点的亮度强，否则该点亮度弱 这样通过观察傅立叶变换后的频谱图，也叫功率图，我们首先就可以看出，图像的能量 分布，如果频谱图中暗的点数更多，那么实际图像是比较柔和的（因为各点与邻域差异 都不大，梯度相对较小），反之，如果频谱图中亮的点数多，那么实际图像一定是尖锐 的，边界分明且边界两边像素差异较大的。对频谱移频到原点以后，可以看出图像的频 率分布是以原点为圆心，对称分布的。将频谱移频到圆心除了可以清晰地看出图像频率 分布以外，还有一个好处，它可以分离出有周期性规律的干扰信号，比如正弦干扰，一 副带有正弦干扰，移频到原点的频谱图上可以看出除了中心以外还存在以某一点为中 心，对称分布的亮点集合，这个集合就是干扰噪音产生的，这时可以很直观的通过在该 位置放置带阻滤波器消除干扰。冲激信号4冲激函数是最基本的函数，其傅里叶变换是系统函数，只要知道系统函数，那么通过这个系统的输出函数并可以确定。在Matlab中产生冲激函