大規(guī)模MIMO系統(tǒng)的改進MCMC檢測算法*

李　蕓，易志強

(杭州電子科技大學 電子信息學院，浙江 杭州310018)

大規(guī)模MIMO系統(tǒng)的改進MCMC檢測算法*

李蕓，易志強

(杭州電子科技大學 電子信息學院，浙江 杭州310018)

大規(guī)模MIMO（多輸入多輸出）技術通過配置大規(guī)模天線陣列提高系統(tǒng)的頻譜和能量效率，接收算法的復雜度是其實現的瓶頸。MCMC（馬爾可夫鏈蒙特卡羅）檢測方法可以較低復雜度獲得接近理論最優(yōu)的性能。提出一種改進的MCMC算法，將超松弛迭代方法應用于MCMC檢測，引入松弛因子加快馬爾可夫鏈收斂速度，降低檢測復雜度。仿真結果表明，該算法能改善系統(tǒng)的誤碼率(BER)性能，解決傳統(tǒng)MCMC算法在高信噪比條件下的“陷入”問題，同時降低運算復雜度。

大規(guī)模MIMO；MCMC檢測；超松弛迭代；Gibbs采樣

0　引言

大規(guī)模MIMO(Large-scale MIMO)又稱Massive MIMO，最早由美國貝爾實驗室的MARZETTA T L提出[1]，該技術在基站配置數十根甚至上百根天線，以獲得更大的空間自由度。文獻[1-4]研究表明，當天線數目趨于無窮大時，瑞利衰落和加性高斯白噪聲等負面影響可以忽略不計，數據傳輸速率能得到極大提高。大規(guī)模天線陣列既帶來了性能增益，也帶來了前所未有的挑戰(zhàn)，如傳輸方案設計、迅速增加的硬件復雜度和計算量等。大規(guī)模MIMO系統(tǒng)中低復雜度、有效的接收算法是該技術從理論到實現的關鍵。

近年來，統(tǒng)計學中的一些方法，如基于蒙特卡羅仿真的 MCMC檢測算法[5-8]被引入到 MIMO系統(tǒng)中，它利用已知符號的后驗條件概率來迭代地求出下一個符號的后驗條件概率，最后得到發(fā)送矢量的后驗概率。MCMC方法的性能主要取決于采樣點數量和迭代次數，與待估計變量維數無關，可以避免算法復雜度隨天線數和調制階數呈指數級增長的問題。傳統(tǒng)的MCMC算法需要經過足夠次數的采樣后，馬爾可夫鏈才能趨于平衡分布。另外，在高信噪比條件下容易出現“陷入”問題(stalling problem)[7]，即采樣“陷入”某一固定狀態(tài)，采樣狀態(tài)減少，導致后驗概率估計誤差。針對上述問題，本文采用超松弛迭代的方法構造馬爾可夫鏈，選擇合適的松弛因子加快馬氏鏈的收斂速度。仿真結果表明，該算法提高了系統(tǒng)檢測性能，降低了運算復雜度。

1　系統(tǒng)模型

本文研究的大規(guī)模MIMO系統(tǒng)由一個天線數為N的基站和K個單天線用戶構成，K≤N，考慮該系統(tǒng)的上行鏈路，基站端的接收向量表示為：

其中 x為 K個用戶的信號發(fā)送向量，x=[x1，x2，…，xK]T；H為N×K維信道矩陣；n是均值為零、協(xié)方差N0IN的加性高斯白噪聲，n=[n1，n2，…，nN]T。

最大似然(Maximum Likelihood，ML)檢測是MIMO系統(tǒng)的最優(yōu)檢測算法，通過遍歷所有可能的發(fā)送信號組合，尋求最優(yōu)的檢測值，即：

式中(χ)K表示發(fā)送向量 x的全部可能取值。隨著收發(fā)天線數及調制階數的增加，ML算法搜索空間呈指數級增長，實際難以實現。

2　大規(guī)模MIMO信號檢測算法

大規(guī)模MIMO基站天線數量可能達到幾百根以上，傳統(tǒng)MIMO的一些準最大似然方法(如基于QR分解的QRM-MLD算法和球形譯碼(SD)算法)在這里難以采用。探索大規(guī)模MIMO系統(tǒng)中高性能、低復雜度的信號檢測算法是要解決的主要問題。文獻[9]將MCMC方法引入大規(guī)模MIMO，并獲得了較好的性能。

2.1MCMC算法

MCMC檢測通過統(tǒng)計抽樣獲得發(fā)送符號矢量，用統(tǒng)計方法估計各符號的后驗概率，其性能和運算量與發(fā)送信號的維數無關，只取決于采樣點數量和迭代次數。

MIMO系統(tǒng)中，多維信號的聯(lián)合概率分布如下：

MCMC算法從條件分布 p(x|y，H)中抽取樣值 x，形成馬爾可夫鏈。MIMO系統(tǒng)中馬爾可夫鏈狀態(tài)數隨x的維數呈指數增加，為了降低采樣復雜度，采用Gibbs采樣構造馬爾可夫鏈。Gibbs采樣[10]是一種基于條件分布的迭代采樣方法，它利用已知符號的后驗條件概率迭代地求出下一符號的后驗條件概率，最后得到發(fā)送信號矢量的后驗概率。第t+1次迭代中第k個符號的Gibbs采樣過程如下：

初始采樣序列 x(0)隨機產生。由于：

其中Ck是歸一化參數，χ是全部發(fā)送信號集合。

2.2MCMC改進算法

傳統(tǒng)MCMC算法需要經過足夠多次迭代才能收斂至平衡分布，提高馬爾可夫鏈收斂速度能夠降低檢測算法的計算復雜度[11]，這對大規(guī)模 MIMO系統(tǒng)尤為重要。因此考慮用超松弛迭代方法(Successive Over-Relaxation，SOR)[12]來提高馬爾可夫鏈收斂速度。

將MIMO迭代檢測器看作一個隨機線性系統(tǒng)，考慮基本線性系統(tǒng)模型如下式：

其中，M=HHH，b=HHy。當K取值很大時，M矩陣求逆運算復雜度很高，考慮采用迭代方法來求解[13]，由于 M 是對稱矩陣，有：

其中，D是對角矩陣，L是嚴格下三角矩陣。由Jacobi迭代得到第 t+1次迭代的解 x(t+1)：

超松弛迭代法(Successive Over-Relaxation，SOR)引入了松弛因子ω，以加快迭代矩陣的收斂速度，令：

當ω∈(1，2)時為超松弛。SOR迭代可以寫作：

其中，w=-HHn，是均值為零、協(xié)方差N0HHH的加性高斯白噪聲。因此可以由以下分布得到采樣值：

其中，Ck是歸一化參數，=N0||hk||2，且：

SOR-MCMC檢測算法實現過程如下：

(1)輸入接收向量y、信道H、迭代次數T；

(2)隨機產生初始向量 x(0)；

(3)while t＜T do

(4)for k=1：K

(5){由式(16)、(15)計算由 SOR迭代獲得的采樣點

(6)由獲得的采樣值，進行對數似然比（LLR）計算，并進行軟判決。

2.3算法復雜度分析

由上面的分析可知，MCMC算法的復雜度主要集中在由迭代采樣構造馬氏鏈，可以先不考慮預條件矩陣處理，即計算和b的運算量。由式(6)可知，傳統(tǒng)MCMC算法獲得采樣值的復雜度為O(N|?|)，完成一次迭代的復雜度為O(KN|?|)；而由式(15)，SOR-MCMC中獲得采樣值的復雜度僅為 O(|?|)，完成一次迭代復雜度為O(K|?|)。當迭代次數增加時，SOR-MCMC復雜度比MCMC算法低更多，特別適用于大規(guī)模MIMO系統(tǒng)。

3　仿真結果及分析

接下來分析上述算法在不同天線方案、不同調制方式下的檢測性能。大規(guī)模MIMO系統(tǒng)上行鏈路收發(fā)天線數分別為N和K，記為K×N，令收發(fā)天線比，要求K≤N。為了便于比較算法性能，采用單用戶匹配濾波器檢測（MFD）算法近似最佳的 MLD算法，即只考慮單個用戶發(fā)送數據，利用MF方法恢復發(fā)送信號。

首先考慮SOR-MCMC算法中松弛因子ω對誤碼率(BER)性能的影響。采用4QAM調制，K=16，信噪比SNR= 10 dB時，迭代次數T取20，基站天線數分別為 8/16/32的仿真結果如圖1。可以看到，β較小時，ω對誤碼率的影響更明顯。β較大時，由于天線分集增益，檢測性能得到了提高。ω=1.4時，SOR-MCMC算法檢測性能最好，后面的仿真中都取該值；ω=1即傳統(tǒng)MCMC算法。

圖1　松弛因子ω對BER性能的影響

圖2是收發(fā)天線數相同時，不同天線配置下幾種算法的性能比較，采用 4QAM調制，ω=1.4，T=20。由圖可知，隨著天線數的增加，MCMC和SOR-MCMC算法的性能都有提高，體現了大規(guī)模MIMO的特性。高信噪比時，傳統(tǒng)MCMC算法由于陷入問題影響了檢測性能，而SOR-MCMC算法克服了這個問題。在16×16/32×32/64× 64等幾種天線配置下，SOR-MCMC的性能都優(yōu)于傳統(tǒng)MCMC方法，且隨著信噪比的增加，不斷逼近單用戶的MFD算法。

圖2　K=N時不同天線方案下算法性能比較

大規(guī)模MIMO的實際應用中，基站天線數一般遠大于用戶數，下面考慮收發(fā)天線數不相等時的算法性能?；咎炀€數N=32、用戶數 K分別為 8/16/32時各算法性能如圖3，仍采用4QAM調制，ω=1.4，T=20。在圖3所示的幾種天線配置下，SOR-MCMC的性能都優(yōu)于傳統(tǒng)MCMC方法，還解決了傳統(tǒng)MCMC方法高信噪比時的陷入問題。β較大時，SOR-MCMC算法性能隨著SNR的增大不斷逼近MFD。由圖3可知，在天線配置為8×32時，SOR-MCMC算法誤碼率達到10-3所需的信噪比與MFD相比只相差0.4 dB。

圖3　K＜N時不同天線方案下算法性能比較

圖4　16QAM時SOR-MCMC算法性能

采用16QAM調制，其余參數不變，SOR-MCMC算法性能仿真結果如圖4?？梢钥吹剑S著信噪比的增加，SOR-MCMC算法檢測性能逼近MFD，即該算法在高階調制下仍然有效。

4　結論

本文主要研究大規(guī)模MIMO中低復雜度、有效的信號檢測技術。MCMC算法因復雜度有限獲得較大關注，但傳統(tǒng)MCMC算法需要經過多次采樣才能趨于平衡分布，且在高信噪比時性能不佳。本文提出了一種改進的MCMC算法，通過超松弛迭代獲得Gibbs采樣，構造平衡分布為目標分布的馬氏鏈，通過選擇合適的松弛因子加快馬氏鏈的收斂速度。仿真結果證明，該算法在不同天線方案、不同調制階數下，均能獲得優(yōu)于傳統(tǒng)MCMC算法的性能，同時計算量更低，非常適用于大規(guī)模MIMO系統(tǒng)。

[1]MARZETTA T L.Noncooperative cellular wireless with unlimited numbers of base station antennas[J].IEEE Trans. Wireless Commun.，2010，9(11)：3590-3600.

[2]HOYDIS J，BRINK S T，DEBBAH M.Massive MIMO in the UL/DL of cellular networks:How many antennas do we need?[J].IEEE Journal on Selected Areas in Communications，2013，21(2)：160-171.

[3]RUSEK F，PERSSON D，LAU B K，et al.Scaling up MIMO：opportunities and challenges with very large arrays[J].Signal Processing Magazine，2013，30(1)：40-60.

[4]LARSSON E G，TUFVESSON F，EDFORS O，et al.Massive MIMO for next generation wireless systems[J].IEEE Commun.Magazine，2014，52(2)：186-195.

[5]CHEN R，LIU J S，WANG X D.Convergence analyses and comparisons of Markov Chain Monte Carlo algorithms in digital communications[J].IEEE Trans.on Signal Processing， 2002，50(2)：255-270.

[6]DOUCET A，WANG X D.Monte carlo methods for signal processing[J].IEEE Signal Processing Magazine，2005，22 (6)：152-170.

[7]BEHROUZ F B，ZHU H D，SHI Z N.Markov chain monte carlo algorithms for CDMA and MIMO communication systems[J].IEEE Trans.On Signal Processing，2006，54(5)：1896-1909.

[8]STEPHEN A L，BEHROUZ F B.Implementation of a markov chain monte carlo based Multiuser/MIMO detector[J].IEEE Trans.on Circuits and Systems，2009，56(1)：246-255.

[9]KUMAR A，CHANDRASEKARAN S，CHOCKALINGAM A，et al.Near-optimal large-MIMO detection using randomized MCMC and randomized search algorithms[C].IEEE ICC，2011：1–5.

[10]SPALL J C.Estimation via markov chain monte carlo[J]. IEEE Control Systems Magazine，2003，23(2)：34-45.

[11]HASSIBI B，HANSEN M，DIMAKIS A G，et al.Optimized markov chain monte carlo for signal detection in MIMO systems：an analysis of the stationary distribution and mixing time[J].IEEE Trans.Signal Processing，2014，62 (17)：4436-4450.

[12]AXELSSON O.Iterative solution methods[M].Cambridge University Press，1994.

[13]GRANT A，SCHLEGEL C.Convergence of linear interference cancellation multiuser receivers[J].IEEE Trans. Commun.，2001，49(10)：1824-1834.

An improved MCMC detection algorithm for massive MIMO systems

Li Yun，Yi Zhiqiang
(Department of Electronic Information，Hangzhou Dianzi University，Hangzhou 310018，China)

Massive MIMO can achieve higher spectral and energy efficiency by large antenna array.The complexity of receive algorithm becomes the bottleneck of implementation.The MCMC detection algorithm can obtain the near optimal performance with low complexity.This paper presents an improved MCMC algorithm which introduces the over-relaxation iterative method to the MCMC detector.A relaxation factor is used to improve the convergence rate and reduce the complexity.Simulation results show that the proposed algorithm ameliorates the BER performance and the computational complexity simultaneously.The stalling problem in high signal-to-noise ratio of conventional MCMC detector is also resolved.

massive MIMO；MCMC detector；over-relexation iteration；Gibbs sampler

TN92

A

10.16157/j.issn.0258-7998.2016.05.028

浙江省自然科學基金資助項目(LQ15F010007)

2016-01-07)

李蕓(1977-)，通信作者，女，博士，講師，主要研究方向：數字系統(tǒng)設計、無線通信，E-mail:liyunr@163.com。

中文引用格式：李蕓，易志強.大規(guī)模 MIMO系統(tǒng)的改進MCMC檢測算法[J].電子技術應用，2016，42(5)：101-103，108.

英文引用格式：Li Yun，Yi Zhiqiang.An improved MCMC detection algorithm for massive MIMO systems[J].Application of Electronic Technique，2016，42(5)：101-103，108.