基于EM算法的OFDM系統(tǒng)聯(lián)合時(shí)變信道估計(jì)與信號檢測技術(shù)研究

任文成

（中國電子科技集團(tuán)公司第五十四研究所 河北 石家莊 050081）

判決指引信道估計(jì)作為一種半盲OFDM信道估計(jì)方法，能夠取得系統(tǒng)頻譜效率與信道估計(jì)性能的折中[1]。然而在時(shí)變衰落信道環(huán)境中，判決指引OFDM信道估計(jì)算法極易受到差錯傳播的影響而惡化系統(tǒng)性能[2]。為此，接收端必須依靠先進(jìn)的信號處理手段提高信道估計(jì)性能，然而時(shí)變信道參數(shù)與未知發(fā)送數(shù)據(jù)是緊密耦合的，彼此之間相互影響，故應(yīng)考慮最優(yōu)的聯(lián)合信道參數(shù)估計(jì)和信號檢測算法。然而最優(yōu)的聯(lián)合估計(jì)方法計(jì)算過于復(fù)雜，難以實(shí)現(xiàn)[3]。期望最大（EM）算法[4]可通過迭代計(jì)算處理，降低上述最優(yōu)聯(lián)合估計(jì)的計(jì)算復(fù)雜度，而且其性能也能夠逼近最優(yōu)，已證明是一種有效的次優(yōu)聯(lián)合估計(jì)算法，而且EM算法可擴(kuò)展性強(qiáng)，目前EM算法已在許多信號處理領(lǐng)域得到了廣泛的應(yīng)用[5]。

基于EM算法的信道估計(jì)[6]通常借助于信道的自回歸（AR）模型[7]實(shí)現(xiàn)快時(shí)變信道跟蹤，但是信道時(shí)變性進(jìn)一步增強(qiáng)時(shí)，AR模型擬合誤差增大，信道跟蹤性能急劇下降，而且在高階調(diào)制下，信號檢測性能對信道估計(jì)誤差極其敏感，導(dǎo)致系統(tǒng)誤碼率性能下降。本文將EM算法與Turbo均衡[8]相結(jié)合，通過在信道估計(jì)與信道譯碼之間相互交換編碼比特的軟信息，提高檢測可靠度的同時(shí)，EM算法本身內(nèi)部迭代過程實(shí)現(xiàn)了對信道參數(shù)和符號的聯(lián)合估計(jì)，從而改善了信號檢測對信道估計(jì)誤差的魯棒性，并基于泰勒級數(shù)展開對信道進(jìn)行建模[9]，在此基礎(chǔ)上配合自適應(yīng)濾波算法實(shí)現(xiàn)了對時(shí)變信道的跟蹤。

1 系統(tǒng)模型

1.1 OFDM系統(tǒng)發(fā)射與接收端模型

OFDM系統(tǒng)發(fā)射端框圖如圖1所示。信息比特流經(jīng)過信道編碼、交織、串并轉(zhuǎn)換以及符號映射后送入IFFT變換模塊，形成時(shí)域發(fā)射信號，添加循環(huán)前綴后，I個(gè)OFDM數(shù)據(jù)符號與IP個(gè)前導(dǎo)符號以時(shí)分復(fù)用的方式送入信道中。每個(gè)OFDM符號持續(xù)時(shí)間為Ts，含有N個(gè)子載波。

OFDM系統(tǒng)接收端框圖如圖2所示。接收端經(jīng)同步、去循環(huán)前綴以及 FFT 變換后，得到頻域接收信號 Y（i），0≤i＜I：

其中 N×1 維向 量 Y（i）=[Y0（i），…，YN-1（i）]T為第 i個(gè)OFDM符號的接收信號，Yn（i）為其中第 n個(gè)子載波接收信號，h（i）=[h0（i），…，hD-1（i）]T為第 i個(gè) OFDM 期間的信道沖激響應(yīng)，hd（i）為第 d 根徑的增益，Y（i）=[Y0（i），…，YN-1（i）]T為零均值加性高斯噪聲，方差為σ2n。對角矩陣X（i）第n個(gè)對角線元素為 Xn（i），是第 n 個(gè)子載波的發(fā)送符號，F(xiàn)H=[f0，…，fN-1]，fd=[1，ej2πn/N，…，ej2πn（D-1）/N]T。

圖1 OFDM系統(tǒng)發(fā)射端框圖Fig.1 The diagram of the OFDM transmitter

圖2 OFDM接收端框圖Fig.2 The diagram of the OFDM receiver

1.2 基于泰勒級數(shù)展開的信道模型

首先基于信道沖激響應(yīng)的泰勒級數(shù)展開，對時(shí)變多徑信道進(jìn)行建模[10]，實(shí)現(xiàn)對信道時(shí)變性的有效擬合，以便能夠?qū)π诺肋M(jìn)行跟蹤。 將hd（i）的u階差分記為h（u）d（i），對信道進(jìn)行U-1階近似，得到第d根徑的過程方程為：

其中，

由式（1）和式（2）得到，第i個(gè)OFDM符號的觀測方程：

E=[e1，eU+1…，e（D-1）U+1]T，其中 eq為 DU×1 維單位向量，只有第q個(gè)元素為1，其余元素為零。

如圖2所示，在幀頭IP個(gè)前導(dǎo)符號期間，進(jìn)行初始信道估計(jì)，在之后第i（i＞IP）個(gè)數(shù)據(jù)OFDM符號期間，基于 EM算法進(jìn)行聯(lián)合信道估計(jì)與符號檢測，迭代L次后生成各個(gè)編碼比特的對數(shù)似然比（LLR）[b（i，n，nB）]，0≤n≤N-1，1≤nB≤NB的后驗(yàn)概率，NB為每個(gè)符號含有的編碼比特?cái)?shù)。

2 基于EM算法的聯(lián)合信道估計(jì)與信號檢測

2.1 EM算法簡介

EM算法能夠解決觀測數(shù)據(jù)（設(shè)為Y）不充分情況下的參數(shù)（設(shè)為θ）最大似然估計(jì)問題。EM算法通過設(shè)置隱藏?cái)?shù)據(jù)（hidden data）i，與觀測數(shù)據(jù)一并組成完整數(shù)據(jù)（complete data），求解完整數(shù)據(jù)下的參數(shù)最大似然估計(jì)。具體來說，EM算法在一次迭代過程中分為兩步，取期望步驟（E步驟）與最大化步驟 （M步驟），E步驟求解一條件對數(shù)似然函數(shù)的期望值，M步則對E步得到的條件對數(shù)似然函數(shù)取最大，得到本次迭代的參數(shù)估計(jì)。上述過程可描述如下：

2.2 聯(lián)合信道估計(jì)與信號檢測

基于E步中獲取的符號軟估計(jì)，M步通過對條件對數(shù)似然函數(shù)期望取最大，得到本次迭代的信道估計(jì)：

其中

其中0＜λ≤1為遺忘因子。

3 計(jì)算機(jī)仿真

3.1 仿真參數(shù)

對上述算法在進(jìn)行仿真，仿真參數(shù)為子載波數(shù)為64，循環(huán)前綴長度為16個(gè)采樣點(diǎn)，調(diào)制方式為QPSK或64QAM，每幀配備1個(gè)前導(dǎo)符號，即IP=1，后跟10個(gè)數(shù)據(jù)OFDM符號，采用1/2碼率的卷積編碼，塊交織器長度為24×16，采用軟輸入軟輸出的維特比信道譯碼算法[14]，歸一化多普勒頻移為fDts=0.05，信道多徑數(shù)為16，采用歸一化的指數(shù)型時(shí)延功率譜和Jakes多普勒功率譜。

3.2 仿真結(jié)果分析

圖3顯示了QPSK調(diào)制和64QAM調(diào)制方式，采用本文基于EM的聯(lián)合信道估計(jì)與檢測算法，系統(tǒng)誤碼率與比特信噪比的關(guān)系，迭代次數(shù)分別為2和4次。由圖中可以看出，相同迭代次數(shù)下，QPSK調(diào)制下算法的性能明顯優(yōu)于64QAM調(diào)制下的性能[15]，而且由于調(diào)制階數(shù)的增加，降低了系統(tǒng)對信道噪聲的抑制能力，在高信噪比區(qū)域，無論QPSK還是64QAM均會有不同程度的誤碼平臺。

圖3 不同調(diào)制方式下的BER性能比較Fig.3 The BER performance under the differentmodulation

圖4顯示了QPSK調(diào)制下本文算法和傳統(tǒng)的MMSE均衡算法性能的比較，其中MMSE均衡采用的信道由幀頭前導(dǎo)符號處的信道估計(jì)提供，由圖中可以看出，當(dāng)信道時(shí)變性增強(qiáng)時(shí)，受限于差錯傳播以及OFDM符號內(nèi)子載波間干擾的影響，MMSE均衡誤碼率性能較差，不能滿足系統(tǒng)需求，而基于EM算法的聯(lián)合信道估計(jì)和檢測則能夠提高系統(tǒng)對時(shí)變信道的魯棒性。

圖4 與傳統(tǒng)MMSE均衡BER性能的比較Fig.4 The comparison with the conventional MMSE equalization

圖5顯示了系統(tǒng)BER性能隨迭代次數(shù)的變化情況，可看出第二次迭代所獲取的增益最大，而三次迭代之后的BER性能則改善不大，算法具有良好的收斂性能，實(shí)際應(yīng)用時(shí)應(yīng)取三次迭代為宜。

圖5 BER性能隨迭代次數(shù)的改善情況Fig.5 The BER performance with the iteration number

4 結(jié) 論

判決指引OFDM信道估計(jì)在時(shí)間選擇性衰落信道中應(yīng)用時(shí)，易受到的差錯傳播影響而惡化系統(tǒng)性能。針對此問題，文中提出的基于EM算法和Turbo均衡的聯(lián)合信道估計(jì)與符號檢測算法提高了信道估計(jì)和信號檢測的可靠度，同時(shí)采用自適應(yīng)濾波跟蹤時(shí)變信道，有效減小了差錯傳播的影響。仿真表明，該算法在快衰落信道環(huán)境中性能明顯優(yōu)于傳統(tǒng)的OFDM信號檢測算法，而且算法具有良好的收斂性能。

[1]趙琦，王海濤，鄒光南.基于LTE的單載波頻分多址技術(shù)研究[J].電子設(shè)計(jì)工程，2011，19（4）:47-54.

[2]李軒，韓笑，關(guān)慶陽，等.基于濾波方法的OFDM信道估計(jì)研究[J].電子設(shè)計(jì)工程，2014，22（12）:145-151.

[3]Hlawatsch F，Matz G.Wireless Communications over Rapidly Time-Varying Channels[M].Elsevier Press， First Edition，2011.

[4]Dempster A P，Larid N M，Rubin D B.Maximum likelihood from incomplete data via the EM algorithm[J].Royal Statistical Society，1977，39（1）:1-38.

[5]Borran M J，Aazhang B.EM-Based Multiuser Detection in Fast Fading Multipath Environments[J].EURASIP Journal on Applied Signal Processing，2002（8）:787-796.

[6]Dogan H，Panayirci E，Poor H V.Low-complexity joint data detection and channel equalization for highlymobile orthogonal frequency division multiplexing systems[J].IET Communications，2010，4（8）:1000-1011.

[7]Wang H，Chang P.On verifying the first-order Markovian assumption for a Rayleigh fading channel model[J].IEEE Trans.Vehicular Technology，1996，45（2）:353-357.

[8]Tuchler M，Singer A C.Turbo equalization:an overview[J].IEEE Transactions on Information Theory，2011，57 （2）:920-952.

[9]LIU Yin-sheng，TAN Zhen-hui.Joint estimation of channel impulse response and carrier frequency offset for OFDM systems[J].IEEE Transactions on Vehicular Technology，2011，60（9）:4645-4649.

[10]王亮，王文策，程斌，等.基于Simulink的發(fā)電機(jī)勵磁系統(tǒng)整流單元的建模與仿真[J].陜西電力，2011（7）:25-28.

[11]施先旺，劉婷婷，李國良.采用有限狀態(tài)機(jī)實(shí)現(xiàn)控制指令的可靠檢測[J].火箭推進(jìn)，2011（5）:63-68.

[12]王坤.基于SVPWM的STATCOM設(shè)計(jì)及實(shí)現(xiàn)[J].陜西電力，2011（1）:11-15.

[13]邱進(jìn).一種應(yīng)用于多斷口光控真空斷路器的同步控制系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J].供用電，2015（5）:69-73.

[14]徐進(jìn)，徐榮森，梅正茂.基于1394b總線仿真設(shè)備的WDM驅(qū)動開發(fā)[J].電子設(shè)計(jì)工程，2015（2）:63-66.

[15]尤裕榮，袁洪濱，王春民，等.PWM高速開關(guān)閥動態(tài)調(diào)節(jié)特性仿真研究[J].火箭推進(jìn)，2013（1）:24-29.