喬旭君 孫澤洲 吳學英
(北京空間飛行器總體設(shè)計部,北京 100094)
深空探測不同于地球軌道航天器,由于星地距離遙遠,接收到的信號很微弱,數(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)注。
當應(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 輸出信號,進入到下行載波中。
遙控信號影響下行載波的第一條途徑,是經(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
遙控信號影響下行載波的第二條途徑,是上行遙控信號經(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)制作用,導致在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
針對遙控信號通過調(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)答機下行載波中的殘留遙控信號,靠近下行載波頻率附近的殘留遙控信號不能得到抑制,可以通過仿真的方法觀察這種方法的改善效果。
針對本文所采用的在鎖相環(huán)VCXO 輸出端增加帶通濾波器濾除下行信號中殘留遙控信號的方案,由于在遙控副載波調(diào)有數(shù)據(jù)后,靠近下行載波頻率附近的殘留遙控信號不能得到抑制,這里通過仿真的方法,可觀察方案的改善效果。
仿真的方法是首先仿真出VCXO 輸出端調(diào)制有殘留遙控的4f0信號,然后將此信號經(jīng)過仿真的帶通濾波器,最后得到濾波后的頻譜圖,對比濾波前后頻譜圖,分析帶通濾波器對殘留遙控信號抑制的效果。
對信號頻譜仿真,一般的思路是按照信號的時域表達式生成信號,然后用傅里葉變換相關(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
本文根據(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]導出到工作空間中。
(3)對VCXO 輸出端信號進行濾波處理。
(4)對濾波后的VCXO 輸出端信號進行傅里葉變換,得到濾波后信號頻譜圖如圖7所示。
圖7 4f0信號濾波后頻譜Fig.7 Spectrum of signal 4f0after the bandpass filter
對比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
在實際工程中,通過在鎖相環(huán)VCXO 輸出端增加晶體帶通濾波器,可以使下行信號中殘留遙控信號比載波低50dB以上,進一步提高了載波的品質(zhì),改善了下行信號的功率分配。
(References)
[1]吳學英,孫澤洲,金蓬嘉.深空探測下行測距調(diào)制度變化的現(xiàn)象和機理分析[J].中國空間科學技術(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].北京:科學出版社,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].北京:清華大學出版社,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].北京:清華大學出版社,2007
Yin Rui,Wan Guolong.Digital signal processing[M].Beijing:Tsinghua University Press,2007(in Chinese)