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('振幅'); grid on n=-10:10; %时间序号n x=5*sin(0.1*pi*n); %产生随时间序列而变化的信号 stem(n,x) title('正弦信号序列') xlabel('时间序号n') ylabel('振幅') axis([-10 10 -5.5 5.5]) grid on