chenmh_lpc.m

来自「Frequently used algorithms for numerical」· M 代码 · 共 7 行

function [A,R,E] = Chenmh_LPC(s,N);
R = conv(s,s(end:-1:1)); 
R = R(fix(length(R)/2)+1:-1:1)/length(s);
% [A,E] = levinson(R(1:N+1),N);
[A,E] = Chenmh_Levinson(R(1:N+1),N);