廖 明,姚 軍,張 偉,王世練
(1.中國工程物理研究院研究生部,四川綿陽621900;2.中國工程物理研究院電子工程研究所,四川綿陽621900; 3.國防科學技術大學電子科學與工程學院,湖南長沙410073)
多速率突發PSK信號全數字無前導字解調器
廖 明1,2,姚 軍2,張 偉2,王世練3
(1.中國工程物理研究院研究生部,四川綿陽621900;2.中國工程物理研究院電子工程研究所,四川綿陽621900; 3.國防科學技術大學電子科學與工程學院,湖南長沙410073)
突發模式通信應用日益廣泛,基于前導字的解調算法降低了系統的數據傳輸效率。提出了一種適于多速率突發信號的全數字無前導字解調方案,通過FIFO操作和采樣率轉換適應多種速率的要求,參數估計、定時同步和載波恢復均采用前向算法,可以在不損失任何數據符號的條件下完成解調,具有較高的數據傳輸效率。算法易于硬件實現,給出了仿真和硬件測試的結果,結果顯示在中高信噪比條件下算法具有較好的性能,驗證了算法的有效性和可靠性。
多速率 突發模式 無前導字
突發通信模式具有良好的短時性、隱蔽性等特點,在衛星通信、短波通信等系統中得到廣泛的應用。突發傳輸中,前后兩個突發數據包相對獨立,其載波頻偏、相偏和定時誤差等參數都不同,故每一個數據包到來時,解調器都要重新估計所需要的參數。文獻[1-3]研究了基于前導字輔助的突發解調算法。其中,文獻[1]通過解耦合處理,可以同時估計出定時誤差、載波頻偏和相位,從而有效減小了前導字的長度;文獻[2]提出一種利用已知前導數據的相關函數來估計載波頻偏的新算法,并對比了幾種不同估計算法的性能。文獻[3]基于Gardner算法提出一種新用法,取定時誤差最小的點作為最佳采樣點。基于前導字的算法性能較好,但是前導字占用了有限的傳輸帶寬,降低了數據傳輸效率,且不適用于盲解調。文獻[4-7]研究了無前導字的突發解調算法。其中,文獻[4]采用相位差分算法估計載波頻偏,所設計的接收機誤碼性能與理論值相比僅下降0.2 dB左右,但是容許的頻偏范圍只有[-5%·Rs,5%Rs];文獻[5]提出一種中頻差分解調方案,不需要突發載波恢復,適合在較低頻率下工作;文獻[6]通過最大化相位差累加矢量的模得到定時誤差,從16個采樣點中選出最佳采樣點,再通過相位平均旋轉度得到頻偏,估計范圍為Rs;文獻[7]給
出了一種適用于大頻偏的解調器,先通過兩級FFT運算,將頻偏縮小到較小的范圍,再通過一個前向的頻率跟蹤環去除殘余的頻偏,其估計范圍可以達到符號率的幾倍。
上述文獻中,突發信號的符號速率都是固定不變并且已知的,而在實際應用的突發通信系統中,每個突發包的符號速率有可能在一定的取值范圍內隨機變化,每個突發包的符號個數不同,相應解調器的結構也要與之適應。針對這種系統,文中提出一種支持多速率突發PSK信號的全數字無前導字解調方案,適用于較大的頻偏范圍,且可以在不損失任何數據符號的情況下完成解調,算法易于工程實現,適合非合作通信方的接收。下面介紹了所設計突發解調器整體結構,論述了各模塊的主要算法,最后給出算法仿真和硬件實現的結果。
1.1 整體結構
多速率突發解調器整體結構如圖1所示,虛線框內是解調器的主體部分。圖1中,中頻信號經ADC采樣轉換為數字信號,采樣時鐘為與符號速率獨立的固定值。突發信號在解調之前先送入一個預處理模塊,預處理模塊給出信號帶寬非常粗略的估計值和載波頻率粗估計值,同時在每一個突發包開始時刻輸出一個標志信號。對PSK信號,由可以得到符號速率的粗估計值。突發解調器接收到數據包后,先根據進行數字下變頻(DDC)將信號搬移到零頻附近,再根據進行多級半帶抽取濾波(HBF)將采樣率限制在符號速率的4~8倍。乒乓FIFO緩存A和B以突發數據包為單位將前后兩個數據包分別存儲,以適應不同的符號速率。低通濾波器(LPF)去除諧波分量后,參數估計模塊計算出當前數據包的載波頻偏和符號速率的精估計值。參數計算的同時,當前突發包進入FIFO緩存C延時,使數據和估計模塊的和同時輸出對齊。頻偏校正模塊進行二次變頻,去除。數據通過采樣率轉換模塊之后采樣率轉換為符號速率的4倍,匹配濾波(MF)后的數據通過定時恢復模塊,恢復出每符號一個最佳采樣點。相偏恢復模塊計算出相偏估計值,然后將信號校正。
不同突發數據包的符號速率不同,后續處理模塊所需的時間也不同。為了適應多速率突發信號,文中設計了一種乒乓FIFO緩存結構,該結構在突發起始時刻標志信號的控制下,將相鄰的兩個突發數據包分離開,以滿足后續模塊不同的延時需要,從而保證不會損失數據符號。

圖1 多速率突發信號無前導字解調器整體結構Fig.1 Structure of burst demodulator without preamble for multiple symbol rates
1.2 符號速率/頻偏估計
PSK信號具有循環平穩特性,由循環平穩有關理論可知[8],信號通過非線性變換之后的頻譜中會出現代表各階循環平穩特性的峰值譜線,譜線的位置一般對應了信號載波頻率和符號速率的線性組合,因此可以通過提取這些譜線來完成信號的參數估計。
文中提出一種基于M次方譜的符號速率/頻偏估計算法,以BPSK為例進行推導,算法可以推廣到MPSK。經過正交下變頻和低通濾波(LPF)后的BPSK信號表示為:

式中,gT(t)為成形濾波器,Ts=1/Rs為符號周期,Δf為載波頻偏,φn為調制相位。對于BPSK而言,φn= {0,π}。對信號進行平方運算,簡化后的結果為:

因為成形濾波函數gT(t)=0,當t?(0,Ts),所以有:
gT(t-nTs)gT(t-mTs)=0,m≠n(3)
結合式(3)可知,僅當m=n時式(2)的乘積項為非零值,故式(2)可化簡為:

由式(4)可知,y2(t)含有(t)以Ts為周期進行延拓與二倍載波頻偏分量的乘積,對其進行FFT運算后,頻譜將在f=2Δf±n(1/Ts)=2Δf±nRs處出現峰值譜線。一般情況下,在n=0,1時,即f=2Δf-Rs,2Δf,2Δf+Rs處譜線最明顯,如圖2所示。

圖2 BPSK信號的平方譜Fig.2 Quadratic spectrum of BPSK signal
圖2(a)是FFT之后的頻譜圖,圖2(b)將零頻搬移到中間,通過峰值搜索求得譜線的位置,從而可以計算出Δf和Rs的估計值。算法的估計范圍較大,允許的符號速率偏差為[-0.2Rs,0.2Rs],載波頻率偏差為[-0.5Rs,0.5Rs]。估計精度由FFT的頻率分辨率決定,頻率分辨率為(1/M0)Rs,M0為參與FFT的符號個數。
1.3 采樣率轉換

式中,mk,μk控制插值的位置,cl(k)為插值濾波器的系數。文獻[10]提出可以采用高效的Farrow結構來實現插值濾波,該結構易于硬件實現,框圖如圖3所示。

圖3 Farrow插值結構Fig.3 Structure of farrow interpolation
1.4 定時恢復
傳統定時同步算法一般是由誤差信號反饋控制NCO,不斷調整ADC的采樣時鐘使其達到最佳采樣位置。對于突發模式下的定時同步,收發兩端工作在穩定獨立的時鐘頻率上,解調器通過前向算法提取信號的定時誤差,然后通過插值算法直接計算出最佳采樣點來,定時恢復的結構如圖4所示。

圖4 突發PSK信號前向定時恢復Fig.4 Feedforward timing recovery for burst PSK signal
定時誤差估計采用O&M算法[11],設經過匹配濾波之后的基帶復信號為r(n),定時誤差為ε,令xn=|r(n)|2,每符號采樣點數為N,L0為統計符號的個數,則O&M算法可以表示為:

因為信號經過采樣率轉換變為4倍符號速率,即N=4,故式(6)可以展開為:

由式(8)可知Xm直接由加減法和累加運算即可求得,復雜度低,適合硬件實現。
計算定時誤差的同時數據送入FIFO緩存中延時,定時控制器根據估計值^ε得到插值控制信號mk和μk,其中mk用于控制FIFO讀出數據的下標index,決定取哪幾個采樣點用于插值,μk用于控制插值濾波器的工作。插值濾波器采用文獻[10]提出的高效Farrow結構,插值算法采用性能較好的立方插值,通過插值計算即可得到最佳采樣點的近似值。
1.5 相偏恢復
經過頻偏校正之后的信號仍有一個非常小的頻偏存在,此極小的頻偏相對于符號速率是一個慢變化量,可以看作是信號的相偏。前向相偏恢復的結構與定時恢復類似,如圖5所示。

圖5 突發PSK信號前向相偏恢復Fig.5 Feedforward phase recovery for burst PSK signal
突發PSK信號的相偏估計采用A.J.Viterbi和A.M.Viterbi提出的V&V算法[12],將MPSK信號表示為r(k)=ρ(k)ejφ(k),則載波相偏的估計值為:

式中,N0為用于估計的符號個數,F[ρ(n)]表示ρ(n)的函數,由文獻[13]的研究可知,為了達到最佳的估計性能,不同條件下F[ρ(n)]的取值如下:

與定時恢復類似,計算相偏的同時數據送入FIFO緩存中延時,^θ送入DDS中生成cos^θ-jsin^θ= e-j^θ,與FIFO的輸出數據進行復數乘法從而去除相偏的影響。
在加性高斯白噪聲信道下對文中的算法進行了仿真,仿真的具體條件為:BPSK調制,Rs=12 Ms/s, Eb/N0=12 dB,計算平方譜的符號長度M0=512,定時估計的符號長度L0=128,相偏估計的符號長度N0=64,初始頻偏為Δf=0.1Rs。經過頻偏校正和采樣率轉換后,信號的采樣率變為f′s=4Rs,殘余的頻偏小于2‰Rs,定時恢復后信號的星座圖見圖6。

圖6 定時恢復后的星座圖Fig.6 Constellation after timing recovery
由圖6看出,因為殘余微小頻偏的影響,信號的星座圖出現了旋轉。極小的頻偏相當于相偏,相偏恢復之后信號的星座圖見圖7,可以看出相偏恢復算法的性能較好,星座圖旋轉現象已經消除。

圖7 相偏恢復后的星座圖Fig.7 Constellation after phase recovery
根據上述算法,采用Xilinx公司V6系列的XC6VSX315T FPGA和TI公司的ADS5474采樣芯片完成了多速率突發PSK信號解調原理樣機的設計與實現。自行設計的突發信號源每秒發送15 000個突發包,實際的信號采用了差分編碼,即DBPSK調制,FPGA的系統工作時鐘為200 MHz,適用的突發信號符號速率范圍為8~50 Ms/s。通過ChipS-cope抓取了硬件調試的實際結果,信號的平方譜如圖8所示。

圖8 ChipScope顯示的平方譜Fig.8 Quadratic spectrum in ChipScope
定時恢復的信號星座圖見圖9,對差分編碼后的信號來說相偏的存在不會影響解調器的正常工作,故硬件實現時省略了相偏恢復模塊,定時恢復后的數據可直接送入差分判決模塊得到輸出結果。

圖9 ChipScope顯示的星座圖Fig.9 Constellation in ChipScope
在實驗室條件下,信號源重復發送同一段突發數據,通過調節模擬噪聲源和信號衰減器來實現不同信噪比,在FPGA內部編程將解調完的結果與信號源進行比對統計差錯符號的個數,得到不同信噪比下系統誤碼率的實際測試值,性能曲線見圖10。

圖10 突發解調器誤碼率性能Fig.10 BER performance of burst demodulator
圖10同時給出了所設計突發解調器誤碼性能的MATLAB仿真結果和理論曲線,算法仿真的性能與理論曲線相差小于0.2 dB。當Eb/N0較小時,硬件測試值與理論值相差較大,隨著Eb/N0增大,測試值逐漸趨近與理論曲線,由此可知在中高信噪比條件下文中所設計的突發解調樣機具有較好的性能。
文中提出一種適用于多速率突發PSK信號的無前導字全數字解調器設計,通過乒乓FIFO結構和符號速率轉換以滿足不同的符號速率,對符號速率和頻偏估計算法進行了詳細推導,算法的估計范圍較大,定時恢復和相偏恢復采用了前向算法,同時通過各級FIFO緩存保證解調過程不會損失任何符號。對算法進行了MATLAB仿真,最后根據設計方案完成了突發信號解調樣機的硬件實現,結果證明了算法的有效性和可靠性。
[1]ANDRONICO M,CASALE S,CORTE A L,et al.A New Algorithm for Fast Synchronization in a Burst Mode PSK Demodulator[C]//1995InternationalConferenceon Communications.Seattle:IEEE,1995:1641-1646.
[2]MENGALI U,MORELLI M.Data-Aided Frequency Estimation for Burst Digital Transmission[J].IEEE Transactions on Communications,1997,45(01):23-25.
[3]雷鳴,梅順良.全數字突發信號接收機實現[J].通信技術,2008,41(10):52-54.
LEI Ming,MEI Shun-liang.An All-digital Receiver for Burst Communication[J].Communications Technology, 2008,41(10):52-54.
[4]朱健軍,姚慶棟.全數字無前導字突發模式PSK接收機[J].電路與系統學報,2001,6(01):80-84.
ZHU Jian-jun,YAO Qing-dong.A Full Digital PreamblelessBurst-mode PSK Receiver[J].Journal of Circuits and Systems,2001,6(01):80-84.
[5]劉乃安,強剛,郭峰.突發擴頻通信中的中頻差分解調器及其性能[J].通信技術,1999(04):57-60.
LIU Nai-an,QIANG Gang,GUO Feng.IF Differential Demodulator and its Performance for Burst SS Communication[J].Communications Technology,1999(04):57-60.
[6]CHUANG J C I,SOLLENBERGER N R.Burst Coherent Demodulation with Combined Symbol Timing,Frequency Offset Estimation,and Diversity selection[J].IEEE Transactions on Communications,1991,39(07):1157-1164.
[7]CHALMERS H,SHENOY A,VERAHRAMI F B.A Digitally Implemented Preambleless Demodulator for Maritime and Mobile Data Communications[C]//Global Telecommunications Conference,1990,and Exhibition.San Diego: IEEE,1990:729-732.
[8]GINI F,GIANNAKIS G B.Frequency Offset and Symbol Timing Recovery in Flat-fading Channels:a Cyclostationary Approach[J].IEEE Transactions on Communications,1998,46(03):400-411.
[9]GARDNER F M.Interpolation in Digital Modems-Part I: Fundamentals[J].IEEE Transactions on Communications,1993,41(03):501-507.
[10]ERUP L,GARDNER F M.Interpolation in Digital Modems -PartII:Implementation and Performance[J].IEEE Transactions on Communications,1993,41(06):998-1008.
[11]OERDER M,MEYR H.Digital Filter and Square Timing Recovery[J].IEEE Transactions on Communications, 1988,36(05):605-612.
[12]VITERBI A J,VITERBI A M.Nonlinear Estimation of PSK-Modulated Carrier Phase with Application to Burst Digital Transmission[J].IEEE Transactions on Information Theory,1983,29(04):543-551.
[13]王軍杰.低信噪比TDMA信號同步技術研究[D].鄭州:解放軍信息工程大學,2010.
WANG Jun-jie.Synchronization Techniques for Low SNR TDMA Systems[D].Zhengzhou:PLA Information Engineering University,2010.

廖 明(1987—),男,碩士研究生,主要研究方向為通信與信息系統;
LIAOMing(1987-),male,graduate student,mainly engaged in communication and information system.
姚 軍(1962—),男,碩士,研究員,主要研究方向為通信信號處理;
YAO Jun(1962-),male,M.Sci.,research fellow,majoring in communication signal processing.
張 偉(1972—),男,碩士,研究員,主要研究方向為數字信號處理技術;
ZHANG Wei(1972-),male,M.Sci.,master,research fellow,mainly working at on digital signal processing.
王世練(1976—),男,博士,副教授,主要研究方向為無線通信。
WANG Shi-lian(1976-),male,Ph.D.,associate professor,mainly engaged in the research of wireless communication.
A Full-digital Burst-mode PSK Demodulator without Preamble for Multiple Symbol Rates
LIAO Ming1,2,YAO Jun2,ZHANG Wei2,WANG Shi-lian3
(1.Graduate School of China Academy of Engineering Physics,Mianyang Sichuan 621900,China; 2.Institute of Electronic Engineering,China Academy of Engineering Physics,Mianyang Sichuan 621900,China; 3.School of Electronic Science and Engineering,National Univ.of Defense Technology,Changsha Hunan 410073,China)
The application of burst-mode communication becomes even popular.The demodulator with preamble may reduce data transmission efficiency of the system.A scheme of full-digital burst-mode demodulator without preamble suitable for multiple symbol rates is proposed.FIFO operation and sampling rate conversion applicable to multiple symbol rates are designed.Feedforward algorithm is suggested for parameter estimation,timing synchronization and carrier recovery,the demodulation and high data transmission efficiency could be done without any loss of data symbols.The algorithm is easy to implement on hardware platform,and the results of simulation and hardware test indicate its effectiveness and reliability in the environment of moderate and high SNR.
multiple symbol rates;burst-mode;without preamble
TN911.72
A
1002-0802(2014)01-0018-06
10.3969/j.issn.1002-0802.2014.01.004