test.m

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

function [xkest] = test(root,c,d,e,i,n) 
      
fmax=-10000000.;
for j=1:3 
	if root(j) >= -1 & root(j) <= 1.
		q=e+2.*root(j)*c+root(j)*root(j)*d;
		f=(i/2.)*log10(1.-root(j)*root(j))-(n/2.)*log10(q/n);
		if f > fmax
			xkest=root(j);
			fmax=f;
		end
	end
end