史德生, 武楠, 王華, 李智信, 匡鏡明
(北京理工大學 信息與電子學院,北京 100081)
為了保證通信鏈路的可用度,衛(wèi)星通信系統(tǒng)必須按照最差信道環(huán)境,預(yù)留一定傳輸裕量下進行設(shè)計,這對衛(wèi)星轉(zhuǎn)發(fā)器功率和容量都存在巨大浪費. 自適應(yīng)傳輸技術(shù)能夠根據(jù)信道質(zhì)量的變化而自適應(yīng)地改變調(diào)制模式和信道編碼保護級別[1],在衛(wèi)星現(xiàn)有轉(zhuǎn)發(fā)器功率和帶寬不變前提下提高傳輸容量和傳輸質(zhì)量. 自適應(yīng)衛(wèi)星傳輸系統(tǒng)對用戶的物理層配置進行逐幀優(yōu)化,使系統(tǒng)在Ku/Ka及更高的頻段以更低的鏈路裕量工作. 因此,系統(tǒng)性能不再受限于最差系統(tǒng)鏈路條件,而是取決于平均系統(tǒng)鏈路條件. 歐洲的DVB-S2標準就將自適應(yīng)編碼調(diào)制模式(ACM)作為其交互式業(yè)務(wù)應(yīng)用的標準模式之一. 本文研究了寬帶自適應(yīng)解調(diào)器中關(guān)鍵的算法,在得到了Matlab性能仿真基礎(chǔ)上,研究了基于FPGA的衛(wèi)星高速寬帶自適應(yīng)調(diào)制解調(diào)器總體結(jié)構(gòu)和硬件實現(xiàn)方案. 研制的調(diào)制解調(diào)器編碼采用LDPC和BCH碼的級聯(lián)編碼,支持BPSK、QPSK、8PSK、16APSK等13種編碼調(diào)制方式. 在誤比特率為10-6時,工作信噪比范圍從1~13 dB,動態(tài)范圍達到12 dB,平均步進約為1 dB,頻譜效率可從0.67 bit/s/Hz(BPSK調(diào)制2/3碼率)~3.5 bit/s/Hz(16APSK調(diào)制8/9碼率).
自適應(yīng)傳輸?shù)幕鶐盘枎Y(jié)構(gòu)需要支持多種調(diào)制方式、多種碼長和碼率的LDPC編碼. 考慮到接收機的實現(xiàn)復(fù)雜度,對寬帶高速自適應(yīng)傳輸中的4種調(diào)制方式(BPSK、QPSK、8PSK和16APSK)設(shè)計了統(tǒng)一的幀結(jié)構(gòu),如圖1所示.
傳統(tǒng)的基于導(dǎo)頻輔助的幀同步檢測器需要在傳輸數(shù)據(jù)幀中插入導(dǎo)頻或者幀頭序列,這種方法占用了一定額外信號帶寬和功率,并且需要大量乘法器實現(xiàn)相關(guān)處理;另外,為了進一步提高幀效率,需要設(shè)計較短幀頭和導(dǎo)頻序列,但幀頭變短又會帶來幀同步性能惡化. 為此,提出了通過使用信道編碼(LDPC碼)的碼字特征的輔助幀同步方法,提升了高速衛(wèi)星自適應(yīng)傳輸系統(tǒng)的頻譜效率.
圖2是幀同步算法結(jié)構(gòu)示意圖. 在幀捕獲的初始階段,將切換控制開關(guān)連接到B,接收到的信號作硬判決后輸入幀同步檢測器,利用LDPC譯碼前信道符號的硬判決計算校驗方程滿足的比例,再根據(jù)門限比較快速判斷是否檢測到幀同步位置. 當幀同步的捕獲完成后,A與切換控制開關(guān)相連,系統(tǒng)進入幀同步跟蹤狀態(tài). 由于在幀同步跟蹤狀態(tài)LDPC譯碼器已經(jīng)開始工作,因此可將其譯碼結(jié)果輸入幀同步檢測器,利用LDPC譯碼信息對幀同步位置進行跟蹤,提高幀同步跟蹤的可靠性.
當采用連續(xù)的M幀進行幀同步捕獲時,則幀同步位置m的估計[2]為
(1)
(2)
式中:N為LDPC碼校驗矩陣H的行數(shù);p0為在非幀同步位置時,幀捕獲檢測器的輸出大于門限的概率質(zhì)量函數(shù);p1為在幀同步位置時,幀捕獲檢測器的輸出大于門限的概率質(zhì)量函數(shù).
在自適應(yīng)傳輸中,每一幀傳輸信號的調(diào)制方式都有可能發(fā)生變化,為每一種模式都設(shè)計相應(yīng)的載波同步算法將顯著增大設(shè)計和實現(xiàn)復(fù)雜度,同時低信噪比下載波同步門限和高階調(diào)制下載波恢復(fù)精度問題都增加了載波同步難度. 載波頻率估計精度最高的方法為數(shù)據(jù)(導(dǎo)頻)輔助載波頻率估計算法. 由于本設(shè)計中的數(shù)據(jù)幀中導(dǎo)頻符號分散在不同位置,不能直接采用傳統(tǒng)方法. 因此,研究給出了一種改進的L&R[3]載波頻率估計算法,其歸一化數(shù)學表達式為
(3)
載波頻率恢復(fù)的結(jié)構(gòu)框圖如圖3所示.
為了選擇與信道狀態(tài)相適應(yīng)的編碼調(diào)制方式,衛(wèi)星發(fā)送端必須獲得信道的衰落特性,從而達到鏈路自適應(yīng)的目的. DVB-S2系統(tǒng)中通過使用基帶信號中的導(dǎo)頻符號累積進行信噪比估計,但需要累積較長時間導(dǎo)頻段才能獲得標準偏差小于0.2 dB的估計精度,降低了信噪比估計的實時性. 本文設(shè)計的編碼輔助APSK信號信噪比估計器對一幀數(shù)據(jù)通過在譯碼器和估計器之間交換軟信息提高參數(shù)估計性能,使估計精度和實時性均優(yōu)于傳統(tǒng)算法[4].
假設(shè)x=[x1x2…xN]T和r=[r1r2…rN]T分別是編碼后發(fā)送和接收到的符號序列,且rk=hxk+nk,k=1,2,…,N,其中h為信道增益,nk的均值為0,方差為σ2的復(fù)高斯白噪聲.
根據(jù)估計理論[5],參數(shù)矢量b的最大似然估計為
(5)
(6)
(7)
(8)
(9)
對于數(shù)據(jù)集y的對數(shù)似然函數(shù)
lnp(y|b)=lnp(r|b)+lnp(x).
(10)
將式(7)、(10)代入式(8),并消去與參數(shù)b無關(guān)項得到
(11)
分別對h和σ2求導(dǎo)并令其等于0,簡化整理后可以得到第i次解調(diào)譯碼迭代后,h和σ2的估計表達式分別為
(12)
(13)
(14)
事實上,基于EM算法的編碼輔助APSK信號信噪比估計器是編碼輔助APSK信號最大似然信噪比估計器的一種實現(xiàn)形式.
(15)
(16)
(17)
(18)
以上兩式最大值相等,都等于用戶可以接受的乒乓切換概率β,且
(19)
本文實現(xiàn)的“門限滯后法”實現(xiàn)框圖如圖4. 門限移位技術(shù)可以有效降低信道估計誤差對系統(tǒng)性能的影響,而門限移位距離的又將直接影響系統(tǒng)的丟包率和傳輸效率. 實際中,作者通過理論計算與Matlab仿真,確定了一組Ka波段下丟包率與傳輸效率折中的雙門限自適應(yīng)切換閾值.
寬帶自適應(yīng)調(diào)制解調(diào)器總體設(shè)計結(jié)構(gòu)如圖5所示. 在寬帶自適應(yīng)調(diào)制器中,ACM控制器根據(jù)自適應(yīng)業(yè)務(wù)調(diào)度器輸入的信道狀態(tài)信息控制信道編碼的碼長和碼率、調(diào)制方式以及組幀的結(jié)構(gòu)等等,實現(xiàn)編碼調(diào)制模式逐幀可變. 在可編程邏輯器件FPGA中實現(xiàn)了所有基帶信號處理模塊,基帶信號經(jīng)過高速(128 MS/s)D/A和濾波器后,通過I/Q調(diào)制器調(diào)制到L波段. 寬帶自適應(yīng)解調(diào)器同樣由I/Q解調(diào)器、高速A/D以及FPGA組成. 接收信號經(jīng)過符號同步和幀同步后將提取當前幀的信令信息,從而控制解調(diào)器的載波同步、比特似然信息計算(解映射)模塊、以及自適應(yīng)信道譯碼器. 解調(diào)器通過信噪比估計模塊對信道狀態(tài)進行實時估計,利用反饋信道提供給自適應(yīng)業(yè)務(wù)調(diào)度器.
根據(jù)式(1)可知,計算歸一化頻率需要計算Rl(m)、累加Rl(m)并且求幅角. 通過觀察看到其中的信號延時以及乘加部分可視為一個N階的復(fù)輸入FIR濾波器. 根據(jù)L&R算法DSP實現(xiàn)結(jié)構(gòu)[3],簡化了改進的L&R算法的結(jié)構(gòu),如圖6所示,進一步減少了硬件實現(xiàn)復(fù)雜度. 由于FIR濾波器采用了流水線結(jié)構(gòu),因此需要對未取共扼路徑上的信號進行延時「N/2?Ts.
在接收機系統(tǒng)解調(diào)后進行LDPC譯碼之前需要解調(diào)符號的軟判決似然信息. 解映射模塊接收到的符號rn,n=1,2…的比特構(gòu)成是bM-1…bk…b1b0,M為調(diào)制階數(shù). 第k個比特bk的對數(shù)似然比表示為
(20)
式中:i=0,…,2M-1;r為接收到的信號矢量;C1為bk=1的星座點集合;C0為對應(yīng)bk=0的星座點集合;S為集合中的具體星座點. 通常適用于APSK信號的軟判決解映射方法由于需要復(fù)雜的對數(shù)操作和指數(shù)操作,硬件實現(xiàn)十分復(fù)雜. 提出簡化了的16APSK軟判決解映射實現(xiàn)結(jié)構(gòu),如圖7中所示. 首先通過劃分星座點集合,改進似然信息計算方法,從數(shù)學表達上把MAX方法[8]的pmax(r|S=Si)=-‖r-Si‖2中的負號提取出來,得到bk的似然信息表達式為
(21)
(22)
經(jīng)過分析發(fā)現(xiàn)16APSK調(diào)制,對每個比特bk最多只要對3個點求距離.
綜合考慮,這種方法比MAX方法節(jié)省了17~20個乘法器,節(jié)省53%~62%;加法比MAX方法節(jié)省了39~42個,節(jié)省80%~85%. 該方法解決了采用高階調(diào)制和迭代譯碼的系統(tǒng)軟判決解映射方法資源消耗過多的問題,同時又不損失系統(tǒng)的整體性能.
本文設(shè)計的自適應(yīng)調(diào)制解調(diào)器發(fā)送端(調(diào)制加編碼)和接收端(解調(diào)加譯碼)均采用了Xilinx公司研制的Kintex-7系列FPGA芯片,具體型號為XC7K325T-1FFG676. 表1給出了ACM發(fā)射機、接收機、編碼器、譯碼器的資源消耗情況. 其中發(fā)射機和接收機Slice資源占用分別為10%和37%.
表1 ISE軟件MAP之后資源占用情況Tab.1 ISE map report of device utilization
在基帶環(huán)路中加入數(shù)字高斯白噪聲,然后通過誤碼儀(AV5233C)對誤比特率進行測量. 不同編碼調(diào)制方式下的誤碼性能測試結(jié)果如圖8、圖9所示,可以看出在誤比特率為10-6時,其硬件實現(xiàn)與理論仿真結(jié)果相差小于0.2 dB.
搭建如圖10所示中頻環(huán)路,利用泰克公司的實時頻譜分析儀RSA6114A觀測調(diào)制器在L波段的頻譜和星座圖,實時測量矢量幅度誤差(EVM)值. 如圖11,設(shè)置的信道符號速率為8.192 MS/ps(symbol per second),采用8PSK調(diào)制,觀測得到的峰值EVM為5.237%. 解調(diào)后的誤比特率與理論值相比,Eb/N0相差小于1.0 dB.
在ACM系統(tǒng)中,信源端根據(jù)接收端估計的信號與噪聲和干擾之比調(diào)整傳輸模式,為了充分利用鏈路自適應(yīng)帶來的系統(tǒng)增益,同時將鏈路中斷概率控制在一定范圍內(nèi),本文對研制的自適應(yīng)傳輸系統(tǒng)在模擬衛(wèi)星信道環(huán)境下進行實時測量和估計. 由圖5中ACM控制器選擇自適應(yīng)編碼調(diào)制控制比特數(shù)據(jù)流的產(chǎn)生. 測試中模擬衛(wèi)星信道的雨衰情況,產(chǎn)生信噪比變化率最大為0.65 dB/s的雨衰情況,符合衛(wèi)星信道最差的雨衰情況[9].
圖12給出了雨衰發(fā)生時真實信噪比與估計信噪比的變化情況. 從圖中可以看出,估計信噪比與真實信噪比的變化趨勢一致. 圖13中給出了在在圖12所示雨衰發(fā)生時,真實信噪比下,使用原始門限To下的ModCod(調(diào)制方式及碼率)、使用移位門限Tth下的ModCod以及系統(tǒng)實際采用的ModCod三者的變化示意圖. 經(jīng)過多次測試,得到該雨衰發(fā)生時共出現(xiàn)5次乒乓切換效應(yīng),平均誤比特率為3.2×10-7,平均頻譜效率損失為0.307 4 bit/s/Hz,. 可見,系統(tǒng)鏈路自適應(yīng)器實現(xiàn)了對雨衰變化的實時有效跟蹤和匹配.
本文研究了衛(wèi)星高速自適應(yīng)調(diào)制解調(diào)器算法與硬件實現(xiàn)方案. 提出的頻偏估計器改進結(jié)構(gòu)和高階調(diào)制軟判決解映射簡化結(jié)構(gòu),實現(xiàn)復(fù)雜度比傳統(tǒng)方法節(jié)省了60%以上,顯著降低了系統(tǒng)的資源開銷. 對研制的自適應(yīng)調(diào)制解調(diào)器實測結(jié)果表明,設(shè)計實現(xiàn)的調(diào)制解調(diào)器基帶性能與理論仿真結(jié)果相差小于0.2 dB. L波段信號在高符號速率高階調(diào)制下誤差矢量幅度小于10%,并且環(huán)路性能與理論值之差小于1 dB. 在衛(wèi)星信道最大變化率達0.65 dB/s下,自適應(yīng)環(huán)路仍可以跟蹤信道變化,平均誤比特率和譜效率損失分別為0.3 bit/s/Hz和10-7,均滿足了設(shè)計指標要求.