czeroflt.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 18 行

function [y] = czeroflt(x,n,ip,a) 
 
y(1)=x(1);
if ip ~= 1
	for i=2:ip
		y(i)=x(i);
		for j=1:i-1
			y(i)=y(i)+a(j)*x(i-j);
		end
	end
end
for i=ip+1:n
	y(i)=x(i);
	for j=1:ip
		y(i)=y(i)+a(j)*x(i-j);
	end
end