喬旭君 孫澤洲 吳學(xué)英

（北京空間飛行器總體設(shè)計部，北京 100094）

1 引言

深空探測不同于地球軌道航天器，由于星地距離遙遠，接收到的信號很微弱，數(shù)據(jù)信號傳輸碼速率低，對測速、測距的精度要求高［1］。為了滿足深空探測的指標要求，通常需要增大發(fā)射EIRP值，提高接收系統(tǒng)的G／T值［2］，同時也需要關(guān)注上行信號對于下行信號的影響。文中關(guān)注的是上行遙控信號對下行載波的影響，這種影響直接關(guān)系到下行載波的品質(zhì)和功率分配。

在統(tǒng)一載波系統(tǒng)中，上行遙控信號進入到下行載波中的殘留遙控信號來源于兩種途徑：一種是上行遙控信號經(jīng)過測距轉(zhuǎn)發(fā)通道對下行載波調(diào)制，進入到下行載波中；另一種是當應(yīng)答機工作在上下行載波相干模式時，上行殘留遙控信號調(diào)制鎖相環(huán)壓控振蕩器（VCXO）的輸出（下行載波基準頻率），經(jīng)倍頻進一步變大后進入到下行載波中。前一種途徑的影響比較大，一般在測距通道打開的情況下都會在鏈路計算中考慮；后一種途徑的影響比較小，但是隨著深空探測的發(fā)展，一些深空探測器根據(jù)探測任務(wù)，已經(jīng)不需要測距轉(zhuǎn)發(fā)，而是只接收上行數(shù)據(jù)信號（相位調(diào)制），下傳下行數(shù)據(jù)信號（二元移相鍵控調(diào)制），在這種情況下，來自第二種途徑的殘留遙控信號就值得關(guān)注。

2 上行遙控信號對下行信號的影響分析

當應(yīng)答機為模擬應(yīng)答機時，應(yīng)答機接收機為模擬鎖相接收機，其簡化鎖相環(huán)［3-4］和遙控副載波對下行載波的影響途徑如圖1所示。當應(yīng)答機工作在相干模式時，在X 頻段中，上下行載波頻率比為749∶880，因此，設(shè)749f0為上行載波頻率，4f0為下行載波頻率的基準頻率，8kHz為遙控副載波頻率。

圖1 殘留遙控信號對下行信號影響的途徑框圖Fig．1 Block diagram of the approach by which the uplink residual telecommand signal affects the downlink signal

由圖1可知，遙控信號可以由兩條途徑影響下行載波，一條是經(jīng)由測距轉(zhuǎn)發(fā)通道對下行載波進行調(diào)制，另一條是通過調(diào)制鎖相環(huán)VCXO 輸出信號，進入到下行載波中。

2．1 遙控信號影響下行載波的第一條途徑

遙控信號影響下行載波的第一條途徑，是經(jīng)由測距轉(zhuǎn)發(fā)通道對下行載波進行調(diào)制。上行測距信號經(jīng)過混頻和濾波之后，在星上應(yīng)答機中經(jīng)過側(cè)音解調(diào)，進入調(diào)制器調(diào)制到下行載波上。這個過程中，由于混頻器和濾波器不能達到理想的狀態(tài)，所以會有一部分噪聲隨測距信號轉(zhuǎn)發(fā)進入到下行信號中，這其中就包括遙控信號。遙控信號經(jīng)過這條途徑在下行載波中產(chǎn)生的影響比較大，所以一般在鏈路預(yù)算中都進行了考慮。

下行殘留遙控信號比殘留載波低的分貝數(shù)為

式中：α為下行信號中殘留遙控信號的調(diào)制指數(shù)；J0（α）為零階貝塞爾函數(shù)；J1（α）為一階貝塞爾函數(shù)；Δ為殘留遙控信號比殘留載波低的分貝數(shù)。

若上行遙控調(diào)制指數(shù)為0．95，下行殘留遙控調(diào)制指數(shù)為0．7，由式（1）可知，在對應(yīng)的下行載波中，殘留遙控信號比載波低約8．5dB。實際下行信號中殘留遙控信號與殘留載波的頻譜如圖2所示，圖中中心頻率處R1點（fc）為下行載波頻率，左右各8kHz（即R2點）是殘留遙控分量。

圖2 經(jīng)測距轉(zhuǎn)發(fā)通道進入下行載波的殘留遙控信號頻譜Fig．2 Spectrum of residual telecommand signal via ranging turn-around channel

2．2 遙控信號影響下行載波的第二條途徑

遙控信號影響下行載波的第二條途徑，是上行遙控信號經(jīng)過混頻、濾波，通過調(diào)制VCXO 輸出信號（下行載波基準頻率）而進入到下行信號中。在應(yīng)答機鎖定后，鎖相環(huán)路隨即將8kHz遙控副載波解調(diào)出來，遙控副載波在鑒相器后低通濾波器前輸出。同時，鑒相器的輸出經(jīng)環(huán)路低通濾波器送至VCXO壓控端。理論上，遙控副載波在經(jīng)過環(huán)路低通濾波器后會出現(xiàn)極大的衰減，但是由于遙控信號不能完全被抑制（根據(jù)一般鎖相環(huán)的經(jīng)驗值，殘留遙控信號被抑制約20dB），所以會有微弱的遙控副載波疊加到VCXO 的壓控端，對VCXO 構(gòu)成調(diào)制作用，導(dǎo)致在VCXO 輸出端4f0信號上存在遙控副載波調(diào)相信號。VCXO 輸出端的信號進一步經(jīng)過倍頻成為下行載波，而這4f0上的殘留遙控信號也會進一步放大，進而在下行載波中產(chǎn)生影響。

這一條途徑的影響，由于在VCXO 輸出端調(diào)制的殘留遙控信號相對于下行載波基準信號4f0較小，所以在測距通道打開的時候，由于第一種途徑的影響存在而沒有得到關(guān)注。但是在一些深空探測任務(wù)中，探測器不再進行測距轉(zhuǎn)發(fā)，在下傳下行數(shù)據(jù)信號時使用BPSK 調(diào)制方式，這時，就不存在第一條途徑的影響，第二條途徑產(chǎn)生的影響就需要得到關(guān)注，因為它直接影響到下行載波的品質(zhì)和功率分配，并影響下行信號的誤碼率。例如，當上行遙控信號調(diào)制指數(shù)為0．95時，接收機鎖相環(huán)VCXO 輸出端4f0信號中，實測4f0信號與8kHz副載波信號幅度差為-77．84dB，經(jīng)過倍頻后，下行信號中，殘留遙控信號比載波低30．57dB，分別如圖3 和圖4 所示。圖3中，4f0頻率是鎖相環(huán)VCXO 輸出端下行載波基準頻率。圖4中，fc是下行載波頻率。

殘留遙控信號經(jīng)過第二條路徑產(chǎn)生的影響，在一些深空探測任務(wù)中成為必須要關(guān)注的問題，我們希望可以進一步減小下行載波中殘留遙控信號的影響，減小其占去的功率，降低信號的誤碼率，從而優(yōu)化深空測控系統(tǒng)。

圖3 經(jīng)殘留遙控信號調(diào)制后的鎖相環(huán)VCXO 輸出信號頻譜Fig．3 Spectrum of output of VCXO modulated by residual telecommand signal

圖4 經(jīng)調(diào)制VCXO 輸出信號進入下行載波的殘留遙控信號頻譜Fig．4 Spectrum of residual telecommand signal via modulating the output of VCXO

3 進一步抑制殘留遙控信號的方法

針對遙控信號通過調(diào)制鎖相環(huán)VCXO 輸出信號（下行載波基準信號）而進入到下行信號中的問題，有兩種解決方案可供考慮。

第一種方案是對圖1 中鎖相環(huán)中的低通濾波器的相關(guān)參數(shù)進行更改，加強其對殘留遙控信號的抑制，使得下行信號中的殘留遙控信號符合指標要求。但是由于這里鎖相環(huán)使用的是無限長單位沖激響應(yīng)（IIR）濾波器，可調(diào)范圍小，調(diào)整后對遙控信號的抑制作用變化不大，而且一般都選用固定的器件進行設(shè)計，一旦改變某一個參數(shù)，鎖相環(huán)其它參數(shù)的選取都會受到影響［5］，所以這種方案并不可行。

第二種方案是在圖1 中調(diào)制器和鎖相環(huán)的VCXO 之間增加帶通濾波器，使用帶通濾波器［6］對VCXO 輸出端4f0信號濾波，濾除4f0頻點之外的信號，只保留中心頻率，這樣可以抑制應(yīng)答機下行信號中的殘留遙控副載波信號，這種方法更易實現(xiàn)。但是當應(yīng)答機上行遙控副載波信號調(diào)有數(shù)據(jù)后，8kHz副載波信號的頻譜會展開，在頻譜上看，不再是單一的8kHz的譜線，而是擴展到了整個帶寬范圍，所以在VCXO 輸出端使用帶通濾波器后，只能部分抑制應(yīng)答機下行載波中的殘留遙控信號，靠近下行載波頻率附近的殘留遙控信號不能得到抑制，可以通過仿真的方法觀察這種方法的改善效果。

4 仿真分析

針對本文所采用的在鎖相環(huán)VCXO 輸出端增加帶通濾波器濾除下行信號中殘留遙控信號的方案，由于在遙控副載波調(diào)有數(shù)據(jù)后，靠近下行載波頻率附近的殘留遙控信號不能得到抑制，這里通過仿真的方法，可觀察方案的改善效果。

仿真的方法是首先仿真出VCXO 輸出端調(diào)制有殘留遙控的4f0信號，然后將此信號經(jīng)過仿真的帶通濾波器，最后得到濾波后的頻譜圖，對比濾波前后頻譜圖，分析帶通濾波器對殘留遙控信號抑制的效果。

4．1 VCXO 輸出端信號仿真

對信號頻譜仿真，一般的思路是按照信號的時域表達式生成信號，然后用傅里葉變換相關(guān)指令生成頻譜圖。

當生成VCXO 輸出端信號時，需要生成遙控碼信號。要使VCXO 輸出端信號里包含1024個遙控碼，若使用循環(huán)嵌套判斷語句則需要完成內(nèi)層十萬次以上，外層上千次的循環(huán)運算。由于MATLAB軟件執(zhí)行循環(huán)語句的效率不高，如果在MATLAB軟件中使用循環(huán)語句生成仿真的信號，將使整個頻譜仿真時間巨大。

雖然MATLAB軟件的循環(huán)運算能力弱，但是其矩陣運算能力強，所以應(yīng)該盡可能地使用矩陣運算，避免使用循環(huán)結(jié)構(gòu)和判斷語句，以期在有限的內(nèi)存空間和盡量短的時間內(nèi)完成運算，這就需要在仿真前做充足的矩陣設(shè)計，而不能利用程序去執(zhí)行循環(huán)、判斷指令。

信號仿真的具體步驟是：

（1）假設(shè)鎖相環(huán)VCXO 輸出端信號中心頻率4f0為40 MHz，以殘留遙控副載波頻率為8kHz，遙控碼速率為1000bit／s為例進行仿真。采樣頻率取中心頻率的4 倍，為160 MHz，是遙控碼速率的160 000倍，即每一個遙控碼被采樣160 000次。當以采樣1024個遙控碼來做仿真時，需要的采樣點數(shù)是1．638 4×108個，這是后期需要處理的總的序列長度。

（2）對于遙控碼，首先隨機產(chǎn)生1024個碼值的矩陣（1024×1），然后使用矩陣變換指令使該矩陣重復(fù)160 000列，生成被采樣后的遙控碼矩陣（1024×160 000）。

（3）對遙控副載波，首先在一個遙控碼周期內(nèi)按采樣頻率采樣，生成一個遙控碼周期內(nèi)的采樣值矩陣（1×160 000），然后將這些采樣值重復(fù)1024行，形成與遙控碼對應(yīng)的遙控副載波的采樣矩陣。

（4）將遙控碼采樣值矩陣與遙控副載波采樣值矩陣對應(yīng)相乘，并按順序重新生成單行的采樣值序列，序列長度是163 840 000。

（5）將遙控副載波的采樣序列加到載波上（中心頻率40 MHz），并增加高斯白噪聲，生成VCXO 輸出端信號。

（6）對VCXO 輸出端信號進行傅里葉相關(guān)運算［7］，得到VCXO 輸出端信號頻譜如圖5所示。

以上程序設(shè)計方法摒棄了傳統(tǒng)的使用循環(huán)語句形成序列的方法，充分利用MATLAB軟件進行矩陣運算，可以在較大數(shù)據(jù)量的處理中有效節(jié)省運算時間。

圖5 4f0信號濾波前頻譜Fig．5 Spectrum of signal 4f0before the band-pass filter

4．2 濾波處理的仿真

本文根據(jù)實際帶通濾波器的參數(shù)進行了濾波器的仿真，然后將4．1中仿真的VCXO 輸出端信號經(jīng)過此濾波器濾波。

濾波仿真的具體步驟是：

（1）使用MATLAB軟件中的fdatool濾波器設(shè)計分析工具進行濾波器仿真［8］。

輸入濾波器指標：3dB 帶寬≥2．8kHz、插入損耗≤6dB、帶內(nèi)波動≤1dB、40dB 阻帶帶寬≤±8kHz。

輸出濾波器的幅頻相應(yīng)曲線，并得到巴特沃斯濾波器設(shè)計系數(shù)（SOS）的矩陣RSOS。濾波器的幅頻響應(yīng)曲線如圖6所示，濾波器設(shè)計系數(shù)的矩陣RSOS如式（2）所示。

圖6 濾波器幅頻響應(yīng)曲線Fig．6 Magnitude response curve of the filter

（2）將SOS矩陣［9-10］導(dǎo)出到工作空間中。

（3）對VCXO 輸出端信號進行濾波處理。

（4）對濾波后的VCXO 輸出端信號進行傅里葉變換，得到濾波后信號頻譜圖如圖7所示。

圖7 4f0信號濾波后頻譜Fig．7 Spectrum of signal 4f0after the bandpass filter

4．3 仿真結(jié)果分析

對比VCXO 輸出端信號在增加帶通濾波器前后的頻譜圖，由圖5 可知，VCXO 輸出端信號濾波前中心頻率4 0MHz左右各8kHz為殘留遙控副載波頻率；由圖7 可知，VCXO 輸出端信號濾波后中心頻率40MHz左右各8kHz的殘留遙控副載波得到了有效抑制。

經(jīng)過對序列的分析，可知f＝3．999 2×107Hz為殘留遙控信號所在頻率，對濾波前后歸一化功率譜進行分析，可知在此頻率處，濾波前殘留遙控信號比中心頻率低31 dB，濾波后低49 dB，說明在VCXO 輸出端增加帶通濾波器后可以抑制殘留遙控信號約18dB。

需要說明的是，實際殘留遙控信號調(diào)制VCXO輸出后，4f0信號與8kHz 副載波信號幅度差為-77．84dB。對VCXO 輸出端信號進行濾波仿真，由于這時在VCXO 輸出端的殘留遙控信號影響已經(jīng)很?。?77．84dB），而MATLAB 軟件存在仿真缺陷，無法用仿真濾波器對這么小的殘留遙控信號進一步抑制。所以本文對VCXO 輸出端信號仿真，是以4f0信號與8kHz 副載波信號幅度差約為-30dB為例進行的。

在實際工程中，可以通過在VCXO 輸出端增加前述的晶體帶通濾波器，使下行信號中的殘留遙控信號進一步被抑制約20dB（從約-78dB抑制到約-98dB），使4f0信號上調(diào)制的殘留遙控信號經(jīng)220倍頻后（增加約47dB），比載波低約51dB，因此，下行信號中殘留的遙控信號滿足指標要求。實際改善效果如圖8所示。

圖8 VCXO 輸出端增加晶體帶通濾波器后效果圖Fig．8 Rendering of the output of VCXO with a crystal filter

5 結(jié)論

在實際工程中，通過在鎖相環(huán)VCXO 輸出端增加晶體帶通濾波器，可以使下行信號中殘留遙控信號比載波低50dB以上，進一步提高了載波的品質(zhì)，改善了下行信號的功率分配。

（References）

［1］吳學(xué)英，孫澤洲，金蓬嘉．深空探測下行測距調(diào)制度變化的現(xiàn)象和機理分析［J］．中國空間科學(xué)技術(shù)，2007，27（6）：48-53

Wu Xueying，Sun Zezhou，Jin pengjia．Discussion about the variance of the downlink ranging modulation index in deep space［J］．Chinese Space Science and Technology，2007，27（6）：48-53（in Chinese）

［2］約瑟夫·袁．深空遠程通信系統(tǒng)工程［M］．吳偉仁，譯．北京：《遙測遙控》編輯部，1988

Yuen J H．Deep space telecommunication system engineering［M］．Wu Weiren，translated．Beijing：Telemetry and Telecommand Editorial Department，1988 （in Chinese）

［3］葉培建，饒煒，孫澤洲，等．嫦娥一號月球探測衛(wèi)星技術(shù)特點分析［J］．航天器工程，2008，17（1）：7-11

Ye Peijian，Rao Wei，Sun Zezhou，et al．Technology characters analysis of Chang’e-1［J］．Spacecraft Engineering，2008，17（1）：7-11（in Chinese）

［4］陳芳允，賈乃華．衛(wèi)星測控手冊［M］．北京：科學(xué)出版社，1992年

Chen Fangyun，Jia Naihua．Handbook of satellite telemetry and telecommand［M］．Beijing：Science Press，1992（in Chinese）

［5］鄭繼禹，萬心平，張厥盛．鎖相環(huán)路原理與應(yīng)用［M］．北京：人民郵電出版社，1976

Zheng Jiyu，Wan Xinping，Zhang Juesheng．Phase lock loop theory and application［M］．Beijing：Post and Telecom Press，1976（in Chinese）

［6］加德納．鎖相環(huán)技術(shù)［M］．姚劍清，譯．北京：人民郵電出版社，2007

Gardner F M．Phaselock techniques［M］．Yao Jianqing，translated．Beijing：Posts ＆ Telecommunications Press，2007（in Chinese）

［7］程佩青．數(shù)字信號處理教程［M］．北京：清華大學(xué)出版社，2008

Cheng Peiqing．Digital signal processing tutorial［M］．Beijing：Tsinghua University Press，2008（in Chinese）

［8］周建興，豈興明，矯津毅，等．MATLAB從入門到精通［M］．北京：人民郵電出版社，2010

Zhou Jianxing，Qi Xingming，Jiao Jinyi，et al．MATLAB：from introduction to mastery［M］．Beijing：Posts＆Telecommunications Press，2010（in Chinese）

［9］王彬，于丹，汪洋．MATLAB數(shù)字信號處理［M］．北京：機械工業(yè)出版社，2010

Wang Bin，Yu Dan，Wang Yang．MATLAB digital signal processing［M］．Beijing：China Machine Press，2010（in Chinese）

［10］殷瑞，萬國龍．數(shù)字信號處理［M］．北京：清華大學(xué)出版社，2007

Yin Rui，Wan Guolong．Digital signal processing［M］．Beijing：Tsinghua University Press，2007（in Chinese）