航天器地面無線測試鏈路中斷現象分析

安金坤，田 林，趙建賀

（中國空間技術研究院 載人航天總體部，北京 100094）

0 引言

航天器電測試是驗證航天器上設備電氣功能和性能的重要手段，在航天器研制流程中起到質量總檢驗的作用，可有效提高航天器的匹配性、兼容性和可靠性，是確保航天器完成任務的重要保證。為檢驗航天器無線傳輸鏈路中收發天線的性能，確定發射場階段整流罩對地面天線接收的影響等，需要對航天器進行無線測試。不同于有線測試，無線測試時信號傳輸鏈路較為復雜，測試現場會存在多徑傳播、阻擋物的陰影效應，以及人員移動、支架展開等引發的微多普勒頻移等。以上這些可能引發測試鏈路接收端信噪比下降，產生誤碼，甚至測試鏈路的中斷，從而影響測試效果和對航天器狀態的 判斷。尤其對于人員移動、支架展開等偶發動態因素，測試現場很難采取及時有效的控制措施。

本文即針對無線測試過程中的偶發動態因素造成測試鏈路中斷現象的作用機理進行著重分析，指出無線鏈路中微小的多普勒頻移即能引起接收端信號幅度的衰減和隨機調頻；為了對抗鏈路中的微多普勒頻移，將天線分集技術引入航天器無線測試領域，以克服信道非理想性。

1 無線信道模型

信道是通信系統中影響傳輸容量的決定因素，無線移動信道是一種隨時間和環境的變化而變化的傳輸媒介，電磁波會受到周圍物體的反射、折射和繞射產生多徑傳播，收發端相對運動會在接收端產生非常復雜的合成干涉信號。通常將電波傳播中的衰落分成大尺度衰落和小尺度衰落[1-4]。前者是指信號在大距離空間上的變化趨勢，主要考慮信號傳播距離、環境遮擋、傳播頻率等因素；后者是指在小距離上接收信號的強度變化情況，主要是由不同到達路徑上信號的相互疊加產生干涉引起接收信號的快速變化。本文主要研究航天器在地面無線測試中的電波傳播情況，傳播距離較短，只考慮后者的影響。本節從時域和頻域的角度介紹多徑傳播和多普勒存在的條件下無線信道的特點，以明確本文研究的著眼點。

1.1 時延色散

時延色散是指由于無線信道存在多條傳輸路徑，在時域發射一個脈沖串后在接收端會收到同一信號、不同延時的多個脈沖。即信道的沖激響應不是δ函數，造成信號的持續時間比發送的信號時長要長。時延色散導致接收信號的頻率選擇性，反映了無線多徑信道的靜態特性。在航天器無線測試中，通過測試前合理選擇接收天線的位置，避開陰影以及信道不理想的空間位置，可以建立較為理想的測試鏈路。因此，時延色散不作為本文的研究對象。

1.2 頻域色散

頻域色散是指由于收發端以及周圍環境的相對變化引入了多普勒效應，造成鏈路的傳輸函數是一個時變函數，使接收信號的頻譜范圍比發射信號的帶寬要寬。頻域色散導致接收信號的時間 選擇性，反映了無線多徑信道的動態特性。在航天器無線測試中，由于現場人員或者物體的短時移動、支架的展開等造成經過移動物體反射或者散射后進入接收機的信號存在多普勒頻移，使得接收端信號靈敏度下降，甚至鏈路中斷。本文即針對微小多普勒頻移[5-7]引發航天器地面無線測試鏈路中斷機理進行分析。

2 鏈路中斷機理分析

微多普勒頻移對接收信號的包絡和頻率（相位）都會產生影響，下面分別進行分析。

2.1 包絡的影響

以圖1所示的兩徑信道模型為例，發射端發射的是載頻為f的單載波信號，接收端接收到來自2 條路徑的信號，其中a 路徑為直射信號，b路徑為經過移動的人反射或者散射后進入接收機的信號。

圖1 無線信道兩徑模型示意圖Fig.1 Two-ray model of wireless channel

b 路徑存在多普勒頻移f1，a 路徑則不存在多普勒頻移。設歸一化a 路徑的信號幅度為1，b路徑與a 路徑信號幅度之比為b；θ1和θ2分別為a 路徑和b 路徑信號的初始相位，則接收端收到的合成信號為

式中：Re 表示取實部。令A=e(jθ1)+be(jθ2)e(j2πf1t)，則式(1)可以簡化為

因此，復信號A就是接收合成信號的等效低通信號，反映接收信號的強度信息，其模值|A|就是合成信號的包絡：

從式(3)可以看到，由于引入了多普勒頻移f1，接收合成信號的包絡產生周期性的起伏，在極小值附近會產生衰落深陷序列，測試無線鏈路可能會中斷。衰落深陷的間隔為多普勒頻移的倒數，即1/f1。實際信號可能會同時存在多條反射路徑，且每條路徑的多普勒頻移都不一樣，因此衰落深陷的概率比兩徑模型的還高，包絡起伏更加劇烈。

若再增加第3 條路徑，假定其為ccos(2πft+ 2πf2t+θ3)，其中，f2為該條路徑的多普勒頻移；c為該路徑與直射路徑幅度比；θ3為該路徑的初始相位。設θ為前2 條路徑的合成初始相位，則可從式(3)遞推出存在3 條路徑時合成信號的包絡為

可見存在3 條路徑時包絡變化的頻率分量更加豐 富，導致接收信號處于衰落深陷的概率增加，同時導致接收信號頻譜的進一步擴展。

以上為對確定信號進行分析的結果，對于實際隨機信號，可采用統計的方法用電平通過率來衡量微多普勒頻移對測試鏈路的影響。電平通過率定義為接收場強在正方向穿過某一電平r的通過率的期望值，可以寫成

式中：r′=dr/dt是r對時間的導數；pdfr,r’是r和r′的聯合概率密度函數；?0和?2分別為多普勒的0階和二階矩陣，n階矩陣定義為

式(6)中，多普勒頻譜的定義為

式中：vmax為最大速度；G(γ)是天線加權模式；γ為入射角度。

2.2 頻率的影響

從接收信號包絡的分析可以定性地看到，含多普勒頻移的多徑分量造成接收合成信號的包絡產生變化，引起接收信號瞬時頻率的展寬。

多徑傳播時接收合成信號r(t)的瞬時頻率定義為

式中：Im 表示取虛部。從式(8)可以看到接收信號的瞬時頻率和信號的包絡|r(t)|成反比，當接收信號的包絡較小時，瞬時頻率會較大。

當接收信號為一隨機信號時，假定接收信號r(t)的同相分量x和正交分量y滿足高斯分布，分別記它們的導數為x′和y′。則一個包含x、y、x′和y′的四維高斯分布的相關矩陣為[2]

通過把變量x、y、x′和y′轉換成極坐標得到包絡r、包絡的導數r′、相位φ和相位的導數（即頻率）φ′的聯合概率密度函數?？紤]到?1=0，并使用雅各比變換得到聯合概率密度分布函數[3]

式(11)對應的累積分布[3]是

式(12)說明接收信號瞬時頻率在-∞到∞之間都可以存在，并不局限于多普勒頻率的范圍內。在指定包絡數值r0條件下瞬時頻率的概率密度函數[3]為

式(13)說明信號的電平越高，方差越小。也就是說當接收信號處于深度衰落的時候，瞬時頻率的方差會很大，很可能在很大的頻率范圍內有分布。這是因為當信號處于深度衰落的時候，一方面較低的電平增加對噪聲的敏感性，另一方面較低的電平會增加隨機調頻的概率。另外需要特別說明的是，由式(10)對φ和φ′的積分就可以得到式(5)中的 pdfr,r'(r,r')。

由以上多普勒頻移對接收信號幅度和頻率的影響分析表明，微多普勒頻移不僅能使接收信號產生衰落，還能使得信號的瞬時頻率發生隨機調制。因此，不論采用相干解調還是非相干解調，對于PSK、FSK 信號的接收信噪比都會產生影響，造成無線測試鏈路的誤碼或中斷。

為了有效地對抗衰落，可以架設多根天 線[8-9]。假定各個天線接收的信號是不相關的，那么多根天線同一時刻都處于衰落的概率相對于單個接收天線將會大大降低，也就能夠大大降低衰落對接收信號的影響。采用空間多天線集中處理的思想如圖2所示，一方面某些接收天線可能收不到含有移動人員反射的信號，如接收機A；另一方面，即使各天線均接收含多普勒頻移的信號，但是多個天線同時位于衰落位置的可能性較小，如接收機B 和接收機C 雖然都接收到含多普勒頻移的信號，但是這2 臺接收機同時處于衰落狀態的概率很小。采用空間多天線集中處理方案的代價是需要增加1 臺信號集中處理器用于對各天線處不相關的接收信號的合成，其常見的合成算法包括互相關法[10]、自適應濾波法[11]、高階統計量法[12]和小波變換法[13]等。

圖2 無線測試空間分集接收示意圖Fig.2 Space diversity receiving scene of a wireless test

假設圖2中單臺接收機A、B 和C 接收信號電平小于門限電平的概率分別為P1、P2和P3，則這個有3 根接收天線的測試現場能夠進行正常測試的概率為P=1-P1P2P3?？梢钥闯?，僅增加2 根天線后，正常測試概率就可以遠大于單根天線正常測試的概率，提高測試可靠性的效果明顯，代價較小。

3 仿真結果

通過計算機仿真驗證以上機理分析的正確性。我們假定發射的信號是頻率為f=300 MHz 的正弦波，采樣速率fs=30 000 MHz，持續時間t=0.002 s。

3.1 包絡仿真

采用圖1所示的兩徑模型，直射路徑幅度為1，反射路徑幅度b=0.5，多普勒頻移fd=1000 Hz。仿真結果顯示的合成信號及其包絡如圖3所示。

由圖3可見，本來發射的恒定包絡為1 的正弦波，由于存在含有多普勒頻移的第2 條路徑，接收合成信號變成調幅信號。信號在包絡極小值處對應的幅度小于1，產生了衰落。相鄰兩個幅度極大值（或極小值）之間的距離就是多普勒頻移的倒數，即信道的相干時間tcoh=1/fd。

3.2 頻率仿真

采用MATLAB 函數rayleighchan(1/fs,fd)[14]，多普勒頻移fd=1000 Hz 的仿真結果如圖4所示?？梢钥吹?，對應于瞬時包絡較?。梢哉J為是衰落發生的時刻）的幾個區域（圖中圈出位置）的瞬時頻率比其他瞬時包絡較大時刻處的瞬時頻率有更大的可能性產生變化。

圖4 多徑信道下接收信號瞬時頻率和瞬時包絡間的關系Fig.4 Received signal of multipath channel∶instantaneous envelope and instantaneous frequency

特別需要說明的是：

1）瞬時頻率的縱軸顯示數值280.65 MHz 與仿真設定的300 MHz 有偏差，這是由于仿真參數設定1 個周期正弦波采集10 個采樣點，相鄰采樣點之間按照直線計算瞬時頻率會導致計算結果與理論值之間有偏差，通過提高采樣率可以降低這種偏差。但本文是對射頻的波形進行仿真，若采用更高的采樣率，則數據量過大，將造成仿真程序無法運行。

2）瞬時頻率的縱軸在不同頻率處顯示的數值相同，都是280.65 MHz，這是由于MATLAB 軟件的數值分辨率不夠高，對于數值較為接近的2 個數在圖形上會顯示為相同的數值。其實圖中瞬時頻率的最大值為max(f)=280.652 979 585 817 5 MHz，最小值為min(f)=280.646 119 628 650 0 MHz。綜合考慮第1）項所述的因素，取以上瞬時頻率最大值和最小值之間的差值Δf=6.859 957 167 506 218× 103Hz。該數值遠大于設定的1000 Hz 多普勒頻移。因此，該仿真證明了接收信號中含多普勒多徑信號時，其瞬時頻率可以遠遠大于多普勒頻移引起的頻率變化范圍。

4 結束語

通過本文分析可見，雖然微多普勒頻移的數值很?。ㄏ鄬τ谳d波頻率幾乎可以忽略），但是它對接收信號的幅度和相位都會產生影響，引起接收信號處于衰落深陷和隨機調頻，從而造成測試鏈路誤碼甚至中斷。因此，在制定航天器無線測試方案時，首先要加強測試現場管理，嚴格避免各類多普勒頻移源的產生；若無法避免微多普勒源時，可以考慮用空間分集技術克服微多普勒頻移的影響，通過對空間多個天線處（一般需超出相關距離）接收信號的集中處理來提高測試鏈路的信號質量。

（References）

[1]Proakis J G.數字通信[M].張力軍, 張宗橙, 鄭寶玉,等, 譯.4 版.北京∶電子工業出版社, 2006∶577-591

[2]Rice S O.Statistical properties of a sine wave plus random noise[J].Bell System Tech J, 1948, 27(1)∶109-157

[3]Jakes W C.Microwave mobile communication[M].New York∶WILEY-IEEE press, 1974∶10-26

[4]Bello P.Characterization of randomly time-variant linear channels[J].IEEE Transactions on Communication System, 1963, 11(4)∶360-393

[5]高紅衛, 謝良貴, 文樹梁, 等.微多普勒的一些工程問題研究[J].系統工程與電子技術, 2008, 30(11)∶2035-2039 Gao Hongwei, Xie Lianggui, Wen Shuliang, et al.Research on several engineering problems of micro- Doppler[J].System Engineering and Electronics, 2008, 30(11)∶2035-2039

[6]陳行勇, 劉永祥, 黎湘, 等.微多普勒分析和參數估計[J].紅外與毫米波學報, 2006, 25(5)∶360-363 Chen Hangyong, Li Yongxiang, Li Xiang, et al.Analysis of micro-Doppler and parameters estimation[J].Journal of Infrared and Millimeter Waves, 2006, 25(5)∶360-363

[7]陳廣鋒, 張林讓, 王純, 等.復合運動目標微多普勒特征的分析和提取[J].西安電子科技大學學報, 2011, 38(3)∶55-62 Chen Guangfeng, Zhang Linrang, Wang Chun, et al.Micro-Doppler feature analysis and extraction of the complex motion target[J].Journal of XIDIAN University, 2011, 38(3)∶55-62

[8]殷勤業, 張瑩, 丁樂, 等.協作分集∶一種新的空域分集技術[J].西安交通大學學報, 2005, 39(6)∶551-557 Yin Qinye, Zhang Ying, Ding Le, et al.Cooperation diversity∶a new spatial diversity technique[J].Journal of Xi’an Jiaotong University, 2005, 39(6)∶551-557

[9]李光球.Nakagami 衰落信道上組合SC/MRC 的性能分析[J].電波科學學報, 2007, 22(2)∶187-190 Li Guangqiu.Performance analysis of combined transmit SC/Receive MRC in Nakagami fadingchannels[J].Chinese Journal of Radio Science, 2007, 22(2)∶187-190

[10]Knapp C H, Cartar G C.The generalized correlation method for estimation of time delay[J].IEEE Transaction on Acoustics, Speech and Signal Processing, 1976, 24(4)∶320-327

[11]甄中建, 陳黎.自適應濾波在擴頻信號分集接收中的應用研究[J].半導體技術, 2007, 32(8)∶700-703 Zhen Zhongjian, Chen Li.Research on diversity reception in DS/SS communication with adaptive filter[J].Semiconductor Technology, 2007, 32(8)∶700-703

[12]宋濤, 許家棟, 梁雙港.一種適用于多發多收數字 信道的盲估計算法[J].西北工業大學學報, 2007, 25(4)∶566-570 Song Tao, Xu Jiadong, Liang Shuanggang.A better blind estimation algorithm for digital MIMO channels[J].Journal of Northwestern Polytechnical University, 2007, 25(4)∶566-570

[13]郭業才, 丁雪潔.正交小波變換和時間分集聯合盲均衡算法[J].南京信息工程大學學報, 2010, 2(1)∶1-5 Guo Yecai, Ding Xuejie.The orthogonal wavelet transform and time diversity combined blind equalization algorithm[J].Journal of Nanjing University of Information Science and Technology, 2010, 2(1)∶1-5

[14]張德豐.MATLAB 在電子信息工程中的應用[M].北京∶電子工業出版社, 2009∶1-25