陳嘉鑫，黎柏春

（中國(guó)民航大學(xué) 航空工程學(xué)院，天津 300300）

0 引言

GPS（global positioning system）系統(tǒng)作為目前世界上應(yīng)用最廣泛的衛(wèi)星導(dǎo)航系統(tǒng)，在國(guó)家經(jīng)濟(jì)建設(shè)和國(guó)防安全方面有著非常重要的地位。在高動(dòng)態(tài)的環(huán)境下，GPS 接收機(jī)載體的速度和加速度會(huì)有較大的變化范圍，再加上會(huì)有其他的干擾源會(huì)散發(fā)干擾信號(hào)，使得GPS 接收機(jī)對(duì)衛(wèi)星信號(hào)的捕捉、跟蹤和定位難度會(huì)進(jìn)一步增大。因此，對(duì)于高動(dòng)態(tài)環(huán)境下GPS性能研究仍然有非常大的研究空缺，需要不斷完善。

目前，已有多種抗干擾技術(shù)，這些技術(shù)包含時(shí)域處理法、頻域處理法、時(shí)頻域處理法、空域處理法和空時(shí)域處理法等［1-3］?，F(xiàn)階段，美國(guó)在研究第3 代GPS，將抗干擾作為核心問(wèn)題，對(duì)其進(jìn)行了大量的研究和總結(jié)，許多研究成果都直接面向?qū)嶋H應(yīng)用。在我國(guó)，抗干擾難題也成為國(guó)內(nèi)許多高校和研究所的研究方向。例如西安電子科技大學(xué)提出的在微分約束下空域最小功率遞推法［4-5］，中國(guó)民航大學(xué)的微分約束最小功率抗干擾算法，國(guó)防科大的GPS 空時(shí)抗干擾的仿真和研究成果等［6-10］。但是國(guó)內(nèi)的抗干擾算法研究還不是太深入，許多成果還是基于國(guó)外提出方法的跟蹤、改進(jìn)和實(shí)現(xiàn)。

1 高動(dòng)態(tài)環(huán)境下陣列接收信號(hào)模型

1.1 高動(dòng)態(tài)環(huán)境的定義

早在1988 年，美國(guó)宇航局噴氣推進(jìn)實(shí)驗(yàn)室在一份關(guān)于高動(dòng)態(tài)GPS 信號(hào)跟蹤技術(shù)的報(bào)告中就明確給出了高動(dòng)態(tài)的2 種定義。第1 種是在1 s 內(nèi)有70 g的加速度斜升，也就是用70g/s 的加加速度持續(xù)1 s；第2 種比第1 種情況還要更劇烈，是0.5 s 內(nèi)有50g的加速度斜升，也就是用100g/s 的加加速度持續(xù)0.5 s，而且50g/s 的加速度持續(xù)2 s，最后在0.5 s 產(chǎn)生50g/s 的加速度斜降。

1.2 高動(dòng)態(tài)下陣列接收信號(hào)模型

設(shè)置M元等距離線陣，各陣元之間距離，其中λ為入射信號(hào)的載波波長(zhǎng)，另設(shè)信號(hào)入射方向方位角為φ，是與陣列法線方向夾角所形成的方位角，φ∈[-π/2，π/2]。設(shè)最左側(cè)的陣元為原點(diǎn)，可以推導(dǎo)出信號(hào)的入射導(dǎo)向矢量為

GPS 接收機(jī)接收到的信號(hào)中包含許多顆衛(wèi)星信號(hào)，還有干擾信號(hào)和噪聲，各信號(hào)及噪聲之間相互獨(dú)立。因?yàn)樾l(wèi)星離地球地面距離較遠(yuǎn)，可近似認(rèn)為相對(duì)于穩(wěn)定的接收機(jī)，接收到的衛(wèi)星信號(hào)來(lái)向不變，而干擾相對(duì)于接收機(jī)距離較近，認(rèn)為干擾信號(hào)來(lái)向會(huì)快速變化。設(shè)有W個(gè)衛(wèi)星信號(hào)、Z個(gè)干擾信號(hào)，陣列接收到的信號(hào)模型表示為

式中：sw(t)表示第w個(gè)衛(wèi)星信號(hào)(w= 1，2，...，W)；φw表示衛(wèi)星信號(hào)的來(lái)向；jz(t)表示第z個(gè)干擾信號(hào)；φz表示干擾信號(hào)的來(lái)向；n(t)表示接收機(jī)熱噪聲矢量，一般是均值為0，方差為σ2n的高斯白噪聲矢量。因?yàn)樾l(wèi)星距離地球遠(yuǎn)，信號(hào)到達(dá)接收機(jī)會(huì)損失很多能量，其電平會(huì)遠(yuǎn)低于噪聲電平，所以可認(rèn)為陣列接受信號(hào)的協(xié)方差矩陣主要由干擾信號(hào)與噪聲信號(hào)決定，設(shè)協(xié)方差矩陣R為

式中：σ2z為第z個(gè)干擾信號(hào)的功率；σ2n為噪聲功率；I為單位矩陣。因?yàn)樽赃m應(yīng)處理要求信號(hào)是穩(wěn)定的，但式（3）中的協(xié)方差矩陣是不斷變化的，所以只有少數(shù)可以用來(lái)估計(jì)自適應(yīng)處理的協(xié)方差矩陣，式（3）轉(zhuǎn)變?yōu)?/p>

式中：φz表示此刻的干擾來(lái)向；?(t)表示此刻快拍數(shù)很少的采樣數(shù)據(jù)，采樣的協(xié)方差矩陣的第m行第n列元素則可以表示為

式中：δmn為Kroneckerδ函數(shù)，即當(dāng)m=n時(shí)，δmn=1；當(dāng)m≠n時(shí)，δmn=0。

2 Mailloux 算法原理及仿真

2.1 Mailloux 算法原理［11］

設(shè)M元等距離線陣，相鄰陣元之間距離為d，接受了P個(gè)遠(yuǎn)距離的窄帶信號(hào)，入射角設(shè)為φP，則陣列接受的快拍數(shù)據(jù)表示為

式中：x(t)=(x1(t)，x2(t)，...xM(t))T，x(t)為M× 1 維陣列數(shù)據(jù)向量；A為陣列流行矩陣，A=(a(φ1)，a(φ2)，...，a(φp))，矩陣中a(φp)表示第p個(gè)信號(hào) 源 的 導(dǎo) 向 矢 量表示為信號(hào)復(fù)包絡(luò)向量，sP(t)是第p個(gè)信號(hào)源的復(fù)包絡(luò)。n(t)=(n1(t)，n2(t)，...nM(t))T為M×1維陣列噪聲向量。陣列的協(xié)方差矩陣定義為

式中：Rs=E[s(t)sH(t)]為信號(hào)的協(xié)方差矩陣；Ι為M維度單位矩陣；σ2n為噪聲的功率。Mailloux 方法是假設(shè)在接收機(jī)接受的干擾信號(hào)源附近，存在K個(gè)相等強(qiáng)度的虛擬干擾，并按照一定的空間間隔排列，用來(lái)構(gòu)建一個(gè)新的協(xié)方差矩陣Rˉ。用新的協(xié)方差矩陣進(jìn)行自適應(yīng)處理時(shí)相比于原協(xié)方差矩陣，會(huì)在接受的干擾方位附近產(chǎn)生較寬的零陷。新的協(xié)方差矩陣可以表示為

式中：R為沒(méi)加虛擬干擾的接收信號(hào)的噪聲協(xié)方差矩陣；符號(hào)⊙表示Hadamard 積；T為擴(kuò)展矩陣［6］，表達(dá)式為

式 中：T(m，n) 表 示 矩 陣T的 第m行 第n列；Δmn=π(xm-xn)β/λ，xm和xn分別為等距線陣上第m，n個(gè)陣元的位置，β=B/(K- 1)，表示虛擬干擾源之間的距離，一般很小，B為虛擬的角度寬度擬干擾源。

根據(jù)式（9）可知，如果等距線陣用Mailloux 方法，不需要知道干擾信號(hào)的來(lái)向。當(dāng)m=n時(shí)，擴(kuò)展矩陣T的元素T(m，n) =K。新的協(xié)方差矩陣會(huì)使原協(xié)方差矩陣的噪聲項(xiàng)功率增強(qiáng)K倍。

2.2 MVDR 算法仿真

仿真時(shí)，為了顯示該方法形成的零陷較寬，考慮需要有對(duì)照方，所以選擇MVDR（minimum variance distortionless response）算法作為常規(guī)方法，即對(duì)照方。設(shè)置陣元數(shù)為32 的等距線陣，快拍數(shù)為100，工作頻率為3 × 109Hz，波長(zhǎng)為光速除以工作頻率，陣元間隔為波長(zhǎng)的一半，干燥比為40 dB，信噪比為0 dB。因?yàn)楦咚拱自肼曅盘?hào)的功率是單位1，所以可以通過(guò)干燥比來(lái)反推得出干擾信號(hào)，通過(guò)信噪比來(lái)反推得出期望信號(hào)。設(shè)置2個(gè)互相獨(dú)立的干擾信號(hào)入射至陣列接收模型，入射角度分別為-30°和30°，干擾頻率 分 別 設(shè) 為5 × 109Hz 和6 × 109Hz，帶 寬 為0.3 GHz，而期望信號(hào)的入射角在0°方向。將期望信號(hào)、干擾信號(hào)進(jìn)行空域處理，將處理后的期望信號(hào)、干擾信號(hào)，和用randn 函數(shù)產(chǎn)生的噪聲信號(hào)疊加，即可得陣列接收信號(hào)x(t)，通過(guò)x(t)根據(jù)式（3）得出協(xié)方差矩陣R。計(jì)算最優(yōu)權(quán)值，進(jìn)行空間角度掃描和計(jì)算空間導(dǎo)向矢量，最后畫出方向圖，得到MVDR 算法的方向圖，如圖1 所示。從圖1 中可以明顯看出MVDR 算法在2 處干擾方向上形成的零陷很窄。

圖1 MVDR 算法的方向圖仿真結(jié)果Fig.1 Simulation results of MVDR algorithm

2.3 Mailloux 算法仿真

Mailloux 方法則是需要在原協(xié)方差矩陣R上乘上擴(kuò)展矩陣T，設(shè)2 個(gè)干擾信號(hào)有方向上的擾動(dòng)，擾動(dòng)的角度最大范圍為6°，即B= 6°，加入7 個(gè)虛擬干擾源。根據(jù)式（9）計(jì)算可得出擴(kuò)展矩陣T，根據(jù)式（8）算出新協(xié)方差矩陣Rˉ。通過(guò)新的協(xié)方差矩陣來(lái)算最優(yōu)權(quán)值，進(jìn)行空間角度掃描和計(jì)算空間導(dǎo)向矢量，最后畫出方向圖，得到Mailloux 算法與MVDR 算法的方向圖對(duì)比，如圖2 所示。從圖2 對(duì)比來(lái)看，Mailloux 方法能夠明顯地加寬干擾方向上形成的零陷。

圖2 Mailloux 算法的方向圖仿真結(jié)果Fig.2 Direction diagram simulation results of Mailloux algorithm

2.4 高動(dòng)態(tài)下算法對(duì)比仿真

2.4.1 加入1 顆衛(wèi)星信號(hào)仿真

以上是沒(méi)有加入GPS 信號(hào)的仿真，本文研究GPS 抗干擾，所以以下是加入了GPS 信號(hào)的仿真。GPS 衛(wèi)星信號(hào)可以用2 種方法產(chǎn)生，第1 種是直接調(diào)用已生成的GPS 數(shù)據(jù)文件，第2 種是用GPS 數(shù)據(jù)模擬器產(chǎn)生所需要的GPS 數(shù)據(jù)文件。本文采用的是第1 種方法，調(diào)用GPS 衛(wèi)星信號(hào)產(chǎn)生程序，該程序調(diào)用已生成的GPS 衛(wèi)星信號(hào)數(shù)據(jù)文件，先調(diào)用了1 顆衛(wèi)星信號(hào)，衛(wèi)星信號(hào)編號(hào)為1，讀取0.1 s 的衛(wèi)星信號(hào)數(shù)據(jù)，設(shè)置采樣率為5 714 000 Hz，所以衛(wèi)星信號(hào)數(shù)據(jù)為1 × 5 714 00 維矩陣。為保證程序能夠運(yùn)行，所以將快拍數(shù)也改為5 714 00。將產(chǎn)生的GPS 衛(wèi)星信號(hào)矩陣替代原仿真程序中的期望信號(hào)。陣列接收信號(hào)模型重新設(shè)置為7 元陣列等距線陣，保留一個(gè)干擾信號(hào)，干擾來(lái)向設(shè)為-30°，干擾信號(hào)和噪聲信號(hào)矩陣改變維度，其他如工作頻率和波長(zhǎng)都不變。

因?yàn)樾l(wèi)星軌道高度在20 000 km 以上，高動(dòng)態(tài)GPS 接收機(jī)距離衛(wèi)星很遠(yuǎn)，所以在短時(shí)間內(nèi)，近似認(rèn)為接收機(jī)接收到的GPS 衛(wèi)星信號(hào)的來(lái)向不變，也就是φw不變，仿真中設(shè)為30°。為了體現(xiàn)高動(dòng)態(tài)，認(rèn)為干擾信號(hào)距離接收機(jī)較近，所以認(rèn)為干擾信號(hào)的波達(dá)方向會(huì)逐漸變化，但變化不大。仿真中為了直觀方便，設(shè)置干擾信號(hào)的前N/2 個(gè)數(shù)據(jù)的波達(dá)方向?yàn)?30°，后N/2 個(gè)數(shù)據(jù)的波達(dá)方向?yàn)?32°。利用干擾信號(hào)的前N/2 個(gè)數(shù)據(jù)求解陣列加權(quán)向量，得到高動(dòng)態(tài)下Mailloux 算法的方向圖，如圖3 所示。

圖3 1 顆衛(wèi)星的Mailloux 算法方向圖仿真結(jié)果Fig.3 Direction diagram of one satellite based on Mailloux algorithm

從圖3 中可看出，零陷明顯加寬。

新建函數(shù)：

式中：WH為利用干擾信號(hào)的前N/2 個(gè)數(shù)據(jù)求解陣列加權(quán)向量；X為后N/2 個(gè)樣本數(shù)據(jù)。

將新建函數(shù)作為輸入?yún)?shù)帶入捕獲程序進(jìn)行捕獲。先計(jì)算每個(gè)C/A 碼包含多少采樣點(diǎn)，取相鄰的2 個(gè)信號(hào)，信號(hào)長(zhǎng)度為1 ms，去直流，計(jì)算采樣時(shí)間間隔，計(jì)算1 ms 信號(hào)內(nèi)的載波相位點(diǎn)，計(jì)算頻率捕獲的個(gè)數(shù)，產(chǎn)生一顆衛(wèi)星的載波偏移初始向量，初始化32 顆衛(wèi)星的載波頻率，對(duì)已產(chǎn)生的當(dāng)前通道的C/A 碼做離散傅里葉變換，掃描頻點(diǎn)，找出最大的值及其對(duì)應(yīng)的頻點(diǎn)，確認(rèn)最大值對(duì)應(yīng)的相位，獲取與最大峰超過(guò)1 個(gè)碼片的次最大峰，求與最大峰距離為1 個(gè)碼片的區(qū)域邊界，確定尋找次最大峰的搜索區(qū)域，保存各通道峰峰比值，判決，如果比值比預(yù)設(shè)閾值大，則表示能接收到該衛(wèi)星信號(hào)。圖4 為MVDR 算法1 號(hào) 衛(wèi)星信號(hào)捕獲圖，圖5 為Mailloux 算法1 號(hào)衛(wèi)星信號(hào)捕獲圖。

圖4 MVDR 算法1 號(hào)衛(wèi)星信號(hào)捕獲圖Fig.4 Signal acquisition diagram of satellite No.1 base on MVDR algorithm

圖5 Mailloux 算法1 號(hào)衛(wèi)星信號(hào)捕獲圖Fig.5 Signal acquisition diagram of satellite No.1 based on Mailloux algorithm

2.4.2 加入4 顆衛(wèi)星信號(hào)仿真

為驗(yàn)證算法可行度，加入多個(gè)衛(wèi)星信號(hào)。調(diào)用4 顆衛(wèi)星信號(hào)，分別是1，2，3，20 號(hào)衛(wèi)星。圖6 是加4顆衛(wèi)星的Mailloux 算法的方向圖，圖7 是MVDR 算法4 顆衛(wèi)星信號(hào)捕獲圖，圖8 是Mailloux 算法4 顆衛(wèi)星信號(hào)捕獲圖。當(dāng)加入多顆衛(wèi)星后，如果捕獲圖8中沒(méi)有出現(xiàn)已加入的衛(wèi)星序號(hào)，只需要將捕獲門限值提高一點(diǎn)即可。從圖7，8 的對(duì)比來(lái)看，Mailloux 算法在高動(dòng)態(tài)條件下能夠捕獲衛(wèi)星信號(hào)。

圖6 4 顆衛(wèi)星的Mailloux 算法的方向圖Fig.6 Direction diagram of four satellites based on Mailloux algorithm

圖7 MVDR 算法4 顆衛(wèi)星信號(hào)捕獲圖Fig.7 Signal acquisition diagram of four satellites based on MVDR algorithm

圖8 Mailloux 算法4 顆衛(wèi)星信號(hào)捕獲圖Fig.8 Signal acquisition diagram of four satellites based on Mailloux algorithm

3 Zatman 算法原理及仿真

3.1 Zatman 算 法原理［12-14］

Zatman 方法是假設(shè)接收機(jī)接受的窄帶干擾信號(hào)具有虛擬帶寬Bw，而且在整條帶寬內(nèi)有平坦的功率譜，用來(lái)構(gòu)建一個(gè)新的協(xié)方差矩陣Rˉ，用新的協(xié)方差矩陣進(jìn)行自適應(yīng)處理時(shí)相比于原協(xié)方差矩陣，會(huì)在接受的干擾方位附近產(chǎn)生較寬的零陷。新的協(xié)方差矩陣可以表示為

式中：T(m，n)表示矩陣T的第m行第n列；xm和xn分別為等距線陣上第m和n個(gè)陣元的位置；φz表示第z個(gè)干擾信號(hào)的來(lái)向；c為光速。

由式（11）可知，Zatman方法需要知道干擾信號(hào)的來(lái)向。當(dāng)m=n時(shí)，擴(kuò)展矩陣T的元素T(m，n) = 1，不會(huì)改變?cè)瓍f(xié)方差矩陣R中的噪聲項(xiàng)功率。

3.2 Zatman 算法仿真

對(duì)照方依舊選擇MVDR 算法。先不加衛(wèi)星信號(hào)，使用期望信號(hào)。陣元數(shù)為32 的等距線陣，快拍數(shù)為100，使用2 個(gè)獨(dú)立的干擾信號(hào)。工作頻率為3 ×109Hz，波長(zhǎng)等于光速除以工作頻率，陣元間隔為波長(zhǎng)的一半，干燥比為40 dB，信噪比為0 dB，設(shè)高斯白噪聲的功率為1。通過(guò)干燥比反推得到干擾信號(hào)，通過(guò)信噪比反推得到期望信號(hào)。將各信號(hào)疊加得到x(t)，剩余的操作與第3 節(jié)中MVDR 算法仿真中的操作一致。MVDR 算法得到的方向圖如圖1 所示。

Zatman 算法是在MVDR 算法的基礎(chǔ)矩陣上再乘上一個(gè)擴(kuò)展矩陣T。假設(shè)干擾信號(hào)的頻帶擴(kuò)展寬度為0.3 GHz，即Bw= 0.3 GHz。根據(jù)式（11）得到擴(kuò)展矩陣，根據(jù)式（8）算出新協(xié)方差矩陣Rˉ。通過(guò)新的協(xié)方差矩陣來(lái)算最優(yōu)權(quán)值，進(jìn)行空間角度掃描和計(jì)算空間導(dǎo)向矢量，最后畫出方向圖，得到Zatman 算法與MVDR 算法的方向圖對(duì)比，如圖9 所示。

圖9 Zatman 算法的方向圖Fig.9 Direction diagram of Zatman algorithm

3.3 高動(dòng)態(tài)下算法對(duì)比仿真

3.3.1 加入1 顆衛(wèi)星信號(hào)仿真

以上是沒(méi)有加衛(wèi)星信號(hào)的仿真，以下是加入衛(wèi)星信號(hào)的仿真。還是用第1 種方法直接調(diào)用已生成的GPS 數(shù)據(jù)文件。先調(diào)用1 號(hào)衛(wèi)星信號(hào)，讀取參數(shù)與第3 節(jié)里一致。

為保證程序能夠運(yùn)行，將快拍數(shù)也改為571 400。將產(chǎn)生的1 號(hào)衛(wèi)星信號(hào)矩陣替代原仿真程序中的期望信號(hào)。設(shè)置7 元陣列等距線陣，保留一個(gè)干擾信號(hào)，干擾來(lái)向設(shè)為-30°。為了體現(xiàn)高動(dòng)態(tài)，也將設(shè)置干擾信號(hào)的前N/2 個(gè)數(shù)據(jù)的波達(dá)方向?yàn)?30°，后N/2個(gè)數(shù)據(jù)的波達(dá)方向?yàn)?32°。用前N/2 個(gè)數(shù)據(jù)求解陣列加權(quán)向量，得到的加權(quán)向量處理后N/2 個(gè)數(shù)據(jù)，得到加寬的方向圖，如圖10 所示。從圖10 中可看出，零陷加寬。

圖10 1 號(hào)衛(wèi)星的Zatman 算法的方向圖Fig.10 Direction diagram of one satellite based on Zatman algorithm

根據(jù)式（10）得到新建函數(shù)，將新建函數(shù)送入捕獲程序中捕獲衛(wèi)星信號(hào)。捕獲程序流程與2.3中所寫一致。MVDR 算法1 號(hào)衛(wèi)星信號(hào)捕獲圖如圖4 所示，圖11 為Zatman 算法1 號(hào)衛(wèi)星信號(hào)捕獲圖。從圖11 中可知，Zatman 算法在高動(dòng)態(tài)條件下可以捕獲到1 號(hào)衛(wèi)星。

圖11 Zatman 算法1 號(hào)衛(wèi)星信號(hào)捕獲圖Fig.11 Signal acquisition diagram of satellite No.1 based on Zatman algorithm

3.3.2 加入4 顆衛(wèi)星信號(hào)仿真

為驗(yàn)證算法有效性，加入多個(gè)衛(wèi)星信號(hào)進(jìn)行仿真。調(diào)用4 顆衛(wèi)星信號(hào)，分別是1，2，3，20 號(hào)衛(wèi)星。圖12 為加4 顆衛(wèi)星的Zatman 算法的方向圖，MVDR算法4 顆衛(wèi)星信號(hào)捕獲圖如圖7 所示，圖13 為Zatman算法4 顆衛(wèi)星信號(hào)捕獲圖。當(dāng)加入多顆衛(wèi)星后，捕獲圖也出現(xiàn)了Mailloux 算法類似的問(wèn)題，同樣只需要將捕獲門限值提高一點(diǎn)即可。從圖7，13 的對(duì)比來(lái)看，Zatman 算法在高動(dòng)態(tài)下能夠捕獲衛(wèi)星信號(hào)。

圖12 4 顆衛(wèi)星的Zatman 算法的方向圖Fig.12 Direction diagram of four satellites based on Zatman algorithm

圖13 Zatman 算法4 顆衛(wèi)星信號(hào)捕獲圖Fig.13 Signal acquisition diagram of four satellites based on Zatman algorithm

4 Mailloux 算 法 和Zatman 算 法 的 優(yōu)缺點(diǎn)對(duì)比分析

根據(jù)2 種算法的推導(dǎo)和仿真結(jié)果的對(duì)比可得，Mailloux 算法運(yùn)算量適中，在線陣模型中不需要估計(jì)干擾來(lái)向，可以直接得到擴(kuò)展矩陣，但在面陣模型中，該算法需提前給出預(yù)估得干擾信號(hào)來(lái)向。此外還有一個(gè)缺點(diǎn)，用T計(jì)算Rˉ時(shí)，協(xié)方差矩陣的噪聲項(xiàng)功率會(huì)增強(qiáng)K倍。該算法在方向圖干擾來(lái)向上形成的零陷比Zatman 算法較寬，比MVDR 算法寬很多。

Zatman 算法運(yùn)算量不大，對(duì)于干擾噪聲的協(xié)方差矩陣的功率沒(méi)有增強(qiáng)，但無(wú)法保證角度的取值，不能進(jìn)行角度的擴(kuò)展，因此該算法不可應(yīng)用于面陣，只可用于線陣。該算法在方向圖干擾來(lái)向上形成的零陷比MVDR 算法寬，但相比Mailloux 算法要窄一些。

通過(guò)2 種算法得對(duì)比可知，Zatman 算法不能用在圓陣以及任意陣模型下，利用Mailloux 算法需要預(yù)估干擾的來(lái)向信息。但是，如果干擾方向隨時(shí)間快速變化，難以得到干擾方向的來(lái)向信息，Mailloux算法便不能解決問(wèn)題，與此同時(shí)也會(huì)增加很大的計(jì)算量。

5 結(jié)束語(yǔ)

本文首先研究了當(dāng)陣列為等距線陣時(shí)接收信號(hào)模型。根據(jù)信號(hào)入射的導(dǎo)向矢量和信號(hào)中不同組成部分，得到陣列接收信號(hào)模型。并由信號(hào)模型求得協(xié)方差矩陣，來(lái)表示各個(gè)維度信號(hào)之間的相關(guān)性。在Matlab 上搭建了陣列接收信號(hào)模型。之后對(duì)Mailloux 算法和Zatman 算法原理進(jìn)行了研究分析。仿真時(shí)在陣列接收信號(hào)模型的基礎(chǔ)上，通過(guò)虛加干擾源和擴(kuò)展頻帶來(lái)獲得擴(kuò)展矩陣，從而得到新的協(xié)方差矩陣。為更直觀驗(yàn)證算法的有效性，選擇MVDR 算法作為對(duì)照。最后為體現(xiàn)高動(dòng)態(tài)的特點(diǎn)，仿真時(shí)改變干擾信號(hào)數(shù)據(jù)后半部分的來(lái)向。用前半部分?jǐn)?shù)據(jù)計(jì)算新的協(xié)方差矩陣，得到的矩陣處理后半部分?jǐn)?shù)據(jù)。

通過(guò)仿真驗(yàn)證可知，在高動(dòng)態(tài)條件下，當(dāng)陣元是等距線陣時(shí)，Mailloux 算法不需要知道干擾的來(lái)向信息，但Zatman 算法需要知道干擾的來(lái)向。相較于MVDR 算法，Mailloux 算法和Zatman 算法都能在方向圖干擾來(lái)向上形成較寬的零陷，能夠有效抑制掉干擾信號(hào)。根據(jù)需要，將3 種算法形成的新建函數(shù)送入捕獲程序后，MVDR 算法不能捕獲到衛(wèi)星信號(hào)，而Mailloux 算法和Zatman 算法能捕獲到衛(wèi)星信號(hào)。

本文只驗(yàn)證了陣列接收信號(hào)模型是在等距線陣時(shí)算法的有效性。由于在飛機(jī)、導(dǎo)彈等一系列高速運(yùn)動(dòng)的載體上，圓陣列的應(yīng)用比較廣泛，后續(xù)將驗(yàn)證圓陣列和其他任意形模型下，Mailloux 算法和Zatman 算法在高動(dòng)態(tài)條件下的抗干擾有效性。