減輕衛(wèi)星信號接收器中干擾的方法和設(shè)備維護

林列書

(1.南京航空航天大學(xué) 繼續(xù)教育學(xué)院，江蘇 南京 211100； 2.廣州民航職業(yè)技術(shù)學(xué)院，廣東 廣州 510470)

1 引言

處于空間軌道當(dāng)中運行的GNSS衛(wèi)星，在地面所接收到的信號非常弱，所以勢必會受到外界各種因素的干擾，所以必須要采取有效的手段抑制這些干擾信號，確保衛(wèi)星接收機工作正常[1]。在GNSS接收機中，抑制干擾的主要方法之一就是進行自適應(yīng)陣列處理[2]。在自適應(yīng)算法的基礎(chǔ)上調(diào)整陣列加權(quán)系數(shù)，并借助空域特征區(qū)分干擾和信號，在信號方位可以直接形成對應(yīng)的波束，所以干擾處可以出現(xiàn)零缺陷[3]。如果已經(jīng)知道了信號導(dǎo)向矢量，在傳統(tǒng)波束的影響下，將會出現(xiàn)最小方差無畸變波束形成算法等來避免形成干擾因素。在GNSS導(dǎo)航的應(yīng)用領(lǐng)域?qū)嶋H上無法正常使用自適應(yīng)波束形成算法的。

估算衛(wèi)星定位的過程可以通過GNSS信號中的歷史信息得以實現(xiàn)[4]。但實際上，如果處于干擾環(huán)境下，是很難獲取歷史數(shù)據(jù)的。整理歸納GNSS接收機當(dāng)中的存儲歷史，從而能夠遞推出相應(yīng)的衛(wèi)星方位，但需要注意的是，無法確保最終結(jié)果的一致性。同時，由于各方面因素的顯著干擾，會造成GNSS接收機平臺出現(xiàn)問題，很難通過衛(wèi)星定位的方式確定GNSS信號導(dǎo)向矢量。信號導(dǎo)向矢量的精確性實際上是自適應(yīng)波束形成算法的關(guān)鍵內(nèi)容，否則將會造成嚴(yán)重的指向偏差，甚至導(dǎo)致信號能力的損失，阻礙接收機的正常運行過程。現(xiàn)階段，已經(jīng)出現(xiàn)了很多不需要信號導(dǎo)向矢量的自適應(yīng)波束形成算法，在無需信號導(dǎo)向矢量時，通過功率倒置法可對強功率干擾進行有效抑制，但在進行弱干擾抑制時，其能力有限[5]。基于此，本文對減輕衛(wèi)星信號接收器中干擾的方法進行了分析研究，并提出對衛(wèi)星接收器設(shè)備進行維護。

2 陣列信號模型

對于單頻信號而言，通過傳播波的波動方程，可得到如公式(1)所示的波的表達(dá)式：

s(r,t)=Aexp[j(ωt-KTr)]

(1)

公式中，波傳播路徑矢量用表示，其中字母A、k分別表示復(fù)振幅和波束矢量，k的單位表示為弧度/米，單位波長的周期數(shù)即為其實際大小。

由于陣列會同時接收多個信號，假設(shè)此時共K顆GNSS衛(wèi)星信號被接收到，有P個干擾被接收到，變頻和采樣后所得到的信號如公式(2)所示。

(2)

3 基于自相關(guān)的波束形成算法

在常見的GNSS接收機結(jié)構(gòu)中，由PLL和DLL等相關(guān)硬件能夠捕獲、跟蹤GNSS信號接收機所收到的信號，并采取解算處理的方式將其實現(xiàn)，GNSS信號信噪比對其捕獲和定位過程將產(chǎn)生極大影響，如果信噪比較低，此時其定位精度將出現(xiàn)混亂現(xiàn)象。信號信噪比損失是干擾抑制過程中需要強調(diào)的重要內(nèi)容。在抑制干擾的過程中，GNSS信號信噪比的最大值出現(xiàn)在輸出過程當(dāng)中。可用公式(3)表示接收信號x(n)的協(xié)方差矩陣：

(3)

(4)

對于PRN1衛(wèi)星發(fā)射的GNSS信號而言，在處理時，GNSS接收機產(chǎn)生的參考信號要與PRN1衛(wèi)星對應(yīng)，其表達(dá)式見公式(5)：

r1(n)=g1(n)cos(ω1n)

(5)

公式中，PRN1衛(wèi)星對應(yīng)的擴頻序列為g1(n)，用公式(6)表示參考信號r1(n)和接收信號x(n)的相關(guān)值c(m)：

(6)

公式中，第i個陣元收到的信號為xi(n)，i=1,2,···,N，積分長度為T。構(gòu)造向量c(m)=[c1(m),c2(m),…,cN(m)]T。使用的擴頻序列，GNSS信號具有非常強的自相關(guān)性，互相關(guān)性較差，因此在相關(guān)值c(m)中，因擴頻增益，PRN1衛(wèi)星信號功率增大。

(7)

(8)

因此，在公式(8)中，最大化輸出信干噪比等價于優(yōu)化目標(biāo)，矩陣對(C,Rx)的最大廣義特征值對應(yīng)的特征矢量為最優(yōu)權(quán)，見公式(9)：

CWopt=λmaxRxWopt

(9)

公式中，矩陣對(C,Rx)的最大廣義特征值為λmax，λmax對應(yīng)的特征矢量為Wopt。在本算法中。GNSS接收機具有并行處理通道的能力，不同的接收機通道與參考信號是相互對應(yīng)的，按照上述內(nèi)容能夠?qū)πl(wèi)星信號干擾抑制的過程進行實現(xiàn)，滿足定位和結(jié)算等相關(guān)過程。

4 算法實施及性能分析

4.1 符號翻轉(zhuǎn)對所提方法性能影響及抑制方法

本研究當(dāng)中存在的算法必須制定與之相適應(yīng)的相關(guān)值矢量的協(xié)方差矩陣Rc,1，從而實現(xiàn)估算的過程，如果協(xié)方差矩陣估計性能損失沒有達(dá)到3 dB時，此時需要增強協(xié)方差矩陣的系統(tǒng)自由度，將其擴展至2倍以上，所以，在對Rc,1進行估算時，首先需要獲得充足的相關(guān)值樣本。導(dǎo)航電文信息包含在GNSS信號中，在計算公式(6)當(dāng)中的相關(guān)值時，如果對導(dǎo)航電文符號進行翻轉(zhuǎn)，此時將會不斷衰減擴頻序列，對最終得到的估算結(jié)果造成顯著影響。將相關(guān)值峰值的有效性充分展現(xiàn)出來，對最終結(jié)果進行估算，如果反轉(zhuǎn)符號，此時可以有效確保協(xié)方差矩陣的正常運行。假設(shè)當(dāng)前已經(jīng)獲得了GNSS信號，最終得出公式(10)：

(10)

公式中，T0表示完整碼序列的長度。在公式(10)當(dāng)中融入公式(2)的接收信號可以得到：

(11)

綜合分析GNSS衛(wèi)星PRN1信號和參考信號的相關(guān)值，在信號接收過程中要重視其他信號的分量，從而簡化其形式，最終能夠得到公式(12)：

(12)

(g1(n)cos(ω1n))*

(13)

因相關(guān)值峰值是參考信號完全與GNSS信號匹配的結(jié)果，在GNSS信號內(nèi)，符號翻轉(zhuǎn)不會發(fā)生，因此可得：

(14)

則協(xié)方差矩陣為：

(15)

因此，使用相關(guān)值峰值對協(xié)方差矩陣進行估計，確保相關(guān)值幅度不會受到符號翻轉(zhuǎn)的影響。在相關(guān)值峰值的基礎(chǔ)上估計相關(guān)矩陣，無需從干擾和噪聲獲得精確相關(guān)值。

4.2 GNSS信號多普勒捕獲誤差對所提方法的影響

因為接收機和衛(wèi)星的相對運動關(guān)系非常顯著，所以如果GNSS信號傳遞到接收機后，此時會產(chǎn)生多普勒頻偏現(xiàn)象。本文是在GNSS信號相關(guān)性和參考信號之上提出的方法，考慮多普勒頻偏會影響所提方法性能，可用公式(16)表示多普勒頻偏對相關(guān)值的影響：

(16)

公式中，相干積分時間為T，多普勒頻率為fd。參考信號和接收GNSS信號的相關(guān)為c(τ)。由公式(16)知，在接收信號有多普勒存在時，通過相關(guān)運算，受sinc函數(shù)調(diào)制，相關(guān)值下降嚴(yán)重。

一般來說，GNSS信號能夠接收到的-10 ～10 kHz的多普勒頻率。多普勒頻率為0時，此時得到最大的GNSS信號和參考信號相關(guān)值。對本地信號的載波頻率進行變換，準(zhǔn)確找出最大相關(guān)峰值，從而對多普勒頻率進行接收，因為搜索波長的限制，再加上載波頻率仍有剩余多普勒，假設(shè)補償后，Δf為剩余多普勒頻率，公式(17)為GNSS信號相關(guān)矩陣：

(17)

多普勒頻偏影響算法主要在于GNSS信號的損耗Rc,1和自相關(guān)值，如果對擴頻增益進行降低處理，此時可得到1 ms的積分時間，為了將損耗控制在1 dB范圍內(nèi)，此時需要確保多普勒剩余頻率不超過500 Hz。

5 仿真與實測結(jié)果

本文所采取的驗證過程為實測數(shù)據(jù)和計算機仿真的手段，同時給出SCORE方法、功率倒置法的結(jié)果進行比較。

5.1 仿真結(jié)果

仿真平臺使用MATLAB軟件，GNSS信號選擇的信號為L1=1 574.31 MHz上的GPSC/A碼，天線是半波長排列的4元方陣，如圖1所示。設(shè)定中頻信號為15.234 MHz，在各仿真情況中，對其余參數(shù)進行具體說明。

圖1 陣列天線模型

假設(shè)有一個窄帶干擾在接收信號當(dāng)中產(chǎn)生作用，此時單頻余弦是以干擾信號的方式存在的，頻率為14.8 MHz，表1為窄帶干擾仿真參數(shù)。

表11 窄帶干擾仿真參數(shù)

對比項信噪比/dB方位角θ俯仰角φ信號-20.5184.5°90.5°干擾25.5129.5°15.5°

圖2為窄帶干擾抑制結(jié)果，由圖2知，相比于另外兩種算法，本研究所使用的方法可以有效增強輸出信息的干燥比，其信號增益接近四元陣列增益上限。通過模歸一的方法處理各算法陣列加權(quán)值，其中信號增益發(fā)揮著至關(guān)重要的作用，分別統(tǒng)計輸入輸出信號，最終得出信號增益。

圖2 窄帶干擾抑制結(jié)果

5.2 實測數(shù)據(jù)分析

GNSS接收機天線陣列是四元方陣非常典型的代表，中頻和采樣頻率分別為7.48 MHz、10.05 MHz。在中頻輸出端得到了-90.5 dBm的GPS信號加噪聲功率，不同的天線單元均有41.5 dB的中頻增益；-50.5 dBm功率的干擾信號，在實驗過程中，設(shè)置單頻余弦的過程即為干擾信號的主要形式，達(dá)到1 574.31 MHz的頻率。在坐標(biāo)系下為方位角θ=0°和俯仰角φ=30°為干擾來波方向。由于鐵桿架設(shè)小天線到GPS陣列陣面的垂直高度已知，根據(jù)相關(guān)數(shù)據(jù)可得出鐵桿水平位置，確定其干擾方向。實測數(shù)據(jù)表明，在多種干擾場景下，所提方法均能干擾抑制能力均良好，相比于SCORE方法、功率倒置法，本研究所提方法在抑制干擾的同時，能保證輸出中GNSS信號信干噪比最大化，可見對陣列誤差穩(wěn)健性很強。

6 衛(wèi)星接收器的維護

衛(wèi)星接收器完成安裝、調(diào)試、固定后，一般在工作的過程中不會有問題出現(xiàn)。但因衛(wèi)星接收器常年都暴露在外，受陽光照射、冷熱變化、酸雨腐蝕、大風(fēng)等天氣影響，會損害衛(wèi)星接收器表面的涂覆，進而出現(xiàn)腐蝕等。因此，在使用衛(wèi)星接收器時，對其定期進行清理、維護、保養(yǎng)等工作。若發(fā)現(xiàn)接收器表面發(fā)生漆層脫落后，要及時對其補修一般需要對接收器進行定期的維護保養(yǎng)，以免造成反射結(jié)果出現(xiàn)較大誤差。定期對傳動系統(tǒng)的軸承等構(gòu)建進行黃油的涂抹，避免出現(xiàn)腐蝕現(xiàn)象。及時清理接收器拋物面中的垃圾，對外露接口件做防水處理；確保天線饋源口薄膜完好。饋源內(nèi)部不能有水珠、水汽等。冬季要及時清除饋源與反射面中的冰凌、積雪等。

7 結(jié)論

本文基于一種對導(dǎo)航信號當(dāng)中產(chǎn)生抑制現(xiàn)象的干擾盲波束形成的算法，以免因為干擾因素的存在而造成造成信息接收過程出現(xiàn)問題，最大化的增益GNSS信號，增強其準(zhǔn)確性和有效性。以該方法為基礎(chǔ)不需要事先確定GNSS信號來波方向，陣列誤差穩(wěn)健。實測數(shù)據(jù)表明，在多種干擾場景下，所提方法均能干擾抑制能力均良好，相比于SCORE方法、功率倒置法，本研究所提方法在抑制干擾的同時，能保證輸出中GNSS信號信干噪比最大化，穩(wěn)健性很強。最后提出對衛(wèi)星接收器的維護措施。