stepdown.m

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

function [aa,rho0,rho] = stepdown(ip,a,sig2)
      
for i=1:ip
	aa(i,ip)=a(i);
end
rho(ip)=sig2;
ip2=ip-1;
for j=1:ip2
	k=ip+1-j;
	den=1.-abs(aa(k,k))^2;
	rho(k-1)=rho(k)/den;
	k1=k-1;
	for i=1:k1
		aa(i,k-1)=(aa(i,k)-aa(k,k)*conj(aa(k-i,k)))/den;
	end
end
rho0=rho(1)/(1.-abs(aa(1,1))^2);