cpoleflt.m

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

function [y] = cpoleflt(x,n,b,iq,y)  
      
y(1)=x(1);
if iq ~= 1
      	for i=2:iq  
      		y(i)=x(i);
      		for j=1:i-1  
			y(i)=y(i)-b(j)*y(i-j);
		end
	end
end
for i=iq+1:n  
      	y(i)=x(i);
      	for j=1:iq  
		y(i)=y(i)-b(j)*y(i-j);
	end
end