music算法

MUSIC算法是空間譜估計測向理論的重要基石。算法原理 ^[2]  如下：

（1） 不管測向天線陣列形狀如何，也不管入射來波入射角的維數如何，假定陣列由M個陣元組成，則陣列輸出模型的矩陣形式都可以表示為：Y(t)=AX(t)+N(t)

陣列方向A矩陣(2張)

其中，Y是觀測到的陣列輸出數據復向量；X是未知的空間信號復向量；N是陣列輸出向量中的加性噪聲；A是陣列的方向矩陣；此處，A矩陣表達式由圖冊表示。

MUSIC算法的處理任務就是設法估計出入射到陣列的空間信號的個數D以及空間信號源的強度及其來波方向。

（2） 在實際處理中，Y得到的數據是有限時間段內的有限次數的樣本（也稱快拍或快攝），在這段時間內，假定來波方向不發生變化，且噪聲為與信號不相關的白噪聲，則定義陣列輸出信號的二階矩：Ry。

Ry矩陣運算及特徵值分解(3張)

（3） MUSIC算法的核心就是對Ry進行特徵值分解，利用特徵向量構建兩個正交的子空間，即信號子空間和噪聲子空間。對Ry進行特徵分解，即是使得圖冊中的公式成立。

（4） U是非負定的厄米特矩陣，所以特徵分解得到的特徵值均為非負實數，有D個大的特徵值和M-D個小的特徵值，大特徵值對應的特徵向量組成的空間Us為信號子空間，小特徵值對應的特徵向量組成的空間Un為噪聲子空間。

（5） 將噪聲特徵向量作為列向量，組成噪聲特徵矩陣 ，並張成M-D維的噪聲子空間Un，噪聲子空間與信號子空間正交。而Us的列空間向量恰與信號子空間重合，所以Us的列向量與噪聲子空間也是正交的，由此，可以構造空間譜函數。

（6） 在空間譜域求取譜函數最大值，其譜峯對應的角度即是來波方向角的估計值。

MUSIC（Multiple Signal Classification多信號分類）算法是1979年由美國人R.O.Schmidt提出的，它標誌着空間譜估計測向進入了繁榮發展的階段。它將“向量空間”的概念引入了空間譜估計領域，經過三十年的發展，可以説其理論已經比較成熟。

自80年代以來，人們對基於特徵分解的超分辨率空間譜估計算法進行了廣泛深入的研究，並提出了一系列高效的處理方法，其中最經典的是多信號分類（MUSIC）算法，這種算法要經過一維搜索才能求出信源的來向，而相對最大似然（ML）和加權子空間擬合（WSF）等多維搜索算法的運算量已經減少了很多。以MUSIC為代表的算法存在一個缺點，即對相干信號處理的不理想。在針對相干信號源的一系列處理方案中，比較經典的是空間平滑技術，如空間平滑（SS）和修正的空間平滑（MSS）算法。然而，空間平滑技術是以損失陣列有效孔徑為代價的，而且只適用於等距均勻線陣（ULA）。

事實上空間譜估計算法都是在已知信號源數目下計算的，而在實際應用中這是不可能的，只能根據觀測數據對源數目進行估計。R.O.Schmidt在他的經典之作中提出了依據陣列協方差矩陣特徵值的分佈來估計信號源的方法。這種方法在理論上是完美的，至少對獨立源和部分相關源是正確的，但實際上由於數據長度有限，很大程度上只能依靠主觀判斷來確定源數。