ex6p3.m

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

C0=0;B1=[2 4;3  1];A1=[1 1 0.9; 1 0.4 -0.4];
B2=[0.5 0.7;1.5  2.5;0.8  1];
A2=[1  -1  0.8;1  0.5  0.5;1 0 -0.5];
[b1,a1]=par2dir(C0,B1,A1)
[b2,a2]=par2dir(C0,B2,A2);
b=conv(b1,b2)
a=conv(a1,a2)
[b0,Bc,Ac]=dir2cas(b,a);
[C0,Bp,Ap]=dir2par(b,a);