北邮数字信号处理MATLAB上机实验报告.docVIP

数字信号处理MATLAB上机实验报告 学 院：信息与通信工程学院 专 业： 通信工程 班 级： 班内序号： 学 号： 姓 名： 任课教师：  设计思路 两张计算方法，通过公式计算理论值，用fft函数计算。将结果误差进行比较，用stem函数绘图。 源代码 Q=0.9+0.3i; %定义Q值 N=32; %定义N值 n=0:N-1; %定义n取值从0到N-1 x=Q.^n; %定义x（n） k=0:N-1; %定义k值 W=exp(-2*pi*i/N); Wk=W.^k; format long ; %显示15位小数 XK1=(1-Q.^N)./(1-Q.*Wk) %计算理论值32点DFT变换XK1 XK2=fft(x,N) %FFT算法计算32点DFT变换XK2 XK3=XK2-XK1 %两种计算方法的误差 subplot(4,1,1),stem(k,x), title(原始数列),xlabel(k),ylabel(x（n）) %画x（n）图 subplot(4,1,2),stem(k,abs(XK1)) title(理论值),xlabel(k),ylabel(|XK1|) %画XK1图 subplot(4,1,3),stem(k,abs(XK2)), title(FFT),xlabel(k),ylabel(|XK2|) %画XK2图 subplot(4,1,4),stem(k,abs(XK3)), title(误差),xlabel(k),ylabel(|XK3|) %画XK3图 频谱图 流程图 结果分析 由XK3可知，理论计算得到的结果和FFT算法得到的结果误差在10的-14次方数量级上，可见FFT具有很高的精度，所以在误差允许的范围内，可以用FFT快速算法计算数列的DFT结果。 设计思路 通过矩形窗函数，取N个点进行DFT变换进行频谱分析，第一项前面系数很小，所以在观察结果的时候要经过放大才可以看到。这里N取1000。 源代码 N1=1000; n=0:N1; xn=0.001*cos(0.45*pi*n)+sin(0.3*pi*n)-cos(0.302*pi*n-0.25*pi); XK=fft(xn,N1); subplot(1,1,1) m=0:1:N1-1; k=m/500 %归一化 subplot(3,1,1),stem(k,abs(XK)), %画0-0.5频谱图 title(频谱),xlabel(x),ylabel(|XK|),axis([0,0.5,0,600]) subplot(3,1,2),stem(k,abs(XK)), %画0.29-0.31频谱图 title(频谱),xlabel(x),ylabel(|XK|),axis([0.29,0.31,0,600]) subplot(3,1,3),stem(k,abs(XK)), %画0.29-0.46频谱图 title(频谱),xlabel(x),ylabel(|XK|),axis([0.29,0.46,0,1]) 频谱图 流程图 结果分析 由上面三张频谱图可以清晰的看到0.3pi、0.302pi、0.45pi的三张谱线。其中，第三张图被放大，因为0.45pi系数为0.001，相对0.3pi、0.302pi太小，所以要放大才能看到。 心得体会 通过这次上机作业，我对MATLAB有了一个初步的了解和认识，掌握了一些函数的功能和使用方法。MATLAB并没有多么复杂的语法，主要是一些功能函数的学习。在写代码调试的过程中也遇到了不少问题，在不断的改正参考书籍以及向同学询问最终完成了整个实验。 MATLAB是个功能很强大的软件，可以在很多学科应用，这里我们通过实践学习并体会了MATLAB在DSP中的应用，借助MATLAB可以快速的计算DFT并且绘制频谱图。同时通过这次上机实验也加深了我对FFT的理解。 开始 设置变量和参数 计算32点DFT 分析频谱并绘图 结束 开始 设置变量和参数 计算DFT 分析频谱并绘图 结束