polefilt.m

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

function [y] = polefilt(x,n,iq,b) 
      
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