limzixiangguan.m

功率谱估计程序

%自相关法，选用电力系统谐波信号作为输入信号。参考书目：现代数字信号处理，皇甫堪等，电子工业出版社。208页。
clc;clear all;close all;
f1=50;ts=1/1000;n=1:256;N=256;fs=1000;delta=0.101043;
x(1:128)=2*sin(2*pi*f1*(1:128)*ts);
x(129:256)=2*sin(2*pi*f1*(129:256)*ts)+0.1*sin(2*pi*350*(129:256)*ts)+0.2*sin(2*pi*400*(129:256)*ts);
p=10
for k=0:p
    r1=0;
    for n=0:N-k-1
        r1=r1+1/N*x(n+1)*x(n+k+1);
    end
    r(k+1)=r1;
end
R=toeplitz(r(1:p));
r2=r(2:p+1)';
a=-inv(R)*r2;

for f=0:fs/2
    s=0;
    for k=1:p
        s=s+a(k)*exp(-i*2*pi*f*k/fs);
    end
    P(f+1)=delta/(abs(1+s))^2;
end
f=0:fs/2;
semilogy(f,P)