劉家昊, 章小梅
(海軍航空大學 青島校區(qū), 山東 青島 266041)
無線通信在實際應用中存在各種干擾,為避免干擾所造成的誤碼影響,使用者通常采用糾錯編碼擴展信源符號集,通過增加冗余度從而降低誤碼率,但是這會導致信道帶寬增加.而網(wǎng)格編碼(TCM)技術將調(diào)制和編碼過程相結合,可在不增加帶寬與功率的情況下,提升誤碼性能,提高通信質量.近年來,網(wǎng)格編碼技術發(fā)展迅速,與很多其他調(diào)制方法相結合形成了TCM-CPM技術、TCM-OFDM技術,并在常規(guī)通信領域外得到擴展與延伸,例如衛(wèi)星通信領域和多用戶檢測領域.
TCM網(wǎng)格編碼調(diào)制是由Ungerboeck在1982年提出的一種將信道編碼與調(diào)制統(tǒng)一考慮的技術[1-3].它將kbit信息分為k1和k2兩部分,其中k1bit被編碼成nbit,k2bit未參與編碼.通過分集映射將nbit對應的2n種組合對應于信號星座中2n個信號子集,編碼所得nbit用來選擇信號子集,未編碼的k2bit選擇各信號子集中2k2個信號點中的一個.TCM可以與mPSK或mQAM等多進制調(diào)制方式相結合,以TCM-8PSK調(diào)制方式為例,4狀態(tài)(3,2,2)網(wǎng)格編碼器結構圖如圖1所示.
圖1TCM-8PSK原理框圖
Fig.1BlockdiagramofTCM-8PSKprinciple
當今電磁環(huán)境日趨復雜,各類軍事訓練任務中通信干擾力度也不斷加大,且美、俄等國針對其所裝備的坦克、戰(zhàn)機、無人機等也進行了相關通信抗干擾研究[5-7].
本文結合某型對空超短波電臺的實際組成結構,并依據(jù)其技術體制和所采用的關鍵技術,構建了與實際裝備貼近的TCM-8PSK調(diào)制的跳頻系統(tǒng)仿真模型,如圖2所示.
圖2 TCM-8PSK調(diào)制方式下的跳頻系統(tǒng)仿真圖Fig.2 Simulation diagram of frequency hopping system in TCM-8PSK modulation mode
跳頻系統(tǒng)可分為6個模塊,分別為信源、編譯碼、調(diào)制解調(diào)、跳頻與解跳、信道和判決分析.其中信源模塊選用Simulink中的Continuous-Time VCO模塊來產(chǎn)生單一頻率的正弦波模擬發(fā)端的音頻信號,頻率可設置為300~3 400 Hz之間的任意值;編譯碼模塊使用PCM編譯碼方式;信道模塊采用加性高斯白噪聲AWGN信道模塊來模擬接收機的內(nèi)部噪聲.下面對調(diào)制解調(diào)、跳頻與解跳及判決分析3個模塊進行詳細說明.
1) 調(diào)制解調(diào)模塊.在本次的仿真模型中調(diào)制解調(diào)模塊采用的是TCM-8PSK模塊.具體結構如圖3所示.信源信號通過約束長度分別為1、3,碼率為2/3的卷積編碼器進入8移相鍵控調(diào)制器進行映射與調(diào)制,經(jīng)過信道后進行解調(diào)與映射.本模塊中使用Simulink模塊庫中Viterbi Decoder模塊進行維特比軟判決譯碼,由于該模塊軟判決譯碼時所需輸入數(shù)據(jù)為0~7之間的整數(shù), 故搭建Zixitong模塊,此模塊功能為轉換信號使得Viterbi Decoder譯碼模塊能正常識別.
圖3 TCM-8PSK模塊Fig.3 TCM-8PSK module
2) 跳頻與解跳模塊.由于m序列相關性較好、序列數(shù)目多等特點,因此,選擇其作為跳頻通信系統(tǒng)的載波信號,偽隨機碼由PN序列模塊產(chǎn)生器來產(chǎn)生[8],信源信號經(jīng)調(diào)制后與PN序列產(chǎn)生器所產(chǎn)生的偽隨機碼相乘完成跳頻功能.解跳模塊中直接使用跳頻模塊中所產(chǎn)生的跳頻信號,再與接收信號相乘完成解跳功能.
3) 判決分析模塊.在本仿真模型中,采用2種判決分析方法來比較跳頻系統(tǒng)的抗干擾能力,分別是誤碼率分析和波形失真度分析.誤碼率是衡量數(shù)據(jù)在規(guī)定時間內(nèi)數(shù)據(jù)傳輸精確性的指標,定義為傳輸中的錯誤碼元數(shù)與所傳輸?shù)目偞a元數(shù)的比,也有將誤碼率定義為用來衡量誤碼出現(xiàn)的頻率[9].誤碼率分析可利用Simulink中的Error Rate Calculation模塊搭建,仿真圖如圖4a所示,Display模塊用于統(tǒng)計數(shù)據(jù)的顯示.
圖4判決分析模塊
Fig.4Decisionanalysismodule
(a)—誤碼統(tǒng)計分析模塊; (b)—失真度分析模塊.
對于跳頻通信系統(tǒng),主要存在以下2種干擾方式,干擾同步系統(tǒng)以及干擾通信信息,這里主要考慮干擾通信信息,以下述3種典型干擾方式為例.
西部礦業(yè)股份有限公司是在2000年12月發(fā)起設立,2007年7月在上海證券交易所掛牌上市的大型礦業(yè)上市公司,總部設在青海省西寧市。公司注冊資本23.83億元,在全國擁有30多家控股公司和子(分)公司,主要業(yè)務是礦產(chǎn)資源綜合開發(fā)利用,主要有銅、鉛、鋅、鐵等有色金屬和黑色金屬的采選、冶煉、貿(mào)易等,分為礦山、冶煉、貿(mào)易、金融四大經(jīng)營板塊。截至2017年底,公司總資產(chǎn)322億元,凈資產(chǎn)115億元,營業(yè)收入278億元。
1) 多音干擾.實施多音干擾的前提是干擾方提前知道跳頻通信載波的一部分頻點,將干擾的功率集中到這些已知的頻率點上進行干擾.對于多音干擾信號的生成,利用多個Sine Wave(正弦波發(fā)生器)模塊來實現(xiàn),仿真模型如圖5a所示[10].不同頻率的信號通過相加器進行組合,就可以產(chǎn)生多音干擾信號.
2) 窄帶干擾.窄帶干擾信號指所占帶寬遠小于信號帶寬的干擾信號[11],仿真模型見圖5b.結合該電臺的實際組成結構,設置模擬濾波器的通帶邊緣頻率范圍為108~150 MHz.
3) 跟蹤干擾.跳頻信號的跟蹤干擾,是指干擾信號能“跟蹤”跳頻信號頻率跳變的一種干擾方式[12].跟蹤干擾要想達到理想的效果,必須滿足2個條件:一是干擾信號與跳頻載波頻率的變化規(guī)律必須一致,二是干擾信號要在載波跳變到下一頻點之前進入收信機.因此在建立跟蹤式干擾仿真模型時,采用的頻率合成器與跳頻載波發(fā)生器中的一致.另外,由于電子設備對信號的分析和處理需要花費一定的時間,所以在頻率合成器后面加入一個延時器,以模擬電子設備的處理時間.仿真模型見圖5c.
圖5 干擾模型Fig.5 Interference module(a)—多音干擾模型; (b)—窄帶干擾模型; (c)—跟蹤干擾模型.
為了分析3種干擾對TCM-8PSK系統(tǒng)的影響,在話音信號頻率(300~3 400 Hz)范圍內(nèi)選取10個頻率點,并在各干擾方式中選取一定強度的干擾變量,對未加干擾及加3種干擾情況下系統(tǒng)進行仿真,且輸出滿足失真度5%要求,仿真結果記錄如表1~表4所示.
通過變換頻率,可以看到隨頻率調(diào)整,TCM-8PSK系統(tǒng)誤碼率的大小與被傳輸信號的頻率大致成線性關系,誤碼率隨頻率的增大而增大.不同干擾情況下該系統(tǒng)的誤碼率變化情況不同,這里使用相關系數(shù)法作為比較判斷2個變量線性關系的方法,相關系數(shù)可以判定變量線性相關關系的強弱[13],如果是線性相關,當相關系數(shù)等于1時,說明變量之間的關系完全是線性關系,若相關系數(shù)等于0,表明變量之間無內(nèi)在關系[14].可用式(1)表示.
表1 TCM-8PSK系統(tǒng)誤碼率數(shù)據(jù)Table 1 Error rate data of TCM-8PSK system
表2 窄帶干擾(方差為1時)情況下系統(tǒng)誤碼率數(shù)據(jù)Table 2 System error rate data in the case of narrowband interference (when variance is 1)
表3 多音干擾(振幅為0.1時)情況下系統(tǒng)誤碼率數(shù)據(jù)Table 3 System error rate data in the case of multi-tone interference (when amplitude is 0.1)
表4 跟蹤干擾(平均值為3時)情況下系統(tǒng)誤碼率數(shù)據(jù)Table 4 System error rate data in the case of tracking interference (when mean value is 3)
表5 跟蹤干擾情況下系統(tǒng)誤碼率、幅值數(shù)據(jù)(t=0.000 01,平均值為3)Table 5 System error rate and magnitude data in the case of tracking interference (when the average value is 3, t=0.000 01)
表6 校正后跟蹤干擾情況下系統(tǒng)誤碼率、幅值數(shù)據(jù)(t=0.000 01,平均值為3)Table 6 System error rate and magnitude data after correction in the case of tracking interference (when the average value is 3, t=0.000 01)
(1)
其中,Cov(X,Y)為X與Y的協(xié)方差,Var[X]為X的方差,Var[Y]為Y的方差.
分別將3種干擾情況與TCM-8PSK系統(tǒng)的誤碼率數(shù)據(jù)進行對比,可知窄帶干擾、多音干擾、跟蹤式干擾的誤碼率數(shù)據(jù)與TCM-8PSK系統(tǒng)的誤碼率數(shù)據(jù)相關系數(shù)分別為0.985 188 693、0.990 130 029、0.962 039 296,可知該系統(tǒng)受跟蹤干擾影響最強,受窄帶干擾影響次之,受多音干擾影響最弱.
本文只考慮對系統(tǒng)影響最強的跟蹤干擾.設置跟蹤干擾的時延t為0.000 01,平均值為3.在話音信號頻率(300~3 400 Hz)范圍內(nèi)選取22個頻率點,且輸出滿足失真度5%要求,記錄仿真結果.選取其中部分數(shù)據(jù)如表5所示.
系統(tǒng)幅頻在兩端頻率處衰落過快,可在系統(tǒng)中增加一均衡模擬LC濾波器[15]對兩端進行補償,彌補幅頻衰落.均衡模擬補償電路如圖6所示(原系統(tǒng)電阻為R0,電容C容抗為XC,電感L的感抗為XL),通過調(diào)整R1、R2、C與L的值得出系統(tǒng)電路的輸出為
其中:VIN、VOUT分別代表校正前后信號波形振幅;F(z)代表所要設計的LC校正電路.
圖6 LC校正電路Fig.6 LC correction circuit
校正后系統(tǒng)仿真結果記錄見表6.
可得二者幅頻特性曲線圖,如圖7所示.
圖7 幅頻特性曲線圖Fig.7 Curve of amplitude-frequency characteristics
此時幅頻特性得到一定修正,由于補償區(qū)域僅為兩端頻率,分別計算補償區(qū)域原系統(tǒng)與校正后系統(tǒng)誤碼率相關系數(shù)即可,400~700 Hz頻段相關系數(shù)分別為0.790 580 741、0.949 772 153, 2 800~3 600 Hz頻段相關系數(shù)分別為0.962 013 687、0.984 357 982.可知隨著幅頻特性的校正,系統(tǒng)誤碼性能較原系統(tǒng)有了一定改善.
在給定仿真條件下,通過改變輸入端信源頻率,得到未加干擾及加3種干擾情況下系統(tǒng)誤碼率數(shù)據(jù),利用相關系數(shù)計算,分別得出窄帶干擾、多音干擾、跟蹤式干擾對TCM-8PSK系統(tǒng)誤碼率數(shù)據(jù)的相關系數(shù),從而可知該系統(tǒng)受3種干擾的影響程度大小,即表明系統(tǒng)抗多音干擾能力最強,抗窄帶干擾能力次之,抗跟蹤干擾能力最弱;并通過分析系統(tǒng)幅頻特性曲線,在跟蹤干擾情況下,比較原系統(tǒng)與校正后系統(tǒng)的誤碼性能,可知隨著幅頻特性得到補償,系統(tǒng)抗跟蹤干擾性能有一定提高.