ms5_4.m

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

nx=-2:10*pi;
x=square(nx/5);       %产生输入信号序列        
subplot(3,1,1);stem(nx,x,'filled','k');            
ylabel('x(n)');title('y(n)=x(n)*h(n)');
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,ny]=convnew(x,nx,hn,nh);                         %卷积运算
subplot(3,1,3);stem(ny,y,'filled','k');
ylabel('y(n)');