華 磊,王亞軍,劉 爽

(江蘇科技大學(xué) 海洋學(xué)院,江蘇 鎮(zhèn)江 212100)

0 引 言

大規(guī)模正交頻分復(fù)用(Orthogonal Frequency Division Multiplexing,OFDM)與多輸入多輸出(Multiple-Input Multiple-Output,MIMO)技術(shù)相結(jié)合可以發(fā)揮兩者的優(yōu)勢(shì),實(shí)現(xiàn)更快的通信速率和更高的頻譜利用率,MIMO-OFDM已經(jīng)成為下一代無(wú)線通信中的核心技術(shù)[1-2]。

但OFDM作為一種多載波技術(shù),當(dāng)各載波頻譜相同或相近時(shí)容易導(dǎo)致信號(hào)的包絡(luò)起伏較大,進(jìn)而造成信號(hào)的峰值平均功率比(Peak-to-Average Power Ratio,PAPR)較高。高峰均比信號(hào)通常需要帶有較大反饋的功率放大器,這就增加了模數(shù)轉(zhuǎn)換器和數(shù)模轉(zhuǎn)換器的復(fù)雜性,同時(shí)對(duì)發(fā)射機(jī)內(nèi)放大器的線性也提出了很高的要求,由此可能導(dǎo)致信號(hào)畸變,使信號(hào)的頻譜發(fā)生變化,導(dǎo)致系統(tǒng)的性能嚴(yán)重惡化[3-4]。而將多用戶大規(guī)模MIMO(Multi-user Massive MIMO,MUM-MIMO)技術(shù)應(yīng)用到工程實(shí)際時(shí),必須要求基站使用低功耗的射頻鏈路和功率高效的射頻功率放大器[5],這就需要發(fā)射信號(hào)具有較低的PAPR。為滿足使用高效的功率放大器和降低多用戶間的干擾這兩個(gè)要求,迫切需要使用更加有效的方法來(lái)降低系統(tǒng)的PAPR。

目前降低OFDM系統(tǒng)PAPR的方法有很多[6-9],盡管可以應(yīng)用于點(diǎn)對(duì)點(diǎn)的MIMO-OFDM系統(tǒng),但由于大規(guī)模MIMO系統(tǒng)中接收端用戶協(xié)同工作極為困難,因此運(yùn)用到大規(guī)模MIMO-OFDM系統(tǒng)中顯得不甚合理。在文獻(xiàn)[10]中,作者給出了基于快速遞推截?cái)嗨惴?Fast Iterative Truncation Algorithm,FITRA)的低PAPR非線性預(yù)編碼設(shè)計(jì)。其設(shè)計(jì)方法聯(lián)合考慮了MUM-MIMO預(yù)編碼、OFDM調(diào)制和PAPR減小,核心思想是通過(guò)利用無(wú)窮范數(shù)代替PAPR,從而得到滿足無(wú)窮范數(shù)約束下的最小化多用戶間干擾(Multi-user Interference,MUI)問(wèn)題。文獻(xiàn)[11]中,作者提出了一種添加干擾信號(hào)降低系統(tǒng)PAPR的方法,并利用交替方向乘子算法(Alternative Direction Method of Multipliers,ADMM)對(duì)其進(jìn)行求解。在該方法中,作者在預(yù)編碼后的頻域信號(hào)中加入一個(gè)干擾信號(hào),同時(shí)對(duì)該干擾信號(hào)做出限制,保證其位于每個(gè)子載波所對(duì)應(yīng)的信道響應(yīng)矩陣的零空間中,因此該方法在獲得好的PAPR性能同時(shí),不會(huì)為整個(gè)系統(tǒng)帶來(lái)額外的信號(hào)失真和帶外輻射。文獻(xiàn)[12]中,作者通過(guò)最速梯度下降方法來(lái)解決多用戶預(yù)編碼和降低PAPR聯(lián)合的問(wèn)題,利用MIMO信道矩陣的高維零空間來(lái)降低信號(hào)的PAPR,但其方法得到的誤符號(hào)率(Symbol Error Rate,SER)性能較差。上述方法都能夠獲得好的PAPR抑制性能,但均存在著算法計(jì)算復(fù)雜度高的缺陷。

針對(duì)以上問(wèn)題,本文利用加速近端梯度算法(Accelerated Proximal Gradient Method,APGM)來(lái)設(shè)計(jì)降低多用戶大規(guī)模MIMO-OFDM系統(tǒng)中發(fā)射信號(hào)PAPR的低復(fù)雜度快速收斂算法,該算法在抑制發(fā)射信號(hào)的PAPR同時(shí)也能獲得好的SER性能。

1 多用戶大規(guī)模MIMO-OFDM系統(tǒng)模型

在如圖1所示的多用戶大規(guī)模MIMO-OFDM系統(tǒng)的下行鏈路中,假定基站配置有M根發(fā)射天線,同時(shí)服務(wù)于K個(gè)單天線移動(dòng)用戶,且滿足K?M。在本文中,假定信道狀態(tài)信息在發(fā)射端完美已知。OFDM載波數(shù)目為N,將載波分為數(shù)據(jù)載波Γ和自由載波Γc兩部分,其中數(shù)據(jù)載波用來(lái)傳輸數(shù)據(jù),而自由載波不傳輸任何數(shù)據(jù),放于頻帶的兩端用作保護(hù)間隔。向量sn∈CK×1為在第n個(gè)子載波上傳輸?shù)陌琄個(gè)終端用戶的數(shù)據(jù),對(duì)于數(shù)據(jù)載波n∈Γ時(shí),sn選自某個(gè)調(diào)制星座字符集;對(duì)于自由載波n∈Γc時(shí),有sn=0K×1。為消除接收端多用戶間的干擾,需對(duì)用戶信息符號(hào)進(jìn)行線性預(yù)編碼:wn=Gnsn,其中,wn∈CM×1為預(yù)編碼后的信號(hào),Gn∈CM×K為相對(duì)于第n個(gè)子載波的預(yù)編碼矩陣。常用的預(yù)編碼方案有匹配濾波(Matched Filtering,MF)、迫零(Zero-forcing,ZF)和最小均方誤差預(yù)編碼等。對(duì)于發(fā)射機(jī)已知信道矩陣Hn∈CK×M,則迫零預(yù)編碼矩陣可以寫(xiě)為,接著將歸一化的頻域信號(hào)wn,?n重新排序后分配到各個(gè)發(fā)射天線上[x1,x2,…,xM]=T[w1,w2,…,wN],T表示圖1中重新排序的映射。其中,向量xm∈CN×1表示第m根天線上的頻域信號(hào)。然后對(duì)頻域信號(hào)xm,?m進(jìn)行離散傅里葉逆變換(Inverse Discrete Fourier Transform,IDFT)得到相應(yīng)的時(shí)域信號(hào)ym,m∈1,2,…,M。完成并行到串行轉(zhuǎn)換之后,通過(guò)對(duì)每根天線上的時(shí)域信號(hào)ym,?m添加循環(huán)前綴來(lái)避免碼間干擾。最后將這些時(shí)域樣本轉(zhuǎn)換成模擬信號(hào),通過(guò)射頻鏈路傳輸?shù)綗o(wú)線信道中。

圖1 大規(guī)模MIMO-OFDM系統(tǒng)框圖

2 峰值平均功率比

OFDM信號(hào)由多個(gè)獨(dú)立的調(diào)制信號(hào)疊加而成,能夠有效地對(duì)抗頻率選擇性衰落,對(duì)于第m根天線上的頻域信號(hào)xm∈CN×1,相應(yīng)的時(shí)域連續(xù)信號(hào)可以寫(xiě)為

式中:T為OFDM符號(hào)時(shí)間;Δf為子載波頻率間隔;N為子載波數(shù)。由于各個(gè)子載波的相位相互獨(dú)立,相互疊加后的OFDM信號(hào)就會(huì)有非常大的動(dòng)態(tài)范圍。通常使用PAPR來(lái)衡量發(fā)射信號(hào)的波動(dòng)。PAPR表示信號(hào)功率的峰值與均值之比,即

在實(shí)際系統(tǒng)中,通常只有發(fā)射信號(hào)ym(t)的一些離散時(shí)間樣本,為了更好地近似連續(xù)時(shí)間信號(hào)的PAPR,在獲得離散時(shí)間樣本時(shí)采用J=4倍過(guò)采樣,即

過(guò)采樣操作也可以表示為ym=FHxm,m=1,2,…,M,其中,FH∈CJN×N表示JN點(diǎn)傅里葉變換矩陣的前N列,J為過(guò)采樣因子。因此具有J倍過(guò)采樣的時(shí)域樣本的PAPR由下式給出:

使用互補(bǔ)累積分布函數(shù)(Complementary Cumulative Distribution Function,CCDF)來(lái)衡量信號(hào)的峰均比,定義為發(fā)射信號(hào)的峰均比超過(guò)門(mén)限值的概率,具體表示為CCDF(γ)=Pr{PAPR＞γ},其中γ表示指定的門(mén)限值。

3 算法設(shè)計(jì)

3.1 加速近端梯度算法

加速近端梯度算法用于求解以下形式的問(wèn)題:

其中:f:Rn→R是連續(xù)可微的凸函數(shù);g:Rn→R∪{+∞}為連續(xù)的凸函數(shù)。算法基本迭代過(guò)程如下:

式中:ωk∈[0,1)是外推參數(shù);γk是步長(zhǎng)因子,通常選取ωk=k/(k+3)。函 數(shù)g(y)的近端算子proxγg(·)定義式如下:

當(dāng)函數(shù)f對(duì)應(yīng)的最小Lipschitz常數(shù)為L(zhǎng)時(shí),使用固定步長(zhǎng)γk=γ∈(0,1/L],該算法以ο(1/k2)的速率收斂于目標(biāo)值。

3.2 低復(fù)雜度APGM算法設(shè)計(jì)

由于發(fā)射天線數(shù)量遠(yuǎn)多于終端用戶數(shù)量,即M?K,將會(huì)有無(wú)限多的預(yù)編碼信號(hào)w?[wH1,wH2,…,wHN]H滿足預(yù)編碼約束,只需找到一組頻域信號(hào)wn,n∈1,2,…,N,使其對(duì)應(yīng)的時(shí)域信號(hào)具有較小的PAPR,同時(shí)能夠消除多用戶間干擾:

式中:∈C(K×|Γ|+M×|Γc|)×1是所有用戶數(shù)據(jù)符號(hào)和自由載波上的0M×1的集合;為塊對(duì)角矩陣,其主對(duì)角塊為Hn,?n∈Γ和IM,?n∈Γc。為了表述方便,將上述將信號(hào)重排序到天線上的過(guò)程寫(xiě)為更緊湊的形式,其中。根據(jù)系統(tǒng)模型,預(yù)編碼后的頻域信號(hào)與原始數(shù)據(jù)及時(shí)域發(fā)射信號(hào)之間的關(guān)系可以用下式表示:

為了簡(jiǎn)化表示,寫(xiě)成更緊湊的形式:

將式(13)轉(zhuǎn)化為可以利用APGM算法求解的形式。將式(13)的放松模型用拉格朗日形式表述為

式中:α是正則化參數(shù)。使用APGM來(lái)求解問(wèn)題(14),通常需要計(jì)算函數(shù)f(y)=‖s-Qy‖22梯度?f(y)的Lipschitz常數(shù),再根據(jù)Lipschitz常數(shù)來(lái)估計(jì)步長(zhǎng)γk。梯度?f(y)的最小Lipschitz常數(shù)對(duì)應(yīng)于L=2σmax(QQH),σmax表示為矩陣的最大奇異值。但是計(jì)算函數(shù)Lipschitz常數(shù)復(fù)雜度太高,不利于實(shí)際工程中的應(yīng)用。可以用一個(gè)線性搜索方法[13]來(lái)確定步長(zhǎng)γk。利用式(8),函數(shù)g(y)=α‖y‖2∞的近端算子可以寫(xiě)成以下形式:

式中:u=xk+1-γ?f(xk+1)。近端算子(15)可以通過(guò)對(duì)文獻(xiàn)[14]中的方法改進(jìn)來(lái)精確獲得(算法1),具體的計(jì)算步驟如下:

3 對(duì)向量b按降序進(jìn)行排序,s=sort(b,‘descending’)

4 fori=1,2,…,JNM

5 令di=∑i k=1sπ(k),其中sπ(k)表示在向量s中第k大的元素;

13 end if

14 end for

15 end for

16 returny

以上就是計(jì)算?2∞近端算子的過(guò)程。根據(jù)公式(6)和(7)分別更新變量xk+1和yk+1,求解模型(14)的APGM算法(算法2)如下:

1 輸入xk,γk-1,β,其中β∈(0,1)

2 令γ=γk-1

3 fork=1,2,…,K

4z=proxγg(xk-γ?f(xk))

5 iff(z)≤^fγ(z,xk),where^fγ(z,w)=f(w)+

6 break

7 end if

8 更新γ=βγ

9 end for

10輸出yk+1=z

3.3 算法復(fù)雜度分析

APGM求解問(wèn)題(14)由算法1和算法2組成。在算法1中,計(jì)算向量b的值需要O(MJN)的復(fù)雜度,對(duì)向量b進(jìn)行排序需O(1/2·MJNlb(JN))的計(jì)算量,計(jì)算j和μ的復(fù)雜度為O(MJN)?！瑉-y‖22需要計(jì)算向量間的乘法運(yùn)算,其復(fù)雜度為O(MJN)。由函數(shù)f(y)帶來(lái)的梯度?f(y)=QH(Qy-s)計(jì)算復(fù)雜度為O(K|Γ|M+1/2·MJNlb(JN))。因此,求解(14)的APGM算法的計(jì)算復(fù)雜度主要來(lái)自近端算子的排序和函數(shù)梯度的計(jì)算,總復(fù)雜度為O(K|Γ|M+MJNlb(JN))。APGM算法與其他算法的算法復(fù)雜度對(duì)比如表1所示,可以看出所提出的APGM復(fù)雜度遠(yuǎn)低于其他方法,同時(shí)算法運(yùn)行時(shí)間也遠(yuǎn)小于其他算法。將各種算法運(yùn)用在大規(guī)模天線陣列中,天線數(shù)目的增加隨之帶來(lái)更高的能量消耗,而APGM具有較低的運(yùn)算復(fù)雜度及更快的運(yùn)行時(shí)間,能夠很好地在實(shí)際系統(tǒng)中應(yīng)用。

表1 算法復(fù)雜度分析

4 仿真分析

將APGM算法與現(xiàn)有的迫零預(yù)編碼(Zero-Forcing(ZF)Precoding) 方法、匹配濾波預(yù)編碼(Matched Filtering(MF) Precoding)方法、PMP方法[10]、PRO-ADMM方法[11]、MU-PP-GDm方法[12]進(jìn)行比較,各種方法的最大迭代次數(shù)如表1所示。為驗(yàn)證APGM的有效性,在大規(guī)模MIMO-OFDM系統(tǒng)下行鏈路中進(jìn)行實(shí)驗(yàn)仿真。在基站處設(shè)計(jì)100根天線,同時(shí)為10個(gè)單天線用戶提供服務(wù)。在OFDM系統(tǒng)中,使用IEEE 802.11n中指定的40 MHz模型,子載波個(gè)數(shù)設(shè)為128,其中Γ=108個(gè)子載波用于數(shù)據(jù)傳輸,在仿真中忽略用作導(dǎo)頻的子載波;Γc=20個(gè)子載波用于頻帶保護(hù),使用互補(bǔ)累積分布函數(shù)來(lái)評(píng)估PAPR降低性能。采用具有格雷映射的16QAM調(diào)制方式,信道為頻率選擇性信道,在時(shí)域具有多抽頭結(jié)構(gòu),且抽頭數(shù)設(shè)為D=8。信道時(shí)域響應(yīng)矩陣{Hd,d=1,2,…,D}的元素服從均值為0,方差為1的復(fù)高斯分布,對(duì)應(yīng)的頻域響應(yīng)矩陣為

圖2給出了APGM與ZF算法、MF算法、PMP方法、PRO-ADMM算法、MU-PP-GDm方法對(duì)PAPR抑制的性能比較。在CCDF=Pr(PAPR＞γ)=10-3時(shí),APGM算法只需30次迭代就可將PAPR降低到2 dB,比ZF方法降低了10 dB以上,比PMP方法的PAPR值小1 dB,比MU-PP-GDm算法的PAPR值小3.2 dB,還比PRO-ADMM算法的PAPR值少3.9 dB;PMP算法可以將PAPR降到3 dB,MU-PP-GDm算法可以將發(fā)射信號(hào)的PAPR降低到5.2 dB,PROADMM算法將信號(hào)的PAPR降低到5.9 dB。APGM只需30次的迭代就可以達(dá)到這樣效果,而PMP方法需進(jìn)行約2 000次的迭代才能將PAPR值降到2.7 dB。從圖中可以看到,APGM的PAPR減少性能優(yōu)于PMP、PRO-ADMM、MU-PP-GDm等算法。

圖2 不同優(yōu)化算法的PAPR值比較

圖3比較了不同算法在加性高斯白噪聲信道上的SER性能,其中信噪比(Sigal-to-Nise Ratio,SNR)定義為SNR=E{‖ym‖}/N0,其中N0表示接收器的噪聲方差?？梢钥吹脚cZF方案相比(在SER=10-2時(shí)),APGM的信噪比僅比理想值多1.3 dB左右,與PMP方法相比少0.7 dB,比PROADMM算法少2.5 dB,與MF方法相比少3.8 dB,同時(shí)還比MU-PP-GDm算法少6 dB以上。

圖3 不同優(yōu)化算法的SER性能分析

圖4給出了不同算法的迭代次數(shù)對(duì)應(yīng)的PAPR對(duì)比。CCDF=Pr(PAPR＞γ)=10-3時(shí)設(shè)計(jì)的APGM經(jīng)過(guò)10次迭代就可以將PAPR降低到3.8 dB,經(jīng)過(guò)20次迭代PAPR降低到2.4 dB,而經(jīng)過(guò)30次迭代PAPR值可以降到2 dB,比迫零預(yù)編碼方法PAPR值減少了10 dB以上;PRO-ADMM方法外迭代200次、內(nèi)迭代2次,總共需400次迭代才能將PAPR值降到5.9 dB;MU-PP-GDm方法經(jīng)過(guò)50次迭代可以將PAPR值降到5.2 dB;而PMP方法在100次迭代后PAPR值僅降到10 dB,PMP方法2 000次迭代后PAPR值降到2.7 dB。

圖4 不同迭代次數(shù)的PAPR值比較

在SER=10-2時(shí),從圖5中可以看到APGM在經(jīng)過(guò)20次迭代之后信噪比值為0.2 dB,與ZF方案相比損失1 dB,30次迭代之后得到的信噪比為0.5 dB。而PMP方法在1 500次迭代后信噪比為1.2 dB,與ZF方法相比損失2 dB,PMP方法經(jīng)過(guò)2 000次迭代之后得到的SNR為1.3 dB。PROADMM方法內(nèi)外循環(huán)迭代共400次之后獲得的信噪比為3 dB,MU-PP-GDm算法經(jīng)過(guò)50次迭代得到的信噪為7 dB。從以上仿真實(shí)驗(yàn)可以看出APGM具有更低的復(fù)雜度和更快的收斂速度,只需少量的迭代就可以獲得更低的SER和更好的PAPR抑制效果,因此APGM比存在的其他方法更高效,更利于工程實(shí)際中的應(yīng)用。

圖5 不同迭代次數(shù)的SER值比較

5 結(jié)束語(yǔ)

本文針對(duì)大規(guī)模MIMO-OFDM系統(tǒng)中較高的PAPR值問(wèn)題提出了一種低復(fù)雜度的APGM,核心是利用發(fā)射天線的額外自由度將MU預(yù)編碼、OFDM調(diào)制和PAPR抑制聯(lián)合成凸優(yōu)化問(wèn)題。仿真結(jié)果表明,APGM只需很少的迭代次數(shù)就可以降低發(fā)送信號(hào)的PAPR,獲得好的SER性能,同時(shí)具有更低的計(jì)算復(fù)雜度,更適合工程實(shí)際中的應(yīng)用。在今后研究中,考慮將APGM用于非完美信道狀態(tài)信息下的系統(tǒng)模型,同時(shí)將發(fā)射信號(hào)的功率及帶外散射作為模型的約束,以進(jìn)一步增強(qiáng)系統(tǒng)的穩(wěn)健性。