ms6_1.m

《数字信号处理实验(MATLAB版)》这

nx=0:8*pi;
x=sin(nx/2)+sin(10*nx)/3;       %产生输入信号序列        
subplot(3,1,1);stem(nx,x,'filled','k');            
ylabel('x(n)');title('求系统的响应');
a=[1,-0.34319,0.60439,-0.20407];
b=[0.1321, 0.3963,0.3963,0.1321];
nh=0:9;
hn=impz(b,a,nh);        %求系统的单位冲激响应
subplot(3,1,2);stem(nh,hn,'filled','k');
ylabel('h(n)');
y=dlsim(b,a,x);         %求系统的响应
subplot(3,1,3);stem(y,'filled','k');
ylabel('y(n)');