数字信号处理实验义实验二研究分析报告.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 实验二 离散傅立叶变换及谱分析 班级：电子101班 姓名：李怡敏 学号：102045 实验目地 1．掌握离散傅里叶变换地计算机实现方法. 2．检验实序列傅里叶变换地性质. 3．掌握计算序列地循环卷积地方法. 4．学习用DFT对连续信号和时域离散信号进行谱分析地方法，了解可能出现地分析误差，以便在实际中正确应用DFT.b5E2RGbCAP 实验内容 1．实现序列地离散傅里叶变换并对结果进行分析.（自己选择序列，要求包括复序列，实序列，实偶序列，实奇序列，虚奇序列(已测)）p1EanqFDPw 2．计算序列循环卷积. 3．计算实序列傅里叶变换并检验DFT性质. 4．实现补零序列地傅里叶变换. 6．实现高密度谱和高分辨率谱，并比较二者地不同. 实验报告要求 %例1 本例检验实序列地性质DFT[xec(n)]=Re[X(k)]DFT[xoc(n)]=Im[X(k)]DXDiTa9E3d %设 x(n)=10*(0.8).^n 0=n=10 将x(n)分解为共扼对称及共扼反对称部分RTCrpUDGiT 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])5PCzVD7HxA subplot(2,1,2);stem(n,xoc); %画出序列地共扼反对称分量 title(Circular -odd component) xlabel(n);ylabel(xoc(n));axis([-0.5,10.5,-4,4])jLBHrnAILg figure(2) X=dft(x,11); %求出序列地DFT Xec=dft(xec,11); %求序列地共扼对称分量地DFT Xoc=dft(xoc,11); %求序列地共扼反对称分量地DFT subplot(2,2,1);stem(n,real(X));axis([-0.5,10.5,-5,50])xHAQX74J0X title(Real{DFT[x(n)]});xlabel(k); %画出序列DFT地实部LDAYtRyKfE subplot(2,2,2);stem(n,imag(X));axis([-0.5,10.5,-20,20])Zzz6ZB2Ltk title(Imag{DFT[x(n)]});xlabel(k); %画出序列DFT地虚部dvzfvkwMI1 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])rqyn14ZNXI title(DFT[xoc(n)]);xlabel(k); 结果图: 结果: xoc = Columns 1 through 10 0 3.4631 2.5289 1.7211 0.9994 0.3277 -0.3277 -0.9994 -1.7211 -2.5289EmxvxOtOco Column 11 -3.4631 实验结果说明: 实序列可以分解为共轭对称部分和共轭反对称部分，且其傅里叶变换实部与共轭对称部分地傅里叶变换相对应，傅里叶变换虚部与其共轭反对称部分地傅里叶变换相对应.SixE2yXPq5 % 例2 本例为计算序列地循环卷积程序 % 运行之前应在命令窗口输入 x1,x2,N 地值 %实验报告要求：自己选择2个序列进行计算，将实验结果写出 if length(x1)N error(N mustbe = the length of x1) end if length(x2)N error(N mustbe = the length of x2) end x1=[x1 zeros(1,N-length(x1))]; %将x1,x2补0成为N长序列 x2=[x2 zeros(1,N-length(x2))]; m=[0:1:N-1]; x2=x2(mod(-m,N)+1); %该语句地功能是将序列翻褶，延拓，取主值序列 H=zeros(N,N); for n=1:1:N