平均周期图 谱估计 matlab版

%%
%%信号源1bartlett法
clear all;
close all;
clc;

L=50;%数据分段段数
M=512;%每段采样数据点数
N=L*M;%采样点个数
MU=0;%噪声均值
D=0.101;%噪声方差
times=100
for i=1:times
v=normrnd（MUsqrt（D）N1）;%随机产生均值为MU方差为D的高斯白噪声
fs=1;%采样频率为1hz
t=1/fs;%采样周期
f1=0.05;
f2=0.40;
f3=0.42;

z=0;
a1=-0.85;

for i=1:N
    z=-1*a1*z+v（i）;
    x（i）=cos（2*pi*f1*i）+cos（2*pi*f2*i）+cos（2*pi*f3*i）+0.5*z;
end
P=0;
for i=1:L
    for j=1:M
    x1（j）=x（j+（i-1）*L）;%X1取每段数据的暂存矢量
    end

 X=fft（x1）;
P=（X.*conj（X））/M+P;
end
P=P/L;
P=10*log10（P）;
n=0:（M/2-1）;
f=n/M;
subplot（222）;
plot（fP（1:M/2））;
hold on;
xlabel（‘频率/hz‘）;
ylabel（‘功率db‘）;
title（‘信号源1的估计功率谱图‘）;
end
%%
%%信号源2
clear all;
close all;
clc;
N=1024;%采样点数
M=32;%每段采样数据点数
L=N/M;%数据段数；
MU=0;
D=1;
v=normrnd（MUsqrt（D）N1）;%随机产生均值为MU方差为D的高斯白噪声
fs=1;%采样频率为1hz
t=1/fs;%采样周期
% f1=0.05;
% f2=0.40;
% f3=0.42;

a1=-1.35;
a2=1.34;
a3=-0.66;
a4=0.24;
x=zeros（N+41）;

for i=5:N+4
x（i）=-a1*x（i-1）-a