ms5r4.m

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

N=10;n=0:N-1;
x=[ones(1,5),zeros(1,N-5)];
a=[1,0,1/3,0];
b=[1/6,1/2,1/2,1/6];
h=impz(b,a,n);
[y,ny]=convnew(x,n,h,n);  %卷积运算
subplot(2,2,1);stem(n,x,'filled','k');            
title('f1(n)');
subplot(2,2,2);stem(n,h,'filled','k');
title('f2(n)');
subplot(2,1,2);stem(ny,y,'filled','k');
title('y(n)');