%%
%%信號(hào)源1bartlett法
clear?all;
close?all;
clc;

L=50;%數(shù)據(jù)分段段數(shù)
M=512;%每段采樣數(shù)據(jù)點(diǎn)數(shù)
N=L*M;%采樣點(diǎn)個(gè)數(shù)
MU=0;%噪聲均值
D=0.101;%噪聲方差
times=100
for?i=1:times
v=normrnd（MUsqrt（D）N1）;%隨機(jī)產(chǎn)生均值為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取每段數(shù)據(jù)的暫存矢量
????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（‘信號(hào)源1的估計(jì)功率譜圖‘）;
end
%%
%%信號(hào)源2
clear?all;
close?all;
clc;
N=1024;%采樣點(diǎn)數(shù)
M=32;%每段采樣數(shù)據(jù)點(diǎn)數(shù)
L=N/M;%數(shù)據(jù)段數(shù)；
MU=0;
D=1;
v=normrnd（MUsqrt（D）N1）;%隨機(jī)產(chǎn)生均值為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