数字信号处理实验实验五

1、数字信号处理实验报告实验名称：应用FFT实现信号频谱分析 学生姓名： z 学生学号： 学生班级： 上课时间： 周二上午 指导老师： 一、 实验目的（1） 能够熟练掌握快速离散傅里叶变换的原理及应用FFT进行频谱分析的基本方法。（2） 对离散傅里叶变换的主要性质及FFT在数字信号处理中的重要作用有进一步的了解。二、 实验原理1、 离散傅里叶变换（DFT）及其主要性质DFT表示离散信号的离散频谱，DFT的主要性质中有奇偶对称性、虚实特性等。通过实验可以加深理解。 实序列的DFT具有偶对称的实部和奇对称的虚部，这可以证明如下：由定义，可得X(k)= =X(N-k)= =所以，X(k)=X*(N-k)

2、实序列DFT的这个特性，在本实验中可以通过实指数序列及三角序列看出来。对于单一频率的三角序列来说，它的DFT谱线也是单一的，这个物理意义可以从实验中得到验证，在理论上可以推导如下：设 其DFT为 X(k)= = =从而，X(0)= X(1)= X(2)=0 X(N-2)=0 X(N-1)=以上这串式中X(0)反映了x(n)的直流分量，X(1)是x(n)的一次谐波，又根据虚实特性X*(N-1)=X(1),而其他分量均为零。当周期减少时显然RN(n)的谱只应该在k=3及k=N-3才有分量，实验者可以通过和上述相同的步骤加以理论证明。由于与相位差，所以它的DFT只包括实部而没有虚部，以上这些性质可在

3、本实验中得到验证。2、利用DFT对信号进行频谱分析DFT的重要应用之一是对时域连续信号的频谱进行分析，称为傅里叶分析，时域连续信号离散傅里叶分析的基本步骤如图5-1所示。 图5-1 时域连续信号离散傅里叶分析的处理步骤其中消混叠低通滤波器LPF（预滤波器）的引入，是为了消除或减少时域连续信号转换成序列时可能出现的频谱混叠的影响。实际工作中，时域离散信号x(n)的时宽是很长的甚至是无限长的（例如语音或音乐信号）。由于DFT需要，必须把x(n)限制在一定时间间隔之内，即进行数据截断。数据的截断相当于加窗处理。因此，在计算机x(n)de DFT之前，用一个时域有限的窗函数(n)加到x(n)上是非常必

4、要的。Xc(t)通过变换器转换成取样序列x(n)时域离散序列。其频谱用表示，它的频率的周期函数，即其中Xc(j)或为xc(t)的频谱。在实际应用中，消混叠低通滤波器的阻带不可能是无限衰减的，故由Xc(j)周期延拓得到的有非零重叠，即出现混叠现象。由于进行DFT的需要，必须对序列x(n)进行加窗处理，即v(n)=x(n)(n)。加窗对频域的影响，用周期卷积表示。 最后是进行DFT运算，加窗后的DFT为 其中假设窗函数长L小于或等于DFT长度N。有限长序列v(n)=x(n)(n)的DFT相当于v(n)傅里叶变换的等间隔取样。V(k)便是sc(t)的离散频率函数。因为DFT频率间隔为，且模拟频率和数

5、字频率间的关系=T，所以离散频率点对应的模拟频率为 显然频率分辨率f为 f=利用DFT计算频谱，只会给频谱或的频率分量，即频率的取样值，而不可能得到连续的频谱函数。如果在两个离散的谱线之间有一个特别大的频谱分量，就无法检测出来了。为了保持原来频谱形状不变的情况下，使谱线加密，即使频域取样点数增加，从而使原来看不到的频谱分量变得可以看到，可以通过在信号数据的末端补加一些零值点，使DFT计算周期内点数增加，但又不改变原有的记录数据的方法实现。3、快速离散傅里叶变换（FFT）快速离散傅里叶变换是计算离散傅里叶变换的一种快速算法，为了提高运算速度，FFT将DFT的计算逐次分解成较小点数的DFT。按时间

6、抽取FFT算法把输入序列x(n)按其n值为偶数或是奇数分解成越来越短的序列。按频域抽取FFT算法是把输出序列X(k)按其k值是偶数或是奇数来分解成越来越短的序列。三、实验内容（1）实验前学生应认真学习数字信号处理中有关章节的内容，掌握DFT的基本理论和应用FFT计算信号频谱的原理与方法；（2）编写一个调用FFT函数的通用程序，可计算下列三种序列的离散频谱；指数序列：v1(n)=(0.9)nu(n);周期为N的余弦序列：v2(n)=cos(n)，且0 n N-1；符合函数序列：v3(n)=0.9sin(n)+0.6sin(n)（3） 计算实指数序列v1(n)的N点离散频谱V1（k），记录N为不同

7、 的2的幂次方时的V1（k）值，并与理论V1(ejk)进行比较。（4） 计算周期为N的余弦序列v2(n)的N点FFT、2N点FFT及（N+2）点FFT，记录结果并分析说明。（5） 已知信号x(t)=0.15sin(2f1t)+sin(2f2t)-0.1sin(2f3t)，其中f1=1Hz，f2=2Hz，f3=3Hz，取样频率为32Hz。编程实现： 32点FFT，画出其幅度谱。 64点FFT，画出其幅度谱，比较两者间的差异，思考实际频率与离散频谱图中横坐标k的对应关系。1,书中例5.1 MATLAB程序如下：clear allN=100;n=0:N-1;xn=2*sin(0.48*pi*n)+c

8、os(0.52*pi*n);XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(x(n);title(x(n) N=100);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|X(n)|);title(X(k) N=100);执行结果如下图形所示： 2，实验内容第（2）小题 .指数序列：v1(n)=(0.9)nu(n)MATLAB程序如下：clear allN=100;n=0:N-1;xn

9、=0.9.n;XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v1(n);title(v1(n) N=100);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V1(n)|);title(V1(k) N=100);执行结果如下图形所示：第（2）小题. 周期为N的余弦序列：v2(n)=cos(n)，且0 n N-1MATLAB程序如下：clear allN=100;n=0:N-1;xn

10、=cos(2*pi/N*n);XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v2(n);title(v2(n) N=100);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V2(n)|);title(V2(k) N=100);执行结果如下图形所示： . 符合函数序列：v3(n)=0.9sin(n)+0.6sin(n)MATLAB程序如下：clear allN=100;n=0:N-

11、1;xn=0.9*sin(2*pi/N*n)+0.6*sin(2*pi/N*3*n);XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v3(n);title(v3(n) N=100);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V3(n)|);title(V3k) N=100);执行结果如下图形所示： 第（3）小题 N=128MATLAB程序如下：clear allN=128;n

12、=0:N-1;xn=0.9.n;XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v1(n);title(v1(n) N=128);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V1(n)|);title(V1(k) N=128);执行结果如下图形所示：N=64 MATLAB程序如下：clear allN=64;n=0:N-1;xn=0.9.n;XK=fft(xn,N);magXK=

13、abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v1(n);title(v1(n) N=64);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V1(n)|);title(V1(k) N=64);执行结果如下图形所示：N=32 MATLAB程序如下：clear allN=32;n=0:N-1;xn=0.9.n; t(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);pl

14、ot(n,xn);xlabel(n);ylabel(v1(n);title(v1(n) N=32);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V1(n)|);title(V1(k) N=32);执行结果如下图形所示： 第（4）小题 . N点FFT MATLAB程序如下：clear allN=100;M=N;n=0:N-1;xn=cos(2*pi/N*n);XK=fft(xn,M);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlab

15、el(n);ylabel(v2(n);title(v2(n) M=N);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V2(n)|);title(V2(k) M=N);执行结果如下图形所示：.2N点FFT MATLAB程序如下：clear allN=100;M=2*N;n=0:N-1;xn=cos(2*pi/N*n);XK=fft(xn,M);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v2(n);ti

16、tle(v2(n) M=2*N);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V2(n)|);title(V2(k) M=2*N);执行结果如下图形所示：.N+2点FFT MATLAB程序如下：clear allN=100;M=N+2;n=0:N-1;xn=cos(2*pi/N*n);XK=fft(xn,M);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v2(n);title(v2(n) M=N+2)

17、;subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V2(n)|);title(V2(k) M=N+2);执行结果如下图形所示： 第（5）小题 .32点FFTMATLAB程序如下：clear allN=32;n=0:N-1;xn=0.15*sin(2*pi/32*n)+sin(2*pi*2/32*n)-0.1*sin(2*pi*3/32*n);XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(n,xn);xlabel(n);ylabel(v3(n);title(v3(n) N=32);subplot(1,2,2);k=0:length(magXK)-1;stem(k,magXK,.);xlabel(k);ylabel(|V3(n)|);title(V3k) N=32);执行结果如下图形所示：.64点FFT MATLAB程序如下：clear allN=64;n=0:N-1;xn=0.15*sin(2*pi/32*n)+sin(2*pi*2/32*n)-0.1*sin(2*pi*3/32*n);XK=fft(xn,N);magXK=abs(XK);phaXK=angle(XK);subplot(1,2,1);plot(