%%%%%%%%%%%%%%
%??task?:?cal?power?spectrum
%??
%%%%%%%%%%%%%%
%%
clear?all;close?all;clc;
Fs?=?1000;
N??=?1024;????????????%?length?of?signals?
Nfft=?1024;
t??=?（0:N-1）/Fs;
xn?=?sin（2*pi*50*t）?+?sin（2*pi*150*t）?+?randn（1N）;
Pxx=?10*log10（?abs（fft（xnNfft））.^2/N?）;??%?power?spectrom
figure;
subplot（211）;
plot（xn）;
xlabel（‘序號‘）;ylabel（‘離散值‘）;title（‘原始信號‘）;
f?=?Fs/Nfft*[0:Nfft/2-1];
subplot（212）;
plot（fPxx（1:Nfft/2））;
xlabel（‘頻率/Hz‘）;ylabel（‘功率譜/dB‘）;title（‘功率譜‘）;
grid?on;
%%
%??上面的程序也可以將xn分成四段，每段256點fft，然后求和平均也可以；即所謂的平均分段法。
%??下面利用