matlab数值信号处理

x=ones(100);t=1:100;%x(n)=1~100,t用于后面绘图；
b=[.001836,.007344,.011016,.007374,.001836];%形成变量b;
a=[1,-3.0544,3.8291,-2.2925,.55075];%形成变量a;
y=filter(b,a,x);%实现式(2.8.3);
subplot(221)
plot(t,x,'g.',t,y,'k-');%将x(n)(绿色)y(n)黑色画在同一图上;
[h,t]=impz(b,a,40),title('单位抽样响应');
subplot(222)
stem(t,h,'.');grid on;
[H,w]=freqz(b,a,256,'whole',1);
Hr=abs(H);
Hphase=angle(H);Hphase=unwrap(Hphase),title('跃阶响应');
subplot(223)
plot(w,Hr),title('幅频响应');
subplot(224)
plot(w,Hphase),title('相频响应');