陣列信號處理—music、capon

1、平面陣中幾種二維DOA估計算法,二維DOA估計算法,接收信號模型二維MUSIC算法二維求根MUSIC算法二維Capon算法二維ESPRIT及二維Unitary-ESPRIT算法二維傳播算子算法PARAFAC技術(shù)各種算法比較,接收信號模型,接收信號模型,平面陣列示意圖,設(shè)平面陣元數(shù)為M×N，信源數(shù)為K。信源的波達(dá)方向為 。第i個陣元與參考陣元之間的波程差為,接收信號模型,設(shè)子陣1沿x軸

2、的方向矩陣為 ，而子陣2的每個陣元相對于參考陣元的波程差就等于子陣1的陣元的波程差加上 ，所以接收信號為協(xié)方差矩陣為,平面陣列示意圖,接收信號模型,其中， 代表有最大的K個特征值構(gòu)成的一個K×K對角陣, 代表有MN-k個較小的特征值構(gòu)成的對角矩陣, 和 分別代表由 和 對應(yīng)的特征值構(gòu)成的特征矢量。,沿x軸的方向矩陣,沿y軸的方向矩陣,二維MUSIC算法,二維MUSIC算法,一維MUS

3、IC算法 陣列協(xié)方差矩陣通過奇異分解，可以劃分為兩個空間，即 ，滿足性質(zhì)：這說明矩陣 中的各個列向量與噪聲子空間正交所以： ，根據(jù)這個性質(zhì)，得到陣列空間譜函數(shù) ，變化角度，通過尋找波峰來估計源信號的角度。,二維MUSIC算法,二維MUSCI算法 類似于一維的情況，此時方向矩陣變?yōu)?所以空間譜函數(shù)為：其中，

4、 分別是 的第k個列向量。,二維MUSIC算法,仿真參數(shù)設(shè)置：快拍數(shù)L=100，目標(biāo)數(shù)K=3，8×8的方陣，假設(shè)源信號的仰角為10°，20°，30°，方位角為15°，25°，35°，信噪比為20dB。仿真結(jié)果如下圖所示：,二維求根MUSIC算法,二維求根MUSIC算法,一維求根MUSIC算法,顧名思義，求根MUSIC方法是MUSI

5、C方法的一種多項式求根形式，我們可以構(gòu)造求根MUSIC的多項式：,式中 ， ，通過求上述多項式的幅度最大的K個根 ，就可以計算出波達(dá)方向估計：,二維求根MUSIC算法,將平面陣的輸出信號進(jìn)行重構(gòu)：,上式可以表示為 ，相當(dāng)于一維DOA估計形式，利用求根MUSIC算法可以求得K個根對應(yīng)的矢

6、量：,進(jìn)一步可求得：,二維求根MUSIC算法,二維求根MUSIC算法,配對之后，可以得到：,式中： 是有噪聲情況下的 估計值,同理，將矩陣看成沿y軸方向的，可以得到：,二維求根MUSIC算法,仿真參數(shù)設(shè)置：快拍數(shù)L=100，目標(biāo)數(shù)K=2，8×8的方陣，假設(shè)源信號的仰角為10°，20°，30°，方位角為15°，25°，35°，信噪比為20dB。仿

7、真結(jié)果如下圖所示：,二維Capon算法,二維Capon算法,一維Capon的空間譜為：,擴展到二維情況：,類似于二維MUSIC算法，我們可以進(jìn)行二維譜峰搜索，尋找譜峰值，即可得到波達(dá)方向。,二維Capon算法,仿真參數(shù)設(shè)置：快拍數(shù)L=100，目標(biāo)數(shù)K=3，8×8的方陣，假設(shè)源信號的仰角為10°，20°，30°，方位角為15°，25°，35°，信噪比為20dB。仿真結(jié)果