%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