魏瑞剛
(中國電子科技集團(tuán)公司第五十四研究所,河北 石家莊 050081)
直升機(jī)經(jīng)過近30年的發(fā)展,其飛行和作戰(zhàn)性能不斷提高,已經(jīng)逐漸成為戰(zhàn)場火力支援、武裝偵察以及人員輸送的主要力量。除此之外,由于我國地震、洪水和森林火災(zāi)等嚴(yán)重自然災(zāi)害頻發(fā),直升機(jī)因其具有快速、機(jī)動(dòng)靈活和不受地理?xiàng)l件限制等特點(diǎn),在搶險(xiǎn)救災(zāi)過程中作為空中救援和偵察平臺(tái),得到了大量應(yīng)用。
隨著直升機(jī)的應(yīng)用范圍越來越廣和任務(wù)復(fù)雜度不斷提高,機(jī)地間的不間斷通信、寬帶數(shù)據(jù)傳輸和多空域的地面協(xié)同指揮的需求也逐漸增長,原有的機(jī)載通信手段難以滿足這種要求。而衛(wèi)星通信具有廣域覆蓋、受地理環(huán)境約束小、支持寬帶組網(wǎng)通信和通信質(zhì)量穩(wěn)定等特點(diǎn),可以很好地解決這一問題,因此研究衛(wèi)星通信在直升機(jī)平臺(tái)上的應(yīng)用具有重要的意義。然而直升機(jī)衛(wèi)星通信信號(hào)的傳輸會(huì)受到直升機(jī)旋翼遮擋的影響,文獻(xiàn)[1]中給出了基于LMS估計(jì)的縫隙預(yù)測方案,該算法能夠完成發(fā)生縫隙時(shí)刻的預(yù)測,但是在工程實(shí)現(xiàn)中檢測的可靠性不高。文獻(xiàn)[2]中給出了基于重發(fā)策略的抗旋翼遮擋方法,該方法的缺點(diǎn)也是傳輸效率較低。文獻(xiàn)[3]中給出了基于信號(hào)檢測的縫隙通信思想,利用縫隙間隔完成信號(hào)的傳輸,但復(fù)雜度較高。本文通過分析直升機(jī)衛(wèi)星通信旋翼遮擋的特點(diǎn),設(shè)計(jì)了一種信道高效傳輸?shù)膸Y(jié)構(gòu),該幀結(jié)構(gòu)的信號(hào)捕獲概率高,虛警概率很小,能夠在不降低傳輸效率的前提下降低旋翼遮擋對信號(hào)捕獲的影響,而且能夠提高頻率估計(jì)的頻差范圍。
由于直升機(jī)物理結(jié)構(gòu)的特點(diǎn),導(dǎo)致直升機(jī)衛(wèi)通天線只能安裝在直升機(jī)旋翼的下方,其安裝示意圖如圖1所示。從圖1可見,直升機(jī)衛(wèi)通天線在收發(fā)衛(wèi)通信號(hào)時(shí),會(huì)受到直升機(jī)旋翼的遮擋,在直升機(jī)飛行的過程中,旋翼周期性地遮擋天線面,造成機(jī)載接收信號(hào)的衰減。
信號(hào)遮擋的周期由旋翼轉(zhuǎn)速和槳葉數(shù)決定,與直升機(jī)飛行航向和衛(wèi)星的仰角和位置等參數(shù)無關(guān)。
T=1/(V×N),
(1)
式中,T為遮擋周期(s);V為旋翼轉(zhuǎn)速(r/s);N為槳葉的數(shù)量。
圖1 直升機(jī)衛(wèi)通天線安裝示意
每個(gè)旋翼的遮擋時(shí)間可以近似依據(jù)式(2)決定,無遮擋時(shí)間近似依據(jù)式(3)決定,文獻(xiàn)[1]中給出在不同的航向夾角β的情況下,旋翼對天線的遮擋時(shí)間長短的計(jì)算公式,可以用γ表示旋翼對天線的遮擋時(shí)間。
(2)
(3)
(4)
λ=h×cotα,
(5)
式中,d為天線與旋翼中心的水平距離;h為天線與旋翼中心的垂直距離;α為天線仰角。
在直升機(jī)旋翼下,用實(shí)時(shí)頻譜儀采集的直升機(jī)旋翼遮擋的實(shí)測數(shù)據(jù)可以得到如圖2所示的信號(hào)衰減模型,信號(hào)衰減周期為旋翼遮擋周期,信號(hào)衰減幅度約為18~20 dB。
圖2 旋翼遮擋信號(hào)衰減
通過上述分析得知,在直升機(jī)旋翼遮擋下,機(jī)載接收機(jī)的接收信號(hào)功率會(huì)按照某種周期衰減,而接收信號(hào)功率的衰減勢必會(huì)影響解調(diào)流程中信號(hào)的同步,而同步信息僅用來完成信道鏈路的同步,不傳輸任何業(yè)務(wù)信息,因此同步信息的長度會(huì)影響信道傳輸?shù)男?,為了便于信道的同步,需要在傳輸?shù)挠行兄芷诓迦胪叫畔ⅲ送饪紤]信道的傳輸效率,要盡量減少同步信息插入的數(shù)量,平衡同步信息的數(shù)量和信道傳輸?shù)男适切盘?hào)幀結(jié)構(gòu)設(shè)計(jì)首先要考慮的問題,另外如果插入的信息能夠被多個(gè)功能利用,同樣會(huì)達(dá)到最優(yōu)設(shè)計(jì)的目的。鑒于以上目的,本文設(shè)計(jì)了一種幀結(jié)構(gòu),其結(jié)構(gòu)示意圖如圖3所示,通過傳輸該幀結(jié)構(gòu)的信息,機(jī)載解調(diào)器能夠在接收信號(hào)被旋翼周期遮擋的條件下正確完成信號(hào)的解調(diào)。
圖3 幀結(jié)構(gòu)示意
圖3中,每個(gè)物理幀由數(shù)據(jù)幀S和導(dǎo)頻幀P組成,其中每個(gè)數(shù)據(jù)幀S由30個(gè)符號(hào)比特構(gòu)成,而導(dǎo)頻幀P是由2種PN碼P1和P2組成,其中P1為周期是128的偽隨機(jī)碼序列,P2為周期為220的偽隨機(jī)碼的截短序列,信道數(shù)據(jù)的捕獲主要由導(dǎo)頻幀P1來完成,頻率同步和相位同步主要由P1和P2共同來完成。假定設(shè)計(jì)的系統(tǒng)編碼幀長度為16 200 Bit,則一幀數(shù)據(jù)中需要插入的導(dǎo)頻比特共有540 Bit,這樣幀效率可達(dá)97%,滿足傳輸系統(tǒng)的要求。
針對以上設(shè)計(jì)的幀結(jié)構(gòu),系統(tǒng)的解調(diào)同步可以按照如下方法進(jìn)行:首先對幀結(jié)構(gòu)中的數(shù)據(jù)信息和導(dǎo)頻信息利用傅里葉變換的方法完成初始的大頻偏估計(jì),對信號(hào)的大頻偏做一校正;然后再利用插入的導(dǎo)頻信息做二次頻率估計(jì),完成對殘余頻偏的頻偏校正,確保最后的頻差能夠滿足載波相位跟蹤的要求;最后利用相位內(nèi)插法完成載波的相位同步。
采用以上幀結(jié)構(gòu)的設(shè)計(jì)方法具有以下優(yōu)點(diǎn):
① 用來捕獲的導(dǎo)頻信息分散地插入數(shù)據(jù)比特中,在低速信號(hào)的傳輸中,降低了旋翼遮擋對信號(hào)捕獲的影響;
② 解調(diào)同步共用了用來捕獲的導(dǎo)頻信息,提高了信道傳輸效率;
③ 采用已知的導(dǎo)頻信息做頻率同步,提高了頻率估計(jì)的頻差范圍。
基于在旋翼遮擋的條件下,信號(hào)捕獲和頻率跟蹤在解調(diào)流程中的重要性,下面重點(diǎn)對信號(hào)的捕獲技術(shù)和頻率跟蹤技術(shù)做詳細(xì)分析。
信號(hào)的捕獲采用基于匹配濾波器的捕獲算法,算法原理框圖如圖4所示。
圖4 信號(hào)捕獲算法原理
接收到的IQ信號(hào)進(jìn)入數(shù)據(jù)抽取模塊,按照每隔30個(gè)符號(hào)抽取一次的規(guī)則進(jìn)行運(yùn)算,抽取后的IQ數(shù)據(jù)分別與本地PN碼產(chǎn)生的偽隨機(jī)序列做相關(guān)運(yùn)算,將相關(guān)運(yùn)算后的數(shù)據(jù)取模并且與設(shè)定的判決門限值做比較,進(jìn)而完成整個(gè)捕獲過程。
由最佳接收的理論可知,對于接收信號(hào)s(t),其匹配濾波器的沖擊響應(yīng)為:
h(t)=s(T-t),
(6)
式中,T為信號(hào)s(t)的持續(xù)時(shí)間,濾波器的輸出為:
(7)
假設(shè)數(shù)字匹配濾波器的輸入信號(hào)的采用速率為fs,濾波器的級數(shù)為N,則式(7)對應(yīng)的離散形式為:
(8)
h(kTs)=PN(NTs-kTs)。
(9)
在實(shí)現(xiàn)過程中,把式(9)中的h(kTs)看作濾波器的系數(shù),可以設(shè)置成PN碼序列{c0,c1,c2,…,cN-2,cN-1}。
在捕獲時(shí),接收到的信號(hào)與本地信號(hào)進(jìn)行相關(guān)運(yùn)算,可以將本地的PN碼序列作為固定的濾波器系數(shù),讓接收到的信號(hào)進(jìn)入濾波器中,每一時(shí)刻都會(huì)得到一個(gè)相關(guān)值,當(dāng)2個(gè)序列的相位一致時(shí),會(huì)出現(xiàn)一個(gè)相關(guān)峰,如果相關(guān)峰的大小超出給定的門限值,則信號(hào)捕獲成功,另外可以在門限判決模塊中加入一種控制策略來減小在旋翼遮擋的條件下信號(hào)捕獲的漏警概率。
下面通過仿真的手段來分析捕獲門限的選取及捕獲虛警的概率。
仿真條件:QPSK調(diào)制,8倍采樣,Es/N0為3 dB,歸一化頻偏為1/32,成型濾波器和匹配濾波器的滾降系數(shù)β=0.35的平方根升余弦濾波器,歸一化時(shí)鐘偏移為1/8,歸一化頻偏為1/32。
門限仿真圖如圖5所示。從圖5可以看到,本地碼與捕獲碼字對齊與不對齊的情況下,相關(guān)值分界分明,門限比較容易設(shè)置。
圖5 捕獲門限仿真
下面分析在不同的頻偏條件下信號(hào)的捕獲概率和虛警概率。假設(shè)如下仿真條件:符號(hào)Es/N0設(shè)為3 dB,判決門限設(shè)為280,旋翼遮擋周期為62.5 ms,遮擋比例為20%,歸一化頻偏以1/27的步進(jìn)在[1/27,1/23]內(nèi)變化,每個(gè)歸一化頻偏樣點(diǎn)蒙特卡洛仿真10 000次,得到的捕獲和虛警概率如表1所示。從表1可以看出,在整個(gè)歸一化頻偏區(qū)間內(nèi),虛警概率不大于0.1%,在工程實(shí)現(xiàn)上,這個(gè)概率在可允許的范圍內(nèi)。
表1 不同頻偏下捕獲概率和虛警概率
歸一化頻偏捕獲概率/%虛警概率/%1/271000.011/261000.031/251000.041/2499.990.031/2399.990.06
下面分析在不同的信噪比條件下,信號(hào)的捕獲概率和虛警概率。假設(shè)如下仿真條件:歸一化頻偏為1/25,判決門限設(shè)為280,旋翼遮擋周期為62.5 ms,遮擋比例為20%,Es/N0以1 dB的步進(jìn)在[0 dB,5 dB]范圍內(nèi)變化,每個(gè)Es/N0樣點(diǎn)蒙特卡洛仿真1 000次,得到相應(yīng)的捕獲和虛警概率如表2所示。從表2可以看出,Es/N0大于3 dB時(shí),捕獲概率可達(dá)100%,虛警概率小于0.1%,滿足實(shí)際的工程需求。
表2 不同符號(hào)信噪比下捕獲概率和虛警概率
Es/N0/dB捕獲概率/%虛警概率/%099.980.08199.980.0621000.0731000.0341000.0351000.02
本文在幀結(jié)構(gòu)的設(shè)計(jì)中,考慮了系統(tǒng)的頻偏估計(jì)問題,頻偏估計(jì)的原理框圖如圖6所示。
圖6 頻率估計(jì)原理
以QPSK調(diào)制方式為例,數(shù)據(jù)處理模塊首先對輸入數(shù)據(jù)作4次方運(yùn)算,并對運(yùn)算后的數(shù)據(jù)進(jìn)行存儲(chǔ)。當(dāng)累計(jì)存儲(chǔ)4096點(diǎn)后用高速時(shí)鐘發(fā)送給FFT模塊做FFT變換。
假設(shè)QPSK調(diào)制信號(hào)為:
s(t)=A×exp(jθM)exp[j(2πΔf+θ)],
(10)
式中,θM為QPSK的調(diào)制相位,θM=π/4×N(N=0,1,2,3);Δf為信號(hào)頻偏;A為調(diào)制信號(hào)幅度。對其進(jìn)行4次方運(yùn)算之后,
s(t)4=A4exp(j×4×θM)×
exp[j×(2π×4×(Δf+θ))]=
A4exp[j×(2π×4×(Δf+θ))]。
(11)
變換之后的信號(hào)為一頻率為4Δf的正弦信號(hào),理論上對其進(jìn)行FFT運(yùn)算后會(huì)在4Δf的頻點(diǎn)出現(xiàn)最大幅度譜線。因此,根據(jù)FFT運(yùn)算后最大值的位置可以計(jì)算出信號(hào)偏移。
假設(shè)如下仿真條件:調(diào)制方式為QPSK,歸一化頻偏為1/64,在門限信噪比Es/N0=3 dB,旋翼遮擋周期為62.5 ms,遮擋比例為20%,F(xiàn)FT長度為4 096條件下通過FFT計(jì)算出的頻偏仿真圖如圖7所示。從圖中可以看出,最大點(diǎn)的位置出現(xiàn)在第257個(gè)點(diǎn)處,則仿真出的頻偏為257/(4 096×4)=1/63.75,其估計(jì)的歸一化精度可達(dá)1/(4 096×4),經(jīng)過前文所述的二次估計(jì)后歸一化精度可達(dá)1/106,滿足工程的需求。
圖7 頻率估計(jì)仿真結(jié)果
在直升機(jī)衛(wèi)星通信系統(tǒng)的應(yīng)用背景下,在幀內(nèi)插入導(dǎo)頻信息解決了直升機(jī)旋翼遮擋對信號(hào)接收的影響,而且滿足所設(shè)計(jì)幀結(jié)構(gòu)的信號(hào)捕獲概率高,錯(cuò)捕概率較低,并且具有較高的頻率估計(jì)精度。由于該算法具有通用性,能夠適應(yīng)不同型號(hào)的直升機(jī),因此該設(shè)計(jì)方法具有極強(qiáng)的實(shí)際使用價(jià)值,給出的原理框圖可以指導(dǎo)工程實(shí)現(xiàn)。
[1] 匡麟玲,倪祖耀,凡明清.一種基于LMS的直升機(jī)旋翼縫隙時(shí)間預(yù)測方法及裝置:CN102521503A[P],2012-06-27.
[2] 董孝東.直升機(jī)衛(wèi)星通信中旋翼遮擋問題及應(yīng)對方法[J].現(xiàn)代電子技術(shù),2014,34(7):42-48.
[3] 姚國義,張德華.一種直升機(jī)旋翼遮擋縫隙檢測方法[J].無線電通信技術(shù),2014,40(3):36-38.
[4] 吳團(tuán)鋒,楊喜根,李際平.一種基于FFT的聯(lián)合幀同步和頻偏估計(jì)算法[J].電子與信息學(xué)報(bào),2006,28(5):840-843.
[5] 趙樹杰.信號(hào)檢測與估計(jì)理論[M].西安:西安電子科技大學(xué)出版社,1998:202.
[6] 潘亞漢,朱自強(qiáng),張志強(qiáng).現(xiàn)代衛(wèi)星通信調(diào)制解調(diào)技術(shù)的發(fā)展[J].微波與衛(wèi)星通信,1997(4):11-15.
[7] 張厥盛,鄭繼禹,萬心平.鎖相技術(shù)[M].西安:西安電子科技大學(xué)出版社,2000:86-89.
[8] 張潤峰.低信噪比條件下的載波恢復(fù)[D].長沙:國防科學(xué)技術(shù)大學(xué),2007:23-24.
[9] 黃振,陸建華,楊士中.衛(wèi)星通信中多普勒頻移的快速捕獲[J].電子學(xué)報(bào),2003,31(7):1052-1057.
[10] XIONG Fuqin.Modem Techniques in Satellite Communications[M].IEEE Journals & Magazines,1994,32(8):84-98.
[11] ERUP L,GARDNER F M,HARRIS R A.Interpolation in Digital Modems-Part II:Implementation and performance[J].IEEE Trans.Comet.,1993,41(6):998-1088.
[12] HOGENAUER E B.An Economical Class of Digital Filter for Decimation and Interpolation-I:Analysis and Design[J].IEEE Trans.on Acoustics,Speech and Signal Processing,1981,ASSP-29(2):155-162.
[13] MITOLA J.The Soft Radio Architecture[J].IEEE Communications Magazine,1995(5):26-39.
[14] 樊昌信,詹道庸,徐炳祥,等.通信原理(第4版)[M].北京:國防工業(yè)出版社,1995.
[15] 潘申富,張麗娜.一種低信噪比解調(diào)的實(shí)現(xiàn)方案及性能仿真[J].無線電通信技術(shù),2011,37(2):55-58.
[16] WOLEJSZA H,CACCIAMANI M.Identification of QPSK Modulation System Phase Ambiguity[J].IEEE Communications Magazine,1996(2):251-259.
[17] 潘申富.基于頻帶共享的衛(wèi)星通信關(guān)鍵技術(shù)研究[D].北京:北京大學(xué),2003:12-28.
[18] 張健.軟件無線電的基本理論構(gòu)架[R].成都:電子科技大學(xué),2000.
[19] MELO A G,SILVA H J A D.Improved Performance QPSK Demodulator for Digital Satellite Radio[C]∥IEEE Comm,2000:94-108.
[20] 樊昌信,曹麗娜,通信原理(第6版)[M].北京:國防工業(yè)出版社,2011:65-72.