梅錦洲

（廣東盛路通信科技股份有限公司，廣東 佛山 528100）

0 引 言

共形天線是一種和物體的外形結(jié)構(gòu)保持一致的天線或天線陣。共形天線可方便地裝備于載體表面而不影響載體的空氣動力學性能，此外，共形天線可采用柱面陣列和錐面陣列等形式，實現(xiàn)了廣角掃描，并在掃描中保持天線的波束形狀和增益。共形天線可使互耦維持在一定范圍，并具有隱形特性[1]。共形天線在衛(wèi)星通信、指揮和控制系統(tǒng)、導彈遙測、武器信管及多普勒雷達等學科占據(jù)重要位置，因此建立了數(shù)學模型，并利用Matlab軟件對傳統(tǒng)算法和子控件類算法進行共形陣列天線DOA的仿真分析[2]。

1 共形陣列天線接收信號模型

為簡化共形陣列天線的分析，假設(shè)入射信號的特性為平穩(wěn)且各態(tài)歷經(jīng)，用時間平均代替統(tǒng)計平均；噪聲為白噪聲，方差為σ2n；忽略共形天線各陣元之間的互耦和極化因素；入射的信號數(shù)小于陣列天線陣元數(shù)，且陣列接收到的所有信號互不相關(guān)、波達方向互不相同；假設(shè)信號為窄帶信號，陣元接收的信號之間無包絡(luò)變化[3]。圓錐共形陣列天線的模型如圖1所示。

其中，每層圓陣中逆時針方向的陣元編號為N（N=0，1，2，…，n-1）；圓陣自上至下的層數(shù)為M（M=0，1，2，…，m-1）；相鄰圓陣之間的間距為d；錐頂角為α。各陣元用直角坐標系可表示為：

圖1 圓錐共形陣列天線示意圖

其中，RMN=Mdtanα。那么t時刻L個遠場信號Si入射到錐面共形陣列天線上時，陣列天線接收信號可表示為：

其中，x(t)為接收信號矢量，H為陣列方向圖矩陣，B為導向矢量矩陣，S(t)為來波信號矢量，N(t)為加性噪聲矢量矩陣。

式（2）是基于建立圓錐共形陣列天線的接收信號的模型，且不考慮極化和互耦的條件下，利用數(shù)學知識推導多個互不相干的遠場信號以不同的方向入射到共形天線各陣元的接收信號。圓錐共形陣列天線在俯仰面上通常使用和圓柱共形陣列天線相同的波束掃描方法，尖銳的圓錐共形陣列天線可為實體提供良好的空氣動力學性能、低雷達散射截面及寬角覆蓋。

2 MUSIC算法

多重信號分類（MUSIC）算法是1979年R.O.Schmidt等人提出。該算法的提出促進了特征結(jié)構(gòu)分類算法的發(fā)展[4]。

假設(shè)陣元方向圖為全向，則窄帶信號的DOA估計的數(shù)學模型為：

其中，x(t)為接收信號矢量，B(θ)為導向矢量矩陣，S(t)為來波信號矢量，N(t)為加性噪聲矢量矩陣。

接收信號的協(xié)方差矩陣為：

其中，B為導向矢量矩陣，RS為入射信號的協(xié)方差矩陣，BH為導向矢量的共軛轉(zhuǎn)置矩陣，σ2為噪聲方差，I為單位矩陣。

由于信號與噪聲是獨立的，所以接收信號協(xié)方差矩陣可分解為信號與噪聲兩個部分，對Rx進行分解可得到：

其中，Us是信號子空間，Un是噪聲子空間。

理想的條件下，噪聲子空間和信號子空間是相互正交的，所以有：

其中，bH(θ)為導向矢量的共軛轉(zhuǎn)置

MUSIC算法是基于上述性質(zhì)提出的，但是接收數(shù)據(jù)協(xié)方差矩陣是用有限采樣數(shù)據(jù)的最大似然估計值：

Music算法是對空間譜函數(shù)PMUSIC進行譜峰值搜索，得到最大值P對應(yīng)的θ值即為來波信號的DOA。具體算法步驟為：由陣列接收的樣本數(shù)據(jù)X(t)得到數(shù)據(jù)協(xié)方差矩陣；對進行特征分解；確定噪聲子空間和信號子空間；根據(jù)信號的參數(shù)范圍進行搜索MUSIC空間譜，找出譜峰；得到譜峰值時所對應(yīng)的角度θ，即為信號的波達方向。

3 MATLAB軟件仿真

模型1：均勻圓陣模型中，陣元M=20，空間中信號L1和L2分別從俯仰角θ1為30°、θ2為70°，方位角φ1為60°、φ2為150°入射，信號頻率f為5 GHz，噪聲方差σ2=0.1，均勻圓陣的半徑r為1.5λ。

在Matlab軟件中采用MUSIC算法對模型1進行DOA仿真，仿真結(jié)果如圖2所示。

圖2 模型1仿真結(jié)果

由圖2可知，MUSIC算法在模型1中完全能識別出波達方向。

模型2：錐面共形陣列天線模型中，陣元M=25，空間中信號L1、L2及L3分別從方位角φ1為60°、φ2為65°、φ3為70°，俯仰角 θ1為35°、θ2為40°、θ3為45°入射，信號頻率f為5 GHz，噪聲方差σ2=0.1，信噪比均為10 dB。圓錐的錐頂角α為20°，兩相鄰均勻圓陣之間的間距d為0.1 m。對模型2進行仿真分析，仿真結(jié)果如圖3所示。

圖3 模型2仿真結(jié)果

由圖3可知，信噪較低的環(huán)境下，MUSIC算法仍能估算出波達方向。

4 結(jié) 論

本文基于已建立的錐面陣列天線模型詳細推導應(yīng)用于DOA估計的MUSIC算法，利用MATLAB軟件基于MUSIC算法對均勻圓陣模型和錐面共形陣列天線模型進行DOA估計仿真。仿真結(jié)果驗證了MUSIC算法在DOA估計中的可行性。