数字信号处理实验(matlab)

《数字信号处理实验(matlab)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、南京研旭电气有限公司介绍dsp知识，为大家提供最新的dsp资讯，更多内容可以去南京研旭电气科技有限公司的官网www.njyxdq.comwww.f28335.com或者官方论坛，嵌嵌dsp论坛www.armdsp.net进行交流学习欢迎大家收听嵌嵌dsp论坛的官方微博http://t.qq.com/qianqiandspDSP之matlab实验1实验题目第1页南京研旭电气有限公司2.程序分析实验一(源代码及分析)A=zeros(1,25)c=complex(0.9,0.3)fori=1:1:25A(i)=c^iendn

2、=0:24k=0:24WN=exp(-j*2*pi/25)第2页南京研旭电气有限公司nk=n'*kWNnk=WN.^nkX=A*WNnk;%X(k)为用定义算法计算的DFT值%下面为32点运算fori=26:1:32A(i)=0endn=0:31k=0:31WN=exp(-j*2*pi/32)nk=n'*kWNnk=WN.^nkX=A*WNnk%下面为32点蝶形运算k1=bin2dec(fliplr(dec2bin([1:32]-1,5)))+1%得到从1到32的倒序顺序B=A(k1)%得到倒序后的A，将其赋给B第3页

3、南京研旭电气有限公司formm=1:5%将DFT做m次基2分解，从左到右，对每次分解作DFT运算Nmr=2^mm;u=1;%旋转因子u初始化WN=exp(-j*2*pi/Nmr);%本次分解的基本DFT因子WN＝exp(-i*2*pi/Nmr)forn=1:Nmr/2%本次跨越间隔内的各次碟形运算fork=n:Nmr:32%本次碟形运算的跨越间隔为Nmr=2^mmkp=k+Nmr/2;%确定碟形运算的对应单元下标(对称性)t=B(kp)*u;%碟形运算的乘积项B(kp)=B(k)-t;%碟形运算的加法项B(k)=B(k

4、)+t;endu=u*WN;%修改旋转因子，多乘一个基本DFT因子WNendend实验二(源代码及分析)第4页南京研旭电气有限公司fori=1:1:1000x(i)=0.001*cos(0.45*i*pi)+sin(0.3*i*pi)-cos(0.302*i*pi-4/pi)endX=fft(x);k=0:500;w=2*pi/1000*k;stem(w/pi,abs(X(1:501)))axis([0.250.35490510])axis([0.30.500.8])3.程序运行结果实验一程序运行结果太长,只附一截图如

5、下:第5页南京研旭电气有限公司实验二结果的有窗与无窗图如下:第6页南京研旭电气有限公司第7页南京研旭电气有限公司4.总结本次实验是第一次使用matlab,实验过程中遇到了不少的问题,经过和同学们的交流以及通过网络查询,最后得到了结果。尽管实验过程很是坎坷,与此同时也收获了很多,不但对原来的是理论知识有了很好的理解,对我们的学习有了一定的帮助,还让我了解到了matlab软件的强大。我想以后再遇到用matlab实验的时候,肯定就不会这么费劲了,这就是我的收获。相信我们在不断地进步下能过做的更好。介绍dsp知识，为大家提供最

6、新的dsp资讯，更多内容可以去南京研旭电气科技有限公司的官网www.njyxdq.comwww.f28335.com或者官方论坛，嵌嵌dsp论坛www.armdsp.net进行交流学习欢迎大家收听嵌嵌dsp论坛的官方微博http://t.qq.com/qianqiandsp第8页