(中國西南電子技術研究所，成都610036)

衛(wèi)星導航(簡稱衛(wèi)導)系統(tǒng)是基于無線電測距的系統(tǒng)，天線是實現(xiàn)無線電信號接收的關鍵部分，通過接收機對不同方向入射的衛(wèi)星信號進行測量，得到不同衛(wèi)星到達接收天線處的距離，從而實現(xiàn)幾何交會定位[1]。在實際中，傳統(tǒng)的測量型衛(wèi)導天線具有穩(wěn)定的相位中心，滿足幾何交會定位的依據(jù)；衛(wèi)星信號從空間到達地面，經(jīng)過路徑衰減后功率位于熱噪聲以下，由于擴頻增益，仍然能夠準確地定位。全球衛(wèi)星導航系統(tǒng)在沒有干擾的環(huán)境下工作良好，但是它仍然極易受到干擾，因此在干擾環(huán)境下，研究人員引入了利用空域分集的陣列信號處理技術[2-3]。陣列天線通過多天線接收，調節(jié)不同陣元信號的幅度和相位進行合成，綜合出期望的天線方向圖形狀，實現(xiàn)期望信號的接收和干擾信號的抑制[4]?，F(xiàn)有的衛(wèi)導陣列處理技術利用了衛(wèi)導信號的負信噪比特性，通過功率約束實現(xiàn)了干擾信號的自適應調零，保證了在強干擾環(huán)境下的衛(wèi)導信號可用性，但是在基于載波相位測量的高精度相對定位中，多天線的各向差異和多射頻前端通道之間的不一致性將導致衛(wèi)導距離測量中產生偏差，違反了幾何交會定位的前提，簡而言之，經(jīng)過陣列處理后不具有穩(wěn)定的相位中心[5]。

為了減小陣列信號處理對衛(wèi)導信號中測距信息的破壞，需要對陣列天線和射頻前端通道進行校準。文獻[6]總結了陣列處理對衛(wèi)導信號的載波相位和偽碼相位產生影響的模型，比較了傳統(tǒng)的功率倒置(Power Inversion,PI)和最小方差無失真約束(Minimum Variance Distortionless Response,MVDR)波束形成器產生的影響，但要對偏差進行約束，需要知道準確的陣列流形矢量。文獻[7]提出了一種基于相關函數(shù)的陣列校準算法，實現(xiàn)天線陣列的空時陣列流形矢量的測量，但是沒有更多地討論校準波束的規(guī)劃。文獻[8]通過設置輔助跟蹤通道，但只測量了陣列天線之間的相對相位，而忽略了相對幅度。文獻[5]總結了衛(wèi)導陣列信號處理最新的進展，指出目前衛(wèi)導陣列信號處理的三大目標是可用性、準確性和完好性，而目前的功率倒置類算法只解決了強干擾環(huán)境下的可用性問題。關于目前衛(wèi)導陣列抗干擾處理的現(xiàn)狀和面臨的挑戰(zhàn)，可以參考文獻[5]。

在衛(wèi)導陣列信號處理需要突破的三大目標中，陣列的現(xiàn)場陣列流形矢量都起到關鍵作用[9]。為此，本文提出一種衛(wèi)導陣列天線的快速現(xiàn)場系統(tǒng)校準算法。文獻[10-11]討論了在相控陣雷達中，通過設置輔助信號源的陣列誤差校正。與現(xiàn)有文獻相比，本文提出的射頻通道和天線分離的系統(tǒng)校準方法，在對天線進行校準時，根據(jù)幾何或者暗室測量的陣列流形矢量來規(guī)劃校準波束，不依賴單獨設置的輔助信號源，實現(xiàn)了衛(wèi)星信號現(xiàn)場陣列流形矢量的快速校準。仿真驗證了系統(tǒng)校準算法的有效性。

2 數(shù)據(jù)模型

含有衛(wèi)導陣列天線的典型陣列接收系統(tǒng)如圖 1所示。衛(wèi)導信號經(jīng)陣列天線接收，進入多通道射頻前端進行信號調理，在數(shù)字模塊中進行模擬數(shù)字轉換器(ADC)采樣、數(shù)字波束形成，再進入多通道數(shù)字接收機進行解調和測量，最后輸出供導航定位算法使用的偽距、載波相位和多普勒等原始距離測量值。為了校準天線和射頻前端對衛(wèi)導信號的影響，首先建立該處理過程的數(shù)據(jù)模型。

圖1 衛(wèi)導陣列抗干擾系統(tǒng)Fig.1 GNSS array anti-jamming receiver

假設第k顆衛(wèi)星信號從(θ,φ)方向入射，到達第i個陣元，其復解析信號可以表示為

exp(j((fRF+fk,D)(t-Δtk)+ξi(θ,φ)))。

(1)

exp(j((fIF+fk,D)(t-Δtk)+ξi(θ,φ)))×

Ciexp(jψi)。

(2)

式中：Ci、ψi分別為通道引入的幅度和相位響應。值得指出的是，通道響應與信號入射方向無關。

將t時刻所有I個陣元的接收信號組成一個矢量，可以得到單顆衛(wèi)星的陣列接收模型(為了簡便，省略了中頻上標)：

exp(j((fIF+fk,D)(t-Δtk)))。

(3)

所有K顆衛(wèi)星信號同時接收，可以得到完整的陣列接收信號模型：

(4)

在衛(wèi)導陣列信號數(shù)字波束形成中，因為衛(wèi)導信號具有負信噪比，常通過功率約束，實現(xiàn)強干擾信號的自適應調零[2]：

minwwHRws.t.wHc=c。

(5)

zk(t)=wHyk(t)=wHa(θk,φk)sk(t)。

(6)

由wHa(θk,φk)產生的相位與信號sk(t)中由衛(wèi)星到接收天線之間距離引起的相位無法區(qū)分，因此沒有經(jīng)過系統(tǒng)校準的陣列流形矢量a(θk,φk)和無約束的波束形成器w將導致衛(wèi)導接收機得到的載波相位測量值無法用于基于載波相位差分的高精度相對定位應用中。為了保證經(jīng)過陣列處理后衛(wèi)導信號中載波相位的可用性，本文提出陣列天線的快速現(xiàn)場系統(tǒng)校準算法，旨在獲得系統(tǒng)使用環(huán)境下的準確陣列流形矢量a(θk,φk)。

3 校準算法

根據(jù)上一節(jié)的分析，陣列流形矢量中需要進行校準的非理想差異，主要包括射頻通道的不一致性和天線對不同方向入射信號的空域響應差異。前者是由射頻通道之間的硬件不一致和硬件隨溫度等環(huán)境的不一致變化引起，通過單獨的通道校準過程進行校準；后者是由于信號入射方向、天線的結構和周圍環(huán)境等因素導致的陣列的接收陣列流形矢量與信號到達波程差不同產生的幾何陣列流形矢量之間存在差異，通過衛(wèi)星信號和數(shù)字接收機的輔助，對陣列天線的現(xiàn)場陣列流形矢量進行測量。

3.1 通道校準

設備自主產生校準所用的單音信號，經(jīng)過數(shù)字模擬轉換器(DAC)變換為模擬信號，再通過1-I功分器送往不同的射頻通道入口。校準時，通過射頻開關選通校準信號，經(jīng)過模擬數(shù)字轉換器(ADC)變換成數(shù)字信號，計算不同通道相對于參考通道之間的相關系數(shù)，獲得不同通道之間的幅相不一致性。

圖2 通道校準處理流程Fig.2 Channel calibration procedure

記歸一化的相關系數(shù)向量為h∈I×1，其各個元素計算方法如下：

(7)

式中：yi為第i個通道的接收信號，yref(t)表示參考通道接收信號，上標*表示復數(shù)共軛，t0表示通道校準的起始時刻，T表示通道校準的持續(xù)時間。通道校準就是獲得歸一化的相關系數(shù)向量的過程。一般在每次開機時運行，當設備需要長時間連續(xù)運行時，可在無任務階段周期地進行通道校準。

3.2 天線校準

從式(6)中可以看出，當控制波束形成器w的幅度和相位時，實際上改變了衛(wèi)導接收機的測量輸入，即等價地影響了衛(wèi)導信號的增益和相位。如果同時設置一系列的已知波束wm,m=1,2,…,M，對同一顆衛(wèi)星信號進行跟蹤，那么

(8)

式中：zk,m為第m個校準波束產生的輸出信號。校準波束wm是自主設定的，zk,m(t)可以通過接收機進行測量，那么未知的現(xiàn)場陣列流形矢量a(θk,φk)可以通過上述觀測方程直接求解：

(9)

式中：?表示矩陣的Moore-Penrose偽逆，α是一個與t時刻信號有關的歸一化標量系數(shù)。通過對式(9)中所得估計向量進行歸一化，可以得到與測量時間無關的第k顆衛(wèi)星入射方向的歸一化現(xiàn)場陣列流形矢量，這就是借助現(xiàn)場衛(wèi)星信號和數(shù)字接收機的現(xiàn)場陣列流形矢量快速校準的基本原理。

從改善式(9)的數(shù)值計算條件，提高數(shù)字接收機的跟蹤測量能力出發(fā)，對快速校準處理的流程進行設計。

Step1 根據(jù)當前第k顆衛(wèi)星在天線坐標系下的入射角θi和φi，計算衛(wèi)星入射方向的幾何陣列流形矢量ageo(θk,φk)：

(10)

式中：pi為第i個陣元在天線坐標系下的位置矢量，作為替代也可以使用暗室測量的陣列流形矢量aac(θk,φk)，用于后續(xù)校準過程，統(tǒng)一記為ac0(θk,φk)。

Step2 規(guī)劃校準波束的權矢量：

wk,m=ac0(θk,φk)*ejrmΔ。

(11)

式中：*表示對應元素相乘的Hadamard乘積；rm為第m個校準波束對應的隨機擾動矢量，每個元素從集合{-1,0,1}中取值；Δ控制相位擾動的大小。

Step4 求解最小二乘問題，計算該方向現(xiàn)場陣列流形矢量：

(12)

根據(jù)式(4)，可將當前通道響應從現(xiàn)場陣列流形矢量中分離，得到凈現(xiàn)場天線陣列流形矢量acal2(θk,φk)：

(13)

式中：hi,i=1,2,…,I為歸一化通道相關系數(shù)向量h的第i個元素。

4 數(shù)值仿真

通過軟件接收機[12]對系統(tǒng)校準算法進行仿真，后文如無特別說明，采樣頻率為典型值62 MHz,中頻為46.42 MHz，陣元數(shù)目為7，中心陣元為參考陣元，蒙特卡洛仿真次數(shù)為500。

4.1 通道校準

在通道校準階段，校準信號由設備自主產生，校準信噪比可控，在仿真中信噪比變化范圍設置為0～10 dB，通道校準中只包含校準信號與熱噪聲，通道平坦度為3 dB，每次蒙特卡洛實驗獨立生成一組通道響應。通道校準仿真結果如圖3和圖4所示。因為陣元4為參考陣元，陣元3和陣元7與其他陣元結果類似，圖中只給出了陣元1、2、5、6的統(tǒng)計結果。

圖3 各通道開機校準幅度誤差Fig.3 Amplitude error of channel calibration

圖4 各通道開機校準相位誤差Fig.4 Phase error of channel calibration

從仿真結果可以看出，校準結果隨著信噪比的提升和快拍數(shù)的增加而變得準確。在信號幅度歸一化為1時，開機校準的誤差在信噪比為0 dB，快拍數(shù)104時，校準精度接近10-2，其中，相位的校準誤差小于1°。通道校準由設備自主產生校準信號，通過接收通道相關，測量幅相差異的校準算法與傳統(tǒng)通道幅相測量算法一致，但是通道校準結果是分離天線現(xiàn)場陣列流形矢量所需，是系統(tǒng)校準不可缺少的一部分。仿真結果給出了工程實踐中信噪比和快拍數(shù)的選擇參考。

4.2 天線校準

仿真場景和通道校準仿真相比，增加了8顆衛(wèi)星信號。分別對可觀測的8顆衛(wèi)星進行在線測量它們的現(xiàn)場陣列流行矢量，每1 ms測量一次。

圖5給出了一組隨機校準波束和本文校準波束衛(wèi)導軟件接收機測量的衛(wèi)導信號載噪比結果，可以看出隨機校準波束對于部分衛(wèi)星信號的跟蹤載噪比低于30 dBHz，不能準確地測量衛(wèi)導信號中的載波相位，無法完成現(xiàn)場陣列流形矢量的測量，而本文校準波束的測量值均大于48 dBHz，能夠準確地測量偽距和載波相位，成功地對天線陣列流形矢量進行校準。

(a)隨機校準波束2

(b)本文校準波束2圖5 本文校準波束和隨機校準波束接收機載噪比測量Fig.5 Carrier-to-noise ratio measurement of the proposed and the random calibration beams

圖6給出了每個歷元的測量誤差，圖7給出了連續(xù)100組測量的滑動平均。從圖6和圖7中可以看出，對于大多數(shù)衛(wèi)星，通過多次在線校準測量值的平均，能夠將現(xiàn)場陣列流形矢量的校準誤差縮小到2°以內，幅度誤差能縮小到0.17 dB以內。

(a)幅度

(b)相位圖6 現(xiàn)場陣列流形矢量單次測量誤差Fig.6 Field array manifold single measurement error

(a)幅度

(b)相位圖7 現(xiàn)場陣列流形矢量滑動平均測量誤差Fig.7 Field array manifold moving average measurement error

從表1中可以看出，隨著滑動平均長度的增加，校準誤差進一步減小。實際應用中，滑動平均長度100可以滿足載波相位差分精密相對定位的需要。與傳統(tǒng)幅相測量添加額外參考信號的方式相比，本文方法在性能上具有類似的精度，但是本文方法不需要假設額外的校準信號和測量儀器就能實現(xiàn)現(xiàn)場陣列流形矢量的快速測量。

表1 不同滑動平均長度下相位校準統(tǒng)計結果Tab.1 Statistics of different moving average length

5 結束語

本文首先建立了衛(wèi)導陣列天線接收信號模型，對天線和射頻前端引入衛(wèi)導接收信號的幅相不一致性進行了建模。在高精度差分相對定位應用中，相位不一致性破壞了載波相位測距的假設，將導致定位解算失敗。分析表明，控制數(shù)字波束形成權矢量的幅度和相位，將直接改變該波束數(shù)據(jù)中衛(wèi)導信號的幅度和相位，據(jù)此，提出了一種通過規(guī)劃校準權矢量，結合多通道數(shù)字接收機對同一顆衛(wèi)星信號進行測量，計算該衛(wèi)星的現(xiàn)場陣列流形矢量的校準方法。同時，結合常規(guī)的通道校準方法，對天線的凈現(xiàn)場陣列流形矢量進行分離，從而實現(xiàn)了一種易于實現(xiàn)的衛(wèi)導陣列天線設備的系統(tǒng)校準方法。與現(xiàn)有方法相比，本文重點研究了校準波束的設置，同等條件下提高了校準成功率。本文方法無需外部特定的校準信號，可以簡化系統(tǒng)結構和降低系統(tǒng)成本，適合工程應用。在進一步的系統(tǒng)實現(xiàn)中，還需要研究陣列天線的整體流形矢量的表示、與數(shù)字接收機的聯(lián)動數(shù)字波束形成、由于多徑導致接收機無法正常測量，以及什么樣的多徑會導致算法失效等問題。

：

[1] KAPLAN E D,HEGARTY C J. Understanding GPS:principles and applications[M].2nd ed. Norwood,MA:Artech House,2006:21-26.

[2] O′BRIEN A J,GUPTA I J,REDDY C J,et al.Space-time adaptive processing for mitigation of platform generated multipath[C]//Proceedings of ION 2010 International Technical Meeting. San Diego,CA:ION,2010:646-656.

[3] COMPTON R T. The relationship between tapped delay-line and FFT processing in adaptive arrays[J].IEEE Transactions on Antennas and Propagation,1988,36(1):15-26.

[4] VAN TREES H L. Optimum array processing:part IV of detection,estimation and modulation theory[M].New York:Wiley-Interscience,2002:1-14.

[5] GUPTA I J,WEISS I M,MORRISON A W. Desired features of adaptive antenna arrays for GNSS receivers[J].Proceedings of the IEEE,2016,104(6):1195-1206.

[6] ZHUANG X B,CUI X W,LU M Q,et al.Analysis of beamforming algorithms for GPS receivers in the jamming environment[J].System Engineering and Electronics,2009,31(3):570-573.

[7] CHURCH C M,GUPTA I J. Calibration of GNSS adaptive antennas[C]//Proceedings of 22nd International Meeting of the Satellite Division of The Institute of Navigation. Savannah,GA:ION,2009:344-350.

[8] XU H,CUI X,SHEN J,et al.A two-step beam-forming method based on carrier phases for GNSS adaptive array anti-jamming[C]//Proceedings of the 2016 International Technical Meeting of The Institute of Navigation. Monterey,CA: ION,2016: 793-804.

[9] 李亞麟,樊訊,胡波,等. 天線校準誤差建模及對開環(huán)波束賦形技術的影響[J].電訊技術,2010,50(3):1-4.

LI Yalin,FAN Xun,HU Bo,et al.Modeling of antenna calibration error and its impact on open-loop beamforming[J].Telecommunication Engneering,2010,50(3):1-4.(in Chinese)

[10] 曹俊鋒,許建文,方云. 接收DBF雷達系統(tǒng)校正技術[J].火控雷達技術,2007,36(4): 63-67.

CAO Junfeng,XU Jianwen,FANG Yun. Correction technology for radar system based on receiving DBF[J].Fire Control Radar Technology,2007,36(4): 63-67.(in Chinese)

[11] 劉宏偉,張守宏. 多饋源空間饋電陣列誤差校正[J].電子學報,1999,27(9):105-107.

LIU Hongwei,ZHANG Shouhong. Error calibration for multi-feed array[J].Acta Electronica Sinica,1999,27(9):105-107.(in Chinese)

[12] BORRE K,AKOS D M,BERTELSEN N. A software-defined GPS and Gallieo receiver:a single-frequency approach[M].Boston:Birkh?user,2007:143-150.