王寧+楊國斌+劉桐辛+蘭婷+周晨+姜春華

摘 要： 為了評估電離層探測系統(tǒng)性能，研究了電離層信道的傳播特性，完成電離層探測系統(tǒng)的仿真建模及后續(xù)的數(shù)據(jù)處理。系統(tǒng)采用偽隨機二相編碼調(diào)制方式，選取Watterson信道模型仿真電離層信道，以超外差接收機和正交解調(diào)方式完成回波信號的采集，最后利用武漢大學(xué)電離層實驗室的雷達數(shù)據(jù)處理軟件對得到的數(shù)據(jù)進行處理，得到高質(zhì)量的雙時響應(yīng)圖和散射函數(shù)圖。仿真結(jié)果驗證電波在電離層中傳播時受到的各項影響，確定系統(tǒng)性能最佳時各模塊的參數(shù)，表明該模型能夠用以指導(dǎo)實際電離層探測系統(tǒng)的設(shè)計和研發(fā)。

關(guān)鍵詞： 電離層探測系統(tǒng)； 仿真建模； 偽隨機二相編碼； Watterson信道； 超外差接收機； 正交解調(diào)

中圖分類號： TN011+.2?34； TN95 文獻標(biāo)識碼： A 文章編號： 1004?373X（2018）06?0001?04

Abstract： To evaluate the performance of ionospheric sounding system and study the propagation feature of ionospheric channels， simulation modeling and subsequent data processing for ionospheric sounding system were achieved. The pseudo?random binary code modulation mode is adopted in the system to select Watterson channel model for simulation of ionospheric channels. Superheterodyne receiver and quadrature demodulation mode are used to complete return signal acquisition. The radar data processing software in the ionosphere laboratory of Wuhan University is used to process the obtained data and obtain the high?quality double?delay response diagram and scattering function diagram. The simulation results verified the influences of radio wave on radio waves propagating in ionosphere， and determined parameters of each module when the system reaches optimum performance， which demonstrates that the model can be used to direct the actual design and development of ionospheric sounding system.

Keywords： ionosphere sounding system； simulation modeling； pseudo?random binary code； Watterson channel； superheterodyne receiver； quadrature demodulation

1925年，美國的G.Breit和M.A.Tuve首次證明了電離層存在，開啟了利用無線電探測技術(shù)研究電離層的時代[1]，此后，各種電離層探測技術(shù)和設(shè)備都得到廣泛發(fā)展。與此同時，隨著計算機技術(shù)的進步，電離層探測系統(tǒng)的仿真研究也日趨重要。通過建模開發(fā)出針對實際系統(tǒng)的軟件仿真模型，能夠?qū)ο到y(tǒng)的硬件設(shè)計和研制提供理論支持，評估系統(tǒng)的各項性能，同時可以從信號級角度上對電波在電離層中傳播的工程化應(yīng)用問題進行仿真和復(fù)現(xiàn)，具有重要的科學(xué)應(yīng)用價值。

目前國內(nèi)系統(tǒng)仿真的主要平臺有Simulink，ADS，SPW，VC++等，兼顧編程難度、通用性、可擴展性、以及人際交互性等因素，Simulink的應(yīng)用較為廣泛，是一個提供動態(tài)系統(tǒng)建模、仿真和綜合分析的集成環(huán)境，具有模塊庫豐富、適應(yīng)面廣、仿真精細等優(yōu)點[2]。本文基于Simulink仿真平臺搭建完整的電離層探測系統(tǒng)模型，并展示了相應(yīng)的仿真結(jié)果及數(shù)據(jù)處理結(jié)果，驗證了仿真模型的正確性。

1 系統(tǒng)仿真建模

1.1 仿真模型

圖1為仿真模型整體框圖。電離層探測系統(tǒng)主要包括發(fā)射模塊、短波信道模塊和接收機模塊，其中接收機模塊又分為模擬部分和數(shù)字部分。發(fā)射模塊完成發(fā)射信號的產(chǎn)生、調(diào)制以及信號功率放大過程，短波信道模塊模擬電離層信道的傳播效應(yīng)，包括多徑效應(yīng)、多普勒頻移效應(yīng)、衰落、噪聲以及各種損耗[3?4]。接收機模塊完成對信號的解調(diào)、接收和數(shù)據(jù)采集，得到的I，Q數(shù)據(jù)流用于后續(xù)的數(shù)據(jù)處理。

1.2 發(fā)射模塊建模

發(fā)射模塊建模要考慮到雷達的最大作用距離和距離分辨力兩個特性，脈沖雷達方程和雷達的距離分辨力如下：

[ Rmax=PtτGtσλ2F2tF2r4π3D0kTsBnL14] （1）

[ΔR=cτ2=c2B] （2）

式中：[Pt]是發(fā)射功率；[τ]為脈沖寬度；[Gt]為發(fā)射天線增益；[Gr]為接收天線增益；[σ]為雷達散射截面；[λ]為信號波長；[Ft]是發(fā)射天線到目標(biāo)的方向圖傳播因子；[Fr]是目標(biāo)到接收天線的方向圖傳播因子；[k]為波爾茲曼常數(shù)；[Ts]為接收系統(tǒng)噪聲溫度；[D0]為可見度系數(shù)；[Bn]為帶寬較正系數(shù)；[L]為損耗因子；[c]為光速；[B]為發(fā)射信號的帶寬[5]；[Ptτ]是發(fā)射脈沖能量。由式（1）可以看出脈寬越大，雷達的發(fā)射功率越大，最大作用距離也越大。但這同時會造成雷達探測距離分辨力的下降。為了解決這一矛盾，系統(tǒng)采用脈沖壓縮技術(shù)，從而使雷達能夠同時獲得較大的作用距離和較高的距離分辨力[6]。脈沖壓縮有相位編碼脈沖和線性調(diào)頻脈沖兩種。相對于線性調(diào)頻脈沖而言，在電離層狀態(tài)變化較平穩(wěn)，多普勒頻移較小的情況下，二相編碼壓縮脈沖具有較高的距離分辨力，且相比于多相編碼，其復(fù)雜性大大降低。因此實際仿真中采用二相編碼脈沖。圖2是發(fā)射模塊的Simulink建模。endprint

首先通過脈沖信號與正弦調(diào)制信號相乘來實現(xiàn)幅度調(diào)制，利用脈沖信號的周期和脈寬設(shè)置發(fā)射信號的碼寬和脈寬。然后采用偽隨機序列進行調(diào)制信號的相位選擇，其中偽隨機m序列的產(chǎn)生可以用Simulink庫的m序列產(chǎn)生器模塊，也可以用延時單元構(gòu)建反饋移位寄存器來實現(xiàn)。在Matlab函數(shù)模塊中編寫的BPSK_a函數(shù)用于選擇正弦載波信號的初始相位。由于偽隨機序列具有良好的自相關(guān)特性，能夠降低系統(tǒng)發(fā)射功率，保證最大主旁瓣峰值比，得到的波形近似一個噪聲調(diào)制的信號，使得雷達能夠兼顧速度和距離上的測試精度和分辨力，增強其抗干擾能力，減小多徑造成的影響[7]。

1.3 Watterson信道建模

電離層信道的傳播特征主要有多徑效應(yīng)、衰落效應(yīng)、多普勒效應(yīng)、噪聲和損耗等。本文采用一種窄帶信道模型（Watterson信道模型）來模擬電離層信道。Watterson信道的數(shù)學(xué)模型的原理是：在窄帶和有限時間的條件下，將時變的非平穩(wěn)電離層信道看作是平穩(wěn)的，從而建立一個靜態(tài)的模型來模擬電離層信道[8]。

圖3是利用Simulink建模得到的Watterson信道模型。首先將信號進行希爾伯特變換，得到解析信號；然后進行數(shù)據(jù)格式轉(zhuǎn)換，將基于采樣的解析信號轉(zhuǎn)化為基于幀格式的信號，轉(zhuǎn)換后的信號進入信道特性模塊。Watterson信道模型是利用抽頭時延線來模擬信道的多徑效應(yīng)，每條路徑可單獨設(shè)置響應(yīng)的時延和增益，增益函數(shù)的幅度服從瑞利分布，相位均勻分布。此特性由模型中的瑞利衰落模塊實現(xiàn)，其中各路徑參數(shù)輸入都以向量的形式輸入，向量中元素的個數(shù)就是信道路徑數(shù)；高斯白噪聲模塊的作用是設(shè)置信號的信噪比，模擬信道的噪聲，信噪比越小，則信道模型產(chǎn)生的噪聲越大，信號受到噪聲的影響也越大。模型中的三個損耗分別是自由空間傳播損耗，地面反射損耗和額外系統(tǒng)損耗（極化損耗，多徑干涉損耗等）。其中自由空間傳播損耗由收發(fā)點的位置確定，地面反射損耗與工作頻率，射線仰角及地面參數(shù)有關(guān)，額外系統(tǒng)損耗是難以確定的參量，通常取其經(jīng)驗值，其值一般與探測時間相關(guān)。

1.4 接收模塊建模

仿真系統(tǒng)的接收機部分采用超外差結(jié)構(gòu)。首先由接收機模擬前端采用二次混頻方案得到中頻信號；然后由數(shù)字中頻模塊經(jīng)過數(shù)字混頻將頻率搬到零中頻；零中頻信號經(jīng)過抽取濾波，從而轉(zhuǎn)化為目的基帶信號。計算機后續(xù)的軟件處理工作都是圍繞著基帶信號展開的。圖4和圖5分別是接收機模擬端建模和數(shù)字中頻模塊建模。

由圖4可知，系統(tǒng)工作在3～30 MHz的短波頻段，接收機模擬端首先通過預(yù)選濾波和亞倍頻濾波器組進行濾波，消除不必要的頻率分量，經(jīng)過直通或高放后輸入第一混頻器單元得到一中頻信號，頻率為41.4 MHz。匹配發(fā)射信號的帶寬，一中頻濾波器通帶設(shè)置為41.36 ～41.44 MHz，該濾波器能夠抑制混頻器產(chǎn)生的不需要的變頻頻率分量，同時也可以抑制鏡像頻率。一中頻經(jīng)過濾波放大后進入二混頻器，二混頻器得到的信號經(jīng)二中頻濾波放大輸出到數(shù)字中頻模塊。二中頻頻率為1.4 MHz。

接收機模擬端輸出的信號經(jīng)過A/D轉(zhuǎn)換以后進入數(shù)字中頻模塊[9]，A/D轉(zhuǎn)換器采樣率設(shè)為10 MHz。圖5是數(shù)字中頻模塊的Simulink模型。為了保證對信號進行較好的接收，下變頻的抽取因子設(shè)置為256，第一級抽取采用積分梳狀濾波器（CIC）實現(xiàn)，CIC的抽取因子設(shè)置為16，第二級抽取由Simulink庫的下采樣模塊實現(xiàn)，抽取因子也是16。數(shù)字壓控振蕩器（NCO）相當(dāng)于一個理想的正余弦函數(shù)發(fā)生器，其作用是實現(xiàn)相干檢波[10]，防止單路檢波產(chǎn)生的過零點檢測盲相。一路信號和余弦信號進行檢波，得到同相支路I，另一路信號和正弦信號進行檢波，得到正交支路Q，通過這兩路信號，能夠同時提取幅度和相位信息。

2 仿真結(jié)果及信號處理

2.1 仿真波形分析

為了測試上述仿真系統(tǒng)的性能，對該仿真系統(tǒng)進行測試，選擇的探測模式及參數(shù)設(shè)置如下：探測頻率為6 MHz；脈沖寬度[τ=25.6 μs]；二相調(diào)制序列為m序列；階數(shù)為6；碼寬為10；脈寬為2；傳播路徑數(shù)為1；路徑時延為1.28×10-3 s；多普勒頻移為2 Hz；信噪比為10 dB；探測次數(shù)為128。圖6是仿真輸出波形。

圖6中：①是發(fā)射波形，發(fā)射信號采用二相編碼調(diào)制方式，占空比為0.2，每個碼元的相位由m序列調(diào)制，這種調(diào)制方式適合單站探測；②是解調(diào)以后的中頻信號，頻率為1.4 MHz，從波形看出，信號經(jīng)過信道傳輸以后信噪比有所下降，可見信號疊加了噪聲的干擾，并且由于窄帶濾波的影響，波形有一定程度的畸變；③是經(jīng)過A/D采樣和正交鑒相后得到的基帶I，Q信號，由波形可以看出，該數(shù)據(jù)流能夠采集到m序列的碼元信息，仿真結(jié)果正確。

2.2 數(shù)據(jù)處理

將仿真得到的I，Q數(shù)據(jù)流進行相關(guān)處理運算得到雙時響應(yīng)圖和散射函數(shù)圖，如圖7所示。

系統(tǒng)距離分辨力的計算公式為[τ?c]。其中[τ]為發(fā)射信號的脈沖寬度，c為光速，由此計算出仿真系統(tǒng)的距離分辨力為7.68 km。由圖7a）可以看出，探測次數(shù)為128次，在群距離為415 km左右處出現(xiàn)回波信號，同時由圖7b）可以看出，信號的多普勒頻移為2 Hz左右。

雙時響應(yīng)圖中的群距離由公式[τs?c]決定，[τs]為信號時延。信號時延來自兩部分：一部分是系統(tǒng)時延，約1×10-4；另一部分是信道時延，由信道參數(shù)知，信道時延為1.28×10-3，信號時延為兩部分時延之和。經(jīng)過計算可知回波應(yīng)出現(xiàn)在415 km左右處，這與圖7a）中的結(jié)果是一致的，通過更改信道時延參數(shù)可以模擬電離層各層回波信號。Watterson信道中設(shè)置的多普勒頻移為2 Hz，這與圖7b）中的結(jié)果基本符合，驗證了仿真系統(tǒng)的正確性。endprint

單次探測的距離響應(yīng)圖和415 km處的多普勒響應(yīng)圖如圖8所示。

從圖8a）可以看出，一次探測時回波能量在群路徑415 km左右處出現(xiàn)較大值，證明回波信號應(yīng)出現(xiàn)在415 km左右處。由圖8b）可以看出，在415 km距離處進行多次探測得到的回波多普勒頻移約為2 Hz，仿真結(jié)果正確。

3 結(jié) 論

基于Matlab/Simulink設(shè)計的電離層探測系統(tǒng)仿真模型能夠準(zhǔn)確模擬實際電離層探測系統(tǒng)的工作過程和電離層信道的傳播特點，仿真結(jié)果正確。在仿真系統(tǒng)搭建和調(diào)試過程中通過不斷更改各硬件仿真模塊參數(shù)，使系統(tǒng)各項性能達到最優(yōu)，對實際探測系統(tǒng)研制過程中遇到的調(diào)制方式選擇，濾波器帶寬設(shè)置，以及解調(diào)方式等都提供了借鑒，且仿真成本小、模型清晰、調(diào)試方便，因此電離層探測系統(tǒng)的仿真研究能夠為實際系統(tǒng)的研發(fā)提供堅實的理論基礎(chǔ)，具有重要的指導(dǎo)意義。

參考文獻

[1] 白亞軍.電離層對高頻電波傳播的影響研究[D].西安：西安電子科技大學(xué)，2014.

BAI Yajun. Study on the effects about the propagation of high?frequency radio waves in ionosphere [D]. Xian： Xidian University， 2014.

[2] 姜成昊，楊進華，張麗娟，等.新型多普勒成像激光雷達原理設(shè)計與仿真[J].紅外與激光工程，2014，43（2）：411?416.

JIANG Chenghao， YANG Jinhua， ZHANG Lijuan， et al. Schematic design and simulation of new Doppler imaging laser radar [J]. Infrared and laser engineering， 2014， 43（2）： 411?416.

[3] 張潔寒.電離層多徑信道均衡方法研究[D].西安：西安電子科技大學(xué)，2014.

ZHANG Jiehan. A study of multipath channel equalization method of ionized stratum [D]. Xian： Xidian University， 2014.

[4] 李雪，柳文，李吉寧，等.電離層相位污染和多模傳播檢測與識別研究[J].電波科學(xué)學(xué)報，2014，29（3）：485?491.

LI Xue， LIU Wen， LI Jining， et al. Detection and recognition of ionospheric phase contamination and multimode propagation [J]. Chinese journal of radio science， 2014， 29（3）： 485?491.

[5] PANG C， ZHANG N， HOU H. Radar equation of high?speed weak target during long time integration [J]. Fire control radar technology， 2015（1）： 7?9.

[6] 余貴水，張潤哲，左雷.脈壓雷達信號處理仿真分析研究[J].海軍工程大學(xué)學(xué)報，2015（3）：95?98.

YU Guishui， ZHANG Runzhe， ZUO Lei. Simulation and analysis of pulse compression radar signal processing [J]. Journal of Naval University of Engineering， 2015（3）： 95?98.

[7] VENKATA K R M. On the periodicity of non?maximal length linear feedback shift register sequences [J]. International journal of engineering & technical research， 2016， 5（5）： 521?525.

[8] 粟根花，夏小霞.寬帶短波信道建模綜述[J].數(shù)字通信，2014，41（3）：48?52.

SU Genhua， XIA Xiaoxia. Summary of broadband shortwave channel modeling [J]. Digital communication， 2014， 41（3）： 48?52.

[9] 陳彥龍，張援農(nóng).電離層多普勒接收機的設(shè)計與實驗[J].空間科學(xué)學(xué)報，2015，35（5）：574?580.

CHEN Yanlong， ZHANG Yuannong. Design and testing of an ionospheric Doppler receiver [J]. Chinese journal of space science， 2015， 35（5）： 574?580.

[10] 林琳.基于FPGA的雷達信號處理板設(shè)計與實現(xiàn)[J].現(xiàn)代電子技術(shù)，2014，37（11）：51?56.

LIN Lin. Design and implementation of radar signal processing board based on FPGA [J]. Modern electronics technique， 2014， 37（11）： 51?56.endprint