ms5_2.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 14 行

N=16;                 
n=0:N-1;
x=[ones(1,N/4),zeros(1,3*N/4)];       %产生输入信号序列        
subplot(2,2,1);stem(n,x,'filled','k');            
title('x(n)');
a=[1,0,1/3,0];
b=[1/6,1/2,1/2,1/6];
hn=impz(b,a,n);                       %求系统的单位冲激响应
subplot(2,2,2);stem(n,hn,'filled','k');
title('h(n)');
y=conv(x,hn);                         %卷积运算
subplot(2,1,2);stem(y,'filled','k');
title('y(n)');