2012年1O月 计算机工程与设计 C()MPUTER ENGINEERING AND DESIGN Oct.2012 第33卷第1O期 Vo1.33 No.10 基于概率和自适应的最大后验概率盲均衡算法 耿天玉,舒 勤,应大力,郑莹莹 (四Jl1大学电气信息学院,四川成都610065) 摘要:为了减小最大后验概率(maximum a posteriori,MAP)盲均衡算法中稳态剩余误差和加快收敛速度,提出了一种 改进的基于概率和自适应 的MAP盲均衡算法。根据概率论的原理,分析了MAP算法中d表示的实际意义,改进算法通 过改变更新方程中的高斯信号方差值 ,即由固定值改为可变值,再运用不等概思想进一步的提高了MAP算法的收敛性 能。仿真结果表明,与MAP盲均衡算法相比,改进型MAP盲均衡算法具有更好的均衡效果和收敛性能。 关键词:盲均衡;最大后验;自适应算法;剩余根均方误差;概率 中图法分类号:TN911.5 文献标识号:A 文章编号:1000—7024(2012)10—3720—05 MAP blind equalization algorithm based on probability and adaptive GENG Tian-yu,SHU Qin,YING Da—li,ZHENG Ying—ying (School of Electrical Engineering and Informantion,Sichuan University,Chengdu 610065,China) Abstract:To reduce the residual mean square error(MSE)and increase the convergence speed in MAP,an improved MAP blind equalization algorithm based on probability and adaptive is presented.According tO the principles of probability theory,the prac— tical significance of the 6 in the MAP algorithm is analyzed.The proposed algorithm changes the Gaussian variance d in renewal equation from fixed value to variable value,which further improves the convergence of MAP blind equalization algorithm using the idea of not equal tO the probability.The simulation results show that the equalization effect and convergence performance of modified MAP blind equalization algorithm are better than those of MAP blind equalization algorithm. Key words:blind equalization;MAP;adaptive algorithm;residual RMSE;probability O引 言 在通信系统中,通信信号在信道系统中传输,会因为 为常数模算法(constant modulus algorithm,CMA)CMA 算法具有很多优点,它的计算量比较小,原理比较清晰, 而且收敛稳定。但同样也存在一些缺点,比如收敛误差较 信道畸变而受到码间干扰(inter symbol interference,ISI) 的作用。使得接收信号存在很高的误码率,能够造成信号 的判决错误从而影响通信质量,为了确保能接收到质量较 高的通信信号,在信号的接收端需要一种有效的信号处理 方法来减小或者消除码间干扰对通信信号产生的影响 。], 一大,收敛速度较慢等_5 ]。后来S Chen将最大后验概率的 思想融入到Bassgang盲均衡算法中,得到了MAP盲均衡 算法_9]。该算法用输出信号的最大后验概率密度函数作为 代价函数的参数,然后根据最陡下降随机梯度方法进行迭 代,从而更新均衡器的权向量。对接收信号进行均衡后具 有稳态误差小,星座图恢复得紧凑等优点,但也存在一系 列问题,该算法中 的选择是靠经验值来选取的,但是根 据概率论的原理,a表示的实际意义是输出信号和相应星座 符号点之间的偏差,因而直接规定其值有失准确性;同时 MAP算法是分阶段进行的,每个阶段的 值是不相同的, 般采用均衡技术来达到这种效果。传统的自适应均衡技 术一般都需要发送训练序列,即在数据传输开始前发送, 发送完后再转到自适应方式开始数据传输,不仅占用了大 量的信道资源,降低了通信的有效通信率,而且在实际应 用中会因为有些情况而受到l4]。盲均衡技术客服了这 方面的缺点,相比之下,它不要发送训练序列,这提高了 把4个阶段的 值规定为相同定值,影响了MAP算法的收 敛性能。本文在MAP算法的基础上对MAP算法进行了改 进,即把原来更新方程中的高斯信号方差值a由固定值改 信道的利用率,同时也较小了系统宽带,因此在目前的通 信领域备受关注。目前在QAM信号中应用最普遍的算法 收稿日期:2011 11—22;修订日期:2012—02—06 作者简介:耿天玉(1986一),男,河北衡水人,硕士研究生,研究方向为现代信号处理;舒勤(1958一),男,重庆人,博士,教授,研究 方向为现代信号处理、计算机网络;应大力(1982一),男,重庆人,硕士,讲师,研究方向为通信信号处理;郑莹莹(1986一),女,吉林 白城人,硕士研究生,研究方向为现代信号处理。E—mail:gengtinayul122@163.corn 第33卷第10期 耿天玉,舒勤,应大力,等:基于概率和自适应的最大后验概率盲均衡算法 .3721. 为自适应值,并且把MAP算法中采用的等概率模式进行不 等概处理,经过明改进算法的收敛速度和稳态剩余 误差均优于MAP算法。 均衡的代价函数可表示为:-,( ( ))一E{J( ( )I ( ))),根据最陡下降的随机梯度方法,更新的均衡器抽 设l『(叫( )f ( ))一alog(f(w(n)j ( )))贝4 MAP盲 1对盲均衡系统模型的描述 图1表示的是盲均衡原理框图[1 。该系统的通信过 程可以描述为:发送端从有限字符集S(S定义为式 头权向量值的式子如下所示 ∑∑e妻妻 ( 一 ) 一 ( )) (1)[1。])先对数据进行等概率选取,构成了发送信号序列 ∑∑e— ly( ) I { (,z)},然后在加性高斯白噪声{,z( )}的影响下通过 某一未知的离散时间信道{h( )},我们把通过信道后的 信号成为盲均衡器的接收信号序列,记为{r( )}; {r ( ))在均衡器{W( ))的作用下,得到输出信号序列 {Y(72)}。 噪声n(n) 输入序列J信道l 、/ 、接收序列J均衡器(反卷积)l恢复的序 {s(n)) 1 h(n)l,\、上,/ {x(n)) l f n)} l (y(n】广 图1盲均衡原理 2对MAP算法的描述 基于MAP盲均衡算法的核心思想就是以均衡器输出信 号的后验概率密度函数为代价函数,对其求导最大化代价 函数,从而得到均衡器抽头权向量的最优值。也就是通过 计算接收的信号是由信号集里的哪个信号点通过信道后而 得到的,即用最低的差错概率将输出信号还原成发送信号。 假设信号星座图为MQAM信号,同时发送端信号序 列{S(n)}是同分布的,星座图中符号点的集合定 义为[123 S一{ 一(2i—Q一1)+ (2k—Q一1),1 i,k Q} (1) 式中:Q一 ,M—QAM信号阶数。 则均衡器的输出信号Y( )可以近似表示为 ( )一z( )+ ( ) (2) 其中的t( )一Iz( 一 ),正整数 表示的是时延 值。7J( )可认为是零均值,方差为a的高斯噪声。这样, Y( )可以用M个高斯信号来近似表示。每个高斯信号的 均值 对应等于星座中的某一符号点s ,即 ( )一 ( )十 ( ) (3) 每个高斯信号的方差即为 ,因此Y( )可以等效表 示为如下的随机过程 ( )一 + ( )以概率P (1三三三 ,z Q) (4) 即 (n)~N( ,d)以概率P (1 志,z三三三Q) (5) 其中P (1≤忌,z≤Q)都相等,等于1/M。 这样,Y( )的后验概率密度函数可以表示为 (( 硼( )J) ( ))一∑∑ 一妻妻 (6) (7) MAP均衡算法是对星座区域进行划分,然后分阶段执 行l1 ,这样可以加快收敛,减小算法的复杂度,每一阶都 只使用对应区域中的4个高斯信号,即等效为4QAM的情 况,迭代式变为 ∑∑ ( 一 ( )) + 一训 ’+ 生土上 ’ (8) MAP算法分段的过程可以简单概括如下口 : 首先是第一阶段,将星座图等效为4QAM模型,对应 一个复平面。该模型运用4个高斯信号,这4个高斯信号 分别对应复平面的4个象限,其象限内所有符号点坐标的 平均值则为4个高斯信号的均值。 第二阶段,采用16个高斯信号,这l6个高斯信号被 分成4组,分别对应复平面的一个象限,而每1组都又含 有4个高斯信号,这样就把象限分割。当第一阶段结束后, 均衡器的输出信号y( )正确的恢复到复平面的4个象限 里,接着执行第二阶段,即在对对应象限进行区域划分, 每个象限划分成4个区域,然后把输出信号划分到对应象 限的某一区域,如果均衡器的输出出现在某一特定的象限, 相应象限的4个分组将用来调整均衡器抽头系数。 同理,MQAM信号会执行H个阶段(H—log M), 其中每个阶段的调整子任务都能非常快而且可靠的完成, 因此均衡器权向量能够更快更可靠地实现收敛。当各个阶 段都执行完以后,便可以得到恢复后的星座图。 以64QAM信号为例,第一阶段,采用均值为{2十 2j},{2--2j},{一2+2j),{一2—2j}的高斯信号,通过 迭代式(8)进行等效的4QAM调整,这一阶段完成后, 信号就能被正确的调整到复平面的某一象限里。使用16个 高斯信号,对应的均值就是16QAM信号星座图中的符号 点Sk/,1Gk,z 4.它们被分成4组,每组对应一个象限。 若是输出信号 ( )属于某一象限,则该象限对应的4个 高斯信号就被用来进行迭代,第二阶段完成后,输出信号 就会被分类到象限的某一区域。第三阶段,使用64个高斯 信号再对星座象限的各个区域进行更细的分割,从而把均 衡器的输出信号恢复到更细的区域中去。于是,通过3个 ・3722・ 计算机工程与设计 2o12正 阶段的均衡调整,信号星座图就成功被恢复了。 咝 3对MAP算法的分析 通过对MAP均衡器抽头权向量式(7)的观察分析可 知,MAP盲均衡算法涉及到主要的参数 和a。迭代步长 的选择对MAP算法的收敛性能的影响很重要,它的选取 合适与否会严重影响到算法的收敛速度和稳态性能,一般 对 取经验值[”]。另一个重要的参数就是对a的选取,由 +。一 + —m=_1 _ 三6孥 ̄1 一 r (13) 当 和P 分别都取固定等值时,式(13)就退化成 式(8)所示的MAP盲均衡迭代式。因此本文提出将更新 方程式中的 由固定值改为自适应值,同时把原来算法中 于MAP是分段进行均衡的,这样对每个阶段的。有比较明 确的区分,如果。取太大,则无法达到星座的理想分割度; 若。取太小,则需要严格控制星座分割的规模l1 。一般对 的选择是靠经验值来确定的,但是根据MAP算法的原 理, 表示的其实是输出信号和星座符号点之间的偏差,因 而直接把它规定为固定值是不准确的,同时算法是分阶段 执行的,每一阶段的 都不一样,如果当系统中其它条件 发生改变时,相应的a值也需要发生相应的改变以便能得 到更好的均衡效果,此时 取固定值了算法的收 敛性能。 再者,MAP盲均衡算法还认为 ( )由区域内的每 个星座符号点通过信道产生的概率是相等的。这个结论在 输出刚开始的时候是成立的,但是当均衡器逐渐达到收敛 的时候,均衡器的输出信号Y( )会越来越靠近发送信 号,这时可以认为离Y( )越近的星座符号点,y( )由 它通过信道而产生的概率越大。当算法最终收敛时,几乎 可以认为Y( )就是由离它最近的星座符号点产生的,因 此采用等概率的模式是不准确地,即这个假设在原理上也 是不准确的。 因此根据上述分析思想,我们可知MAP盲均衡算法是 分阶段执行的,前H一1阶段所使用的高斯信号对应均值都 不是发送信号星座中的符号点,只有最后一阶段采用的高 斯信号对应的均值才是发送信号星座中的符号点。因此设 Dn={ ,m=1,2,3,4),Dn为相应执行阶段中划分 出的第n个区域,Era为该区域对应的高斯信号的均值。则 Y(n)~N( , )以概率P 。 当均衡器的输出信号Y( )∈ 时,则Y( )的后 验概率密度函数可局部近似为 ( ( )l ( ))一∑ (9) m 1 V 7盯m 为了确保快速收敛,防止算法发散,设 =mim (1O) J (议^( )I ( ))一O'olog( (议^( )l ( ))) (11) 则算法的代价函数可表示为 J (ws(n))一E{.厂 (姒( )I ( ))) (12) 然后对代价函数利用最陡下降的随机梯度方法进行更 新,如下式所示 采用的等概率模式进行不等概处理,即把 也进行自适应 的调整以提高MAP盲均衡算法的收敛性能。 4改进型MAP算法 4.1自适应仃的MAP盲均衡算法 根据MAP盲均衡算法的原理,随着均衡的进行,无论 Y( )最有可能由哪个星座符号点 得来,其偏离这个准 直的幅度都会越来越小,当达到收敛的时候,几乎可以认 为Y( )就是由离它最近的星座符号点得来的了。而根据 式(5),o的实际意义是均衡器输出Y( )偏离它对应准 值的方差值。即 一E{i ( )一e埘I。} (14) 设 ^ 一l ( )一 } (15) ^ 可以证明 是 的无偏估计。根据无偏估计的定 义有 ^ E{ 一 )一E{I ( )一£ l。一E{I ( )一 I }} 一E{l ( )一 l }一E{E{l ( )一em I。}} A —E{ }一E{l ( )一 I } ^ 一E{d }一 (16) ^ 用 来作为 的估计值,就得到了自适应a改进算 ^ 法的思想。 采用无偏估计子 进行估计,则0"0一rain A d ,Pm都相等,等于1/4。 故式(13)简化为 三m=三1 二 一 ( +1)一 ( )+ 4.2基于概率和自适应仃的MAP盲均衡 由前面的分析可知,假设MAP盲均衡算法输出信号Y ( )由区域内每个星座符号点通过信道产生的概率是相等 的,在原理上是不准确的,所以应该对 进行自适 应调节。 算法的基本思想:根据Y(n)离区域内各个星座符号 第33卷第10期 耿天玉,舒勤,应大力,等:基于概率和自适应的最大后验概率盲均衡算法 .3723. 点的距离来确定它们的概率,距离越近概率越大,当距离 为零时概率为1,也就是说概率随距离单调递减。用符号表 示如下 头权向量初始化为:中心抽头值为1,其余为0,Rum 50 (即每5O个点取一次RMSE),信道为典型的电话信道_】 ] H( )一0.005+0.009z- 一0.024z- -4- 0.854z-。+一0.218z-4+ Dis 一l (n)一£州l 则 应满足以下3个条件 4 (18) 0.049z- 一0.O16z- (25) f∑P 一1 图2采用64QAM信号,根据前面分析可知均衡算法 分3个阶段执行,为比较MAP算法中不同d值的选取对 {P >P fi D/s <D/s (19) MAP算法也有影响,我们分别取出3种情况来进行比较, 【Pm一1 fi Dis ≈0 其中m,nE[1,2,3,4],对于第三个条件,只要 当Dism比较小的时候我们就可以认为对应的P 为1。仿真 时取D西 4o.05。 故设计P 和Dis 的关系为 r1 if iDs 0.05 l『 0 fi iDs >0.05&minDis 0.05 P 一 m l1 ̄Dis s l∑(1//)/ (n+1)一 ( )+ —∑4一 P丁 n 2( 一 ( )) +。一 + — 二=_— —一 r ( ) 5仿真结果及其分析 本文采用的评价标准是算法的收敛速度和稳态剩余误 差。文献E14—15]中采用的性能指标是MSE,定义为 1 14q ̄.n ——MS )一 I(7z)一OEy( ]l ( ) + 式中:Q[ ( )]—— (n)的判决值,定义如下 Q[ ( )]一arg arin{l (,z)一Ski I } (23) sk/c。 本文对式(22)稍作修改,得到根均方误差(root mean square error,RMSE)指标,如式(24)所示,算法 的收敛速度和稳态剩余误差通过RMSE曲线来分析,因此 本文中的稳态剩余误差代表的实际意义是稳态剩余根均 方误差 RMSE(是) √ / 1 kNum I 卜Q[ )]I2 + (24) 本文通过matlab仿真来比较算法的收敛性能。仿真 中,发送信号取自64QAM信号集,取发送信号的长度 一 5000,均衡器的阶数K一11,信噪比SNR=40,均衡器抽 其中3种情况均取 =O.0003,第一种情况MAP1中的a按 文献E13]来选取,即第一阶段6=0.8,第二阶段d=I.6, 第三阶段d—o.6。第二种情况MAP2取第一阶段d=0.6, 第二阶段口一0.6,第三阶段a=0.6,第三种情况MAP3取 第一阶段d—I.6,第二阶段d=1.6,第三阶段d=1.6。由 图2可知,直接把自适应a取固定值是不准确的,不同值 对MAP算法的收敛性能有影响。 迭代步数 —————・———一MAP1—————{ ———一MAP2——MAP3 图2各阶段不同d值对MAP算法的影响 图3中为基于概率和自适应d的MAP算法与MAP算 法和基于自适应d的MAP算法的比较,我们可以看出基于 自适应。的算法的剩余根均方误差小于MAP盲均衡算法, 但是收敛速度慢于MAP盲均衡算法,原因是Y(n)在向 产生它的星座符号点靠近的过程中,当二者的偏差 不是 特别小的时候,区域内其它星座符号点对收敛的影响比较 大,所以算法从没收敛到完全收敛的这段过程,RMSE会 大于MAP盲均衡算法。而基于概率和自适应d的MAP算 迭代步数 一基于概率和自适应(r的MAP +MAP算法 一基于自适应 的MAP ̄t法 图3 3种算法的比较 ・3724・ 计算机工程与设计 2012正 法的剩余稳态误差和收敛速度均为3种算法最佳。 [-72 Chen&Semi—blind fast equalization of QAM channels using concurrent gradient-Newton CMA and soft decision-directed 6结束语 MAP盲均衡算法中的高斯信号方差值。的实际意义是 均衡器输出Y(n)偏离它对应准值的方差值,它是随着均 衡的进行而不断减小的,本文将其由固定值改进为可变值, scheme[J].Int J Adaptive Control and Signal Processing, 2010,24(6):467—476. r8]Chen S,Yao W,Hanzo L.CAM and soft decision-directed scheme for semi—blind beamforming of QAM systems[c]. Calgary,Canada:Proc VTC2008一Fall,2008:21 24. 仿真结果表明稳态剩余误差进一步的减小,但使得收敛速 度减慢,同时注意到概率随距离单调递减,所以运用不等 [9]Chen S.Low complexity concurrent constant modulus algo rithm and soft decision directed scheme for blind equatisation 概思想,对算法进一步的改进,得到了基于概率和自适应 的MAP盲均衡算法,改进算法比MAP算法和单纯的自适 应的MAP算法具有更好的均衡效果和收敛性能。 参考文献: E12 YANG Huibiao,LUO Daisheng,YU Yanmei.An improved constant modulus blind equalization algorithm FJI.Computer &Digital Engineering,2011,39(5):7-9(in Chinese).[杨 慧彪,罗代升,余艳梅.一种改进的恒模盲均衡算法[J3.计 算机与数字工程,2011,39(5):7-9.] E2]YU G,COWAN C F N.A self-structured adaptive decision feedback equalizer[J].IEEE Signal Processing Letters, 2006,13(3):169—172. E3]Wei X SH,David G M,Mulgrew C R A unified approach to dynamic length algorithms for adaptive linear equalizers EJ]. IEEE Trans on Commun,2007,55(3):908—920. E42 ZHOU Jun,LIU Peng,CAO Suzhi,et a1.Improved variable step-size dithered signed error constant modulus algorithm EJ]. Chinese Journal of Scientific Instrument,2008,29(6):699— 704(in Chinese).[周俊,刘鹏,曹素芝,等.改进的变步长 抖动符号误差常数模算法[J].仪器仪表学报,2011,32 (3):701—706.] L5]ZHANG Tian—yu.Research of OFDM system and Simulink simulation[J].Journal of Changchun University of Technolo— gY(Natural Science Edition),2008,29(6):699—704(in Chinese).[张天瑜.OFDM系统研究及其Simulink仿真[J]. 长春工业大学学报;自然科学版,2008,29(6):699—704.] [6]Icart S,Comon P,Rota L Blind parauniary equalization[J]. Signal Processing,2009,89(3):283-290. [J].IEE Proceedings of the Vision,Image and Signal Proce- ssing,2003,150(5):312—320. -[10]Hu Hengyun.XIE Ning.WANG Hui.Concurrent dithered signed—error constant modulus algorithm(DSE-CMA)with soft decision directed(SDD)scheme rC].International Con ference on future Computer and Communication,2010: 557—560. -[11]Shafayat Abrar.Evaluation of dispersion constants for signe ̄ derror blind equalization algorithms[c].International Con ference on Emerging Technologies,2005:188—192. [12]SU Jing,SHU Qin,WANG Ling wei.An improved CMA blind equalization algorithm using MAP estimation[J].Com— puter Engineering and Design,2010,31(21):4572—4575 (in Chinese).[苏静,舒勤,王灵伟.基于最大后验估计的 恒模盲均衡算法_J].计算机工程与设计,2010,31(21): 4572—4575.] [13]Chen S,Cook T B,Anderson I C.A comparative study of two blind FIR equalizers EJ].Digital Signal Processing,2004 (14):18—36. [142 GUO Ye-cai,ZHANG Yan ping.Hyperbolic tangent error function based variable step size blind equalization algorithm [J].Ship Science and Technology,2007,29(6):59—61(in Chinese).[郭业才,张艳萍.基于双曲正切误差函数的变步 长盲均衡算法[J].舰船科学技术,2007,29(6):59—61.] [15]ZHANG Cheng-yu.Research on the blind equalization algo rithm based on PDF and the modified CMA rD].Dalian: Dalian Maritime University,2010(in Chinese). [张成榆. 基于PDF盲均衡算法和修正恒模算法的研究[D].大连:大 连海事大学,2010.]