latcfilt.m

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

function [y]=latcfilt(K,x)
%格型滤波器的实现
%y为输出序列；K为格型滤波器反射系数数组；x输入序列
Nx=length(x)-1;
x=K(1)*x;
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,:);