功率谱计算matlab程序

%%%%%%%%%%%%%%
%  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，然后求和平均也可以；即所谓的平均分段法。
%  下面利用