潘 兵
(北京圣非凡電子系統(tǒng)技術(shù)開發(fā)有限公司 北京 102209)
CRM通信調(diào)制解調(diào)技術(shù)研究?
潘 兵
(北京圣非凡電子系統(tǒng)技術(shù)開發(fā)有限公司 北京 102209)
論文采用了匹配濾波器和分?jǐn)?shù)階傅里葉變換對CRM(chirp-rate modulation)通信信號進行解調(diào),推導(dǎo)了其理論解調(diào)誤碼率公式,將匹配濾波器解調(diào)和分?jǐn)?shù)階傅里葉變換解調(diào)的誤碼性能進行了對比分析;得到了影響解調(diào)性能的解析表達式,CRM通信信號抗多普勒頻移能力較強,但對碼同步誤差較為敏感。通過仿真分析,驗證了上述理論推導(dǎo)的正確性,分?jǐn)?shù)階傅里葉變換解調(diào)較匹配濾波解調(diào)性能差了約0.6dB~1.6dB,但其可避免載波恢復(fù),大大簡化了系統(tǒng)的同步設(shè)計。
CRM通信;匹配濾波器;分?jǐn)?shù)階傅里葉變換
目前用于chirp調(diào)制技術(shù)主要有二進制正交鍵控(Binary Orthogonal Keying,BOK)調(diào)制和直接調(diào)制(Direct Modulation,DM)兩種[1],在Chirp-BOK調(diào)制中,chirp信號被用于表示調(diào)制后的符號,采用調(diào)頻率攜帶信息,也即調(diào)頻率調(diào)制(Chirp Rate Modulation,CRM);而在Chirp-DM調(diào)制中,chirp信號僅用于擴展己調(diào)信號的頻譜,采用相位攜帶信息。由于CRM通信信號通過匹配濾波器的輸出具有明顯的壓縮脈沖[2],能量被集中在很短的時間內(nèi),因此可以用匹配濾波器來進行解調(diào);而chirp信號在分?jǐn)?shù)階傅里葉變換(Fractional Fourier transform,F(xiàn)rFT)域有很好的能量聚集性[3],因此還可用FrFT來檢測接收CRM通信信號能量聚集的峰值,以進行取樣判決,恢復(fù)出數(shù)據(jù)。本文對CRM通信信號的上述兩種解調(diào)方法進行研究,并對其進行性能仿真與分析。
2.1 chirp-rate調(diào)制原理
利用chirp信號進行數(shù)據(jù)傳輸?shù)幕舅枷胧前褦?shù)據(jù)信息調(diào)制到chirp信號的各參數(shù)上,主要參數(shù)包括線性調(diào)頻率、相位的正反極性、幅度的包絡(luò)等。chirp信號的表達式為[4]
s(t)=a(t)cos(2πf0t+μπt2) -T/2≤t≤T/2 (1)
式(1)中 ,a(t)為 chirp信 號 的 包 絡(luò) ,a(t)=1,||t≤T/2、f0為chirp信號的中心頻率、μ為調(diào)頻斜率、T為chirp脈沖寬度,μ>0時為up-chirp信號,μ<0時則為down-chirp信號。如圖1所示,chirp-rate調(diào)制通信時用up-chirp信號(正調(diào)頻斜率)表示信碼“1”,down-chirp信號(負(fù)調(diào)頻斜率)表示信碼“0”。圖2為chirp-rate調(diào)制原理。
圖1 up-chirp和down-chirp
圖2 chirp-rate調(diào)制原理
2.2 chirp-rate調(diào)制仿真
表1為本文chirp-rate調(diào)制參數(shù),圖3為chirp-rate信號的時域和頻域波形。
圖3 chirp-rate調(diào)制信號的時域和頻域波形
表1 chirp-rate調(diào)制參數(shù)
3.1 匹配濾波解調(diào)
3.1.1 匹配濾波解調(diào)原理
up-chirp信號和down-chirp信號具有尖銳的自相關(guān)特性,其相關(guān)函數(shù)為Sinc函數(shù)形式,up-chirp信號通過down-chirp匹配濾波器的輸出脈沖壓縮表達式為[5]
顯然最佳的抽樣判決時刻是在t=0時刻,包絡(luò)取最大幅值為 TB。從式(2)可看出,信號的波形具有Sinc函數(shù)的特點,如圖4(a)所示,主瓣時寬為2 B,其幅度為 TB;如圖4(b)所示,壓縮后峰值功率較未壓縮時有 TB倍增益,該特性非常利于解調(diào)恢復(fù)信號。
圖5比較了chirp信號通過匹配濾波器和非匹配濾波器的情況??梢钥闯觯?jīng)匹配濾波的輸出是明顯的壓縮脈沖,能量被集中在很短的時間內(nèi)。而經(jīng)非匹配濾波的輸出沒有被壓縮,能量則均勻分布在時域上,因此抽樣判決器可把對應(yīng)的信號辨別出來。
圖4 匹配濾波輸出的波形
圖5 chirp信號經(jīng)匹配濾波輸出波形
圖6所示為chirp-rate調(diào)制解調(diào)系統(tǒng)框圖,匹配濾波器通過比較輸出峰值的大小判定接收信號是up-chirp信號或down-chirp信號,進而恢復(fù)出原始信息。
圖6 chirp-rate調(diào)制解調(diào)系統(tǒng)框圖
3.1.2 匹配濾波解調(diào)誤碼率
如圖2所示,當(dāng)發(fā)送信號為s1(t)時,經(jīng)高斯白噪聲信道后s(t)=s1(t)+n(t)(n(t)均值為0,方差為n02的高斯白噪聲)。圖6上支路匹配濾波器輸出為
所以在高斯白噪聲信道下,chirp-rate調(diào)制在匹配濾波接收時的平均誤碼率計算公式為
從式(7)可看出,chirp-rate調(diào)制系統(tǒng)的性能與up-chirp和down-chirp的互相關(guān)系數(shù) ρ有關(guān)。當(dāng)ρ=0時,兩信號完全正交,系統(tǒng)誤碼率與chirp-rate相干解調(diào)誤碼率相同;當(dāng)ρ=0.5時,系統(tǒng)誤碼率與傳統(tǒng)BOK相干解調(diào)誤碼率相同[6]。在此本文給出互相關(guān)系數(shù) ρ的表達式[7]為 ρ=C(D)/D ,(C(x)為Fresnel積分,D=TB為時寬帶寬積)。表明互相關(guān)系數(shù)只與時寬帶寬積D有關(guān)。由表1參數(shù)知D=T*B=240,代入互相關(guān)系數(shù) ρ表達式可得ρ=0.032 接近0,即 s1(t)和 s2(t)準(zhǔn)正交。圖7(a)是互相關(guān)系數(shù) ρ與時寬帶寬積D的關(guān)系圖,圖7(b)是不同ρ值的誤碼率曲線圖。
圖7 (a) 互相關(guān)系數(shù) ρ與時寬帶寬積D關(guān)系
圖7 (b) 不同 ρ 值的誤碼率
可見,ρ與 D有關(guān),ρ越小,up-chirp和down-chirp信號越接近正交,系統(tǒng)性能也越好。從圖7(b)可以看出,當(dāng) ρ=0.032時,BER曲線基本與ρ=0曲線重合,基本可認(rèn)為是兩路正交信號。在BER=10-4時,所需信噪比相差僅0.15dB左右。
3.2 分?jǐn)?shù)階傅里葉變換解調(diào)
3.2.1 分?jǐn)?shù)階傅里葉變換解調(diào)原理
分?jǐn)?shù)階傅里葉變換(FrFT)的基函數(shù)是分?jǐn)?shù)階域上一組正交的chirp基,一個chirp信號在特定的分?jǐn)?shù)階傅里葉域中將表現(xiàn)為一個沖擊響應(yīng)[8]。如圖8所示,分?jǐn)?shù)階傅里葉變換在某個分?jǐn)?shù)階傅里葉域中對給定的chirp信號具有非常好的能量聚斂特性。
圖8 (b)down-chirp信號任一階次分?jǐn)?shù)域波形
實chirp信號的分?jǐn)?shù)階Fourier譜與相同參數(shù)復(fù)chirp信號的分?jǐn)?shù)階Fourier譜是一樣的,只是能量降低一半,多了一個對稱普。在實際通信系統(tǒng)中,復(fù)信號可通過I,Q兩路實信號來產(chǎn)生[9]。從分析信號的分?jǐn)?shù)階Fourier譜的角度出發(fā),可用復(fù)信號模型來替代實信號模型。設(shè)基帶chirp信號x(t)=Aexp(j2πf0t+jπkt2+jφ),則其分?jǐn)?shù)階 Fourier變換為Xα(μ)
對于給定的基帶信號形式,可以事先計算出“1”、“0”碼元相應(yīng)峰值的分?jǐn)?shù)階Fourier域采樣位置值um如圖9所示,從而在接收端直接在該點進行采樣判決。
從圖9中可以看出無論是復(fù)chirp信號還是實chirp信號在最佳階數(shù)的FrFT均產(chǎn)生尖峰脈沖,但尖峰脈沖的位置不同,因此可通過這一區(qū)別實現(xiàn)基于FrFT的解調(diào)。
根據(jù)實chirp信號和圖9所示的實chirp信號的分?jǐn)?shù)階傅里葉譜特點,可得到本文的FrFT解調(diào)方式,如圖10所示。
1)將接收后的信號混頻到基帶,再將實基帶信號按每個碼元周期做 p階FrFT變換處理;
2)對變換結(jié)果取模平方后,按照預(yù)先確定的采樣位置um1、um0進行采樣判決。
當(dāng) k=-cotα、f0=ucscα?xí)r,得到 x(t)的分?jǐn)?shù)階Fourier幅度譜峰值輸出為
圖9 碼元1和0對應(yīng)Chirp信號的FrFT
圖10 FrFT非相干解調(diào)
3.2.2 分?jǐn)?shù)階傅里葉變換解調(diào)誤碼率
在加性高斯白噪聲信道條件下,解調(diào)器輸入端的接收信號為 y(t)=x(t)+n(t)。x(t)為發(fā)送信號,n(t)是方差為的高斯白噪聲,y(t)服從高斯分布。由于FrFT為線性變換,因此對信號作 p階FrFT等價為將信號通過參數(shù)為 p的線性濾波器,高斯過程通過線性系統(tǒng)后仍服從高斯分布[10]。當(dāng)發(fā)送“1”碼時,有:
式(14)即為本文推導(dǎo)出的FrFT解調(diào)的理論誤碼率計算公式。
4.1 誤碼率性能分析
圖11所示為CRM通信的匹配濾波解調(diào)和FrFT解調(diào)的誤碼率仿真和理論曲線。
圖11 高斯白噪聲下CRM通信誤碼率性能
圖11結(jié)果表明,式(14)表示的理論曲線與實際仿真結(jié)果基本吻合,驗證了FrFT解調(diào)理論公式的正確性。另外,從圖中可以看出在BER=10-5~10-2時,F(xiàn)rFT解調(diào)比匹配濾波解調(diào)有約0.6~1.6dB的性能損失。由于FrFT解調(diào)屬于非相干解調(diào)方式,該結(jié)論是在預(yù)料之中的。(注:為使FrFT仿真解調(diào)性能接近于FrFT理論解調(diào)性能,需合理設(shè)計CRM通信參數(shù)使已調(diào)信號相位連續(xù),需設(shè)置chirp-rate信號中心頻率與信號時寬乘積為正整數(shù))。
本文所采用的FrFT算法是目前復(fù)雜度最小的DFrFT算法-Close-form算法[12]。
4.2 相位誤差對CRM通信解調(diào)性能的影響
4.2.1 相位誤差對匹配濾波解調(diào)性能的影響
若存在碼同步誤差τ,輸入chirp信號變?yōu)閤(t)=Aexp(j2πf0(t-τ)+jkπ(t-τ)2)。
經(jīng)過匹配濾波后,計算得出存在碼同步誤差的
普勒頻移效果,對移動通信中抗多普勒效應(yīng)具有良好的抑制效果。
圖12 衰減系數(shù)與頻移比的關(guān)系
對上述解調(diào)方式進行了簡單總結(jié),得到表2如下。
本文對CRM通信的兩種解調(diào)方式進行理論研究與仿真分析,F(xiàn)rFT解調(diào)算法比匹配濾波法性能稍差,但匹配濾波(相干解調(diào))需要嚴(yán)格估計載波的頻偏和相位,同步系統(tǒng)設(shè)計繁瑣,運算復(fù)雜度大。FrFT解調(diào)雖然比匹配濾波法相差0.6dB~1.6dB,但它對相位誤差影響小,且多普勒頻移與碼同步誤差的影響近似,F(xiàn)rFT解調(diào)屬于非相干解調(diào),因此可避免載波恢復(fù),這使得系統(tǒng)的同步設(shè)計大大簡化。由于CRM通信本身屬于一種擴頻通信體制(如本文擴頻增益約24dB),因此CRM通信可應(yīng)用于短波極低速(如本文信息速率為12.5bps)抗干擾隱蔽通信方面,具有很強的抗偵察抗截獲能力。
表2 FrFT解調(diào)與匹配濾波解調(diào)的比較
[1]Zhuang Weihua,Shen Xuemin,Bi Qi. Ultral wideband wireles communications[J].Wireless Communications and Mobile Computing,2003(3):663-685.
[2]Ji W L,Zheng G X,Bao M Q,etAl.Ultra-narrowband wireless communication technology based on QVMCK modulation[C]∕Proc. of China-Japan Joint Microwave Conference,2008:164-166.
[3]陶然,鄧兵,王越.分?jǐn)?shù)階傅里葉變換及其應(yīng)用[M].北京:清華大學(xué)出版社.2009:32-33.
[4]鄭森,鄭繼禹,仇洪冰,王玫.現(xiàn)代擴頻通信[M].西安:西安電子科技大學(xué)出版社.2013:127-128.
[5]賀鵬飛.基于Chirp-BOK調(diào)制的超寬帶無線通信系統(tǒng)研究[J].南京郵電大學(xué)學(xué)報,2006,26(2):21-25.
[6]Proakis J G.On the probability of error for multichannel reception of binary signals[J].IEEE Trans. Commun.Technol,1968,COM-16(1):68-71.
[7]潘兵.基于變時寬的CSS擴頻通信新技術(shù)[J].現(xiàn)代軍事通信,2016,24(1):15.
[8]Pei S C,Yeh M H,Tseng C C.Discrete Fractional Fourier Transform Based on Orthogonal Projections[J].IEEE Transactions on Signal Processing,1999,47(5):1335-1348.
[9]劉鋒.分?jǐn)?shù)階Fourier變換中量綱歸一化因子的選?。跩].系統(tǒng)工程與電子技術(shù),2011,33(2):237-239.
[10]劉建成,劉忠,王雪松.高斯白噪聲背景下的LFM信號的分?jǐn)?shù)階Fourier域信噪比分析[J].電子信息學(xué)報,2007,29(10):2337-2340.
[11]曹志剛,錢亞生.現(xiàn)代通信原理[M].北京:清華大學(xué)出版社,1992:276-278.
[12]Soo-Chang Pei,Jian-jiun Ding. Closed-form discrete fractional and affine Fourier transform[J].IEEE Trans.on SP,2000,48(5):1338-1353.
Research on CRM Communication Modulation and Demodulation Technology
PAN Bing
(Beijing Sheng Fei Fan Electronic System Technology Development Co.,Ltd,Beijing 102209)
This paper adopted the matched filter and fractional Fourier transform of the CRM(chirp rate modulation)communication signal demodulation,and its theory of demodulation error rate formula is deduced,the matched filter demodulation and the ber performance of fractional Fourier transform demodulation is analyzed.The analytical expressions that affect demodulation performance are obtained,and the communication signal of CRM is better than that of doppler shift,but it is sensitive to code synchronization error.Through the simulation analysis,validate the correctness of the theoretical derivation,the fractional Fourier transform demodulation is matched filtering demodulation performance difference is about 0.6dB~1.6dB,it can avoid carrier recovery,greatly simplifies the system design of synchronization.
CRM communication,matched filter,fractional Fourier transform
TN91
10.3969∕j.issn.1672-9730.2017.10.013
Class Number TN91
2017年4月8日,
2017年5月26日
潘兵,男,碩士研究生,工程師,研究方向:通信信號處理。