第32卷第6期 字 航 学 报 Vo1.32 N0.6 2011年6月 Journal of Astronautics June 2011 卫星高速宽带数传系统的变步长LMS均衡 算法研究 李海滨,韩延光,庄述燕,关新平 (燕山大学工业计算机控制工程河北省重点实验室,秦皇岛066004) 摘要:对码间干扰实现快速有效的均衡是卫星高速宽带数字传输的关键之一,针对s型函数(Sigmoid函数 和双曲正切函数)变步长LMS均衡算法在误差e(n)接近零处 (n)变化大的缺陷及算法复杂度高的问题,提出了 一种改进的变步长LMS自适应均衡算法。算法通过简化步长更新式的复杂度,增加误差信息量,优化改进的s型 函数结构,建立了新的步长因子 ( )与误差函数e( )非线性关系,使算法误差接近零处步长因子变化平滑缓慢, 消除不相关噪声序列的影响,且算法复杂度比其它同类改进算法小。仿真结果表明新算法性能优于其它改进算法。 文中给出了卫星高速宽带数字传输系统的仿真模型,将新算法应用到系统模型的自适应均衡器中,系统误码性能 得到较大改善。 关键词:卫星高速宽带数字传输;自适应均衡;变步长;LMS算法 中图分类号:TN911.72 文献标识码:A 文章编号:1000—1328(2011)06—1423—05 DOI:10.3873/j.issn.1000—1328.2011.06.033 A Modified Variable Step LMS Adaptive Filtering Algorithm for Satellite High-Speed Broadband Digital Transmission System LI Hai—bin,HAN Yan—guang,ZHUANG Shu—yan,GUAN Xin—ping (Key Lab of Industrial Computer Control Engineering of Hebei Province,Yanshan University,Qinhuangdao 066004,China) Abstract:Equalizing inter—symbol interference(ISI)quickly and effectively is one of keys for satellite high—speed broadband digital transmission.A kind of S-type function(Sigmoid function and hyperbolic tangent function)variable step LMS algorithms is discussed in this paper.These algorithms have flaws that (n)changes a lot when error e(n)is close to zero and high complexity.Therefore,an improved variable step LMS adaptive equalization algorithm is proposed.Through simplifying the complexity of step formula update,increasing the amount of error information,optimizing the structure of modiifed S-type functions and building new nonlinear relations between step factor (n)and error function e(n),the step of new algorithm changes not only smoothly and slowly when error e(n)is close to zero,but also eliminates the impact of uncorrelated noise sequences.The complexity of the new algorithm is smaller than other similar modiifed algorithms.Simulation results show that the performance of new algorithm is superior to other modiifed algorithms.The satellite high—speed broadband digital transmission system simulation model is also given in this paper.When the new algorithm is applied to the adaptive equalizer of system model, the system bit error rate(BER)performance is improved greatly. Key words:Satellite high—speed broadband digital transmission;Adaptive equalization;Variable tap—length;LMS algorithm 言T脯孵壮ubeAmph fier, ̄喜 在卫星高速宽带数字传输系统中,星上转发器 特性,当它工作在饱和区或接近饱和区时,其群时延 收稿日期:2010—06-07; 修回日期:2OlO一07—14 1424 宇航学报 第32卷 曲线的不理想,会引起信道的严重非线性失真。采 用均衡器补偿信道失真产生的码问干扰是一种重要 的方法。实际的卫星高速宽带数字传输系统要求均 衡器能够基于对信道特性的测量随时调整自身的系 数,以适应信道特性的变化,而寻求高速的白适应均 衡算法是自适应均衡器实现的关键。 传统的定步长LMS白适应均衡算法在收敛速 度、时变系统跟踪速度与稳态精度方面对步长因子 的要求是相互矛盾的。为此,文献[1]提出了一种 Sigmoid函数变步长的LMS算法(SVSLMS),并给出 了步长因子的调整原则。然而,文献[2]分析指出 此算法的Sigmoid函数过于复杂,且在误差e( )接 近零处变化太大,不具有缓慢变化的特性,使得 SVSLMS算法在自适应稳态阶段仍有较大的步长变 化。文献[3,4]通过调整更新式,提高了算法性能, 但缺点是复杂度提高,工程实现困难。基于此,本文 提出一种改进的变步长LMS自适应均衡算法,利用 步长因子和误差函数新的非线性关系,不仅误差接 近零处变化缓慢,而且步长更新式更简单,大大减少 了算法的运算量。与文献[2—4]算法进行仿真比 较,本文算法具有较小的算法复杂度,且收敛速度和 跟踪速度最快。应用到卫星高速宽带数字通信均衡 系统中,有效改善了系统的误码性能。 1本文变步长LMS算法 最小均方误差(LMS)均衡算法 是一种以期 望响应d( )和输出信号Y( )之间误差的均方值最 小为准则的。输入信号 依次经过m一1个延迟单 元,在n时刻构成了一个信号矢量X(n),输入信号 矢量为: X(/7,)=[ ( ), (n一1),…, (n一,n+1)] (1) 加权矢量为: W( )=[W W …,W ] (2) 输出信号为: y(n)=∑ ̄1)ni (n—i+1) (3) 误差信号为: e(n)=d(n)一Y(n) =d(n)一W(n)‘X(n) (4) 基本LMS均衡算法的迭代方程 为: W(n+1)=W(n)+ ( )e( ) (/7,) (5) 其中 (n)为步长因子,当 (n)为常数时,算法为 固定步长LMS算法;当肛(n)为一个变数时,算法 为变步长LMS算法 。为了保证算法在迭代过 程中收敛,步长因子需满足0< (g/,)<1/A…,其 中A…为输入信号矩阵X( )的最大特征值。 本文提出更新式如下: (n)= (n)exp(O/}e(n)1) (6) 如果 (n) …(n),则有 (n)= (n)。其 中 …(n)为步长调整值的最大值。参数 >0和 >0用于控制肛(n)与e(n)关系曲线的形状。 当 :0.1、 …(n)=0.2时, (n)与e(n)的 关系曲线随ot的变化情况如图1所示。当 =5、 …(n)=0.2时, (n)与e(n)的关系曲线随卢的 变化情况如图2所示。由图1和图2可知:在收敛范 围内,初始阶段J e(n)l值较大,对应的步长值 ( ) 较大。当自适应均衡算法达到稳态阶段,此时 le(n)l值变小,相应的步长值变小,且具有缓慢变 化的特性。这里对肛(n)与e(n)关系曲线的平滑度 作如下定义: =a/b,其中0为当 (n)=tx…(n) 时e(n)跨越的区问长度,b为当IX( )= 时e(凡) 跨越的区间长度, 。为曲线接近零时步长选取的参 考极小值。易得,当 越小时曲线接近零处的平滑性 越好,反之,曲线越尖锐,平滑性越差。参数 、 和 …(n)的选择取决于实际系统的具体条件,且需要 满足下列条件:0< 。 (n)<1/A…。以下参数选取 的讨论均是在满足收敛条件下进行的,将不再重复 说明。 图1 固定, (n)与e(n)的关系曲线图 Fig.1 Relation diagram of ( )and e(n) when卢is fixed 由图1可以看出,当 =0.1和 =0.1时,曲 线开口较大,导致自适应均衡算法在迭代初始阶段 I e(n)l较大,步长值 (n)却相对较小,起不到快速 收敛和快速跟踪的作用,且l e(n)I接近零时,曲线 第6期 李海滨等:卫星高速宽带数传系统的变步长LMS均衡算法研究 1425 表现相对尖锐,说明步长值变化和稳态误差较大。当 =100和 =0.1时,曲线开口较小,算法在初始阶 2均衡算法-眭能分析 段始终以最大的步长值更新权向量,从而实现快速 收敛和对时变系统的快速跟踪。均衡算法步长值经 收敛性能、时变信道跟踪性能和稳态性能是检 验均衡算法性能优劣的重要指标。下面分析参数OL 和卢对改进的变步长LMS自适应均衡算法性能的 过一段骤然下降的过程,该过程同样反映算法的收 敛速度和跟踪速度较快,在l e(n)l接近零时, (n)曲线变化平滑缓慢地接近于零。易得,参数oL 对l e(n)l与 (n)关系曲线的开口度影响更大。 影响。仍然采用传统s型函数改进算法” 中的模 拟条件:自适应均衡器阶数L=2;不妨先将卫星高 速宽带数字传输系统模拟为FIR滤波器,其系数为 0 2 0・15 一 0.05 0 图2 固定, ( )与e(n)的关系曲线图 Fig.2 Relation diagram of/z(n)and e(n) when is fixed 由图2可以看出,当O/=5和 =0.001时,曲线 底部较大,同样导致自适应均衡算法在迭代的初始阶 段J e(n)J较大,步长值 (n)却相对较小,起不到快 速收敛和快速跟踪的效果,且l e(n)l还相对较大 (如l e(n)J:0.5)时,曲线就过早地进入平滑缓慢 的变化状态,说明步长值过早地变小,导致收敛速度 和跟踪速度变慢。当 =5和 =10时,曲线底部宽度 较小,均衡算法在初始阶段同样以最大的步长值更新 权向量,从而实现快速收敛和对时变系统的快速跟 踪。算法步长值经过一段平滑下降的过程,该过程同 样反映算法的收敛速度和跟踪速度较快,在l e(n)I 接近零时, (n)曲线变化相对(Ot:100和/3=0.1) 尖锐地接近于零。易得,参数 对l e(n)I与/x(n)关 系曲线的底部斜率有更大的影响。 综上分析,Ot越大,曲线上部的开口度越小,说 明算法在初始阶段收敛速度越快,但如果O/过大,则 会导致平滑度k越尖锐,使e(n)接近零时步长值变 化大。 越大,虽然曲线底部斜率变大,但如果此时 选择最佳的Ot值,使e(n)接近零时曲线平滑度k最 小,则可以得到最佳的曲线形状。实验得出最佳参数 为Ot=1000和/3=150。 =[0.8,0.5] ,在第500个迭代次数时刻卫星 高速宽带数字传输系统发生时变,系数矢量变为 =[0.4,0.2] ;参考输入信号 (n)是零均值, 方差为1的高斯白噪声; (n)是与 (n)不相关的高 斯白噪声,其均值为零,方差为0.04。分别做200次 的仿真,迭代次数为1000,然后求其统计平均, 得出学习曲线。 0 _2 —4 譬.6 s 至一1。 .12 —14 当/3:0.O01,/x…(n)=0.08时,均衡算法收 敛曲线随 的变化情况如图3所示。可以看出,当O/ =5时,白适应均衡算法收敛曲线的性能很差。算法 迭代400次后才趋于收敛,稳态误差约为一12dB,空 间信道发生突变后,算法跟踪性能很差。当Ot=30 时,算法性能得到很大改善。算法迭代100次后收 敛,且稳态精度高,约为一14dR,卫星高速宽带数字 传输系统发生突变后,算法跟踪速度较快。当 = 1000时,均衡算法性能进一步提高,算法迭代70次 后收敛,稳态精度变化不大。卫星高速宽带传输系 统发生突变后,算法跟踪速度进一步加快。图中参 数 在5到30之间变化时,算法收敛曲线变化较明 显,性能改善较快,说明对算法性能的影响较大,在 30到1000之间变化时,算法收敛曲线变化较小,性 能提高较慢,说明对算法性能的影响较小,随着Ot的 1426 O 4 Erp ∞ 氤一0 罐m 宇航学报 第32卷 进一步增大,算法收敛曲线变化不明显,性能几乎 没有提高。 另一类s型函数(双曲正切函数)算法存在的缺陷 进行的改进,该算法获得了同其它s型函数改进算 法类似的性能,但算法复杂度高,运算量仍相对较 大。本文均衡算法运算量与文献[2]、文献[3]和文 献[4]的运算量比较如表1所示,其中N为均衡器 一.12 —14 一l6 当 :1000, …(几)=0.08时,自适应均衡算 法收敛曲线随|8的变化情况如图4所示。可以看出,当 口=0.0001时,自适应均衡算法收敛曲线的性能很 差。算法迭代300次后才趋于收敛,稳态误差约为 一13dB,卫星高速宽带数字传输系统发生突变后,算 法跟踪性能很差。当 =0.1时,算法性能得到很大 改善。算法迭代60次后收敛,且稳态精度高,约为 14dB,未知系统发生突变后,算法跟踪速度较快。当 =150时,算法性能进一步提高,算法迭代40次后收 敛,稳态精度变化不大。卫星高速宽带传输系统发生 突变后,算法跟踪速度进一步加快。图中参数IB在 0.0001到0.1之间变化时,算法收敛曲线变化较明 显,性能改善较快,说明对算法性能的影响较大,在 0.1到150之间变化时,算法收敛曲线变化较小,性能 提高较慢,说明对算法性能的影响较小,随着 的进 一步增大,算法收敛性能会有所改善,但不明显。综上 对参数 和 的分析,实验得出使算法收敛速度快且 稳态性能好的最佳参数为O/=1000和 =150。 3 仿真研究 3.1本文算法与其它算法性能比较 文献[2]提出的均衡算法克服了Sigmoid函数 在自适应稳态阶段步长调整过程中的不足,并指出 Sigmoid函数过于复杂。但该文算法复杂度同样相 对较大。文献[3]提出的均衡算法同样是针对 Sigmoid函数存在不足进行的改进,但步长值更新式 只是对Sigmoid函数的简单修改,复杂度比Sigmoid 函数还高,其实质还是Sigmoid形式。文献[4]针对 阶数 表1均衡算法运算量 Table 1 Equalization algorithm operations 由表1可以看出,本文均衡算法同其它三种均 衡算法相比,本文算法比文献[2]和[4]算法乘法运 算增加一次,但加减运算分别减少2次和3次。由 于DSP芯片(如TMS320系列)具有专用的硬件乘 法器,乘法可在一个指令周期内完成,故本文算法的 运算量相对较小,复杂度低。 0 一矗 兽一6 ∞ 矍.8 2 10 一l2 一l4 .16 0 200 400 600 800 1000 迭代次数 图5本文算法与文献[2]、[3]和[4]算法性能比较图 Fig.5 Performance comparison diagram among algorithms in this paper,Ref.[2],Ref.[3]and Ref.[4] 当本文算法选取参数 …(n)=0.08、Ot= 1000和I8=150,文献[2]算法参数 =300和 = 0.2,文献[3]算法参数 :10000和 =0.2,文献 [4]算法参数Ot=80、 =0.2和h:2000。分别做 300次仿真,本文算法与文献[2]、[3]和[4]算 法性能比较图如图5所示。由图中可以看出,在四 种均衡算法稳态性能相同的条件下,本文算法具有 更快的收敛速度、跟踪速度和较好的稳态性能。因 此,本文均衡算法比文献[2]、[3]和[4]提出的均衡 算法性能更优越。 第6期 李海滨等:卫星高速宽带数传系统的变步长LMS均衡算法研究 1427 3.2本文算法在卫星高速宽带数传系统中的应用 卫星高速宽带数传系统的仿真模型如图6所示。 原始基带 信号m)厂—————1 ——— _ 调制器 — r——一发射机 卫星高速宽 带传输信道 二 H 图6卫星高速宽带数传系统仿真模型 Fig.6 Simulation model of satellite high—speed broadband digital transmission system 在此选取某卫星高速宽带数字传输系统参数,原 始基带信号的符号速率375M符号/秒,卫星信道群时 延失真在2倍奈奎斯特带宽内为6ns,调制方式为 QPSK。数据量10 个,将本文算法应用到该系统的自 适应均衡器中,得到如图7所示的误码率比较图。 懈 露 媸 图7误码率比较图 Fig.7 Bit error rate(BER)comparison diagram 由图7可以看出,发送信号经过卫星高速宽带 数字传输信道后,产生了严重的码间干扰,致使系 统均衡前误码率很大。当引入基于本文算法的自适 应均衡器后,系统误码性能得到较大改善。误码率为 尸 =10 时,信噪比损失约为3.9dB。 4 结 论 本文提出了改进的变步长LMS自适应均衡算 法,研究了参数 和 对算法的影响,分析了一类S 型函数(Sigmoid函数和双曲正切函数)算法及其改 进算法的缺陷,重点说明了新算法比其它改进算法 复杂度低、运算量小,在保证算法稳态性能的同时, 提高了算法收敛和跟踪性能。仿真结果验证了新算 法的优越性。在卫星高速宽带数字传输系统中具有 很好的应用前景。 参 考 文 献 [1] 覃景繁,欧阳景正.一种新的变步长自适应滤波算法[J].数 据采集与处理,1997,12(3):171—174.[Tan Jing.fan, Ouyang Jing—zheng.A noval variable step size LMS adaptive ifltering algorithm based on sigmoid function[J].Journal of Data Acqusition&Processing,1997,12(3):171—174.] [2] 高鹰,谢胜利.一种变步长LMS自适应滤波算法及分析[J]. 电子学报,2001,29(8):1094—1097.[Gao Ying,Xie Sheng- li.A variable step size LMS adaptive filtering algorithm and its analysis[J].Acta Electronica Sinica,2001,29(8):1094— 1097.] [3] 罗小东,贾振红,王强.一种新的变步长LMS自适应滤波算 法[J].电子学报,2006,34(6):1123—1126.[Lno Xiao. dong,Jia Zhen—hong,Wang Qiang.A new variable step size LMS adaptive filtering algorithm[J] Acta Electronica Sinica,2006, 34(6):1123一I126.] [4] 张中华,张瑞金.一种新的变步长LMS自适应滤波算法及性 能分析[J].系统工程与电子技术,2009,31(9):2238— 224 1.[Zhang Zhong—hua,Zhang Rui—jin.New variable step size LMS adaptive filtering algorithm and its performance analysis[J]. Systems Engineering and Electronics,2009,3 1(9):2238— 2241.] [5] 邱天爽,魏东兴.通信中的自适应信号处理[M].北京:电子工 业出版社,2005. [6] 张贤达,保铮.通信信号处理[M].北京:国防工业出版社, 2000. [7] Kwong R H,Johnston E W.A variable step size LMS algorithm [J].IEEE Transactions on Signal Processing,1992,40(7): 1636—1642. [8] Aboulnasr T,Mayyas K.A robust variable step-size LMS—type algorithm:analysis and simulations[J].IEEE Transactions on Signal Processing,1997,45(3):631—639. 作者简介:李海滨(1978一),男,副教授,硕士生导师,主要从 事航天测控技术、机器视觉等领域的研究。 通信地址:秦皇岛燕山大学电气工程学院(066004) 电话:(0335)8387556 E-mail:hbli@ysu edu.cn (编辑:余未)