關(guān)慶陽 趙洪林 郭慶

(哈爾濱工業(yè)大學(xué)通信研究所,黑龍江哈爾濱 150001)

作為未來通信的多址接入技術(shù),正交頻分多址接入(OFDMA)技術(shù)已經(jīng)得到更廣泛的應(yīng)用.由于正交頻分復(fù)用(OFDM)的子載波的正交性為多用戶的接入提供了更便利的條件,來自不同用戶的傳輸信息在時(shí)域內(nèi)占用互相重疊的子載波,因而,OFDMA系統(tǒng)被IEEE 802.16E選為多址接入技術(shù)的物理層標(biāo)準(zhǔn)之一[1].與OFDM系統(tǒng)相比,OFDMA系統(tǒng)對(duì)子載波頻率偏移更為敏感,OFDMA系統(tǒng)需要更多的子載波和更窄的載波間隔,同時(shí)每個(gè)接入用戶的子載波都有不同的頻率偏移[2],因而導(dǎo)致整個(gè)系統(tǒng)的干擾消除更為困難.與OFDMA下行鏈路相比,多用戶頻率偏移對(duì)上行鏈路的影響更大些[3].為克服這種多載波頻率偏移的影響,最小二乘法(LS)被提出[4],但是其需要獲知頻率偏移的先驗(yàn)知識(shí).Huang等[5]提出一種頻域循環(huán)卷積方法,在接收端進(jìn)行載波頻率偏移(CFO)補(bǔ)償以及多用戶消除,但是這種算法僅限于載波頻率偏移的波動(dòng)間隔很大的情況下,補(bǔ)償效果較好.文獻(xiàn)[6]提出一種SAGE(Space-Alternating Generalized Expectation-Maximization)算法來消除OFDMA上行鏈路非同步引入的干擾,但是算法的復(fù)雜度高.文獻(xiàn)[7]提出最大似然估計(jì)算法來進(jìn)行干擾估計(jì),雖然降低了復(fù)雜度,但是難以實(shí)現(xiàn).文獻(xiàn)[8]提出最小均方誤差(MMSE)均衡算法來消除上行鏈路干擾,但是隨著子載波數(shù)目的增大,算法的復(fù)雜度也增加.文獻(xiàn)[9]指出串行干擾消除方法(SIC)的性能與干擾消除順序有關(guān).文獻(xiàn)[10]提出按照功率大小進(jìn)行子載波排序,然后采用基于最小均方準(zhǔn)則的串行干擾消除算法(MMSE-SIC)消除頻率偏移干擾.文獻(xiàn)[11]提出按照信干比順序進(jìn)行子載波排序.文獻(xiàn)[9-11]提出的SIC算法受第一級(jí)判決的約束,如果第一級(jí)判決錯(cuò)誤,積累的判決誤差將引起后面各級(jí)性能的下降,嚴(yán)重降低整個(gè)系統(tǒng)的性能.

文中提出一種時(shí)域自適應(yīng)最小均方誤差算法(TD-LMS)在時(shí)域內(nèi)消除載波頻率偏移引起的干擾.首先對(duì)接收的頻率偏移信號(hào)進(jìn)行預(yù)估計(jì),利用估計(jì)量對(duì)頻率偏移值進(jìn)行調(diào)整,然后采用TD-LMS算法在時(shí)域內(nèi)消除殘留頻率偏移的干擾.

1 系統(tǒng)模型

假設(shè)OFDMA系統(tǒng)上行鏈路共有T個(gè)用戶與基站進(jìn)行通信;并且每個(gè)用戶已經(jīng)預(yù)先被分配了特定的子載波.用戶i共有N個(gè)子載波,用戶i所發(fā)送的一個(gè)數(shù)據(jù)幀的信息為Xi(0),Xi(1),…,Xi(N-1),然后通過反傅里葉變換調(diào)制送上各自獨(dú)立的子信道.發(fā)送前,需要在每個(gè)數(shù)據(jù)幀前添加循環(huán)前綴形成保護(hù)間隔來防止信道干擾.設(shè)M(i)表示用戶i對(duì)子載波的占用情況,如果用戶i占用了第k個(gè)子載波,此時(shí)M(i)可以表示為M(i)(k,k)=1.用戶i發(fā)送的第n個(gè)信號(hào)經(jīng)調(diào)制可以表示為

用戶i信號(hào)經(jīng)過信道后可以表示為

式中:*表示線性卷積;hi(n)為用戶i的信道響應(yīng).引入頻率偏移后,用戶i的時(shí)域信號(hào)為

式中:ξi為用戶i引入的載波頻率偏移;g(n)表示均值為零,方差為 σ2的加性高斯白噪聲.將接收的所有用戶信號(hào)y(n)進(jìn)行離散傅里葉變換(DFT)解調(diào)后可得

式中:Si(k)為用戶i接收信號(hào)的頻域表示;G(k)為高斯白噪聲的頻域表示;P(ξi)為頻率偏移引起用戶i的頻域N×N干擾矩陣,可以表示為

矩陣元素p(l,ξi)可以表示為

式中:l為子載波索引號(hào)的間隔距離.令m表示用戶i中除載波k外的其余載波,u表示OFDMA系統(tǒng)中用戶j占用的子載波,其載波數(shù)設(shè)為θ,則pj(u-k,ξj)為用戶j的由載波頻率偏移引起的對(duì)用戶i的干擾,將頻率偏移引入的干擾分離出來,式(4)可化簡為

公式(7)第1項(xiàng)表示用戶i接收的無載波干擾信號(hào),第2項(xiàng)表示用戶i中,除載波k外的其余載波m的干擾,即載波間干擾,第 3項(xiàng)表示用戶 j的所有載波對(duì)用戶i的多址接入干擾,第 4項(xiàng)為加性高斯白噪聲的頻域表示.而文中所采用的TD-LMS算法就是在時(shí)域內(nèi),降低頻率偏移ξi對(duì)OFDMA系統(tǒng)上行鏈路的影響.

2 OFDMA上行鏈路頻率偏移消除

為實(shí)現(xiàn)OFDMA上行鏈路載波頻率偏移的消除,首先在時(shí)域內(nèi)對(duì)帶頻率偏移的信號(hào)進(jìn)行估計(jì),并采用估計(jì)量對(duì)載波頻率偏移值進(jìn)行預(yù)消除,然后通過TD-LMS算法消除系統(tǒng)殘留載波頻率偏移的干擾.這是因?yàn)門D-LMS算法直接處理較大的頻率偏移干擾時(shí),性能有限,所以需要在前端進(jìn)行頻率偏移干擾預(yù)消除,以提高消除性能.

2.1 頻率偏移時(shí)域估計(jì)及預(yù)消除

令頻率偏移的均值 ˉζ為用戶 i頻率偏移 ξi的估計(jì)值,r(n)為對(duì)應(yīng)于yi(n)、收發(fā)端已知的無頻率偏移理想信號(hào),令其所占用的載波數(shù)為L,并且L?N.利用已知信號(hào)r(n)進(jìn)行頻率偏移粗估計(jì)及預(yù)消除,同時(shí)可以將已知信號(hào)作為參考信號(hào)來完成TD-LMS算法的自適應(yīng)學(xué)習(xí).由此,頻率偏移可以表示為:

公式(3)代入公式(10)后,可得

式中:ρi為預(yù)消除頻率偏移后的殘留頻率偏移,滿足.后端采用的時(shí)域TD-LMS算法就是對(duì)殘留頻率偏移ρi進(jìn)行自適應(yīng)消除.

2.2 用TD-LMS算法消除殘留頻率偏移

經(jīng)預(yù)消除后,殘留頻率偏移 ρi數(shù)值較小,可以通過時(shí)域自適應(yīng)TD-LMS算法消除,具體的消除過程如圖1所示.

圖1 TD-LMS算法頻率偏移消除原理Fig.1 Elimination principle for frequency offset with TD-LMS algorithm

TD-LMS算法的消除過程如下:載波頻率偏移預(yù)消除后,帶有殘留頻率偏移ρi,將信號(hào)zi(n)作為TD-LMS濾波器的輸入,x(n)為經(jīng)過TD-LMS算法處理后的輸出信號(hào).當(dāng)n≤L時(shí),r(n)作為收發(fā)端預(yù)先共知的參考信號(hào),經(jīng)過自適應(yīng)學(xué)習(xí),TD-LMS算法達(dá)到收斂狀態(tài);當(dāng)L＜n＜N-1時(shí),濾波器判決后的輸出信號(hào)r(n)作為參考信號(hào)來完成對(duì)剩余偏差信號(hào)zi(n)的消除.定義時(shí)域代價(jià)函數(shù)為

式中:E為對(duì)信號(hào)求期望值運(yùn)算,H表示向量的共軛轉(zhuǎn)置.代價(jià)函數(shù)對(duì)權(quán)值向量w(n)求梯度可得

經(jīng)運(yùn)算可得

式中:ΔJ(n)為代價(jià)函數(shù)的梯度表示;Rz,z(n)為接收的第n個(gè)時(shí)域信號(hào)的自相關(guān)函數(shù),Rz,z(n)=E[zi(n)?zi(n)];Rz,r為接收第n個(gè)時(shí)域信號(hào)與參考信號(hào)的互相關(guān)函數(shù),Rz,r=E[zi(n)rH(n)].通過權(quán)值向量的遞推關(guān)系可以得到

式中:μ為自適應(yīng)步長.更新權(quán)值向量,可以得到

將頻率偏移干擾代入公式(16),誤差信號(hào)e(n)可以表示為

并且

式中:D為信號(hào)向量經(jīng)運(yùn)算后得到的實(shí)常數(shù),D= r(n)rH(n).

更新權(quán)值矩陣w(n),可以得到

取極限可得

輸出信號(hào)可以表示為

式中:Θ(ρi)為算法消除后的干擾因子.從公式(21)可得,TD-LMS算法可以獲得將帶有頻率偏移的信號(hào)zi(n)與消除乘積的效果,并且可以提取出無頻率偏移干擾理想信號(hào)si(n);事實(shí)上,當(dāng)ρi≠0時(shí), Θ(ρi)乘積因子會(huì)改變TD-LMS算法的輸出幅度,并會(huì)帶來固定相位旋轉(zhuǎn).由于TD-LMS算法所針對(duì)的是經(jīng)過預(yù)消除后的殘留頻率偏移,頻率偏移值很小,所以幅度以及相位影響可以忽略.

3 仿真結(jié)果及分析

3.1 仿真條件

采用Matlab軟件驗(yàn)證文中提出的TD-LMS算法.假定OFDMA上行鏈路系統(tǒng)的子載波總數(shù)為256,用戶數(shù)為 4,DFT窗口為 256,傳輸帶寬設(shè)定為30MHz.設(shè)定每個(gè)用戶都有獨(dú)立的子信道,基帶信號(hào)映射方式為正交相移鍵控(QPSK),頻率偏移的范圍設(shè)定為均勻分布,令分布范圍為(-ξmax,ξmax),這里ξmax是每個(gè)子用戶的最大頻率偏移范圍.

3.2 星座圖仿真

通過對(duì)星座圖的描繪能更加直觀地看到 TD-LMS算法載波頻率偏移的消除能力.本研究所采用的映射方式為QPSK.為了更清晰地看出基帶信號(hào)星座圖的變化情況,令信噪比SNR=10dB,取每數(shù)據(jù)幀OFDM數(shù)據(jù)符號(hào)為2000個(gè),并且設(shè)單用戶的固定頻率偏移值為 0.15,圖2示出了帶有載波頻率偏移的子載波信號(hào)星座圖和經(jīng)過TD-LMS算法消除后的子載波信號(hào)星座圖.從星座圖可以清晰看到TD-LMS算法能夠消除由于頻率偏移所引入的干擾.

3.3 誤碼率分析

圖3示出了對(duì)于所有用戶,ξmax=0.1,0.2,在不同信噪比下的平均誤碼率曲線,這包括無頻率偏移干擾的情況、頻率偏移干擾無消除的情況、文獻(xiàn)[8]提出的MMSE修正算法、文獻(xiàn)[4]提出的LS修正算法以及文中提出的TD-LMS自適應(yīng)消除算法.本次仿真所采用的自適應(yīng)TD-LMS算法的步長μ=0.075,權(quán)矩陣w取單位矩陣.圖 3的條件是假定各子信道都是獨(dú)立的.從圖3中的曲線可以看到,與存在頻率偏移干擾的情況相比,MMSE算法在一定程度上抑制加性白噪聲的干擾,性能優(yōu)于LS算法,文中提出的TD-LMS算法對(duì)頻率偏移消除的誤碼率性能明顯提高,對(duì)由載波頻率偏移引起的多用戶干擾消除效果明顯.當(dāng)ξmax= 0.2時(shí),與LS算法相比,文中提出的TD-LMS算法在相同的誤碼率條件下,平均可獲得3dB的增益.但是隨著頻率偏移范圍的擴(kuò)大,文中提出的算法的曲線與理論上的無頻率偏移干擾的曲線偏離較大,這主要是由于較大的載波頻率偏移范圍引起了干擾波動(dòng).

圖3 載波頻率偏移分布于(-0.1,0.1)、(-0.2,0.2)時(shí)的平均信噪比的誤碼率曲線Fig.3 Bit-error-rate curves under different average signal-tonoise ratio when carrier frequency offset distributing in (-0.1,0.1),(-0.2,0.2)

3.4 TD-LMS算法自適應(yīng)學(xué)習(xí)曲線

圖4 TD-LMS算法自適應(yīng)學(xué)習(xí)曲線Fig.4 Adaptive learning curves for TD-LMS algorithm

采用TD-LMS算法的均方誤差的學(xué)習(xí)曲線見圖4.從圖 4中可以看出,經(jīng)過一些數(shù)據(jù)包的加權(quán)運(yùn)算,均方誤差明顯減小,而且經(jīng)過 130個(gè)數(shù)據(jù)包的學(xué)習(xí),均方誤差可以降低到 10-4,這也可以證明TDLMS算法在時(shí)域內(nèi)可以完成對(duì)OFDMA系統(tǒng)頻率偏移的消除.

4 結(jié)語

文中提出了一種基于時(shí)域自適應(yīng)TD-LMS算法的補(bǔ)償方案來消除OFDMA上行鏈路的各載波頻率偏移引起的干擾,TD-LMS算法自適應(yīng)學(xué)習(xí)的能力較強(qiáng),結(jié)合時(shí)域頻率偏移估計(jì)及預(yù)消除過程,可以自適應(yīng)地有效消除不同載波頻率偏移引起的干擾.仿真結(jié)果表明,該補(bǔ)償算法對(duì)由載波頻率偏移引起的干擾補(bǔ)償效果明顯.

[1] The Institute of Electrical and Electronics Engineering. IEEE 802.16E-2005.IEEE standard for local and metropolitan area networks part 16:air interface for fixed and mobile broadband wireless access systems amendment 2: physicaland medium access control layers for combined fixed and mobile operation in licensed bands[S].New York:The Institute of Electrical and Electronics Engineering,2005:126-129.

[2] Morelli M.Tim ing and frequency synchronization for the uplink of an OFDMA system[J].IEEE Transactions on Communications,2004,52(2):296-306.

[3] Kivanc D,Li Guoqing,Liu Hui.Computationally efficient bandwidth allocation and power control for OFDMA[J]. IEEE Transactions on Wireless Communications,2003,2 (6):1150-1158.

[4] Cao Zhongren,TureliU,Yao Yu-dong.Analysis of two receiver schemes for interleaved OFDMA uplink[C]∥Signals,Systems and Computers,Conference Record of the Thirty-Sixth Asilomar Conference.Asilomar:IEEE,2002: 1818-1821.

[5] Huang D,Letaief K B.An interference cancellation scheme for carrier frequency offsets correction in OFDMA systems[J].IEEE Transactions on Communications, 2005,53(7):1155-1165.

[6] Lee Jong-ho,Kim Seong-cheol.Time and frequency synchronization for OFDMA uplink system using the SAGE algorithm[J].IEEE Transactions on W ireless Communications,2007,6(4):1176-1181.

[7] Pun O,Morelli M,Kuo J.Maximum-likelihood synchronization and channel estimation for OFDMA uplink transmissions[J].IEEE Transactions on Communications, 2006,54(4):726-736.

[8] Cao Zhongren,Tureli U,Yao Yu-dong,et al.Frequency synchronization for generalized OFDMA uplink[C]∥Global Telecommunications Conference,2004.Dallas: IEEE,2004:1071-1075.

[9] Trajkovic V D,Rapajic P B,Kennedy R A.Adaptive ordering for imperfect successive decision feedback multiuser detection[J].IEEE Transactions on Communications, 2008,56(2):173-176.

[10] Tevfik Y,Huseyin A.Carrier frequency offset compensation with successive cancellation in uplink OFDMA systems[J].IEEE Transactions on Wireless Communications,2007,6(10):3546-3551.

[11] Hou Shengwei,Ko CC.Intercarrier interference supp ression for OFDMA up link in time and frequency selective Rayleigh fading channels[C]∥VTC Spring 2008.Hannover:IEEE,2008:1438-1442.