[数字信号处理实验14.docxVIP

厦门理工学院实验报告书课程名称：数字信号处理实验实验名称：离散系统的零极点分析专业班级：姓名：学号：日期：地点：成绩：教师签名：光电与通信工程学院一.实验目的(1)加深对快速傅里叶变换（FFT）基本理论的理解。(2)了解使用快速傅里叶变换（FFT）计算有限长序列和无限长序列信号频谱的方法。(3)掌握用MATLAB语言进行快速傅里叶变换时常用的子函数。二．实验原理1．用MATLAB提供的子函数进行快速傅里叶 从理论学习可知，DFT是唯一在时域和频域均为离散序列的变换方法，它适用于有限长序列。尽管这种变换方法是可以用于数值计算的，但如果只是简单的按照定义进行数据处理，当序列长度很大时，则将占用很大的内存空间，运算时间将很长。 快速傅里叶变换是用于DFT运算的高效运算方法的统称，FFT只是其中的一种。FFT主要有时域抽取算法和频域抽取算法，基本思想是将一个长度为N的序列分解成多个短序列，如基2算法、基4算法等，大大缩短了运算的时间。 MATLAB中提供了进行快速傅里叶变换(FFT)的子函数，用fft计算DFT，用ifft计算IDFT。2．用FFT计算有限长序列的频谱 （1）基本概念 一个序号从n1到n2的时域有限长序列x(n)，它的频谱X(ejw)定义为它的离散傅里叶变换，且在奈奎斯特(Nyquist)频率范围内有界并连续。序列的长度为N，则N＝n2－n1＋1。计算x(n)的离散傅里叶变换(DFT)得到的是X(ejw)的N个样本点X(ejwk)。其中数字频率为 式中：dw为数字频率的分辨率；k取对应－(N－1)/2到(N－1)/2区间的整数。 在实际使用中，往往要求计算出信号以模拟频率为横坐标的频谱，此时对应的模拟频率为 式中：D为模拟频率的分辨率或频率间隔；Ts为采样信号的周期，Ts＝1/Fs；定义信号时域长度L＝NTs。在使用FFT进行DFT的高效运算时，一般不直接用n从n1到n2的x(n)，而是取x（n） 的主值区间(n＝0，1，…，N－1)的数据，经FFT将产生N个数据，定位在k＝0，1，…，N－1的数字频率点上，即对应［0，2p］。如果要显示［－p，p］范围的频谱，则可以使用fftshift(X)进行位移。三．实验任务已知有限长序列x(n)＝［1，0.5，0，0.5，1，1，0.5，0］，要求：用FFT算法求该时域序列的DFT、IDFT的图形； 解： MATLAB程序如下：xn=[1,0.5,0,0.5,1,1,0.5,0]; N=length(xn); n=0:N-1;k=0:N-1; Xk=fft(xn,N); subplot(2,1,1);stem(k,abs(Xk)); title(Xk=DFT(x(n))); xn1=ifft(Xk,N); subplot(2,1,2);stem(n,xn1); title(x(n)=IDFT(Xk));假定采样频率Fs＝20 Hz，序列长度N分别取8、32和64，使用FFT来计算其幅度频谱和相位频谱。其matlab程序如下：Fs=20;C=[8,32,64]; for r=0:2; N=C(r+1); xn=[1,0.5,0,0.5,1,1,0.5,0]; D=2*pi*Fs/N; k=floor(-(N-1)/2:(N-1)/2); X=fftshift(fft(xn,N)); subplot(3,2,2*r+1);plot(k*D,abs(X)); title(幅度频谱);xlabel(rad/s); subplot(3,2,2*r+2);stairs(k*D,angle(X)); title(相位频谱);xlabel(rad/s); end2 .已知一个无限长序列x(n)＝0.5n(n≥0)，采样周期Ts＝0.2 s，要求序列长度N分别取8、32和64，用FFT求其频谱。其matlab程序如下：Ts=0.2;C=[8,32,64]; for r=0:2; N=C(r+1); n=0:N-1; xn=exp(0.5*n); D=2*pi/(N*Ts); k=floor(-(N-1)/2:(N-1)/2); X=fftshift(fft(xn,N)); subplot(3,2,2*r+1);plot(k*D,abs(X)); title(幅度频谱);xlabel(rad/s);subplot(3,2,2*r+2);stairs(k*D,angle(X)); title(相位频谱);xlabel(rad/s);end四．实验涉及的MATLAB子函数（1）fft功能：一维快速傅里叶变换(FFT)。调用格式：y=fft（x）；利用FFT算法计算矢量x的离散傅里叶变换，当x为矩阵时，y为矩阵x每一列的FFT。当x的长度为2