数字信号处理上机实验期末考试.docxVIP

数字信号处理实验期末考试 P31 第16题 本页由陈恕提供 P31 16 ys=0;B1=1;A1=[1,-0.6,0.08]; xn=[1,zeros(1,30)]; xi=filtic(B1,A1,ys); hn1=filter(B1,A1,xn,xi); n=0:length(hn1)-1; subplot(2,2,1);stem(n,hn1,.) title((a)系统(1)的单位脉冲响应);xlabel(n);ylabel(h(n)); xn=ones(1,30); sn1=filter(B1,A1,xn,xi); n=0:length(sn1)-1; subplot(2,2,2);stem(n,sn1,.) title((b)系统(1)的单位阶跃响应);xlabel(n);ylabel(s(n)) ys=0;B2=[2,0,-1];A2=[1,-0.7,0.1]; xn=[1,zeros(1,30)]; xi=filtic(B2,A2,ys); hn2=filter(B2,A2,xn,xi); n=0:length(hn2)-1; subplot(2,2,3);stem(n,hn2,.) title((c)系统(2)的单位脉冲响应);xlabel(n);ylabel(h(n)) xn=ones(1,30); sn2=filter(B2,A2,xn,xi); n=0:length(sn2)-1; subplot(2,2,4);stem(n,sn2,.) title((d)系统(2)的单位阶跃响应);xlabel(n);ylabel(s(n)) 数字信号处理实验期末考试 P31第16题(图)、第17题 本页由符德山提供 17 B=0.866;A=[1,-0.8,0.64]; xn=[1,zeros(1,48)]; hn=filter(B1,A1,xn); n=0:length(hn)-1; subplot(2,2,1);stem(n,hn,.) title((a)系统的单位脉冲响应);xlabel(n);ylabel(h(n)); xn=ones(1,100); sn=filter(B,A,xn); n=0:length(sn)-1; subplot(2,2,2);stem(n,sn,.);axis([0,30,0,2]) title((b)系统的零状态单位阶跃响应);xlabel(n);ylabel(s(n)) 数字信号处理实验期末考试 P109第23题 本页由杨迎洲提供 P109 23 xn=[1 2 3 3 2 1]; %输入时域序列向量x(n) N=32;M=1024; Xjw=fft(xn,M); %计算xn的1024点DFT,近似表示序列的傅里叶变换 Xk32=fft(xn,N); %计算xn的32点DFT xn32=ifft(Xk32,N); %计算Kk32的32点IDFT %以下为绘图部分 k=0:M-1;wk=2*k/M; %产生M点DFT对应的采样点频率(关于π归一化值) subplot(3,2,1);plot(wk,abs(Xjw)); %绘制M点DFT的幅频特性图 title((a) FT[x(n)]的幅频特性图);xlabel(ω/π);ylabel(幅度) subplot(3,2,3);plot(wk,angle(Xjw)); %绘制x(n)的相频特性图 line([0,2],[0,0]) %画横坐标轴线 title((b)FT[x(n)]的相频特性图); xlabel(ω/π);ylabel(相位);%axis([0,2,-3.5,3.5]) k=0:N-1; subplot(3,2,2);stem(k,abs(Xk32),.); %绘制64点DFT的幅频特性图 title((c)32点DFT的幅频特性图); xlabel(k);ylabel(幅度);axis([0,32,0,15]) subplot(3,2,4);stem(k,angle(Xk32),.); %绘制64点DFT的相频特性图 line([0,32],[0,0]) %画横坐标轴线 title((d)32点DFT的相频特性图) xlabel(k);ylabel(相位);axis([0,32,-3.5,3.5]) figure(2) k=0:M-1;wk=2*k/M; %产生M点DFT对应的采样点频率(关于π归一化值) subplot(3,2,1);plot(wk,abs(Xjw)); %绘制M点DFT的幅频特性图 title((e) FT[x(n)]和32点DFT[x(n)]的幅频特性);xlabel(ω/