cas2dir.m

来自「<Matlab 7辅助信号处理技术与应用》的源码」· M 代码 · 共 19 行

function [b,a]=cas2dir(b0,B,A);
%变级联形式为直接形式
%[b,a]=dir2cas(b0,B,A)
%b0=增益系数
%B=包含各因子系数bk的K行3列矩阵
%A=包含各因子系数ak的K行3列矩阵
%a=直接型分子多项式系数
%b=直接型分母多项式系数
%初始化矩阵
[K,L]=size(B);
b=[1];
a=[1];
%将多项式项乘所得的系数即为所求
for i=1:K
   b=conv(b,B(i,:));
   a=conv(a,A(i,:));
end
b=b*b0;