[实验二预作结果序列的DFT及谱分析.docVIP

实验二 离散傅立叶变换及谱分析预做实验 实验目的 1．掌握离散傅里叶变换的计算机实现方法。 2．检验实序列傅里叶变换的性质。 3．掌握计算序列的圆周卷积的方法。 4．熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 5．学习用DFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差，以便在实际中正确应用DFT。 实验内容 1．实现离散傅里叶变换。 2．计算序列圆周卷积。 3．计算实序列傅里叶变换并检验DFT性质。 4．实现连续信号傅里叶变换以及由不同采样频率采样得到的离散信号的傅里叶变换。 5．实现补零序列的傅里叶变换。 6．实现高密度谱和高分辨率谱，并比较二者的不同。 三、实验结果 例1．计算离散傅里叶变换 设序列为 x(n)={ 1 2 3 4 5 6 7 8},对其进行DFT n=0:10; x=10*(0.8).^n; [xec,xoc]=circevod(x); subplot(2,1,1);stem(n,xec); title(Circular -even component) xlabel(n);ylabel(xec(n)); axis([-0.5,10.5,-1,11]) subplot(2,1,2);stem(n,xoc); title (Circular -odd component) 　 xlabel(n);ylabel(xoc(n)); axis([-0.5,10.5,-4,4]) figure(2) 该图中上图为序列的共轭对称分量；下图为共轭反对量。 X=dft(x,11); Xec=dft(xec,11); Xoc=dft(xoc,11); subplot(2,2,1);stem(n,real(X));axis([-0.5,10.5,-5,50]) title (Real{DFT[x(n)]});xlabel(k); subplot(2,2,2);stem(n,imag(X));axis([-0.5,10.5,-20,20]) title (Imag{DFT[x(n)]});xlabel(k); subplot(2,2,3);stem(n,Xec);axis([-0.5,10.5,-5,50]) title (DFT[xec(n)]);xlabel(k); subplot(2,2,4);stem(n,imag(Xoc));axis([-0.5,10.5,-20,20]) title(DFT[xoc(n)]);xlabel(k); 例2．计算实序列x(n)=10*(0.8).^n 0=n=10 的DFT并检验DFT的性质 DFT[xec(n)]=Re[X(k)] DFT[xoc(n)]=Im[X(k)] 程序： n=0:10; x=10*(0.8).^n; [xec,xoc]=circevod(x); subplot(2,1,1);stem(n,xec); title(Circular -even component) xlabel(n);ylabel(xec(n));axis([-0.5,10.5,-1,11]) subplot(2,1,2);stem(n,xoc); title(Circular -odd component) 　 xlabel(n);ylabel(xoc(n));axis([-0.5,10.5,-4,4]) figure(2) X=dft(x,11); 可以看出 DFT[xec(n)]=Re[X(k)] DFT[xoc(n)]=Im[X(k) Xec=dft(xec,11); Xoc=dft(xoc,11); subplot(2,2,1);stem(n,real(X));axis([-0.5,10.5,-5,50]) title(Real{DFT[x(n)]});xlabel(k); subplot(2,2,2);stem(n,imag(X));axis([-0.5,10.5,-20,20]) title(Imag{DFT[x(n)]});xlabel(k); subplot(2,2,3);stem(n,Xec);axis([-0.5,10.5,-5,50]) title(DFT[xec(n)]);xlabel(k