最近很糾結於burg算法，levinse算法還有shur算法，這些用來求線性預測係數的方法。
Burg算法 2009-01-06 22:24 參考張賢達現代信號處理第二版編寫 改正了昨天程序中的小bug function a = burg(x,ORDER) % MADE BY H.X.GUAN,2008/12/18 % BURG ALGORITHEM,AR model %%%%% Initial %%%%%% N = length(x); f = x; % 前向 g = x; % 後向 % 假設x爲行向量 P0 = x*x'/N; PM = P0; m = 2; if ORDER>N disp('階數超過向量長度無法計算'); a = 0; return; end a = zeros(1,ORDER+1); am = a; a(1) = 1; am(1) = 1; %%%%%%%%%%%%%%%%%%%% while ( m <= ORDER+1 ) sum1 = 0; sum2 = 0; for n=m:N sum1 = sum1 + f(n)*g(n-1); sum2 = sum2 + f(n)*f(n) + g(n-1)*g(n-1); end Km = -1.0*sum1 / (sum2*0.5); for i=2:m am(i) = a(i) + Km*a(m-i+1); end % am(m) = Km; for i=2:m a(i) = am(i); end PM = (1-Km*Km)*PM; for n=2:N fm(n) = f(n) + Km*g(n-1); gm(n) = Km*f(n) + g(n-1); end f = fm; g = gm; m = m + 1; end