朱亞立，鄭為民，童 力，童鋒賢，張 娟，劉 磊

(1. 中國科學院上海天文臺，上海 200030；2. 中國科學院射電天文重點實驗室，江蘇 南京 210008； 3. 上海市導航定位重點實驗室，上海 200030；4. 中國科學院大學，北京 100049)

甚長基線干涉測量[1](Very Long Baseline Interferometry, VLBI)是20世紀60年代后期出現的一種射電天文技術。它具有空間分辨率高、測量精度高的特點，現已廣泛應用于天體物理、天體測量、天文地球動力學和深空探測領域[2]。我國于2007年首次在探月工程中使用VLBI技術，建立了中國VLBI網(Chinese VLBI Network, CVN)，包括上海佘山站、北京密云站、云南昆明站、烏魯木齊南山站和上海VLBI數據處理中心。中國VLBI網已參加了探月工程嫦娥系列探測器的多次任務，對關鍵弧段的快速精密測定軌發(fā)揮了不可替代的重要作用[3-5]。

VLBI測站由天線、接收機、數據采集終端、時間頻率系統等組成，工作時，各測站接收機天線匯集的射頻信號放大后，經獨立本振進行混頻得到中頻信號，然后傳送至數字化終端，高速模數采樣后作數字基帶轉換，得到所需頻率通道的基頻信號，最后經過編碼，記錄在磁盤或者直接通過網絡傳送到VLBI中心[5]。

為研究VLBI中心對不同探測器的測量方法，評估分析中心各功能模塊的能力及整個系統的正確性、可靠性，目前常規(guī)做法是在每次飛行任務前開展大量實際目標聯測。對單目標探測器的實際觀測演練可評估測站與系統的可靠性。但對于那些全新的深空探測任務，在正式任務前有時難以找到合適的探測器作為試驗對象。比如在后續(xù)探月工程中，將要開展針對月球軌道的交會對接，屆時需要對月球距離的兩個動態(tài)目標進行同波束測量。顯然，目前沒有類似的實際探測器可作為動態(tài)雙目標的觀測對象。若采用探測器VLBI信號仿真方法，則可根據探測器的數量、軌道、信標等特點，模擬生成VLBI終端接收的不同探測器信號，提供月球軌道雙目標同波束VLBI測量數據，用于研究和改進VLBI處理模塊。通過仿真與實際演練兩種方法結合，可較好地評估VLBI中心與整個系統對雙(多)目標探測器的測量能力。因此，研究VLBI信號仿真技術有重要意義與迫切工程需求。

在探月工程中，探測器下行VLBI射頻信號為S/X頻段，采用調相方式，調制信號在形式上主要包括點頻信號和有限帶寬信號兩類。在仿真信號時，沒有直接仿真射頻信號，而是先將信號從射頻變換至數字中頻，然后再對中頻信號采樣，提取通道信號，最終生成用于驗證VLBI數據處理系統的仿真數據，這樣既完整地保留了調制信號頻譜的主要特性，又避免了直接處理需極高采樣率的射頻信號。

1 探測器信號結構

在探月工程中，探測器VLBI下行射頻信號采用調相方式，包括遙測信號、單向差分測距(Differential One-way Ranging, DOR)信號、測距信號、殘留遙控信號(圖1)。其中，下行載波、測距信號(測距主音和測距次音)和單向差分測距信號為點頻信號；遙測信號和遙控信號為經過相位調制的有限帶寬信號。這些信號統一調制在下行信號的主載波上，數學模型[5]為

(1)

其中，PT為總功率；fc為載波頻率；d1(t)為遙測信號；θ1為遙測信號的調制指數；d2(t)為遙控信號；θ2為遙控信號的調制指數；f1為測距信號；φ1為測距信號的調制指數；f2為單向差分測距信號；φ2為單向差分測距信號的調制指數。

圖1 X波段下行信號上邊帶頻譜示意圖
Fig.1 The spectrum of X band downlink signal

2 仿真原理與實現方法

2.1 探測器信號時延模型構造

VLBI測站與探測器之間存在相對運動，因此在仿真時要考慮相對運動效應引起的測站接收探測器信號的時間延遲(時延)及其變化[5]。時延包含幾何時延、引力時延、天線時延修正、源結構時延修正、中性大氣與電離層時延修正等[6]，其中，幾何時延是理論時延的主要部分，隨目標的運動而變化。

設地心O、測站A及探測器B(圖2)，g0(t)表示t時刻參考點(地心)接收到信號波前的光行時；g1(t)表示t時刻參考點(地心)信號波前超前到達測站的光行時；g2(t)表示t時刻參考點(測站)接收到信號波前的光行時。若t2時刻從探測器B發(fā)出信號，同一波前到達地心O、測站A的時刻分別為t0,t1，則根據定義，g0(t0)，g1(t0)及g2(t1)分別為

g0(t0)=t0-t2，

g1(t0)=t1-t0，

g2(t1)=t1-t2，

(2)

顯然，

g2(t1)=g0(t0)+g1(t0).

(3)

在實際應用中，每個觀測弧段的時延模型都以五次時延多項式表示：

g(t)=b0+b1t+b2t2+b3t3+b4t4+b5t5，

(4)

其中，b0,b1, ...,b5為時延模型多項式系數。

根據實際時延模型g0(t),g1(t)，用最小二乘法作五次多項式擬合，得到時延模型g2(t)。

2.2 時延模型線性近似

為減少計算量，在短時間內假設時延模型隨時間線性變化，即

(5)

其中，系數b0,b1在短時間內為恒常數，忽略時延模型g2(t)高次項。

圖2 VLBI時延示意圖
Fig.2 The diagram of VLBI delay

在短時間內用(5)式表示線性化時延模型存在一定誤差，誤差大小取決于用線性近似的時間長度和五次時延多項式中高次項系數。相關處理機[7]使用皮秒量級的時延模型精度足以滿足應用需求，故誤差值在0.01 ps(1 e-14 s)以內的模型誤差即可視為滿足仿真精度要求。在較短的觀測時間內，時延模型g2(t)的二次項是線性近似的主要誤差項，其系數b2一般為e-11量級，要想使得該項引起的時延誤差小于0.01 ps，選定的線性化時間長度應小于0.01 s。其他高次項的系數b3,b4,b5，一般在e-14量級，由此引起的誤差可忽略。

2.3 構造中頻信號

探測器下行載波信號的類型主要包括點頻信號和有限帶寬信號。其中點頻信號包括：載波、測距信號和單向差分測距信號。有限帶寬信號包括遙測信號和遙控信號，是以正弦波為載波的調相信號，調制方式為PCM/PSK/PM，一般可表示為sin(2πft+θ)，用相位θ的變化傳遞數字信息。對于二進制相移鍵控信號，θ以0/π交替變化表示二進制位0/1， 1對應相位0，-1對應相位π。因此遙測和遙控信號又可表示為Asin(2πft)，A隨相位θ的變化而變化。

為完整地保留調制信號及其邊帶頻譜，選擇中頻信號約為60 MHz。有限帶寬信號的碼速率最大為2 048 bit/s，碼元變化周期是采樣周期的兩萬九千多倍。在一個碼元周期內，有限帶寬信號保持不變，對其采樣可仿真遙測和遙控有限帶寬信號。因此，需選擇合適的時間尺度對點頻信號和有限帶寬信號采樣，仿真中頻信號。

為簡約起見，根據(1)式探測器轉發(fā)的下行信號的一般形式，分析其中典型的點頻信號、有限帶寬信號多普勒效應影響。信號表達式如下：

f(t)=sin[2πfct+θ1Asin(2πf3t)+φ2sin(2πf2t)]，

(6)

其中，fc,f2分別為載波和點頻單向差分測距信號；f3為有限帶寬遙測信號；A為遙測信號的幅值；θ1,φ2為調制系數。

(7)

將(5)式和(6)式代入(7)式，得到：

(8)

根據(8)式，探測器下行信號在頻域上表現為關于載波左右對稱，其頻譜主要分布在[(1-b1)(fc-fmax), (1-b1)(fc+fmax)]，其中fmax為max(f3,f2)，b1?1。

(9)

其中，fd為信號起始頻點(1-b1)(fc-fmax)。

對比(9)式與(8)式，只是載波頻率發(fā)生變化，其余信號不變。即通過降低載波頻率，將信號從射頻搬移到中頻，且信號f3,f2前后保持不變。

綜上，VLBI測站接收信號是添加時延的射頻信號，由(1)式添加時延并降低載波頻率，得到包含多普勒效應的中頻信號，形式如下：

(10)

其中，f4為遙控信號；A′為遙控信號幅值。

2.4 VLBI終端信號仿真

圖3顯示了仿真VLBI終端接收信號實現流程。先仿真帶有多普勒效應的探測器VLBI下行中頻信號，再濾波得到某個通道的頻率區(qū)間，通過希爾伯特變換保留正頻率部分，再通過數字變頻移頻至基帶，經希爾伯特逆變換取實信號，最后通過重采樣降低采樣率得到通道實信號。通道寬度為2 MHz時，采樣率取4 MHz。

3 數值試驗

3.1 仿真探測器下行中頻信號

根據上述方法采用某月球探測器信標實際參數仿真了探測器VLBI下行中頻信號，其中包含測距信號(500 kHz主音頻率)、遙測信號(副載波65.536 kHz且碼速率256 bit/s)和兩對單向差分測距信號。高頻信號起始頻點fd為8 460 MHz，降低載波頻率后生成探測器下行中頻信號，采樣率117 MHz，快速傅里葉變換點165 000。

圖4顯示帶寬58.5 MHz的探測器下行中頻信號頻譜。圖中，箭頭3為載波信號，2和4是距載波約3.85 MHz的第1對單向差分測距信號(DOR1)；1和5是距載波約19.27 MHz的第2對單向差分測距信號(DOR2)，載波兩側500 kHz處是測距主音信號。另外，有限帶寬的遙測信號(帶寬約0.5 kHz)靠近載波。

圖3 信號仿真流程圖
Fig.3 Steps of signal simulation

測距信號、遙測信號、單向差分測距信號關于載波對稱，調制系數決定各信號的功率分配，其中載波、測距信號和單向差分測距信號較為明顯。由于快速傅里葉變換的柵欄效應使圖中載波兩側對應頻點幅值產生細微差別。其他信號諧波項能量較小，幅度基本在-180 dB以下。

圖4 探測器下行中頻信號頻譜圖
Fig.4 The spectrum of downlink probe quasi-intermediate-frequency signal

3.2 可行性驗證

通過仿真信號與VLBI終端接收信號的對比以及對VLBI相關處理結果的分析，驗證了探測器VLBI信號仿真方法的可行性。

圖5(a)為實際信號單通道自功率譜，從中可觀察到0.96 MHz處的載波信號；距離載波500 kHz、65.536 kHz的測距主音信號和遙測信號；其他為信號諧波項。

圖5(b)為仿真VLBI終端接收信號自功率譜(采樣率4 MHz，快速傅里葉變換點262144，采樣1 s)，包含載波信號、測距信號(500 kHz主音頻率)、遙測信號(副載波65.536 kHz且碼速率1 024 bit/s)、高斯白噪聲信號。頻率0.968 MHz處為載波信號，測距主音信號和有限帶寬的遙測信號分別距離載波500 kHz、65.536 kHz，關于載波對稱。距離測距主音65.536 kHz、131.072 kHz且關于測距主音對稱的為遙測信號諧波項。

圖5 單通道自功率譜圖 (a)實際信號；(b)仿真信號
Fig.5 The auto-power spectrum of channel signal (a) The auto-power spectrum of real signal; (b) The auto-power spectrum of simulation signal

仿真信號與對應的實際信號對比表明，載波信號、測距信號、遙測信號的頻率和功率分配與實際信號一致，噪聲信號的功率和分布也與實際信號相同，驗證了仿真方法的可行性。

探測器VLBI下行中頻信號分別添加不同時延模型，仿真兩臺站終端接收信號(采樣率4 MHz，快速傅里葉變換點512，采樣1 s)。經過VLBI相關處理，可得互功率譜幅度和相位圖(圖6)。圖6(a)的幅度譜顯示，載波信號和遙測信號在[0.8 MHz，1.2 MHz]頻段，遙測信號的諧波項和測距信號在[0.4 MHz, 0.6 MHz]，[1.4 MHz, 1.6 MHz]頻段。圖6(b)的相位譜顯示，在載波信號、遙測信號、測距信號分布的主要頻段，相位譜有明顯相關條紋，其它頻段的點由于信號較弱，相位分布彌散，條紋較弱，驗證了不同臺站接收的探測器VLBI下行仿真信號的相關性。

圖6 仿真兩站信號互功率譜圖 (a)幅度譜；(b)相位譜與條紋
Fig.6 The cross-power spectrum of two station simulation signal. (a) The amplitude spectrum; (b) The phase spectrum and strip

4 結 論

本文在研究探測器信號結構的基礎上，分析了包含探測器下行點頻信號和有限帶寬信號的VLBI終端信號仿真方法，構造了中頻信號仿真VLBI數字終端接收信號。采用最小二乘法可構造出滿足精度的時延模型，用于模擬多普勒效應。通過仿真數字終端信號與實際信號的比較以及對不同臺站的仿真信號作相關處理，驗證了探測器VLBI信號仿真方法的可行性。后續(xù)建立仿真各模塊誤差分析模型，完善仿真系統，為月球探測器VLBI信號仿真的應用，乃至火星探測器VLBI信號仿真的研究打下良好的基礎。

[1] Thompson A R, Moran J M, Swenson Jr G W. Interformetry and synthesis in radio astronomy[M]. Weinheim: Wiley Vch & Sons, 2004: 33-37.

[2] 魏杰昌, 鄭為民, 童力. 帶有時延的VLBI數字基帶信號仿真[J]. 中國科學院上海天文臺年刊, 2015(36): 124-135.

Wei Jiechang, Zheng Weimin, Tong Li. Simulation of VLBI digital baseband signal with time delay[J]. Annals of Shanghai Astronomical Observatory Chinese Academy of Sciences, 2015(36): 124-135.

[3] 江悟, 沈志強, 舒逢春, 等. 利用DiFX處理CVN天文觀測數據進展[J]. 天文學進展, 2014, 32(4): 516-524.

Jiang Wu, Shen Zhiqiang, Shu Fengchun, et al. The progress of adapting DiFX correlator for Chinese VLBI network[J]. Progress in Astronomy, 2014, 32(4): 516-524.

[4] 鄭鑫, 劉慶會, 吳亞軍, 等. 基于同波束VLBI差分相延遲的 “玉兔” 月球車動作監(jiān)視分析[J]. 中國科學: 物理學 力學 天文學, 2014, 44(8): 872-878.

Zheng Xin, Liu Qinghui, Wu Yajun, et al. Motion monitoring and analysis of Chang′E-3 rover based on same-beam VLBI differential phase delay[J]. Science Sinica: Physica, Mechanica & Astronomica, 2014, 44(8): 872-878.

[5] 錢志瀚, 李金嶺. 甚長基線干涉測量技術在深空探測中的應用[M]. 北京: 中國科學技術出版社, 2012: 1-14.

[6] 童鋒賢. 深空探測器VLBI相位參考測量關鍵技術研究[D]. 北京: 中國科學院大學, 2016.

[7] 陳中, 鄭為民. VLBI軟件相關處理機現狀和發(fā)展趨勢[J]. 天文學進展, 2015, 33(4): 489-505.

Chen Zhong, Zheng Weimin. Current status and development of VLBI software correlator[J]. Progress in Astronomy, 2015, 33(4): 489-505.