资源描述：
《2012数字信号处理实验讲义 - 实验二报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验二离散傅立叶变换及谱分析班级：电子101班姓名：李怡敏学号：102045一、实验目的1．掌握离散傅里叶变换的计算机实现方法。2．检验实序列傅里叶变换的性质。3．掌握计算序列的循环卷积的方法。4．学习用DFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差，以便在实际中正确应用DFT。二、实验内容1．实现序列的离散傅里叶变换并对结果进行分析。（自己选择序列，要求包括复序列，实序列，实偶序列，实奇序列，虚奇序列(已测)）2．计算序列循环卷积。3．计算实序列傅里叶变换并检验DFT性质。4．实现补零序列的傅里叶变换。6．实现
2、高密度谱和高分辨率谱，并比较二者的不同。三、实验报告要求%例1本例检验实序列的性质DFT[xec(n)]=Re[X(k)]DFT[xoc(n)]=Im[X(k)]%设x(n)=10*(0.8).^n0<=n<=10将x(n)分解为共扼对称及共扼反对称部分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)');ax
3、is([-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)X=dft(x,11);%求出序列的DFTXec=dft(xec,11);%求序列的共扼对称分量的DFTXoc=dft(xoc,11);%求序列的共扼反对称分量的DFTsubplot(2,2,1);stem(n,real(X));axis
4、([-0.5,10.5,-5,50])title('Real{DFT[x(n)]}');xlabel('k');%画出序列DFT的实部subplot(2,2,2);stem(n,imag(X));axis([-0.5,10.5,-20,20])title('Imag{DFT[x(n)]}');xlabel('k');%画出序列DFT的虚部subplot(2,2,3);stem(n,Xec);axis([-0.5,10.5,-5,50])title('DFT[xec(n)]');xlabel('k');subplot(2,2,4);ste
5、m(n,imag(Xoc));axis([-0.5,10.5,-20,20])title('DFT[xoc(n)]');xlabel('k');(1)结果图:结果:xoc=Columns1through1003.46312.52891.72110.99940.3277-0.3277-0.9994-1.7211-2.5289Column11-3.4631(1)实验结果说明:实序列可以分解为共轭对称部分和共轭反对称部分，且其傅里叶变换实部与共轭对称部分的傅里叶变换相对应，傅里叶变换虚部与其共轭反对称部分的傅里叶变换相对应。%例2本例为计算序
6、列的循环卷积程序%运行之前应在命令窗口输入x1,x2,N的值%实验报告要求：自己选择2个序列进行计算，将实验结果写出iflength(x1)>Nerror('Nmustbe>=thelengthofx1')endiflength(x2)>Nerror('Nmustbe>=thelengthofx2')endx1=[x1zeros(1,N-length(x1))];%将x1,x2补0成为N长序列x2=[x2zeros(1,N-length(x2))];m=[0:1:N-1];x2=x2(mod(-m,N)+1);%该语句的功能是将序列翻褶
7、，延拓，取主值序列H=zeros(N,N);forn=1:1:N%该for循环的功能是得到x2序列的循环移位矩阵H(n,:)=cirshftt(x2,n-1,N);%和我们手工计算循环卷积得到的表是一致的endy=x1*H'%用矩阵相乘的方法得到结果选择序列及结果如下：x1=[-1,1,2,-3,4,5,6,0,8,10];>>x2=[1,1,1,1,1,1,1,1,1,1];>>y=xhjj(x1,x2,10);y=32323232323232323232%例4本例说明补零序列的离散傅立叶变换%序列，已给出序列的傅立叶变换程序和将原序
8、列补零到10长序列的DFTn=0:4;x=[ones(1,5)];%产生矩形序列k=0:999;w=(pi/500)*k;X=x*(exp(-j*pi/500)).^(n'*k);%计算离散时间傅立叶变换X
