example5_2.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 8 行

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)%直接型的输出