parfiltr.m

来自「数字信号处理的MATLAB仿真加实验指导」· M 代码 · 共 10 行

function y=parfiltr(C,B,A,x);   %IIR并联型实现
[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);