DFT在信号频谱分析中的应用

1、设计一 DFT在信号频谱分析中的应用一、设计题目 DFT在信号频谱分析中的应用二、设计目的1. 熟悉DFT的性质。2. 加深理解信号频谱的概念及性质。 3. 了解高密度谱与高分辨率频谱的区别。三、设计原理所谓信号的频谱分析就是计算信号的傅里叶变换。连续信号与系统的傅里叶分析显然不便于直接用计算机进行计算，使其应用受到限制，而DFT是一种时域和频域均离散化的变换，适合数值运算，成为分析离散信号和系统的有力工具。 工程实际中，经常遇到的连续信号Xa(t),其频谱函数Xa(jW)也是连续函数。数字计算机难于处理，因而我们采用DFT来对连续时间信号的傅里叶变换进行逼近，进而分析连续时间信号的频谱。四、

2、实现方法离散傅里叶变换是有限长序列的傅里叶变换，它相当于把信号的傅里叶变换进行等频率间隔采样，并且有限长序列的离散傅里叶变换和周期序列的离散傅里叶级数本质是一样的。快速傅里叶变换（FFT）并不是一种新的变换，它是离散傅里叶变换的一种快速算法，并且主要是基于这样的思路而发展起来的：（1）把长度为N的序列的DFT逐次分解成长度较短的序列的DFT来计算。（2）利用WN(nk)的周期性和对称性，在DFT运算中适当的分类，以提高运算速度。（对称性，;周期性，r为任意整数） 离散傅里叶变换的推导：离散傅里叶级数定义为 （1-1）将上式两端乘以并对n在0N-1求和可得 因为 所以 这样用k代替m得（1-2）

3、令则（1-2）成为DFS （1-3）（1-1）成为IDFS （1-4）式（1-3）、（1-4）式构成周期序列傅里叶级数变换关系。其中都是周期为N的周期序列，DFS·表示离散傅里叶级数正变换，IDFS·表示离散傅里叶级数反变换。习惯上，对于长为N的周期序列，把0nN-1区间称为主值区，把称为的主值序列，同样也称为的主值序列。由于，对于周期序列仅有N个独立样值，对于任何一个周期进行研究就可以得到它的全部信息。在主值区研究与是等价的，因此在主值区计算DFS和DFT是相等的，所以DFT计算公式形式与DFS基本相同。其关系为 所以离散傅里叶正变换 0kN-1离散傅里叶变换（DFT）定

4、义：设有限长序列x (n) 长为N（0nN-1），其离散傅里叶变换是一个长为N的频率有限长序列（0kN-1），其正变换为 0kN-1 （）离散傅里叶变换的实质是：把有限长序列当做周期序列的主值序列进行DFS变换，x(n)、X(k)的长度均为N，都是N个独立值，因此二者具有的信息量是相等的。已知x(n)可以唯一确定X(k),已知X(k)可以唯一确定x(n)。虽然离散傅里叶变换是两个有限长序列之间的变化，但它们是利用DFS关系推导出来的，因而隐含着周期性。构造离散傅里叶变换的Matlab实现程序如下： functionXk=dft(xn,N) n=0:1:N-1; k=n; WN=exp(-j*2

5、*pi/N); nk=n'*k; WNnk=WN.nk; Xk=xn*WNnk快速傅里叶变换（FFT）并不是与DFT不同的另外一种变换，而是为了减少DFT计算次数的一种快速有效的算法 共轭对称性：设有限长序列的长度为N，以N为周期的周期延拓列为 周期序列的共轭对称分量和共轭反对称分量分别为 (1-5) (1-6)同样可以证明，它们满足 （1-7） （1-8） 则有限长序列的圆周共轭对称分量和圆周共轭反对称分量分别定义为： （1-9） （1-10）由于满足 故 （1-11）显然，长度为N的有限长序列可以分解为圆周共轭对称分量和圆周共轭反对称分量之和，和的长度皆为N。利用有限长序列与周期序

6、列的共轭对称分量和反对称分量的关系式（1-9）和式（1-10），以及式（1-11）可以推导出DFT的一系列的对称性质（1） DFT 式中表示的共轭复序列。证明：DFT 又因为 所以DFT（2） 复序列实部的DFT等于DFT的圆周共轭对称部分，即 DFT证明：DFTDFT=DFT+DFT=利用DFT的对称性可求得的DFT:设 则DFT因为 所以DFTDFT=五、设计内容1. 用MATLAB语言编写计算序列x(n)的N点DFT的m函数文件dft.m。并与MATLAB中的内部函数文件fft.m作比较。程序如下：function Xk=dft(xn,N)if length(xn)<N xn=xn

7、,zeros(1,N-length(xn);endn=0:N-1;for k=0:N-1 Xk(1,k+1)=sum(xn.*exp(-1)*j*n*k*(2*pi/N);End运算量估计：对于N=点序列进行时间抽选奇偶分解FFT计算，需分M级，每级计算N/2个蝶。每一级需N/2次复乘、N次复加，因此总共需要进行：复乘： 复加：直接计算N点的DFT，需要次复乘、N(N-1)次复加。N值越大，时间抽选奇偶分解FFT算法越优越。例如当N=2048点时，时间抽选奇偶分解FFT算法比直接计算DFT速度快300多倍可以用一下Matlab程序比较DFT和FFT的运算时间N=2048;M=11;x=1:M,

8、zeros(1,N-M);t=cputime;y1=fft(x,N);Time_fft=cputime-tt1=cputime;y2=dft(x,N);Time_dft=cputime-t1t2=cputime;运行结果：Time_fft = 0.0469Time_dft = 15.2031由此可见FFT算法比直接计算DFT速度快得多2. 对离散确定信号 作如下谱分析：(1) 截取使成为有限长序列N()，写程序计算出的N点DFT ,画出时域序列图xnn和相应的幅频图。程序如下： （假设N取15，即0n14 时, 编写程序，计算出X(n)的15点DFT Xk）n = 0:14;xn=cos(0.

9、48*pi*n)+cos(0.52*pi*n);Xk = fft (xn, 15);subplot(2,1,1); stem(n, xn); grid;subplot(2,1,2); stem(n, abs(Xk); grid; (2) 将 (1)中补零加长至M点，长度M自己选,（为了比较补零长短的影响，M可以取两次值，一次取较小的整数，一次取较大的整数），编写程序计算的M点DFT, 画出时域序列图和两次补零后相应的DFT幅频图。程序如下： （假设M取20和M取65，即分别补5个0和50个0，得补零后20点的序列xn1和65点的序列xn2，编写程序，计算出xn1的20点DFT Xk1和 xn2

10、的65点DFT Xk2）n = 0:14; xn=cos(0.48*pi*n)+cos(0.52*pi*n);n1 = 0:19; xn1 = xn, zeros(1,5);n2= 0:64; xn2 = xn, zeros(1,50);Xk1 = fft(xn1, 20);Xk2 = fft(xn2, 65);subplot(3,1,1); stem(n, xn); grid;subplot(3,1,2); stem(n1, abs(Xk1); grid;subplot(3,1,3); stem(n2, abs(Xk2); grid;(2) 用补零DFT计算 (1)中N点有限长序列频谱并画出

11、相应的幅频图。程序如下：（假设M取150）n = 0:14; xn=cos(0.48*pi*n)+cos(0.52*pi*n);n1 = 0:149; xn3 = xn, zeros(1,135);Xk3 = fft(xn3, 150);plot(n1, abs(Xk3); grid;3. 研究高密度谱与高分辨率频谱。对连续确定信号以采样频率fs=32kHz对信号采样得离散信号，分析下列三种情况的幅频特性。(1)采集数据长度取N=16点，编写程序计算出的16点DFT,并画出相应的幅频图。 (2) 采集数据长度N=16点，补零加长至45点，利用补零DFT计算 的频谱并画出相应的幅频图。(3) 采

12、集数据长度取为45点，编写程序计算出M点采集数据的的频谱并画出相应的幅频图。程序如下：T=1/(32*103);t=(0:15);xn=cos(2*pi*6.5*103*t*T)+cos(2*pi*7*103*t*T)+cos(2*pi*9*103*t*T);Xk=fft(xn,16);subplot(2,1,1);stem(t,xn);grid;subplot(2,1,2);stem(t,abs(Xk);grid;T=1/(32*103);t=(0:15);xn=cos(2*pi*6.5*103*t*T)+cos(2*pi*7*103*t*T)+cos(2*pi*9*103*t*T);n1=

13、0:45; xn1=xn,zeros(1,30);Xk1=fft(xn1,46);subplot(2,1,1);stem(n1,xn1);grid;subplot(2,1,2);plot(n1,abs(Xk1);grid;T=1/(32*103);t=0:45;xn=cos(2*pi*6.5*103*t*T)+cos(2*pi*7*103*t*T)+cos(2*pi*9*103*t*T);Xk2=fft(xn,46);subplot(2,1,1);stem(t,xn);grid;subplot(2,1,2);plot(t,abs(Xk2);grid;六、设计报告要求1. 简述设计目的及原理。2. 完成设计中