limmusic.m

功率谱估计程序

%Music算法。选用电力系统谐波信号作为输入信号。参考书目：近代信号处理理论与方法，张明友等著，国防工业出版社。131-132页
clear all;close all;
f1=50;ts=1/1000;n=1:256;N=256;fs=1000;
x=sin(2*pi*f1*n*ts)+0.1*sin(2*pi*350*n*ts)+0.1*sin(2*pi*360*n*ts);D=0;p=30;
y=x';
noise=randn(1,256);
noise=noise';
dB=30;
y=sigmerge(y,noise,dB);
x=y';
for m=0:N-1
    s=0;
    for n=0:N-1-m
        s=s+x(n+1)*x(n+m+1)/N;
    end
    r(m+1)=s;
end
r=toeplitz(r);
[a,b,c]=svd(r);
for f=0:fs/2
    for n=0:N-1
        e(n+1)=exp(i*2*pi*f*n*ts);
    end

    s=0;
    for k=p+1:N
    s=s+(abs(e*c(:,k)))^2;
    end
    P(f+1)=1/s;
end
f=0:fs/2;
semilogy(f,P);