exp2

a1=[1,1.5,0.9]; %设置filter函数的A向量
b1=[1,-0.8]; n=0:100;   %设置filter函数的B向量
x1=[1 zeros(1,100)];    %输入单位样本
y1filter=filter(b1,a1,x1);  %y[n]+1.5*y[n-1]+0.9*y[n-2]=x[n]-0.8*x[n-1]
stem(n,y1filter);
title('信号y'); 
xlabel('时间序号n'); 
ylabel('振幅');