毕业论文-CMA盲均衡算法仿真研究

word文档可自由复制编辑CMA盲均衡算法仿真研究摘要盲均衡是一种新兴的自适应均衡技术，它不需要参考输入的训练序列来维持正常工作，仅依据接收序列本身的先验信息来均衡信道特性。自它出现后，就得到广泛的关注，并在许多领域中得到应用。本文系统地分析研究和归纳总结了盲均衡的基本理论。重点分析了Bussgang类盲均衡算法中的恒模(CMA,ConstantModulusAlgorithm）盲均衡算法。分析了传统CMA盲均衡算法的收敛性能，由于采用固定步长，使得收敛速度和收敛精度之间相互制约，其应用受到很大的限制。为了解决这一矛盾，本文提出了一种基于均方误差（MSE,MeanSquareError）的CMA盲均衡算法，这是一种利用时变步长来代替固定步长的自适应变步长CMA盲均衡算法，并进行了计算机仿真。结果表明改进算法相对于CMA算法收敛性能有一定的提高。关键字:盲均衡,恒模算法,变步长,均方误差

CMAblindequalizationalgorithmSimulationABSTRACTThispaperanalyzedsystematicallystudiesandsummariedtheblindbalancedelementarytheory.AnalysisfocusedontheBussgangtypeblindequalizationofconstantmodulusalgorithm(CMA,ConstantModulusAlgorithm)algorithmforblindequalization.ThispaperanalyzesofthetraditionalCMAblindequalizationalgorithmperformance,asaresultoftheuseoffixed-step,makingconvergencespeedandresidualerrorbecomeacontradiction,whichmakestheapplicationfieldsofCMAalgorithmlimited.Inordertosolvethecontradiction,thispaperderivesanimprovedCMAblindequalizationalgorithmutilizingthevaryofMSE.Thisisanadaptivevariablestep-sizeCMAblindequalizationalgorithm,whichusesatime-varyingstepsizetoreplacethefixedstepsize.ThesimulationwithcomputershowstheimprovedalgorithmshavethebetterconvergenceperformancethanCMAalgorithm.Keywords:blindequalization,ConstantModulusAlgorithm,variablestep-size,MeanSquareError目录摘要（中文） I摘要（外文） II1绪论 11.1研究盲均衡的目的和意义 11.2盲均衡的研究现状 21.3衡量算法收敛性能的指标 32恒模算法 42.1盲均衡的基本结构 42.2Bussgang类盲均衡算法 62.2.1决策指向算法 72.2.2Sato算法 72.2.3Godard算法 82.3恒模算法的提出 82.4恒模算法的理论推导 92.5步长因子对恒模算法收敛性能的影响 113基于剩余误差的变步长恒模盲均衡算法 173.1恒模算法中剩余误差的分析 173.2基于MSE的变步长恒模盲均衡算法 183.2.1基于MSE的变步长恒模盲均衡算法的表达形式 183.2.2算法性能分析 183.3基于MSE的变步长恒模算法的MATLAB实现 19结论 24参考文献 25附录 26致谢 321绪论盲均衡是一种新兴的自适应均衡技术，它不需要参考输入的训练序列来维持正常工作，仅依据接收序列本身的先验信息来均衡信道特性。因此，在数据通信系统中不必发送训练序列，可以提高信道效率，同时盲均衡技术还可以获得更好的均衡性能。盲均衡技术优越的性能使它受到更加广泛的关注，并在许多领域中得到应用。盲均衡技术可有效地应用于数字通信、雷达、地震和图像处理等系统。盲均衡技术己成为数字通信领域中热点研究的课题之一。在盲均衡的几种算法中，又以CMA（ConstantModulusAlgorithm）恒模算法的研究最为广泛。1.1研究盲均衡的目的和意义在数字通信系统中，带限发射、接收滤波器、放大器、时延与多径效应、发射机与接收机之间的相对运动、祸合效应和多址干扰等因素综合作用会使信号序列在传递过程中产生码间干扰和信道间干扰.为了降低误码率，必须对码间干扰进行适当的补偿。传统的克服码间干扰的方法是在接收端加均衡器，使均衡器的特性正好与信道的特性相反，使之能够准确补偿传输信道的特性，从而消除码间干扰。有些应用场合如无线移动通信中信道是时变的，为了准确地补偿信道的特性，均衡器应有及时调整参数、动态跟踪信道变化的能力，具有这种“智能特性”的均衡器称之为自适应均衡器。这种均衡器在数据传输之前，通常需要预先发送一段收端和发端都已知的训练序列。接收机测量出该序列通过信道后产生的变化或误差，并依据该误差信息对均衡器参数进行调整，最终使均衡器正好补偿信道特性，从而使接收机能够从均衡器输出中得到几乎无错的发送信号，保证数据的可靠传输。这段过程被称为训练，此时均衡器被称为工作在训练模式。训练过程结束后，数据传输开始，此时发送信号是未知的，为了动态跟踪信道特性可能发生的变化，接收机将均衡器输出的判决信号作为参考信号，用来测量信号通过信道后产生的误差，对均衡器输出的信号继续进行调整，此时均衡器工作在判决(DecisionDirected)模式。根据自适应滤波理论，均衡器在判决修正模式下能正常工作的条件是输入信号的眼图预先张开到一定程度(判决结果的错误率极低)，以保证均衡器可靠地收敛。如果这个条件不满足，就要由发端发送一个收端已知的训练序列对均衡器进行训练，使之收敛。因而训练过程也被称为均衡器的学习过程，对一般通信系统来讲是不可缺少的阶段。然而训练序列的使用有如下几点缺陷:(1)由于训练序列的传输占用了部分时间，有效的信息速率降低了。(2)对于严重的衰落信道，训练序列必须频繁发送。(3)当通信发生短时中断时，每一次新的通信开始之前必须发送训练序来初始化接收机。(4)在某些特殊应用场合，接收机无法得到训练信号(如在破译截获的敌方信号时)。由于自适应均衡器具有上述缺陷，使之不能适应现代数字通信系统高速度、大容量的发展趋势。因此，近年来人们致力于研究不借助训练序列，仅仅根据接收到的信号序列本身进行自适应均衡的技术盲均衡。与普通均衡器相比，盲均衡器具有收敛域大、应用范围广等特点。1.2盲均衡的研究现状1975年，日本学者Y·Sato在对传统的自适应均衡的均方误差函数进行了简单改进后，第一次提出应用于多幅度调制数据传输中的自恢复均衡的概念，后称之为盲均衡。自此以后，许多专家学者都投入到盲均衡的研究中，从不同方面采用各种代价函数和优化方法，得出许多应用于不同场合的盲均衡算法。目前，盲均衡的研究主要分为以下几类：基于高阶谱的盲均衡一般情况下，基于二阶统计量的盲均衡算法只能解决最小或最大相位信道的均衡问题，对非最小相位信道则无能为力。但是系统输出序列的高阶统计量既能反映信道传递函数的幅度信息和相位信息，又能有效抑制信道中的加性高斯噪声，从而能用于各种信道辨识与参数估计。基于神经网络的盲均衡信道均衡也可以看作为分类问题，把均衡器看成判决器，从而尽量精确地恢复发送序列。因此有很强分类功能的神经网络就很适合做均衡器。神经网络为非线性动态系统，它具有很大规模并行处理、高度的鲁棒性等特征，尤其适于处理复杂的非线性问题。基于信号检测的盲均衡有些文献将基于信号检测理论的盲均衡算法从原理上分为最大似然序列估计盲均衡算法，贝叶斯估计盲均衡算法，以及最小错误概率盲均衡算法等。Bussgang类盲均衡Bussgang类盲均衡以横向滤波器为结构，利用信号的物理特征选用合适的代价函数和误差控制函数来调节均衡器抽头，使得恢复信号接近于源信号。此类算法是以一种迭代方式进行盲均衡，并在均衡器输出端对输出信号作无记忆非线性变换。由于它是在传统自适应滤波的基础上发展而来，因此保留了传统自适应算法的简单性，复杂度低，运算量小，概念清楚，易于实现。但这类算法的缺点是算法收敛时间长，手电后稳态剩余误差大，对非线性或存在零点的信道均衡效果不好等。目前桥位经典的Bussgang类算法由Sato算法、决策指向算法、BGR算法、StopandGo算法、Godard算法等。1.3衡量算法收敛性能的指标衡量算法收敛性能的指标主要有收敛速度、误码特性、运算复杂度、跟踪时变信道的能力和抗干扰能力等。收敛速度均衡器开始工作后，需要一个收敛过程才能使均衡器的抽头系数由初值逐渐过渡到最优值，收敛速度越快，收敛过程所需时间越短，通信初期的误码数越少。误码特性在不增加算法计算复杂度和收敛速度满足要求的前提下，降低均衡器的误比特率（BER）具有重要意义。运算复杂度许多均衡算法尽管收敛速度快，但计算量太大，因而对硬件和软件要求很高，使其实际应用受到很大的限制。因此，在误码率满足要求的前提下，应降低均衡算法的计算复杂度。跟踪时变信道的能力算法跟踪时变信道的能力，主要体现在信道发生时变的情况下，算法能否收敛和稳定的问题。算法的跟踪能力受其原理和参数的制约。抗干扰能力抗干扰能力是算法对信道中叠加的噪声，尤其是突发强噪声干扰的抵抗能力。抗干扰能力差的算法遇到强噪声干扰时收敛性能变差甚至无法收敛。2恒模算法2.1盲均衡的基本结构图2-1为盲均衡原理框图。其中是发送序列，是未知信号的冲激响应（包含了发射滤波器、传播媒介和接受滤波器的综合作用），为系统接收序列，同时也是盲均衡器的输入序列，为噪声信号，为均衡器的冲激响应，为被均衡器恢复的信号，为判决输出信号。图2-1盲均衡系统输入序列假设为独立同分布序列，通过一未知时变离散时间传输信道，考虑加性信道噪声，得到均衡器接收序列可表示为：=(2-1)可知,是由和卷积而成,要想从中获得,就需要对进行反卷积或解卷积运算,或等价辨识传输信道的逆信道.当和已知时,可以获得。均衡器的训练就属于此种情况但当未知时,即3个参数中只有一个是已知,求解就相当困难,这就是盲均衡或盲解积。均衡器是线性自适应滤波器系统，它的输出为(2-2)若不考虑信道噪声的影响，则由信道输入端到均衡器输出端的冲激响应等于(2-3)因此均衡器输出可以写成(2-4)盲均衡的目的是通过算法调节均衡器权值使均衡器输出序列逼近于信道输入序列，这就要考虑到代价函数的选取以及采用的优化算法。如果通过以上的选取获得了一个理想均衡器，也即一个理想的逆滤波器，令表示理想均衡器的冲激响应，则它与信道冲激响应之间满足“理想逆关系”，表达如下（2-5）式中，为Kronecker函数。目前的盲均衡算法一般采用有限长抽头式横向滤波器，其结构如图2-2所示。图2-2横向滤波器的结构图其中，横向滤波器的长度为L，横向滤波器的输入为（2-6）滤波器的抽头系数为（2-7）则横向滤波器的输出可表示为==（2-8）理想的滤波器是无限长的，图2-2所示滤波器是截断的有限长滤波器，它是理想滤波器的近似模型，这就必然带来剩余码间干扰，滤波器的输出仅仅是源信号的估计值。因此误差信号为==（2-9）训练过程的任务是求出一组抽头系数，使均衡器能最有效地消除码间干扰，这组抽头系数称为最佳抽头系数。为了使均衡器获得最佳抽头系数，需要根据不同应用场合选用不同的优化算法，盲均衡算法用对均衡器输出信号的无记忆非线性变换来代替自适应算法中的期望信号。2.2Bussgang类盲均衡算法图2-3Bussgang盲均衡器的原理图图2-3为Bussgang类盲均衡器原理图。Bussgang类盲均衡算法作为盲均衡算法的一个分支，是在传统的自适应滤波器的基础上发展起来的。早期的盲均衡器以横向滤波器为基本结构，利用信号的物理特征选择合适的代价函数和误差控制函数来调节均衡器的权系数。这类算法是以一种迭代方式进行盲均衡，并在均衡器的输出端对数据进行非线性变换，当算法以平均值达到收敛时，被均衡的序列表现为Bussgang统计量。因此，此类算法称为Bussgang类盲均衡算法。Bussgang类盲均衡算法的显著特点是算法思路保持了传统自适应均衡的简单性，物理概念清楚，没有增加计算复杂度，运算量较小，便于实时实现。缺点是算法的收敛时间较长，收敛后剩余误差较大，没有解决均衡过程中的局部收敛问题，对非线性信道和存在零点的信道均衡效果不佳。Bussgang类盲均衡器采用一个非线性估计函数g(·)，使，用近似代替。如果一个随机过程满足下式条件时：(2-10)则该过程叫做Bussgang过程。式(2-10)揭示出，Bussgang过程应具有下述特性:均衡器输出序列的自相关函数等于用该输出序列作变元的无记忆非线性函数g(·)与输出序列之间的互相关函数。1952年了J.J.Bussgang第一个发现任何相关的高斯过程均具有上述性质。1955年，J.F.Barrett和D.G.Lampard进一步证明了所有具有指数衰减自相关函数的随机过程均具有这一性质，进一步推广了Bussgang的结论。不同的Bussgang类盲均衡算法具有不同的无记忆非线性函数g(·)，但都必须满足式(2-10)。归纳起来，Bussgang类盲均衡算法主要由以下两个公式表述，其中，式(2-11)为均衡器输出，式(2-12)为抽头系数迭代公式。(2-11)-2(2-12)式中，2L+1为均衡器长度，，µ为迭代步长因子。关于Bussgang算法的收敛性，有以下重要结论:若输入序列是亚高斯的，并且的二阶倒数为负值，则Bussgang算法是收敛的。Bussgang算法有三个非常有名的特例—(DD)决策指向算法、Sato算法、Godard算法。下面再分别介绍一下。2.2.1决策指向算法当Bussgang算法收敛，并且眼图“张开”时，均衡器便以决策指向模式工作，均衡器横向滤波器的抽头系数的最小均方误差即可以象自适应均衡器一样进行控制。图2-4决策指向均衡器的方框图决策指向(Decision-Directed)模式使用的无记忆非线性函数是一“阀值决策装置”。给定横向滤波器输出信号,阂值决策装置根据发射信号的字符集，对做出决策判断，使判断结果与最接近，例如，在二进制等概率数据序列的简单情况下，数据和决策取值分别为(2-13)将决策指向算法与Bussgang算法作一比较，可见决策指向算法是取g(.)=sgn(.)的Bussgang算法。2.2.2Sato算法M进制PAM(脉冲幅度调制)系统的盲均衡最早是Sato于1975年提出的。在Sato算法里，将代价函数定义为:=E(2-14)式中，为常数，定义为=.很显然，Sato算法是Bussgang算法取g(.)=sgn(.)时的一个特例。2.2.3Godard算法D.N.Godard[2]于1980年提出了一种可用于二维数据通信系统的盲均衡算法，它最大的特点是将幅度的均衡和相位恢复独立进行，互不干扰，因而允许灵活采用载波同步方案，这对载波偏移较大的系统特别有用。Godard在算法中应用了一种新的代价函数(2-15)式中，为一常数定义为(2-16)将式(2-15)两边对均衡器权向量求导可得代价函数对的梯度(2-17)去掉上式中的数学期望操作即为Godard迭代算法中的随机梯度，因此，均衡器抽头系数的更新公式为:(2-18)由上式可知，Godard算法是Bussgang算法中的无记忆非线性函数(2-19)2.3恒模算法的提出Godard最早提出了恒模算法(CMA)，它是Bussgang类盲均衡算法中最常用的一种。Godard算法无记忆非线性函数。表达式g(·)如下:(2-20)式中，p=1,2,当p=2时，Godard算法就是CMA算法。它通过调节线性均衡器的抽头增益来达到使代价函数减小的目的。CMA以其计算复杂度低、易于实时实现等优点，成为通信系统中广泛应用的盲均衡技术。恒模盲均衡算法适用于所有具有恒定包络（简称恒模）和一部分非恒包络（如QAM）的发射信号的均衡。CMA算法无记忆非线性函数g(.)为:(2-21)式中，是常数。根据信号传输理论和图2-1可知:均衡器的输入为:==(2-22)均衡器的输出为:===(2-23)CMA算法的权值迭代公式为(2-24)式中，为迭代步长因子，通常取足够小的正常数，它决定收敛的速度。2.4恒模算法的理论推导CMA算法的代价函数为:(2-25)选取这个代价函数的合理性在于，发送信号的功率应该是恒定的，均衡器输出信号的功率也应该是恒定的。按照最速下降法的迭代公式:(2-26)有:(2-27)因为=，故有:==2=2(2-28)于是:=4E(2-29)用随机梯度代替梯度的期望值，得到算法公式:4a(2-30)现在进一步考虑应该取什么值才是合理的。对均衡器的要求是:当达到理想均衡时，必须有:=0(2-31)所谓达到理想均衡，就是均衡器输出序城n)是发送序列x(n)的一个延时版本，即:=(2-32)其中，是一个固定的相位。由=0和式(3-9)得到:(2-33)也就是对应元素相等i=0,(2-34)注意到均衡器输入序列可以一般地写成:(2-35)式中，包括发送滤波器、信道和接收机前端(不含均衡器)的复合信道冲激响应;是频率偏移和相位抖动引起的时变相位移。各个序列统计独立，随机相位与发送序列互不相关。在向量中的元只有满足的项对和有贡献。这时显然有:=kE(2-36)以及:E=kE(2-37)式中，k是信道引入的确定性贡献。既然要求: (2-38)则对取值的要求就是:=(2-39)表2-1给出了Godard算法或常数模算法小结。表2-1中，CMA是对常数膜性能曲面进行随机梯度最小化运算的。与经过训练的均衡器的单峰MSE性能曲面相比，盲均衡器的常数模性能曲面是多峰的。误差曲面的多模式性和缺少期望响应信号大大影响了CMA的收敛性能。CMA在初始化、收敛速率与超量MSE等方面有它自己的特点。表2-1Godard算法或常数模算法小结运算等式均衡器误差更新Godard常数(1)初始化由于CMS误差曲面是非凸的，算法可能会收敛于一个非期望的最小值，这就说明了初始化过程的重要性。在实际中，所有的均衡器都用选择中心方法来初始化，即除了中心（参考）系数设定为大于某一常数外，所有其他的系数都设为零。(2)收敛速率经过训练的LMS算法有一个有界的收敛速率，因为二次误差曲面的Hessian矩阵（它决定了曲率）是恒定的。由于常熟模准则的误差曲面是多峰的，并且包含鞍点，所以CMA的收敛速率在鞍点附近较低，它与在一个局部最小值附近经过训练的LMS收敛速率相当。(3)超量MSE在经过训练的LMS算法中，超量MSE由步长、MMSE、滤波器系数的数量和输入信号的功率决定，并且CMA的超量MSE也取决于原信号的峭度。2.5步长因子对恒模算法收敛性能的影响实验一：用MATLAB对CMA算法进行了仿真，输入信号采用4QAM调制方式，信噪比为20dB,滤波器阶数为11,信道采用典型电话信道。步长分别为0.01、0.005、0.001，仿真实验运行总次数为3000次。=0.005+0.009-0.024+0.854-0.218+0.049-0.016(2-40)(a)收敛曲线(b)4QAM信号的星座图(c)均衡器输入星座图(d)步长0.01对应的均衡器输出星座(e)步长0.005对应的均衡器输出星座(f)步长0.001对应的均衡器输出星座图2-5不同步长CMA算法仿真图2-5(a)为4QAM信号通过典型电话信道采用不同步长值对应的收敛曲线比较。图2-5(b)为4QAM信号的星座图。图2-5(c)~(f)为4QAM信号通过典型电话信道采用不同步长值对应的均衡前后的星座图。图2-5(a)的仿真结果证实，采用大步长，能够加快收敛速度，但同时会带来大的稳态剩余误差和误码率。为了减小算法收敛后的稳态剩余误差和误码率应采用小步长，但这样会使算法收敛速度变慢。从图2-5(b)~(f)中可以看出，算法均衡后的星座更加集中、清晰，具有更小的稳态剩余误差和误码率。实验二：用MATLAB对CMA算法进行了仿真，输入信号采用4QAM调制方式，信噪比为15dB,滤波器阶数为7,信道采用普通信道。步长分别取0.01、0.001，仿真实验运行总次数为3000次。=1+0.3-0.3+0.1-0.1(2-41)(a)收敛曲线(b)4QAM信号的星座图(c)均衡器输入星座图(d)步长0.01对应的均衡器输出星座(e)步长0.001对应的均衡器输出星座图2-6不同步长CMA算法仿真图2-6(a)为4QAM信号通过普通信道采用不同步长值对应的收敛曲线比较。图2-6(b)为4QAM信号的星座图。图2-6(c)~(e)为4QAM信号通过普通信道采用不同步长值对应的均衡前后的星座图。图2-6(a)的仿真结果证实，采用大步长，能够加快收敛速度，但同时会带来大的稳态剩余误差和误码率。为了减小算法收敛后的稳态剩余误差和误码率应采用小步长，但会使算法收敛速度变慢。从图2-6(b)~(e)中可以看出，算法均衡后的星座更加集中、清晰，具有更小的稳态剩余误差和误码率。由实验一和实验二得知在两种不同信道下，迭代步长值越大,收敛速度就越快，但收敛后的稳态误差也就越大；减小步长值可以降低收敛后的稳态误差，但是会导致收敛速度的降低。Bussgang类盲均衡算法的一般格式是，先建立一个代价函数，使理想系统对应于代价函数的极小值点，然后采用某种自适应算法一步一步调整均衡器的抽头系数来寻找代价函数的极值点，当代价函数达到极值点后，抽头系数也达到了最优值。步长在算法收敛过程中起着非常重要的作用，采用大步长，每次调整抽头系数的幅度就大，体现到收敛性能上就是算法收敛速度和跟踪速度快，当均衡器抽头系数接近最优值时，抽头系数将在最优值附近一个较大的范围内来回抖动而无法进一步收敛，因而会有较大的稳态剩余误差。反之，采用小步长，每次调整抽头系数的幅度就小，算法收敛速度和跟踪速度慢，但当均衡器抽头系数接近最优值时，抽头系数将在最优值附近一个较小的范围内来回抖动而无法进一步收敛，因而稳态剩余误差较小。恒模算法采用固定步长，算法在收敛速度和收敛精度方面对调整步长的要求是相矛盾的，因而制约了恒模算法收敛性能的进一步提高。解决这一矛盾的最好方法是将自适应均衡中的变步长思想应用于恒模算法。在算法收敛期加大步长，提高收敛速度。算法收敛后降低步长，提高收敛精度。目前，变步长自适应均衡算法的主要研究成果有，用MSE作为控制步长变化的参量、用剩余误差的非线性变换作为控制步长变化的参量、用剩余误差的自相关函数作为控制步长变化的参量、用剩余误差的峰度作为控制步长变化的参量、用剩余误差和均衡器输入信号的互相关作为控制步长变化的参量，用梯度自适应变步长的方法来控制步长的变化，还有用误差信号的范数来控制步长的变化。后续章节将研究将变步长思想应用于恒模算法，来克服恒模算法采用固定步长所存在的缺陷，提高恒模算法的收敛性能。3基于剩余误差的变步长恒模盲均衡算法将变步长思想应用于恒模算法就是在算法收敛初期加大步长，以加快收敛速度，当算法收敛后，减小步长，以减小稳态剩余误差。在本章中，提出了基于剩余误差的变步长恒模盲均衡算法，分析了剩余误差的变化规律，指出将剩余误差直接用于步长控制的不足之处，提出将剩余误差的一种变换MSE、作为控制步长的参量，形成一种基于剩余误差的变步长恒模盲均衡算法，并通过计算机仿真实验验证了改进算法的收敛性能。3.1恒模算法中剩余误差的分析假设均衡器的时变最优权矢量为：(3-1)则有: (3-2)式中为零均值，独立同分布的干扰信号。将式(3-2)代入剩余误差的表达式，可得:=-===(3-3)式中，称为权误差矢量。在算法收敛过程中，由于逐渐向靠近，所以权误差矢量呈逐渐减小趋势，最后趋于零，所以式(3-3)中第一项也逐渐减小，最后趋于零。第二项为干扰信号。以上理论分析表明，剩余误差信号的变化趋势是由大到小，在算法开始时，均衡器权矢量距离最优权矢量最远，剩余误差最大，在算法收敛过程中剩余误差逐渐减小，算法收敛后达到最小。从以上分析可见，剩余误差的变换规律与变步长思想对步长变化规律的要求基本一致，但将剩余误差直接用于步长控制存在一些缺陷。首先变步长算法在收敛之前应一直采用较大步长才能真正起到加快收敛速度的作用，用剩余误差作步长控制往往是开始时步长较大，收敛速度也快，但剩余误差迅速下降，步长随之很快变小，收敛速度变慢，总体来看收敛速度得不到提高。其次，从式(3-3)可以看到剩余误差对干扰信号敏感，尤其是算法收敛后，如果信道中有突发的强干扰信号时，会很大，随之产生的大步长会引起误调，严重时可能会使算法发散。为更适合于步长控制，本章提出将剩余误差进行适当变换后，再来控制步长的变化。3.2基于MSE的变步长恒模盲均衡算法MSE的含义为=E=E，是剩余误差平方的期望值，本小节分析了用MSE来控制步长的合理性，提出了基于MSE的变步长恒模盲均衡算法。3.2.1基于MSE的变步长恒模盲均衡算法的表达形式参照图2-1盲均衡系统的原理框图可知，均衡器的输入为:==(3-4)均衡器的输出为:===(3-5)改进算法中抽头系数的迭代采用下式:(3-6)式中,为可变步长，其参数变化由下式来控制。==aE=aE(3-7)公式(3-4)~(3-7)就构成了基于MSE的变步长恒模盲均衡算法。在实际应用中，得到的方法是先对取平方，然后使之通过长度为L的移动矩形窗，再取平均来得到其估计值。a为比例因子，用于控制步长的取值范围。3.2.2算法性能分析(1)步长变化特性分析将式(3-3)代入式(3-7)得:=aE=aE{+2=aE+aE(3-8)在算法收敛过程中，由于逐渐向靠近，所以权误差矢量呈逐渐减小的趋势，最后趋于零，所以式(3-8)中第一项也是逐渐减小，最后趋于零。第二项为噪声的平均功率。以上理论分析证实，步长因子随着算法的收敛逐渐减小。用MSE控制步长的优势在于，如果信道中有突发强干扰信号时，变大，但由于经过加窗取平均，则可以削弱干扰信号的影响，使得MSE变化不大，这样可以减小因步长变化太大而引起的误调。(2)比例因子a的确定原则为了确保算法收敛，还必须合理选择式(3-7)中的参数a.a用于控制的取值范围，使得的最大值小于步长上界。由下式确定.=2/3tr(R)(3-9)式中，R为均衡器输入信号的自相关矩阵，tr(R)为R的迹。a具体取值应满足上述限定条件，并在仿真实验中调整确定。(3)矩形窗函数的长度L对算法性能的影响在盲均衡算法的实际应用中，信道可能是时变的甚至是突变的，还可能随机产生强噪声。在设计算法时必须考虑这些因素。信道的突变和强噪声都会引起剩余误差急剧增加，由此而引起的均方误差变化大小取决于矩形窗函数的长度L。由于均方误差的估计值是由L项剩余误差的平方再平均得到的，所以L越大引起的均方误差变化越小，步长变化也越小。因此，L的选择对算法的影响非常大，L取值越小，步长对信道突变和突发噪声越敏感，即对信道时变的跟踪能力越强。但步长对突发噪声敏感，就有可能对均衡器造成大的误调，二者是互相矛盾的。在选择参数L时，要根据实际应用场合具体确定。在信道时变严重，而强噪声干扰较少的环境下应选择较小的L值，以提高均衡器对信道时变的跟踪能力;在信道比较稳定，而干扰噪声较强的环境中应选择较大的L值，以减小强噪声引起的误调;在信道时变严重，干扰噪声也较强的环境下，只能折中处理。3.3基于MSE的变步长恒模算法的MATLAB实现实验一：输入信号分别采用4QAM调制方式，信噪比为20dB波器阶数为11。信道采用典型电话信道,仿真实验运行总次数为2000次。典型电话信道:(3-10)(a)收敛曲线(b)4QAM信号的星座图(c)均衡器输入星座图(d)CMA均衡器输出星座(e)改进CMA均衡器输出星座图3-1两种算法仿真图图3-1(a)给出了4QAM信号通过典型电话信道后改进算法和恒模算法的收敛曲线.图3-1(b)是4QAM信号的星座图。图3-1(c)~(e)给出了4QAM信号通过典型电话信道后恒模算法和改进算法均衡前后的星座图。从图3-1(a)中可以看出，4QAM信号通过典型电话信道时，改进算法均具有较快的收敛速度和较小的稳态剩余误差。从图3-1(b)~(e)中可以看出，算法经过均衡后改进算法的星座更加集中、清晰，即改进算法具有更小的稳态剩余误差和误码率。实验二：普通信道输入信号分别采用4QAM调制方式，信噪比为20dB,波器阶数为11。信道采用普通信道,仿真实验运行总次数为2000次。普通信道:.(3-11)(a)收敛曲线(b)4QAM信号的星座图(c)均衡器输入星座图(d)CMA均衡器输出星座(e)改进CMA均衡器输出星座图3-2两种算法仿真图图3-2(a)给4QAM信号通过普通信道后改进算法和恒模算法的收敛曲线。图3-2(b)给出4QAM信号的星座图。图3-2(c)~(e)给出了4QAM信号通过普通信道后改进算法和恒模算法均衡前后的星座图。从图3-2(a)中可以看出，4QAM信号通过普通信道时，改进算法具有较快的收敛速度和较小的稳态剩余误差。从图3-2(b)~(e)中可以看出，信号经过均衡后改进算法的星座更加集中、清晰，即改进算法具有更小的稳态剩余误差和误码率。由实验一和实验二得知在两种不同信道下改进CMA算法收敛速度都明显快于CMA算法而且有较小的稳态剩余误差。由此可见改进CMA算法的性能要优于CMA算法。盲均衡优于传统自适应均衡之处在于：不需要用训练序列，发信端无需任何改动，仅在接收端改变算法即可大幅度提高通信系统可靠性，因此可以很好的运用于多点通信系统和广播系统中的均衡问题。此项技术的实际应用，对于提高接收信号的质量、保证信息的准确可靠，具有十分重要的意义。本文介绍了CMA算法的特点，给出了CMA迭代算法的具体步骤。然后，采用CMA盲算法对自适应滤波器进行均衡，并对CMA盲均衡算法的性能进行了研究并提出了一种基于MSE的变步长恒模算法。计算机仿真结果表明，盲均衡算法迭代步长的选取会影响到算法的收敛速度和稳态误差的大小，在具体选择迭代步长时，在保证盲均衡算法收敛的范围内，可根据实际需要在二者之间作出折中选择;均衡器的阶数在比较高时，对盲均衡算法收敛性能的影响己经非常小，因此，对于收敛性能相当的均衡器，可以选择其中阶数较小的。还表明了，改进算法性能优于恒模算法。另外，随着信噪比的增加，CMA盲均衡算法的收敛性能也相应提高。结论盲均衡是一种新的自适应均衡技术，它不再需要参考输入的训练序列来维持正常工作，仅依据接收序列本身的先验信息来均衡信道特性。因此，在数据通信系统中不必发送训练序列，这样就提高了信道效率。同时盲均衡技术还可以获得更好的均衡性能。盲均衡技术优越的性能使它受到更加广泛的关注，并在许多领域中得到应用。盲均衡技术可有效地应用于数字通信、雷达、地震和图像处理等系统，现已成为数字通信领域中热点研究的课题之一，相信未来的盲均衡技术将会得到更大的发展和更广的应用。本文所做的主要工作有:(1)分析了传统CMA算法的收敛性能，并进行计算机仿真实验加以验证。传统CMA算法采用固定步长造成收敛速度与收敛精度的矛盾，使其应用受到很大的局限。为了提高CMA算法的收敛性能，将自适应均衡算法中变步长的思想引入到CMA算法中。(2)提出了一种基于MSE变换的CMA盲均衡算法。对改进算法进行了理论分析，研究了改进算法中参数的选取原则，计算机仿真表明改进算法相对于CMA算法收敛性能有一定的提高。盲均衡是一项涉及许多知识领域的新兴的综合技术。特别是随着通信技术的飞速发展，使盲均衡技术的应用领域更加广泛。本文关于盲均衡技术所做的研究工作是十分有限的，有关盲均衡技术需要做的研究工作还很多。(1)Bussgang类盲均衡算法代价函数的非凸性，使之容易产生局部收敛。因此，Bussgang类盲均衡算法代价函数的凸性，是一研究的热点。(2)可以将更加先进的数学理论方法引入到盲均衡技术的研究中，提高盲均衡算法的性能。(3)对有盲均衡技术的理论知识加以实现、应用。由于作者的知识水平有限，文中还存在许多不足的地方，敬请各位老师批评指正，提出宝贵意见。参考文献[1]陈怀琛.MATLAB及其在理工课程中的应用指南[M].西安：西安电子科技大学出版社，2004[2]周炯槃，庞沁华，续大我，杨鸿文.通信原理(第三版)[M].北京：北京邮电大学出版社，2008.8:148-156.[3]高鑫.无线通信系统中Bussgang族盲均衡算法的研究[D].大连，大连海事大学通信与信息系统专业，2007[4]郭丽华.自适应盲均衡算法在通信系统中应用的研究[D].哈尔滨，哈尔滨工程大学信号与信息处理专业，2003[5]张雄.基于Bussgang技术盲均衡算法的研究[D].太原，太原理工大学电路与系统专业，2003[6]刘术平.恒模算法及其在盲均衡中的应用[D].南京，南京邮电大学通信与信息系统专业，2006[7]郭晓宇.改进型恒模盲均衡算法的研究[D].太原，太原理工大学电路与系统专业，2005[8]丁志高.常模量算法(CMA)及在通信信道均衡中的研究[D].南京，东南大学信号与信息处理专业，2004[9]赵宝峰.变步长盲均衡算法的研究[D].太原，太原理工大学电路与系统专业，2004[10]陈怀琛.数字信号处理教程—MATLAB释义与实现[M].北京：电子工业出版社，2008[11]梁虹.信号与系统分析及MATLAB实现[M].北京：电子工业出版社，2002[12]AlbertBenveniste,MauriceGoursat.BlindEqualizers[J].IEEETrans.oncommunications,1984(32):871~879.[13]PhilipSchniter,C.RichardJohnson.DitheredSigned-ErrorCMA:Robust,Computation-lly,EfficientBlindAdaptiveEqualization[J].IEEETrans.onsignalprocessing,1999,(47):1959-1597. 附录第2章实验一程序：clearallM=4;n=5000;m=3000;h1=[0.0050.009-0.0240.854-0.2180.049-0.016];L=11;u1=0.01;u2=0.005;u3=0.001;mse_av1=zeros(1,n-L+1);mse_av2=zeros(1,n-L+1);mse_av3=mse_av2;snr=20;forj=1:ms=randint(1,n,M);s1=qammod(s,M);m1=sum(abs(s1).^4);m2=sum(abs(s1).^2);R2=m1/m2;s2=filter(h1,1,s1);x=awgn(s2,snr,'measured');w1=[00000100000];w2=w1;w3=w1;fori=1:n-L+1y=x(i+L-1:-1:i);z1(i