parfiltr.m

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

function y=parfiltr(C,B,A,x);
%滤波器并联实现
%[Y]=parfiltr(C,B,A,x);
%C=当分子多项式阶数大于分母多项式阶数时产生的多项式
%B=K列3行bk系数矩阵
%A=K列3行ak系数矩阵
%x=输入序列

%初始化矩阵
 [K,L]=size(B);
N=length(x);
w=zeros(K+1,N);
w(1,:)=filter(C,1,x);
for i=1:K
%对于每个子系统而言，与直接型机构计算一样
   w(i+1,:)=filter(B(i,:),A(i,:),x);
end
y=sum(w);