%Samp9_8
clf;Fs=1000;   %采樣頻率
N=1024;Nfft=256;n=0:N-1;t=n/Fs; %數據長度、時間序列
window=hanning(256);   %選用的窗口
noverlap=128;   %分段序列重疊的采樣點數（長度）
dflag='none';  %不做趨勢處理
randn('state',0);   %設置產生隨機數的初始狀態
xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);   %帶噪聲原始信號
Pxx=psd(xn,Nfft,Fs,window,noverlap,dflag);   %功率譜估計
f=(0:Nfft/2)*Fs/Nfft;   %求得對應的頻率向量
plot(f,10*log10(Pxx));  %繪制功率譜
xlabel('頻率/Hz');ylabel('功率譜/dB');
title('PSD—Welch方法'); grid on