latcfilt_m.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 13 行

function [y]=latcfiltm(K,x);
%y=输出序列
%K=格型结构滤波器的系数矩阵
%x=输入序列
Nx=length(x)-1;
M=length(K)-1;K=K(2:M+1);
fg=[x;[0 x(1:Nx)]];
for m=1:M
   fg=[1,K(m);K(m),1]*fg;
   fg(2,:)=[0 fg(2,1:Nx)];
end
y=fg(1,:);