%Samp9_13
clf;Fs=1000;   %采樣頻率
N=1024;Nfft=256;n=0:N-1;t=n/Fs;   %數據長度、分段數據和時間序列
randn('state',0)  %設置產生隨機數的初始狀態
xn=sin(2*pi*100*t)+2*sin(2*pi*200*t)+randn(1,N); %原始信號
pmusic(xn,[7,1.1],Nfft,Fs,32,16);   %采用多信號分類法估計功率譜
xlabel('頻率/Hz'); ylabel('功率譜/dB')
title('通過MUSIC法估計的偽譜')