数字信号处理实验(matlab).docVIP

北京邮电大学 DSP之matlab实验 1实验题目 2. 程序分析 实验一(源代码及分析) A=zeros(1,25) c=complex(0.9,0.3) for i=1:1:25 A(i)=c^i end n=0:24 k=0:24 WN=exp(-j*2*pi/25) nk=n*k WNnk=WN.^nk X=A*WNnk;%X(k)为用定义算法计算的DFT值 %下面为32点运算 for i=26:1:32 A(i)=0 end n=0:31 k=0:31 WN=exp(-j*2*pi/32) nk=n*k WNnk=WN.^nk X=A*WNnk %下面为32点蝶形运算 k1=bin2dec(fliplr(dec2bin([1:32]-1,5)))+1%得到从1到32的倒序顺序 B=A(k1)%得到倒序后的A，将其赋给B for mm=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) for n=1:Nmr/2 %本次跨越间隔内的各次碟形运算 for k=n:Nmr:32 %本次碟形运算的跨越间隔为Nmr=2^mm kp=k+Nmr/2; %确定碟形运算的对应单元下标(对称性) t=B(kp)*u; %碟形运算的乘积项 B(kp)=B(k)-t; %碟形运算的加法项 B(k)=B(k)+t; end u=u*WN; %修改旋转因子，多乘一个基本DFT因子WN end end 实验二(源代码及分析) for i=1:1:1000 x(i)=0.001*cos(0.45*i*pi)+sin(0.3*i*pi)-cos(0.302*i*pi-4/pi) end X=fft(x); k=0:500; w=2*pi/1000*k; stem(w/pi,abs(X(1:501))) axis([0.25 0.35 490 510]) axis([0.3 0.5 0 0.8]) 3. 程序运行结果 实验一程序运行结果太长,只附一截图如下: 实验二结果的有窗与无窗图如下: 4. 总结 本次实验是第一次使用matlab,实验过程中遇到了不少的问题,经过和同学们的交流以及通过网络查询,最后得到了结果。尽管实验过程很是坎坷,与此同时也收获了很多,不但对原来的是理论知识有了很好的理解,对我们的学习有了一定的帮助,还让我了解到了matlab软件的强大。我想以后再遇到用matlab实验的时候,肯定就不会这么费劲了,这就是我的收获。相信我们在不断地进步下能过做的更好。北京邮电大学信息与通信工程学院 investigations by person/in accordance with credit rules, personnel of customer credit investigation and analysis. 1) active marketing, expanding customers; 2) arrange interviews with our customers; 3) qualifying first instance; 4) accepts an application for the borrower; 5) collect, collate, review customer information; 6) due investigations by person/in accordance with credit rules, personnel of customer credit investigation and analysis. 1) active marketing, expanding customers; 2) arrange interviews with our customers; 3) qualifying first instance; 4) accepts an application for the borrower; 5) collect, collate, review customer information; 6) due第7页 北京邮电大学电信工程学院 investigations by person/in accordance with credit rules, personnel of customer credit investigation and analysis. 1) active marketing, expanding c