casfiltr.m

新版Matlab 7辅助信号处理技术与应用

function y=casfiltr(b0,B,A,x);
%滤波器的级联实现
%[y]=casfiltr(b0,B,A,x);
%b0=增益系数
%B=包含各因子系数bk的K行3列矩阵
%A=包含各因子系数ak的K行3列矩阵
%x=输入序列

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