本文首发在我的个人博客,http://www.weekreport.cn/archives/440,欢迎来踩。
1. M U S I C MUSIC MUSIC算法简介
M U S I C MUSIC MUSIC算法由R.O.Schmidt在1979年提出,该算法利用了信号的频率向量和噪声子空间的正交性构造空间谱函数,通过谱峰搜索来估计信号频率。
M U S I C MUSIC MUSIC算法需要构造扫描函数如下: P ^ M U S I C ( w ) = 1 a H ( w ) G G H a ( w ) = 1 ∑ i = K + 1 M ∣ a H ( w ) u i ∣ 2 , w ϵ [ − π , π ] \hat{P}_{MUSIC}\left( w \right ) = \frac{1}{a^{H}\left(w \right)GG^{H}a\left( w \right )} = \frac{1}{\sum_{i=K+1}^{M}\left |a^{H}\left( w \right )u_{i} \right |^{2}},\quad w\epsilon \left [ -\pi,\pi \right ] P^MUSIC(w)=aH(w)GGHa(w)1=∑i=K+1M∣aH(w)ui∣21,wϵ[−π,π]
上式中, M M M代表信号自相关矩阵的阶数, K K K代表频率源个数, a w = [ 1 e − j w ⋮ e − j ( M − 1 ) w ] a_{w}=\begin{bmatrix}1\\e^{-jw}\\ \vdots \\e^{-j\left(M-1 \right )w} \end{bmatrix} aw=⎣⎢⎢⎢⎡1e−jw⋮e−j(M−1)w⎦⎥⎥⎥⎤是信号频率向量; u i u_{i} ui是信号自相关矩阵的特征向量,可以看到我们只取后 ( M − K ) \left( M-K \right) (M