psd_burg.m

用Burg法进行功率谱预测的函数

% 离散随机信号处理 Project I
% 功率谱估计
% 冯伟  2005310367  [e-mail]
% 功能: Burg法进行谱估计   2007年5月18日

function S=PSD_Burg(x,p,M)     
N=length(x);
pi=3.1415926;
S0 = zeros(1,M+1);

%初始值
P1=1/N*sum(x*x');
f1=x(2:N);
b1=x(1:N-1);
a1=zeros(1,p);

% 递推---AR模型参数估计
for m = 1:p
    f0 = f1;
    b0 = b1;
    a0 = a1;
    P0 = P1;
    km_up = 0;
    km_dn = 0;
    for nn = m:N-1
        km_up = km_up - 2*conj(f0(nn))*b0(nn);
        km_dn = km_dn + abs(f0(nn))^2 + abs(b0(nn))^2;
    end
    km = km_up/km_dn;
    P1 = (1-abs(km)^2)*P0;
    if m>1
        for jj = 1:m-1
            a1(jj) = a0(jj) + km*conj(a0(m-jj));
        end
    end    
    a1(m) = km;
    f1(m+1:N-1) = f0(m+1:N-1)+conj(km)*b0(m+1:N-1);
    b1(m+1:N-1) = b0(m:N-2)+km*f0(m:N-2);
end
a=a1;
Pv=P1;

% 功率谱估计
for n=1:M+1
    k=1:p;
    S0(n)=10*log10(Pv/abs(1+a*exp(-j*pi/M*(n-1)*k).')^2);
end
S = S0;