parfiltr.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 12 行

function y=parfiltr(C,B,A,x)
%IIR滤波器的并联实现；
%C为当M>=N时的多项式部分；B级联型传输函数的分子；A级联型传输函数的分母；
%y为输出序列；x为输入序列
[K,L]=size(B);
N=length(x);
w=zeros(K+1,N);
w(1,:)=filter(C,1,x);
for i=1:1:K
    w(i+1,:)=filter(B(i,:),A(i,:),x);
end
y=sum(w);