%music算法以自相關矩陣分析為基礎，將系統自相關矩陣分為兩個子空間：信號子空間和噪聲子空間。
%為求功率譜估計，music算法計算信號空間和噪聲空間的
%特征值向量函數，使得在正弦信號頻率處功率譜現峰值，而其他地方函數最小。
例子：用music法估計信號x(t)=sin(2pif1t)+2sin(2pif2t)+w(t)的功率譜。其中f1=50Hz，
f2=120Hz,w(t)為白噪聲，采樣頻率Fs=1000Hz.
程序如下：
N=1000;
n=0:N-1;
Fs=1024;
t=n/Fs;
xn=sin(2*pi*100*t)+2*sin(2*pi*200*t)＋randn(1,N);
subplot(211);
pmusic(xn,7)