example5_3.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

b=[1    -3   11    -27   18];
a=[16   12    2    -4    -1];
[C,B,A]=dir2par(b,a)
format long;
delta=impseq(0,0,7);
hpaar=parfiltr(C,B,A,delta)
hdir=filter(b,a,delta)