信号变换.docxVIP

第四章 信号的变换 4.3 快速傅里叶变换 傅里叶变换能将满足一定条件的某个函数表示成三角函数（正弦和/或余弦函数）或者它们的积分的线性组合。在不同的研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。最初傅里叶分析是作为热过程的解析分析的工具被提出的。但随着傅里叶变换的丰富和发展，极大地促进了信息科学的丰富和发展。现代的信息科学和技术离不开傅里叶变换的理论和方法。 4.3.1傅里叶级数和傅里叶变换 描述周期现象最简单的周期函数是物理学上所说的谐波函数，他是由正弦或者余弦函数来表示 (4-1) 利用三角公式，上式可以写成 (4-2) 由于φ是常数，令a=Acosφ，b=-Asinφ，那么可以得到 (4-3) 其中公式中 (4-4) 从这里可以看出：一个带初相位的余弦函数可以看成是有一个不带初相位的正弦函数与一个不带初相位的余弦函数的合成。 谐波函数是周期函数中最简单的函数，它描述的也是最简单的周期现象，在实际中所碰到的周期现象往往要比它复杂很多。但这些复杂的函数都可以近似分解成不同频率的正弦函数和余弦函数。下面就介绍一种复杂的函数分解为一系列不同频率的正弦函数和余弦函数的方法。 法国数学家傅里叶发现，任何周期函数都可以用正弦函数和余弦函数构成的无穷级数来表示（选择正弦函数与余弦函数作为基函数是因为它们是正交的），后世称为傅里叶级数（法文：série de Fourier，或译为傅里叶级数）一种特殊的三角级数。 在高等代数中有这样一个问题，怎么将一个周期为2l函数分解成傅里叶级数，给出的解答式为 (4-5) 其中 (4-6) 如果f(x)是奇函数，积分上下限相互对称，则此时f(x)cosnπxl项成为奇函数，可以知道an均为零，得到的傅里叶正弦级数 (4-7) 上式中bn的积分可以简写为 (4-8) 如果f(x)是偶函数，同样因为积分上下限相互对称，这时f(x)sinnπxl为奇函数，故bn均 为零，得到的傅里叶级数是余弦级数 (4-9) 式中，an可以简写为 (4-10) 现在我们将上述公式应用于离散的傅里叶级数中。 （1）在信号处理中，我们遇到的常常不是一个函数，而是一个离散的的数列，举一个例子，等间隔时间取样的时间序列{x0，x1，x2，…， xN-1}，在这里数据的个数是N，一般取N为偶数，如果取2的对数对应的偶数能够加快计算速度。下面对取值范围进行改造，首先，我们得到的数字信号只能在正的时间段取值，在负的时间段不能取值。但由于我们取的是无限长的周期序列，周期为2l，因此，我们把取值范围(-l，l)修改为（0,2l），这样就可以避免在负的时间段取值。 （2）由于处理的是离散的数据序列，因此不能再用积分，而应用积分的离散形式，用求和来表示，即 (4-11) (3)我们在（0,2l）里等间隔取了N个取值点，取样时间间隔为dx→?t，其中l=N?t2。 有了上述的改正，我们可以得到 (4-12) (4-13) 所以式(4-5)的离散形式为 (4-14) 式中 (4-15) (4-16) (4-17) 在实际数据处理中，k一般取N2，此时波的周期最小，获得的频率范围最大，所以想要获得高频率的信号，就需要缩短取样间隔。 4.3.1快速傅里叶变换的MATLAB实现 有限长序列可以通过 (DFT）将其频域也离散化成有限长序列。但其计算量太大，很难实时地处理问题，离散傅里叶变换因此引出了快速傅里叶变换(FFT). 1965年，Cooley和Tukey提出了计算离散傅里叶变换（DFT）的快速算法，将DFT的运算量减少了几个数量级。从此，对快速傅里叶变换（FFT）算法的研究便不断深入，数字信号处理这门新兴学科也随FFT的出现和发展而迅速发展。根据对序列分解与选取方法的不同而产生了FFT的多种算法，基本算法是基2DIT和基2DIF。FFT