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);