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

一、 实验目的 学握离散傅里叶变换的计算机实现方法。 检验实序列傅里叶变换的性质。 掌握计算序列的圆周卷积的方法。 熟悉连续信号经理想采样前后的频谱变化关系，加深对吋域釆样定理的理解。 学习用DFT对连续信号和时域离散信号进行谱分析的方法，了解町能出现的分析误差, 以便在实际中正确应用DFT。 二、 实验内容 实现离散傅里叶变换。 计算序列圆周卷积。 计算实序列傅里叶变换并检验DFT性质。 实现连续信号傅里叶变换以及由不同釆样频率采样得到的离散信号的傅里叶变换。 实现补零序列的傅里叶变换。 实现高密度谱和高分辨率谱，并比较二者的不同。 三、实验结果 例1 ?计算离散傅里叶变换 Circular ?e\en component该图中上图为序列的共轨对称分量;下图为共轨反对量。设序列为 x(n) = { 1 2 3 4 5 6 7 8},对 其进行DFT Circular ?e\en component 该图中上图为序列的共轨对称分量;下图为共轨反对量。 n=0:10; x二 10*(0. 8). n; [xec, xoc]=circevod(x); subplot (2, 1, 1); stem(n, xec); title(，Circular -even componentJ) xlabel ( n) ; ylabcl ( xec (n)); axis([-0. 5, 10. 5,-1, 11]) subpl ot (2, 1,2); stem(n, xoc); title (，Circular -odd component) xlabel ( n) ; ylabcl ( 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 C DET[xec(n)]) ;xlabel ( k‘)； subplot (2, 2, 4) ; stem(n, imag(Xoc)) ;axis ([-0. 5, 10. 5, -20, 20]) titleC DFT[xoc (n) Y ) ; xlabel ( k‘); 例2.计算实序列x(n)二10*(0. 8).、 0=n=10的DFT并检验DFT的性质 DFT [xec (n) ] =Re [X (k) ] DFT [xoc (n)]=Im[X (k)] 程序： n=0:10; 50403020104030100RealDFTlx(n)]10ImagDFQx(n)]0 5 1010DFTlxoc(n)]0 5 10x二 10*(0. 8). *n; 50 40 30 20 10 40 30 10 0 RealDFTlx(n)] 10 ImagDFQx(n)] 0 5 10 10 DFTlxoc(n)] 0 5 10 [xec, xoc]=circevod(x); subplot (2, 1, 1); stem(n, xec); ????ti tle(* Circular -even component) ???? xlabel ( n) ; ylabel ( xec (n)) ;axis ([-0. 5, 10. 5, ~1, 11 DFT{xec(n))1) DFT{xec(n)) subplot (2, 1, 2); stem (n, xoc); title(，Circular -odd component) xlabel ( n) ;y label C 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); subpl ot (2, 2, 1); stem(n, real (X)) ;axi s ([-0. 5, 10. 5, -5, 50]) tit