cas2dir.m

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

function y=cas2dir(b0,B,A);
[K,L]=size(B);
b=[1];
a=[1];
for i=1:1:K
b=conv(b,B(i,:));
a=conv(a,A(i,:));
end
b=b*b0;