Burg算法的MATLAB实现

clc
clear all
clf

N=128;  
n=1:N;       
w（n）=randn（1N）;     %产生零均值方差为1的高斯白噪声
xn=sqrt（20）*sin（2*pi*0.2*n）+sqrt（2）*sin（2*pi*0.213*n）+w（n）;
a=zeros（N-1N-1）;
p0=sum（xn.^2）/N;
f0=xn;   %%%
g0=xn;   %%%
Km=0;
P=zeros（1N-1）;
f=zeros（N-1N-1）;
g=zeros（N-1N-1）;
zi=0;
mu=0;
% w=0:2*pi;
m=1;
for n=m+1:N
            zi=zi+f0（n）*g0（n-1）;
            mu=mu+f0（n）*f0（n）+g0（n-1）*g0（n-1）;
end
        Km=-2*zi/mu;
%         a（10）=1;
        a（m1）=Km;
        P（m）=（1-Km^2）*p0;      %%%
        for n=2:N
        f（1n）=f0（n）+Km*g0（n-1）;
        g（1n）=conj（Km）*f0（n）+g0（n-1）;
        end

fo