第二次DSP上机实验报告MATLAB实验报告

《第二次DSP上机实验报告MATLAB实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、DSP-MATLAB实验报告实验名称：学生姓名：班级：班内序号：学号：日期：1.实验要求与目的要求：(1)用DFT分析频谱，得到高分辨率频谱清楚的谱线；⑵利用基2时间抽选的FFT算法，计算有限点的DFT,并与理论值相比较。目的：掌握(a)用傅立叶变换进行信号分析吋基本参数的选择。(b)经过离散时间傅立叶变换(DTFT)和有限长度离散傅立叶变换(DFT)后信号频谱上的区别，前者DTFT时间域是离散信号，频率域还是连续的，而DFT在两个域中都是离散的。(c)离散傅立叶变换的基本原理、特性，以及经典的快速算法(基2时间抽选法)，体会快速算法的效率。(d)获得一个高密度频谱和高分辨率频谱的

2、概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析(例如小波)的学习和研究打下基础。2关键算法分析代码的重点是利用fft函数计算离散序列的DFT,难点是矩阵的加减乘除及乘方中所涉及的matlab语句的特点。实验代码如下：(1)N=1000;n=(0:N・l);y=0.001*cos(0・45*n*pi)+sin(0.3*n*pi)・cos(0.302*n*pi・pi/4);y=0・001*cos(0.45*n*pi)+sin(0・3*n*pi)・cos(0.302*n*pi・pi/4);XK=fft(y,N);-fft函数两个参数分別是进行DFT变换的序列，及D

3、FT的长度subplot(3,1,1);stem(n,y);title(,时域波形');xlabel(h);ylabel('y');subplot(3,l,2);stem(abs(XK)/.');axis([0,500,0,600]);-限定x轴与y轴的范围titleC频域前500点波形,);xlabelCk,);ylabel(,XK,);subplot(3,l,3);stem(abs(XK),T);axis([800,900,0,600]);—x轴的范围很小，目的是能看到两条很近的谱线titleC频域后500点部分波形’);xlabel(k);ylabeK'XK);(2)N=25

4、;Q=0・9+j*0・3;n=(0:24);x=Q.An;WN=exp(-j*2*pi/N);k=(0:24);WK=WN.Ak;XK=(l-Q.AN)./(l-Q*WK);stem(XK);—理论值的计算XKl=fft(x,32);-32点基2时间抽选法的计算subplot(3,1,1);stem(n,x);titleC时域波形')；xlabel(,n,);ylabel(,x,);subplot(3,1,2);stem(XK);title(f频域理论值,);xlabel(,k,);ylabel(,XK,);subplot(3,l,3)；stem(XKl);title(f频域波形'

5、);xlabel(k);ylabel('XKf);3.程序运行结果(1)(2)4.问题分析这次实验的题目相对简单，编程过程中遇到的问题主要是矩阵的加减乘除及乘方运算所对应的matlab语句的规则和fft函数的使用，fft函数的使用方法通过help语句便可查得，其两个参数：第一个是指进行fft运算的序列，第二个值fft运算的长度。而矩阵的运算规则需要上网查看相关资料或与同学讨论解决。总的来说，这次实验比第一次编程时语句的运用要熟练，而但编程比较顺利5.心得与总结第一次用matlab编程时遇到了很多困难，其实现在想一下，用matlab编程时不愁没有吋间编，就愁拿出了好多吋间，但真到编吋

6、总是报错，乂找不到错在哪里。这是让人很纠结的。这次编程我放弃了闭门造车，一个人冥思苦想的方式，而是在编的过程中与同学讨论，互相纠错，最后相对顺利的把程序编了出来。另一个比较有感触的是熟能生巧。第一次编程时什么都不懂，遇到不会编的语句就到网上去查，查到语句后又不知道那些语句中的参数代表什么，通过老师课上的讲解，我学会了用help语句，而且也积累了一些基本的函数，像axis,subplot,stem等，所以显示波形的操作做起来就比较得心应手了，多练习总是有好处的。