您的当前位置:首页正文

基于后验概率估计的软输出MMSEMIMO排序串行干扰抵消检测算法

来源:画鸵萌宠网
第15卷第2期 2010年4月 文章编号:1007—0249(2010)02.0059.07 电路与系统学报 JOURNAL OF CIRCUITS AND SYSTEMS VO1.15 NO.2 April, 2010 基于后验概率估计的软输出M MS E M I M O 排序串行干扰抵消检测算法 王军舟, 王振伟, 李少谦 (电子科技大学通信抗干扰技术国家级重点实验室,四川I成都611731) 摘要t该文研究了信道编码无线MIMO通信系统中的软输出MMSE排序串行干扰抵消检测算法。利用MMSE滤 波器输出的高斯近似表示,提出了基于后验概率估计的判决误差传播抑制算法,以提高检测算法输出的编码比特对数 似然比的可靠性;利用检测器的矩阵结构和矩阵求逆引理,进一步提出了降低复杂度的实现方法。仿真结果表明:提 出的检测算法优于现有算法,当BER为10一,在所给的仿真条件下有1~3dB的性能增益。 关键词t MIMO MMSE;排序串行干扰抵消;复杂度 中图分类号,TN92 文献标识码t A 1 引言 对于采用信道编码的无线MIMO通信系统,MIMO检测器应当输出编码比特的对数似然比 (1og—likelihood ratio,LLR)以获得理想的错误比特率(bit error rate,BER)性能¨J。众所周知,最小 BER意义下的最优检测器是最大似然(maximum likelihood,ML)检测器,但ML检测器的复杂度随着 发射天线数目和调制星座阶数成指数增长,当天线数目和调制星座阶数较大时,难以实际应用。另一 方面,线性最小均方误差(minimum mean square error,MMSE)检测器虽然复杂度较低,但由于不能 充分利用MIMO系统提供的分集增益【2 J,其性能通常难以满足实际应用的要求。 为了在可以实现的复杂度 获得满意的性能,本文将基于MMSE的排序顺序干扰抵消(ordering sequential interference cancellation,OSIC)检测算法应用到信道编码的MIMO无线通信系统。但传统的 OSIC检测算法由于假设 想下扰抵消  ̄4J,使得检测器输出的编码比特的LLR不可靠。为此,需要抑 制不理想干扰抵消导致的饯 十扰。文献[5,6]提出采用星座分解和误符号率估计的方法估计非理想干 扰抵消下残留干扰的统计特性,但由于只考虑了与判决符号相邻的调制星座点,并且为与判决符号具 有相同欧氏距离的星座点分配相同的符号概率,导致其性能受限。为此,本文利用MMSE滤波器输出 的高斯近似,通过估计 检测符号的后验概率估计干扰抵消残留误差的统计特性,从而有效避免现有 算法存在的上述问题。为了降低检测器的复杂度,本文进一步利用检测器的矩阵结构,采用矩阵求逆 引理,提出了降低复杂度的实现方法。仿真结果表明:本文提出的算法优于现有算法,当BER为10~, 在所给的仿真条件下有1-3dB的性能增益。 本文安排如下:第2部分介绍本文研究的MIMO系统的模型和基于硬判决的MIMO MMSE OSIC 检测算法;第3部分推导基于后验概率估计的软输出MIMO MMSE OSIC检测算法;第4部分分析比 较各种算法的复杂度;第5部分给出仿真结果并进行讨论;第6部分是全文的总结。 2 系统模型与硬判决MlMO MMSE OSlC检测算法 2.1系统模型 对于具有 根发射天线和Ⅳ 根接收天线的无线MIMO通信系统,其接收信号可以用下式描述 ・收稿日期・ 2007—12・11 修订日期:2008—03-21 基金项且t “863”课题(2009AA011801,2009AA012002);“973”课题(2009CB32O4O5);国家重大专项课题(2008ZX03005.001 2009ZX03007.004,2009ZX03005.002,2009ZX03005—004);国家基础科研项目(A1420080150) 通信作者・ Email:junwang@ueste.edu.cn 电路与系统学报 』V 第l5卷 Y=Hs+n=∑hi_一 Isf+, l i=1 (1) 其中,Y=[ ,…, ] 是 ×1的接收符号矢量;,l=[ 一, ] 是独立零均值加性高斯白噪声(AWGN) 矢量,满足E[nn ]= ;H=[ 一,h蜥]是 × 的MIMO信道系数矩阵,其元素 ,表示从第J根 ,发射天线到第f根接收天线之间的平坦Rayleigh衰落信道系数。本文假设空间独立的MIMO信道模型, 并且仅有接收机知道精确的信道信息,即理想信道估计。 =[ 一,SN 是M根天线同时发射的复IT TMIMO符号矢量,其中,S ∈A(i=1,…,N )是由编码比特矢量bi从 阶复调制星座 映射得到的调 制符号(例如,M.PSK或M—QAM)。本文假设所有发射天线采用相同的调制方式,并且所有发射天 线采用等功率分配,即E{ssH}=E sIN 。2.2硬判决MIMO MMSE OSCI检测算法 MIMO MMSE OsIc算法的基本思想是【4】:基于最大检测后信干噪比(signa1.to—interference and noise ratio,SINR)排序准则,得到第i次检测的数据符号索引ki,ki∈{1,2,…,Nr};在进行第i次检测时, 首先抵消已检测出的i一1个符号的硬判决估计,在假设理想干扰抵消的前提下,利用MMSE滤波器抑 制余下的Nr—i个未检测符号和噪声的干扰。基本算法步骤如下[ 1: ● 初始化: f÷_1 Y HH +鲁lN kl argmax(SlNR ̄。 、 ● 递归: = 畸Yf ; =.Q(rk,) Y 1=Yi一 日f(H,H +譬luR] ki+1 argmax(SINR!Ⅲ ) f÷-f+1 在上面的算法描述中,符号wi表示第f次检测时的MMSE滤波器, 表示 的第 行;符号Q( .) 表示对检测得到的符号rk,作硬判决;SINR ̄。表示第f次检测时对应第 路发送数据符号的SINR,其定 义为: SINR 墨 :; (2) W t(Wi'J,l +E ∑WjhP 3 软输出MMSE MIMO 0SIC检测算法 上一节描述的硬判决MIMO MMSE OSIC检测算法(为了简单起见,后文中简述为传统OSIC检 测算法)总是假设理想判决和理想干扰抵消,因此,在有信道编码的MIMO系统中将会导致输出的编 码比特的LLR不可靠,导致系统性能降低。为此,本节在考虑不理想干扰抵消的影响的基础上,推导 软输出MMSE MIMO OSIC检测算法。该算法核心思想是:在计算编码比特的LLR值时,通过计算残 留干扰抵消误差协方差矩阵,利用MMSE滤波器抑制判决误差传播的影响;并利用残留干扰抵消误差 协方差矩阵的结构特点,通过矩阵求逆引理降低计算复杂度。 第2期 王军等:基于后验概率估计的软输出MMSE MIMO排序串行干扰抵消检测算法 61 3.1 MMSE滤波器与LLR计算 采用与传统OSIC检测算法相似的思想,令ki∈{1,2,…, }是采用最大检测后SINR排序得到的第i 次检测的发射符号索引。于是,对应的干扰抵消后的接收信号矢量可以表示为 kNT kil Yf=hk +∑^, ,+∑ ,( ,一; )+,l (3) {: ,(: 1U ID 上式中,,,,是来自没有检测的数据符号的干扰; 。是由于已经检测的发射符号的判决误差造成的不理 想干扰抵消导致的干扰; 表示 的估计。根据(3)式,用于抑制Iv、,。和,l的MMSE滤波器为 ,, . 一2 、一0 =^ I H 竹 嚣 + 1 R ir。+ J J 4) 其中,H =[ ,…,h 蜥]是去掉H的kI尼 ,…,E kH列后得到的子信道矩阵;协方差矩阵 E 为 =研 J ]=Hkl:kiQj : 。 (5) 1 这里,Q 是(j一1)x(i—1)的残留干扰抵消误差协方差矩阵,其元素 瓦 由下式计算 . ^ ==后=  ]●●●●●J ● 一 lI 一 于是,对应 的近似高斯信道表示为 ,k =w l yi ks k+ k7) iii 式中, = h , CN(O,0. 2=( 一 ) )。这样,利用对数和近似(1og—sum approximation) , 编码比特 的LLR可以根据(7)式采用下式近似计算 '6】: 1 ( min l ̄ 一 卜 8 上式中,Ao和 分别表示第 个比特为0和1的调制星座符号的集合。 3.2残留干扰抵消误差协方差矩阵计算 与传统OSIC检测算法相比,软输出MIMO MMSE OSIC检测算法需要计算(5)式中的残留干扰 抵消误差协方差矩阵Q 。为此,文献【5,6]利用(7)式的高斯近似,采用AWGN下的错误符号概率公 式估计S 的错误判决概率,并利用该错误判决概率计算g:。该方法的一个重要不足之处是:认为调 , 制星座中与判决星座点间的欧式距离相同的所有星座点具有相同的发送概率,显然与实际情况不符。 事实上,这一概率应当取决于这些星座点与判决前的接收符号瓦之间的欧式距离,而不是将瓦.判决 后的距离。为此,本文利用后验概率估计的方法来计算在接收到符号瓦的情况下,所有星座点的发射 概率,并以此来计算g 基于(7)式中的瓦,发射符号 = ,∈A的后验概率为 p : I 。瓦 ): 三pL 七.  (9) 考虑每个信道编码码字足够长并且理想交织,每个调制星座符号可以认为等概率发送。于是,(9)式 可以简化为 p =Or )∞exP(-J瓦。一['lkiat l/ ) (10) 于是,根据约束条件:∑/s = 『l瓦 J:l,可以方便地采用下述方法计算得到 = fl瓦,j ,∈A: 62 电路与系统学报 第15卷 p = ,I瓦 )= i三 c t 同时,令 岛=∑ ,p(s = ,I瓦 ),即采用软干扰抵消。可以得到 E 一 l  lJ-∑l ,一 = ,I ) (12) 注意到E【( 一 )l J=0,即Q 是对角矩阵。另外,如果令(5)式中的Q 为全零矩阵,则认为干扰 3.3排序规则 根据(7)式,容易得到第ki路检测的数据符号的检测后SINR为 SNIR : (13) “ 1一/-lk 贝0, m ax{SINR )==m ax{ z }==max( ̄ h } (14) ..因此,检测时选择矩阵日¨H蜥[Hki:kNTHkn: +击 ,D+善J j Hki:kNT最大对角元素对应的 作为待检测 符号的索引。于是,相比(2)式,采用(14)式确定检测符号顺序具有较低的计算复杂度。 3.4降1氐复杂度的实现 对于每个MIMO矢量符号 ,为了求出所有 个发射符号对应的编码比特的LLR值,前面描述 的算法需要求Ⅳ 个Ⅳ ×Ⅳ 的复矩阵的逆。对于复矩阵求逆,数值计算比较稳定的方法是采用奇异值 分解(singular value decomposition,SVD)的方法[ ,其复杂度为:13N3+ 次复数乘/除法和13 一 次复数加/减法【引,复杂度很高。为此,下面进一步利用(4)式中的矩阵结构给出一种降低实现复杂 度的实现方法。 降低矩阵求逆复杂度的一个常用的方法是利用矩阵求逆引理[ 。在3.2节中已经分析说明Q 是对 角矩阵,因此,可以将(4)式重新写成如下形式: w h kmN T 击 铂D. 1 H手 = f (15) =k m=k . 注意到:q , =g,√i+1,因此,可以得到 Af=Afl+ 岛h n—.(g ,f/E 一1),i=2,…,Nr 删H+鲁 6 利用矩阵求逆引理[9】, 可以表示为 一 /(g 等鬻 2,i一 )+^ |Il ,…, … 由于 :(AT,J_ ̄hk )H,因此仅需要计算 二 。因此,在检测每个MIMO矢量符号对应的 个发送 第2期 王军等:基于后验概率估计的软输出MMSE MIMO排序串行干扰抵消检测算法 63 符号时,仅需要一次矩阵求逆运算。这样,在计算 ,i=2,…,Ⅳ,时,仅需要3Jv;+Ⅳ +1次复数乘/除法 和2Ⅳ:+1次复数DHI减法。由于1次复数乘/除需要6 次浮点运算(floating—point operations,FLOPS),1 次复数Dril减需要2次FLOPS J,当发射和接收天线数 目相同时,每检测一个MIMO发射符号矢量,不同天 线数目下,采用上述方法相对于第3节的直接实现方 法,FLOPS减少的比例如图1所示。从图中可以看出, 采用提出的降低复杂度的实现方法,可以极大地降低 实现的复杂度。例如,当N =NR=4时,每检测一个 MIMO发射矢量符号,所需FL0PS可以降低大约 2 4 6 8 10 12 14 16 天线数目 7】%。 图l 检测每个MIMO符号矢量复杂度降低的比例 3.5 算法总结 与第2节描述的传统OSIC检测算法类似,本文所提算法可以总结为如下递归过程: ● 初始化: f 1 丑g肇遨一∞ 0 一 礤 0 0 0 0 0 J,l卜J, l,tv1=/-/I4 + 根据(14)式确定检测符号序号kl ● 递归: 根据(7)、(8)式计算映射为S 的编码比特 的LLR值; 根据(3)式计算YⅢ 根据(7)、(11)、(12)、(15)~(17)式计算 + 根据(14)式确定检测符号序号k 4 复杂度比较 这部分给出本文提出的算法、文献[4,5]的算法、传统OSIC检测算法和MMSE算法的复杂度比较。 由于在所有这些检测算法中,主要的复杂度均来自于计算MMSE滤波器需要的矩阵求逆运算,其它部 分(如计算LLR)的计算过程基本相同,且复杂度相对较小,因此,在下面的比较中,为了简单起见, 忽略这些部分的计算复杂度。 相对于文献[5,6]的算法,本文提出的算法的主要区别在于:利用估计的后验概率获得残留干扰抵 消误差相关矩阵Q ,并且采用软干扰抵消,使得Q 为对角矩阵。而文献[5,6]的算法采用硬判决抵消, 导致Q 矩阵非对角元素不为零。为了降低复杂度,文献[5,6]直接忽略非对角元素,以0.3~0.5dB的性 能损失为代价,强制Q 为对角矩阵。在这种情况下,文献[5,6]的算法也可以采用本文第4节提出的降 低复杂度的实现方法,此时,两者的计算复杂度相似。每检测一个MIMO矢量符号,两者的复杂度均 为:13N3+(3N 一2)N2+D(JV )次复数乘/除法、13N3+(2 一3)Ⅳ +D(Ⅳ )次复数DrJl减法。 根据第4节的分析,相对于传统OSIC检测算法,本文提出的算法增加的复杂度主要来自于计算 (g /Es一1),i=2,…,JV )以及(16)式中Ai引入的矩阵相加运算。而 在检测第一路信号计算 对应的MMSE滤波器时已经算出,因此,检测一个MIMO矢量符号需要增加的复杂度为:( 一1)Ⅳ; 次复数乘法、2(Nr—1)Ⅳ 次复数加法。另一方面,由于MMSE检测算法仅需要1次复矩阵求逆,根据 电路与系统学报 第15卷 第4节的分析,检测一个MIMO矢量符号的复杂度为:13N3+ 次复数乘/除法和13N.'一 次复数an/ 减法。 根据上述分析,表l给出了上述几种算法的复杂度比较。从表中可以看出,本文提出的算法与传 统OSIC检测算法复杂度相当;相对MMSE检测算法,复杂度有所增加,但由于采用了本文第4节提 出的降低复杂度的实现方法,增加的复杂度并不显著。例如,当Nr=NR=4时,相对于简单的MMSE 检测算法,本文所提算法所需的复数乘/除次数约增加16%,所需复数)Jn/减次数约增加12%。 表l 不同检测算法检测1个MIMO矢量符号的复杂度比较 5 仿真结果与分析 仿真参数如下:MIMO系统的发射机和接收 天线数均为4,即Nr=N =4。MIMO信道为空 间独立、平坦Rayleigh衰落。仿真中采用了码率 为1/2、约束长度为7的卷积码,其8进制表示 的生成多项式为(171,133) ,每个信息序列的长度 为8000bits,对应编码后的码字长为16000bits, 每个码字采用160x100的块交织器,信道译码采 用APP译码算法;采用Gray映射的QPSK和 16一QAM调制。同时仿真了本文提出的算法、文 献[4,5]的算法、传统OSIC检测算法和简单的 ¥NR(dB) 图2不同检测算法的BER性能比较(QPSK调制) MMSE检测算法的性能。其中,如前文所述,传 统OSIC检测算法认为干扰理想抵消、忽略判决 误差带来的误差传播效应;MMSE检测算法则仅 采用一次MMSE检测,不作干扰抵消。 图2和图3分别给出了采用QPSK和 16.QAM调制时不同检测器算法的性能比较。从 图中可以看出,本文所提算法明显优于现有算法。 相对与文献[4,5】的算法,在雎尺=1O4时,采用 QPSK调制时本文所提算法有约2dB的性能增 益;采用16.QAM调制时,本文所提算法约有1dB 的增益。相对于传统OSIC检测算法和MMSE检 测算法,这一增益分别超过了3dB。本文提出的 SN R《dB) 算法与文献[4,5】的算法相比,16.QAM相对QPSK 图3不同检测算法的BER性能比较(16.QAM调制) 性能增益下降的主要原因是由于:当能量归一化后,高阶调制时不同星座点到接收符号之间的距离差 比低阶调制小,因此文献【4,5】的算法引入的误差相对减小。值得一提的是:传统OSIC检测算法由于 忽略了误传传播的影响,导致输出的编码比特的LLR值不可靠,其性能甚至比简单的MMSE检测还 差。 6 总结 本文研究了信道编码无线MIMO通信系统中的软输出MMSE OSIC检测算法。提出了基于后验符 第2期 王军等:基于后验概率估计的软输出MMSE MIMO排序串行干扰抵消检测算法 65 号概率估计的方法抑制不理想干扰抵消导致的误差传播,并进一步给出了降低复杂度的实现方法。仿 真结果表明,与现有的算法相比,在所给的仿真条件下,本文提出的算法在BER为l04时,有1-3dB 的性能增益。 参考文献: [1] Miiller—Weinfutner S.H.Coding approaches for multiple antenna transmission in fast fading and OFDM[J].1EEE Trans.Signal Process, 2002,59(1 o】:2442—2450. [2]Paulraj A,Nabar R,Gore D.Introduction to Space-Time Wireless Communications[M】.U.K.:Cambridge Univl_Press,2003. [3] Wolniansky P W,Foshchini G J,Golden G D,et a1.V—BLAST:An arhitecture for realizing very high data rates over the rich—scattering wireless channel[A].in Proc.URSI Int.Symp.Signals,Syst.,Electron.[C].Italy,1998.259-300. [4]Debbah M,Muquet B,Courville M d,et a1.A MMSE successive interference cancellation scheme for a new adjustable hybrid spread OFDM system[A].In Proc.VTC 2000・Spring[C】.Tokyo,Japan,2000.745—749. [5]Lee H,Lee B,Lee I.Iterative detection and decoding with an improved V—BLAST for MIMO-OFDM systems[J].IEEE J.Se1.Areas Commun.,2006,24(3):504-513. [6】Lee H,Lee I.New approach for error compensation in coded layered V—BLAST OFDM systems[J].IEEE Trans.Commun.,2007,55(2): 345.355. 【7】 Wang X,Poor V,Wireless Communications:Advanced Techniques for Signal Reception【M].Prentice Hall PTR,2003—03. [8] Bennesty J,Huang Y A,Chen J.A fast recursive algorithm for optimum sequential signal detection in a BLAST system[J].IEEE Trans. Signal Process.,2003,51(7):1722-1730. [9] Golub G H,Van Loan C F.Matrix Computations[M].third ed.Baltimore,MD:Johns Hopkins Univ.Press,1996. 作者简介。王军(1974.),男,副教授,目前研究方向:无线与移动通信中的信号处理技术;王振伟(1986一),男, 硕士生,目前的研究方向是:认知无线网络中的频谱感知技术;李少谦(1957一),男,教授,博士生导师,目前研究 方向:扩跳频抗干扰通信技术、无线与移动通信技术。 A posteriori probability based soft—output MMSE MIMO ordering sequential interference cancellation detection WANG Jun, WANG Zhen—wei,LI Shao—qian f National Key Laboratory of Science and Technology on C0mmunications,UESTC,Chengdu 6 1 1 73 1,China) Abstract:For channel coded wireless multiple—input multiple-output(MIMO)communication systems,soft-output minimum mean square error(MMSE)ordering sequential interference cancellation(OSlC)detection is investigated in this paper.By Gaussian approximation of the output of MMSE filter,a posteriori symbol probability estimation based scheme is proposed to suppress the impact of error propagation due to imperfect interference cancellation so that the reliability of the log—likelihood ratio of each coded bit can be improved.According to the matrix structure of the detector,complexity-reduced implementation is also proposed by utilizing the matrix inverse lemma.The results show that the proposed detector outperforms the existing detectors with gains of 1-3dB at BER=10一. Key words:Multiple-input multiple-output(MIMO);minimum mean square error(MMSE);ordering sequential interference cancellation(OSlC);complexity 

因篇幅问题不能全部显示,请点此查看更多更全内容

Top