%  exa3-1_filter.m, for example 3-1
% to test filter.m and to obtain the step response.

clear;
x=ones(1,100);
t=1:100;
b=[0.00183,0.00734,0.01101,0.00737,0.00183];
a=[1,-3.054,3.829,-2.292,.551];
% 
y=filter(b,a,x);
% 求所給系統的輸出，本例實際上是求所給系統的階躍響應；
plot(t,x,'r.',t,y,'k-');
grid on;
ylabel('x(n) 與 y(n)')
xlabel('n')