MUSIC演算法

MUSIC（英語：MUltiple SIgnal Classification，全稱：多訊號分類）是一種用於頻率估計^[1]和無線電測向^[2] 的算法。

發展歷史

現實生活中有許多根據測量結果估計接收訊號所依賴的一組常數參數訊號處理問題。已經有幾種解決這些問題的方法，包括Capon（1969）的最大似然（ML）方法和Burg的最大熵（ME）方法。雖然這些方法經常成功並且被廣泛使用，但是它們具有某些基本限制（尤其是參數估計中的偏差和靈敏度），主要是因為它們使用測量的不正確模型（例如自迴歸模型而不是特殊ARMA模型）。

Pisarenko（1973）是最早利用數據模型結構的人之一，在使用協方差方法估計加性噪聲中復雜平面波的參數的情況下這樣做。 Schmidt（1977）在Northrop Grumman工作時獨立（1979）是第一個在任意形式的傳感器陣列情況下正確利用測量模型的人。特別是施密特通過首先在沒有噪聲的情況下推導出完整的幾何解決方案來實現這一目標，然後巧妙地擴展幾何概念以在存在噪聲的情況下獲得合理的近似解。得到的算法稱為MUSIC（多重訊號分類）並且已被廣泛研究。

在基於數千次類比的詳細評估中，麻省理工學院的林肯實驗室得出結論，在目前公認的高解析度算法中，MUSIC是最有前途的，也是進一步研究和實際硬體實現的主要候選者^[3]。然而，雖然MUSIC的性能優勢很大，但它們是以計算（搜索參數空間）和存儲（陣列校準數據）為代價的。^[4]

頻率估算的應用

MUSIC使用特徵值和特徵向量方法估計訊號或Autocorrelation matrix的頻率內容。該方法假設訊號${\displaystyle x(n)}$由存在高斯白噪聲的 ${\displaystyle p}$複指數組成。給定${\displaystyle M\times M}$次M自相關矩陣，${\displaystyle \mathbf {R} _{x}}$，如果特徵值按降序排序，對應於 ${\displaystyle p}$個最大特徵值（即最大可變性的方向）的特徵向量跨越訊號子空間。其餘的 ${\displaystyle M-p}$ 特徵向量跨越正交空間，其中只有噪聲。注意，對於 ${\displaystyle M=p+1}$，MUSIC與Pisarenko harmonic decomposition相同。一般的想法是使用平均來改善Pisarenko估計的性能。

MUSIC的頻率估算功能是

${\displaystyle {\hat {P}}_{MU}(e^{j\omega })={\frac {1}{\sum _{i=p+1}^{M}|\mathbf {e} ^{H}\mathbf {v} _{i}|^{2}}},}$

其中， ${\displaystyle \mathbf {v} _{i}}$ 是噪聲的特徵向量，

${\displaystyle \mathbf {e} ={\begin{bmatrix}1&e^{j\omega }&e^{j2\omega }&\cdots &e^{j(M-1)\omega }\end{bmatrix}}^{T}}$

被稱為導向矢量，在這種情況下表示為均勻陣列。

估計函數的 ${\displaystyle p}$個最大峰值的位置給出了 ${\displaystyle p}$訊號分量的頻率估計。

MUSIC是Pisarenko harmonic decomposition的推廣和計算。在該方法中，僅使用單個特徵向量並將其視為一組自回歸係數，其零值可以通過分析或多項式根尋找算法找到。相比之下，MUSIC假設已將幾個此類函數添加在一起，因此可能不存在零。相反地，存在可利用計算，搜尋估計函數的局部最小值。

與其他方法比較

當預先知道組件的數量時，MUSIC優於簡單的方法，例如:存在噪聲的情況下選取DFT譜的峰值，因為它利用該數字的知識來忽略其最終報告中的噪聲。

與DFT不同，它能夠以高於一個樣本的精度估計頻率，因為可以針對任何頻率評估其估計函數，而不僅僅是DFT頻段的頻率。這是一種超解析度成像的形式。

它的主要缺點是它需要事先知道組件的數量，因此原始方法不能用於更一般的情況。存在用於僅從自相關矩陣的統計特性估計源分量的數量的方法。例如， ^[5] 此外，MUSIC假設共存源不相關，這限制了其實際應用。

最近的迭代半參數方法提供了強大的超解析度成像，儘管有高度相關的來源，例如:SAMV (algorithm) ^[6][7]

其他應用

MUSIC的修改版本，表示為"時間反轉-MUSIC"（TR-MUSIC），最近已應用於計算時間反轉成像。^[8][9]。 MUSIC算法也被以 C 語言庫的形式實現用於以快速檢測DTMF頻率（雙音多頻信令） - libmusic ^[10]。

參見

• 譜密度
• 週期圖法
• Welch's method
• Bartlett's method
• SAMV演算法
• 無線電測向
• 基音檢測算法
• 超高解析度顯微鏡學