资源描述：
《实验二预作结果 序列的DFT及谱分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验二离散傅立叶变换及谱分析预做实验一、实验目的1．掌握离散傅里叶变换的计算机实现方法。2．检验实序列傅里叶变换的性质。3．掌握计算序列的圆周卷积的方法。4．熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。5．学习用DFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差，以便在实际中正确应用DFT。二、实验内容1．实现离散傅里叶变换。2．计算序列圆周卷积。3．计算实序列傅里叶变换并检验DFT性质。4．实现连续信号傅里叶变换以及由不同采样频率采样得到的离散信号的傅里叶变换。5．实现补零序列的傅里叶变换。6．实现高密度谱和高分辨率谱，并比较二者的不同
2、。三、实验结果1．计算离散傅里叶变换设序列为x(n)={12345678},对其进行DFT程序：x=[12345678]；Xk=dft(x,8);结果：Xk=Columns1through336.0000-4.0000+9.6569i-4.0000+4.0000iColumns4through6-4.0000+1.6569i-4.0000-4.0000-1.6569iColumns7through8-4.0000-4.0000i-4.0000-9.6569i2．计算实序列x(n)=10*(0.8).^n0<=n<=10的DFT并检验DFT的性质DFT[xec(n)]=Re[X(k)
3、]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-evencomponent')xlabel('n');ylabel('xec(n)');axis([-0.5,10.5,-1,11])subplot(2,1,2);stem(n,xoc);title('Circular-oddcomponent') xlabel('n');ylabel('xoc(n)');axis([-0.5,10.5,-4,4])figure(2
4、)X=dft(x,11);9Xec=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(
5、n)]');xlabel('k');subplot(2,2,4);stem(n,imag(Xoc));axis([-0.5,10.5,-20,20])title('DFT[xoc(n)]');xlabel('k');该图中上图为序列的共轭对称分量；下图为共轭反对称分量。9该图中左上角的图为序列DFT的实部；右上角图为DFT的虚部左下角的图为序列的共轭对称部分的DFT；右下角的图为序列的共轭反对称部分的DFT可以看出DFT[xec(n)]=Re[X(k)]DFT[xoc(n)]=Im[X(k)3．计算序列的循环卷积程序：iflength(x1)>Nerror('Nmustbe>=th
6、elengthofx1')endiflength(x2)>Nerror('Nmustbe>=thelengthofx2')endx1=[x1zeros(1,N-length(x1))];x2=[x2zeros(1,N-length(x2))];m=[0:1:N-1];x2=x2(mod(-m,N)+1);H=zeros(N,N);forn=1:1:NH(n,:)=cirshftt(x2,n-1,N);endy=x1*H'将以上程序存入文件a2.m，然后在matlabcommand窗口输入：x1=[352768];x2=[14637];N=6;则运行结果为：y=10613210611
7、082115若改变循环卷积长度：N=10则运行结果为：y=3174054821151031156656计算两个序列的线性卷积，即z=conv(x1,x2)则运行结果为：z=3174054821151031156656通过以上结果可以看出：当循环卷积的长度大于等于线性卷积的长度时，二者结果相等；当循环卷积的长度小于线性卷积的长度时，二者是不等的。4．实现连续信号傅里叶变换以及由不同采样频率采样得到的离散信号的傅里叶变换令9，绘制其幅频特性。对其用不同采样频率进行采样，绘
