coeff.m

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

function [p,q,r,c,d,e] = coeff(i,ip,s,n,aa,p,q,r,c,d,e) 
      
if i ~= 1
	i1=i-1;
	i2=i+1;
	ap(1)=1.;
	for j=1:i1 
		ap(j+1)=aa(j,i1);
	end
	for j=1:i
		b(j)=ap(i+1-j);
	end
	c=0.;
	d=0.;
	for l=1:i 
		for k=1:i 
			c=c+ap(k)*b(l)*s(k,l+1);
			d=d+b(k)*b(l)*s(k+1,l+1);
		end
	end
else
	c=s(1,2);
	d=s(2,2);
	e=s(1,1);
end
p=((n-2.*i)/(n-i))*(c/d);
q=(i*e+n*d)/(d*(i-n));
r=n*c/(d*(i-n));