clear all
close all
clc
x=ones(1,100);
t=1:100;
b=[0.00183,0.00734,0.01101,0.00737,0.00183];
a=[1,-3.054,3.825,-2.292,0.551];
[C,B,A]=dir2par(b,a)
[b,a]=par2dir(C,B,A)
y=parfiltr(C,B,A,x);
%y=filter(b,a,x);
plot(t,x,'r',t,y,'k-');
grid on;
ylabel('x(n)與y(n)');
xlabel('n');