psd_zxfc.m

来自「周期图法进行功率谱估计的函数」· M 代码 · 共 29 行

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

function S=PSD_zxfc(x,p,M)     
N=length(x);
pi=3.1415926;
xf=x(p+1:N)';
A_H=zeros(p,N-p);
S0 = zeros(1,M+1);
for ii=1:p
    for jj=1:N-p
        A_H(ii,jj)=x(p-ii+jj);
    end
end

% AR模型参数估计
a=-(A_H*A_H')^-1*A_H*xf;
Pv=1/(N-p)*(xf'*xf+a'*A_H*xf);

% 功率谱估计
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;