《离散傅里叶变换及性质c.docVIP

实验报告格式 院系： 专业： 班级： 实验名称：离散傅立叶变换及性质 课程名称： 数字信号处理 实验室： 实验C楼201 成绩 实验人 学号： 指导教师：赵发勇 实验目的： 1、进一步加深离散傅立叶变换DFT算法原理和基本性质的理解 2、掌握离散傅立叶变换DFT的MATLAB实现，理解DFT的频谱泄漏和分辨率问题。 实验内容： xn=input(xn=); N=length(xn); subplot(211);stem([0:N-1],xn); grid; X=zeros(1,N); for k=0:N-1 for n=0:N-1 X(k+1)=X(k+1)+xn(n+1)*exp(-j*2*pi*n*k/N);%按定义式计算序列的DFT end; end; subplot(212);stem([0:N-1],abs(X)); grid; 实验结果分析 1.1 xn=[1+[0:3],8-[4:7]] 1.1.1 xn=[1+[0:3],8-[4:7],0*[8:15]] 或xn=[1+[0:3],8-[4:7],zeros(1,8)] 1.1.2 xn=[4-[0:3],[4:7]-3,0*[8:15]] 或xn=[4-[0:3],[4:7]-3,zeros(1,8)] 1.2.1 xn=[4-[0:3],[4:7]-3] 1.2.2 xn=sin(pi*[0:7]/8) 1.2.3 xn=sin(pi*[0:15]/8) 1.3.1 xn=cos(8*pi*[0:15]/64)+cos(16*pi*[0:15]/64)+cos(20*pi*[0:15]/64) 1.3.2 xn=cos(8*pi*[0:31]/64)+cos(16*pi*[0:31]/64)+cos(20*pi*[0:31]/64) 1.3.3 xn=cos(8*pi*[0:63]/64)+cos(16*pi*[0:63]/64)+cos(20*pi*[0:63]/64) 1.4.1令x(n)=x2(n)+x3(n)，计算其 8离散傅里叶变换， X(k)=DFT［x(n)］ xn=x2(n)+x3(n)=sin(pi*[0:7]/8)+[4-[0:3],[4:7]-3] 1.4.2 令x(n)=x2(n)+x3(n)，计算其 8离散傅里叶变换， X(k)=DFT［x(n)］ xn=sin(pi*[0:15]/8)+[4-[0:3],[4:7]-3,0*[8:15]] 2.1.1给定正弦信号。现对其进行抽样，设抽样点数N=16，给定抽样频率. xn=sin(2*pi*50*[0:15]/150) 2.1.2 给定正弦信号取时，观察16点DFT的频谱。 xn=sin(2*pi*50*[0:15]/200) 2.2给定正弦信号, 取fs=200HZ,N=32时观察32点DFT的频谱。 xn=sin(2*pi*50*[0:31]/200) 2.(2)xn=sin(2*pi*50*[[0:15],0*[16:31]]/200) 取时，在抽样点后再补N个零，观察32点DFT的频谱。 思考题解答 （1）实验1中序列1和序列2在N=8时的幅频特性是否相同，为什么？N=16呢？ 由以上图像可知序列1和序列2在N=8时的幅频特性相同，因为周期延拓，它们的图像有一部分相同，当抽样段不同他们的幅频特性就可以相同。当N=16时是周期的整数倍，此时它们的幅频特性就不相同，如1.1.1图和1.1.2图所示。 （2）结合实验中所得的给定序列的幅频特性曲线，与理论结果比较，并分析说明误差产生的原因及用FFT作谱分析时有关参数的选择方法。 误差产生的原因：1、信号做频谱分析时的“频率分辨率”影响。2、加窗对频域分析带来的频谱的“泄漏”，边瓣越大，且衰减的越慢，泄漏越严重；在频谱分析中，泄漏往往会模糊原来真正谱的形状，窗函数过大的边瓣有可能产生虚假的峰值。 用FFT作谱分析时有关参数的选择方法：1、若已知信号的最高频率fc,为了防止混叠，选定采样频率fs：fs=2fc; 2、根据实际需要，选定频率分辨变化量f,一旦选定后，即可确定FFT所需的点数N,N=fs/(变化量f）;变化量f越小越好，但变化量f越小，N越大，计算量、存储量也随之增大。一般取N为2的整次幂，变用FFT计算，若已给N，可用补零方法使N为2的整次幂； 3、fs和N确定后，即可确定所需相应模拟信号x(t)的长度T=N/fs.分辨率变化量f反比于T，而不是N，在给定的T情况下，靠减小T