chenmh_lpc.m

Frequently used algorithms for numerical analysis and signal processing

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);