1 程序清单及波形显示: clc;close all;clear all; %======内容1:调用filter 解差分方程,由系统对u(n)的响应判断稳定性====== A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B 和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号 x1(n)=R8(n) x2n=ones(1,128); %产生信号 x2(n)=u(n) y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) n=0:length(y1n)-1; subplot(2,2,1);stem(n,y1n,'.'); title('(a) 系统对R8(n)的响应y1(n)');xlabel('n');ylabel('y1(n)'); y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) n=0:length(y2n)-1; subplot(2,2,2);stem(n,y2n,'.'); title('(b) 系统对u(n)的响应y2(n)');xlabel('n');ylabel('y2(n)'); hn=impz(B,A,58); %求系统单位脉冲响应h(n) n=0:length(hn)-1; subplot(2,2,3);y=hn;stem(n,hn,'.'); title('(c) 系统单位脉冲响应h(n)');xlabel('n');ylabel('h(n)'); 020406000.20.40.60.8(a) 系统对R8(n) 的响应y1(n)ny1(n)05010015000.51(b) 系统对u(n) 的响应y2(n)ny2(n)020406000.050.1(c) 系统单位 脉 冲 响应h(n)nh(n) %===内容2:调用conv 函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号 x1(n)=R8(n) h1n=ones(1,10); h2n=[1 2.5 2.5 1 ]; y21n=conv(h1n,x1n); y22n=conv(h2n,x1n); figure(2) n=0:length(h1n)-1; subplot(2,2,1);stem(n,h1n); title('(d) 系统单位脉冲响应h1n');xlabel('n');ylabel('h1(n)'); n=0:length(y21n)-1; 2 su bplot(2,2,2);stem(n,y 21n); title('(e) h1(n)与 R8(n)的卷积y 21n');x label('n');y label('y 21(n)'); n=0:length(h2n)-1; su bplot(2,2,3);stem(n,h2n); title('(f) 系统单位脉冲响应h2n');x label('n');y label('h2(n)'); n=0:length(y 22n)-1; su bplot(2,2,4);stem(n,y 22n); title('(g) h2(n)与 R8(n)的卷积y 22n');x label('n');y label('y 22(n)'); 051000.51(d) 系统单位脉冲响应h1nnh1(n)0510152002468(e) h1(n)与 R8(n)的 卷 积 y21nny21(n)01230123(f) 系统单位脉冲响应h2nnh2(n)051002468(g) h2(n)与 R8(n)的 卷 积 y22nny22(n) %=========内容 3:谐振器分析======================== u n=ones(1,256); %产生信号 u (n) n=0:255; x sin=sin(0.014*n)+sin(0.4*n); %产生正弦...