经典谱估计BT法和周期图法 matlab及实验报告

%BT法求功率谱
N=1024;%采样点数
n=0:N-1;
fs=20000;%采样频率
fz1=200;%正弦信号频率
fz2=400;
f1=fz1./fs;
f2=fz2./fs;
x=sin（2*pi*f1*n）+3*sin（2*pi*f2*n）;%正弦信号
figure（1）
subplot（221）
plot（x）; 
xlabel（‘n‘）;
ylabel（‘幅值‘）;
title（‘原始正弦信号‘）;
u=randn（1N）;
y=x+u   %y为x加白噪声的信号
subplot（222）
plot（y）;       
title（‘混有高斯白噪声的正弦图形‘）;
xlabel（‘n‘）;
ylabel（‘幅值‘）;
z=xcorr（y‘unbiased‘）;  %对y做自相关
subplot（223）;
plot（z）;       
title（‘加白噪声后的信号一次自相关‘）;
subplot（224）
Z=fft（zN）;
Zmax=max（Z）;
Z=Z/Zmax;%幅度归一化
f=（0:N-1）*fs/N;
plot（f20*log（abs（Z）））;
title（‘BT法求 功率谱‘）;
%周期图法求功率谱
figure（2）
subplot（221）
plot（x）; 
xlabel（‘n‘）;
ylabel（‘幅值‘）;
title（‘原始正弦信号‘）;
su