example5_2.m

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

b=[1    -3    11    -27    18];
a=[16    12    2    -4    -1];
[bo,B,A]=dir2cas(b,a)%直接型转化为级联

delta=impseq(0,0,7)
format long
hcas=casfiltr(bo,B,A,delta)%级联的输出
hdir=filter(b,a,delta)%直接型的输出