ex6p10.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 9 行

b=[1,2,2,1];
a=[1,13/24,5/8,1/3];
[K,C]=dir21adr(b,a)
[x,n]=impseq(0,0,7)
format long
hdirect=filter(b,a,x)
hladder=ladrfilt(K,C,x)