arupt.m

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

function [aa,e,pow] = arupt(aa,i,s,n,c,d,e) 
      
if i ~= 1
	i1=i-1; 
	for j=1:i1 
		aa(j,i)=aa(j,i1)+aa(i,i)*aa(i-j,i1);
	end
end
e=e+2.*c*aa(i,i)+d*aa(i,i)^2;
pow=e/n;