递归最小二乘RLS自适应均衡算法

1、第三章递归最小二乘算法 第三章递归最小二乘(RLS)自适应均衡算法引言在自适应滤波系统中，最陡梯度( LMS)法由于其简单获得了广泛的应用。但各种 LMS 算法均有收敛速度较慢(收敛所需码元数多)， 对非平稳信号的适应性差( 且其中有些调整延时较大) 的缺点。究其原因主要是LMS 算法只是用以各时刻的抽头参量等作该时刻数据块估计时平方误差均最小的准则，而未用现时刻的抽头参量等来对以往各时刻的数据块均作重新估计后的累积平方误差最小的原则( 即所谓的最小平方( LS) 准则 ) 。为了克服收敛速度慢，信号非平稳适应性差的缺点，根据上述内容，可采用新的准则，即在每时刻对所有已输入信号而言重估的平方误

2、差和最小的准则( 即 LS准则 ) 。从物理概念上可见，这是个在现有的约束条件下利用了最多可利用信息的准则，即在一定意义上最有效，信号非平稳的适应性能也应最好的准则。这样建立起来的迭代方法就是递归最小二乘( RLS： Recursive Least Square)算法，又称为广义 Kalman 自适应算法。用矩阵的形式表示RLS 算法非常方便，因此我们首先定义一些向量和矩阵。假定在时刻t，均衡器的输入信号为rt，线性均衡器对于信息符号的估计可以表示为KI?(t)cj (t 1)rt j式 (3-1)jK让 cj(t 1)的下标 j从 j 0到 j N 1 ，同时定义y(t) vt K，则 I?

3、(t)变为N1I?(t)cj (t 1)y(t j)j0CN(t 1)YN (t)式 (3-2)其中CN(t 1) 和YN (t) 分别为均衡器系数cj(t 1) ， j 0,1, ,N 1 和输入信号y(t j)， j 0,1, ,N 1的列向量。类似的， 在 DFE 均衡器结构中，均衡器系数cj (t)，j 0,1, ,N 1 的前 K1 1个系数为前向滤波器系数，剩下的K2N K1 1 为反馈滤波器系数。用来预测I?(t)的数据为rtK1 ,rt1,It1,ItK2， 其中Itj ,1 j K 2为判决器先前作出判决的数据。这里，我们忽略判决器判错的情况，因而It j I t j ,1

4、jK 2。同时为方便起见定义vt K1 j (0 j K1)y(t j)式 (3-3)I t K1 j (K1 j N 1)YN(t) y(t), y(t 1), ,y(t N 1)rt K1, ,rt 1,rt,It 1, ,It K2式 (3-4)RLS 自适应算法RLS 算法对于I?(t) 的估计可以从下面的式子得到。假定我们的观测向量为YN(n)，n 0,1, t，我们期望得到均衡器的系数向量CN(t) 使得均方误差的加权平方和t(n) wt n | eN (n, t) |2式 (3-5)n0最小。其中误差定义为eN(n,t) I(n) CN(t)YN(n)式 (3-6)w 代表遗忘因

5、子，0 w1。 这样我们对过去的数据引入了一个指数权，这对于信道特性为时变的情况非常合适。关于权向量CN (t) 的 (n) 最小化便得到下面的线性方程RN(t)CN(t)DN (t)式 (3-7)其中RN(t)为信号的自相关矩阵，定义为tRN (t)wt nY*N (n)Yn(n)式 (3-8)n0DN(t)为互相关向量tDN(t)wt nI(n)Y*N(n)式 (3-9)n0式 (3-7)的解为著名的Wiener Hopf 方程CN(t)RN1(t) DN(t)式 (3-10)为了避免复杂的求逆运算，引入一N N 矩阵PN (t)CNPN(t) PN(t1) KN (t)YN(t)PN(t

6、 1)式(3-28)(t)由式 (3-8)有t tnRN (t)wt nYN*(n)YN (n)n0wRN (t 1) YN* (t)YN (t)又由矩阵求逆引理有：RN1(t)1 RN1(t 1)RN1(t 1)YN*(t)Y1N(t)RN1*(t 1)N w Nw YN (t)RN1(t 1)YN* (t)在上式中定义PN (t)RN1 (t)，令N(t) YN (t)PN (t 1)YN* (n)KN(t)PN(t 1)YN* (t)w N (t)N (t)为一标量，KN (t)为一N 维矢量，称为Kalman 增益向量。则PN(t)1 PN (t 1) KN (t)YN(t)PN(t

7、1)w假定我们在式(3-16)两边右乘以YN* (t)，1PN(t)YN*(t)PN(t1)YN*(t)KN(t)YN(t)PN(t1)YN*(t)w1w N(t)KN(t) KN(t) N(t) wKN (t)Kalman 增益向量可以被定义为PN (t)YN* (t)由于CN (t) PN (t)DN (t)DN (t) wDN (t 1)I(t)YN* (t)我们得到1CN(t)PN(t 1) KN(t)YN(t)PN(t 1)wDN(t 1)I(t)YN*(t)w1*PN(t 1)DN(t 1)I(t)PN(t 1)YN*(t)w式 (3-11)式 (3-12)式 (3-13)式 (3

8、-14)式 (3-15)式 (3-16)式 (3-17)式 (3-18)KN(t)YN (t)PN (t 1)DN (t 1)1*I (t)KN(t)YN (t)PN(t 1)YN*(t) wCN(t 1)KN(t)I(t) YN (t)CN(t 1)式 (3-19)YN (t)CN (t 1)为均衡器在t时刻的输出，也就是I?(t)YN(t)CN (t 1)式 (3-20)而eN(t,t1) I(t)I?(t)eN(t)式 (3-21)为期望信号与估计信号之间的误差。因此，CN (t) 可以根据下式来递推更新CN(t)CN(t 1)KN(t)eN(t)式 (3-22)式 ( 3-22)表明：

9、 t时刻最佳的CN(t) 值可由 t 1时刻的最佳CN (t 1) 值加一修正量得到。这就是递推最小二乘算法或Kalman 算法。将上述在推导过程中出现的各式予以整理，可得到正规RLS 算法的计算步骤。由于此算法为迭代型，故应在已得迭代式组外，还注意在计算的初始部分设置合理的初始值组。根据经验设定则一般可得到较快的收敛效果。由于矩阵RN (t) 类似于统计自相关矩阵，而向量DN (t)近似于互相关向量。应该注意到RN(t)不是一个Toeplitz 矩阵，对于较小的t，RN(t) 可能处于病态条件；因而通常初始时需要在RN(t)上加上一个IN， 为一个1 的正常数，IN 为单位阵。由于对于过去的

10、信号引入了指数权，加上I N 的作用将随着时间增加而减弱。正规 RLS 算法的计算步骤如下：步骤 1： 初始化：令 CN (0) YN(0) 0，式 (3-23)PN (t)I N ( 一般取1 的正数)， n 0式 (3-24)步骤2：更新t t 1eN (t)I (t) YN (t)CN (t 1)式 (3-25)N (t)YN (t)PN (t 1)YN* (t)式 (3-26) 加 (减 )法 2.5N2 1.5N 次。KN(t)PN (t 1)YN (n)w N (t)式 (3-27)我们看到均衡器系数随时间的变化量等于预测误差乘以Kalman 增益向量。由于 KN (t)为 N 维

11、，KN (t)的每一个元素有效地控制着均衡器每一个系数，因而能够得到快速的收敛性质。相反，最陡梯度算法( steepest-descent algorithm)均衡器系数的更新可表示为CN (t) CN(t 1)YN*(t)eN (t)式 (3-30)唯一变化的参数为步长。图 3.1 给出了这两个算法初始收敛速度的比较，信道选自3,具有固定参数f0 0.26，f1 0.93， f2 0.26。信道的特征值为max/ min 11 。均衡器的所有系数在初始迭代时置为0。最陡梯度算法的步长选为0.020。与最陡梯度算法相比 RLS 算法具有较快的跟踪性能和收敛性能。这对于时变信道来说极为重要。例如

12、，短波( HF)信道变化非常快，用梯度算法无法对信号进行均衡。而Kalman算法就能够足够快地跟踪这种变化。3.1 Kalman 算法与梯度算法性能比较3.3 几种改进型快速跟踪的RLS 算法指数遗忘的加窗RLS 算法和Reset-RLS算法RLS 算法广泛的应用于自适应滤波，系统辨识与信号预测。该算法只有在方程误差为0 均值的高斯白噪声以及系统模型非时变时才能保证渐进趋于真值。该算法的另一个显著特点是，为了减小预测中的噪声影响，当参数慢慢趋向于真值时，增益向量便接近于0。因此，RLS 算法就有可能跟踪不上信道参数的变化。为了解决这一问题，在实践当中，人们提出了许多改进的RLS 算法。例如指数

13、遗忘的加窗RLS 算法，避免了增益向量变成0。这一算法的优点是它对于信道参数的变化总是能够起到预防的作用；然而也因为非0 的增益向量使得该算法对信道的扰动和噪声都非常敏感。另外一个方法是一旦检测到信道的变化，就重新初始化迭代协方差矩阵P(t)，如何检测信道参数的变化就成为该算法的关键。在实际操作中，我们可以通过设置适当的门限来检测信道的突跳。一旦迭代误差超过该门限，RLS 算法便被重新初始化。我们称此方法为复位RLS( Reset RLS) 算法。SPRLS算法在文献 4 中 Park 和 Jun 提出了 SPRLS 算法 ( Self Perturbing Least SquaresAlgo

14、rithm ) 。它是一种基于前向预测误差来调整RLS 算法的迭代矩阵，前向预测误差越大，意味着信道发生变化的可能性越大。我们再次给出正规RLS 算法：?(t)?(t 1) k(t)e(t)式 (3-31)e(t)y(t) ?T (t 1) (t)式 (3-32)P(t 1) (t)k(t)PT (t 1) (t)式 (3-33)1 T (t)P(t 1) (t)P(t) P(t 1) k(t) T(t)P(t 1)式 (3-34)P(0) I /式 (3-35)?(0) 0式 (3-36)其中， t为离散采样时刻，?(t) 要预测的滤波器权向量，y(t)是期望的输出信T号， y(t) (t)

15、 n(t) ， 为信道的真值向量，n(t) 为量测噪声，e(t)为前向预测误差，k (t)为自适应增益向量。SPRLS 算法在 (3-34) P(t)的迭代式中加入一项与预测误差有关的项，即P(t) P(t 1) k(t) T(t)P(t 1) N I NTe2(t 1) I式 (3-37)其中 为常数，e(t)为前向预测误差，定义为e(t) y(t) ?T (t) (t)式 (3-38)NINT 函数为四舍五入取整函数，为灵敏增益系数，根据系统量测噪声的大小调整其大小，I 为单位阵。在迭代矩阵中加上一项与前向预测误差有关的项，我们称之为自扰动项。当信道发生变化时, 迭代矩阵受到前向预测误差的

16、扰动，变成非0 从而跟踪信道的变化。该误差平方在SPRLS算法中起着关键的作用，当e2 (t) 0.5时，这一扰动项为0，因而控制着启动该扰动项的最小误差限。当信道参数发生突跳时，该算法比其它普通的RLS 算法性能要好。这一算法在信道参数变化很大而噪声很小时确实很有效，然而当噪声电平很高时上述算法几乎不能工作，原因在于信道参数的变化以及噪声的影响都可能导致较大的前向预测误差。针对上述问题，J.Jiang和R.Cook在文献5中提出了一个新的方案,也就是我们下面将要介绍的MRLS 算法。MRLS 算法在 MRLS 算法当中，迭代矩阵的更新基于一个观测向量与预测向量相关矩阵的函数。由于相关性，这种

17、算法不仅在噪声存在的情况下非常稳健，而且能够快速的跟踪信道参数的变化。正规 RLS 算法的迭代公式已经给出。我们可以看出当迭代时间很长时，增益向量 k(t)已经趋于0，即使信道参数已经改变，(t) 和 (t 1) 也非常接近，也就是说 RLS算法失去了跟踪能力。因此J.Jiang和 R.Cook提出只要信道参数变化就在 P(t) 的迭代式中加上一项正的对角阵，即P(t) P(t 1) k(t) (t)P(t 1) NINT | fy,y?, n,M,(t 1)| I 式 (3-39)1 t1 21 t12f y, y?, n ,M ,(t 1)y (i)y(i)y?(i) n式 (3-40)M

18、 it1MM it1M、 定义如同SPRLS算法。M 代表了运行窗口的大小，它取决于噪声的强度。一旦该扰动项被激活，它便一直起作用直到f 函数的返回值接近于0。从式(3-40)可以看出，f 函数是包括y(t) 的自相关函数，y(t) 与 y?(t) 的互相关函数，以及噪声方差的代数和。由于n(t) 是高斯白噪声，与(t) 不相关，很容易看出， 不管噪声强度多大，当预测参数已经收敛到真值时，一个相当大的M 就使得 f函数的返回值为0。当然，M 的选择应当正比于噪声的方差。另一方面，在任何噪声电平下，当?不等于*时， f 函数的值非0。也就是说f 函数对于信道参数的变化很敏感而对于噪声不敏感。由于

19、采用了运行窗口M ，该算法对于信道参数变化的检测有一个时延。M 越大，收敛的速度越慢。事实上这一算法也是在检测时延和抗噪声之间的一种折衷。结果表明，在噪声很大的情况下，这种改进的RLS 算法可以很好的跟踪时变信道。因为该算法要调整M 的大小，为方便起见我们称这种改进算法为MRLS 算法。ISPRLS 算法虽然 MRLS 算法的性能要优于SPRLS算法，但其灵敏参数以及 M 窗口的大小要根据噪声的强度进行选择，仍然很不方便。于是Kwang Seop Eom等在文献 6中提出另一种基于Kalman 滤波的改进型SPRLS 算法。我们称之为ISPRLS算法 ( Improved Self-Pertu

20、rbing Recursive Least Squares Algorithm) 。 它不需要在不同的噪声条件下改变灵敏参数及窗口大小，能够很好的克服噪声的影响，快速的跟踪信道。我们同样假定P(t)为 N N 协方差矩阵，(t)为输入向量，y(t)为系统输出向量，?(t) 为滤波器的权向量，e(t)为前向预测误差，Q(t)为一自扰动项。ISPRLS 算法如下：式 (3-41)k(t) P(t 1) (t)1 T (t)P(t 1) (t)e(t) y(t)?(t1) (t)式 (3-42)?(t)?(t 1) k(t)e(t)式 (3-43)P(t)Ik(t) T(t)P(t 1) Q(t)式

21、 (3-44)Q(t) 对于决定该算法的性能是否良好是一个关键的因素。因此Q(t) 的设计是主要要解决问题。为了避免根据噪声的大小来确定参数，该算法基于Kalman 滤波的原理设计出一个新的自扰动项。假定在某一时刻t m，系统的真值参数从(0)0变到 (t) m，时变系统可以表示为(t) (t 1) m (t) (t m)式 (3-45)y(t) T (t) (t) n(t)式 (3-46)其中 (n m)为 Kronec ker delta 函数。 n(t) 为方差为2的量测噪声。由Kalman滤波原理可以得到式(3-44)中的Q(t) 为22Q(t) N I NT | Ee (t)2| I

22、式 (3-47)2其中 Ee2(t) 可以由 e2 (t)代替，e2 (t)e2 (t 1) (1)e2 (t)式 (3-48)为实常数，01 。在式(3-47)和(3-48)中，和 不需要根据噪声的大小进行选择。因为分母上的2 会自动根据噪声的大小来调整Q(t) 。前面介绍的四种改进型的RLS 算法应用的背景均为信道估值( 辨识 ) ，鉴于自适应辨识与自适应均衡都属于自适应滤波系统，我们也期望这些改进型的RLS 算 法应用于均衡器时可以获得较好的性能。在下面一节，我们将给出各种算法的计 算机仿真以及结果分析。3.4 计算机模拟及结果分析本节通过计算机模拟来比较各种RLS 算法的跟踪能力和收敛

23、性能。在上述SPRLS、 MRLS、 ISPRLS 三种具有相同思想的改进型RLS 算法中，由于SPRLS对于噪声过于敏感，相比MRLS、 ISPRLS算法性能较差，故我们选取正规RLS，Reset-RLS， MRLS， ISPRLS这几种算法分别用在信道估值和信道均衡中，并进行 了比较。 3.4.1 计算机仿真仿真 1，我们选择文献4中的例子，在不同的信噪比下比较Reset-RLS 和ISPRLS的均方误差性能。假定未知信道的参数在第70次迭代时由0.1 0.2 0.3 0.40.5 0.6 0.7 0.8 0.9 1.0 变为0.8 0.9 1.3 1.1 -0.7 0.1 0.9 1.5

24、 -0.2 -0.3， Reset-RLS在信道变化时重新启动迭代，而ISPRLS在不同信噪比情况下参数不用改变。3.2 均方误差曲线3.3 均方误差曲线第三章递归最小二乘算法 递归最小二乘算法 可以看出，Reset-RLS与 ISPRLS对于突变信道的跟踪能力大体相当。但这时Reset-RLS需要对于信道参数的变化实时地检测，门限的选取便成为该算法的关键而 ISPRLS算法只是在迭代矩阵中加进一项扰动项，该扰动项会对信道参数仿真2，我们选择正规RLS， MRLS， ISPRLS 三种算法进行比较，信道条件1，仍在第70 次迭代时发生突变。3.4 参数误差曲线3.5 参数误差曲线图中给出的是信

25、道参数估计的误差曲线，参数预测误差被定义为*? ，我们取2范数，即xx12 x22xn2 。可以看出，在第70 次迭代信道发生变化时， MRLS 算法和 ISPRLS算法可以迅速的跟上信道的变化，具有较快的收敛能力。 MRLS 和 ISPRLS 算法相比正规RLS 对于信道的突跳具有快速的跟踪能力和收敛性能。MRLS 和 ISPRLS相比，由于MRLS 的运行窗口M 的限制，在刚开始迭代时，它对于信道的跟踪有一段时延。而且在信噪比较低时( 10dB、 0dB)，由于在 MRLS 和 ISPRLS 算法中多加的扰动项受到噪声的影响，在信道突跳前的一段时间，两者收敛时的参数误差要大于正规RLS 算

26、法。因而较好的选择应该将正规RLS 与 MRLS 或 ISPRLS结合起来，在信道不变时，用正规RLS 就能达到较好的收敛效果，一旦信道发生变化，启动MRLS 或 ISPRLS快速的跟踪上信道的变化，可见 MRLS 和 ISPRLS 算法对于时变信道快速的跟踪能力。我们将它们用在均衡系统中，下面的仿真给出了算法的模拟结果。仿真 3，将 MRLS， ISPRLS 算法应用于信道均衡，采用线性均衡器，信号为2PSK， 信 道 采 用 复 值 信 道 ， 在 第 70 次 迭 代 时 信 道 参 数 从 1+0j 0+0j 0.5*(0.996+0.087j) 0+0j 0.3*(0.985+0.1

27、74j) 变 为 1+0j 0.3*(0.996+0.087j) 0+0j 0.2*(0.985+0.174j) 0+0j。3.6 均方误差曲线图 3.7 均方误差曲线从曲线可以看出，在高信噪比情况下( 30dB、 20dB) ， ISPRLS 这种改进型算法，在信道发生突跳时具有较快的跟踪能力和较好的收敛性能，MRLS 算法略优于 RLS 算法， 但二者的跟踪能力明显不如ISPRLS算法。 在低信噪比情况下( 10dB、0dB), MRLS 和 ISPRLS算法却显示不出它们对于正规RLS 算法性能的改善。信道的变化早已淹没在很强的噪声之中，区分不开。这与在仿真2 当中 MRLS 和 ISP

28、RLS算法的优良性能似乎还有一定的差距。我们期望的改进型RLS 算法用与均衡系统时，在低信噪比下仍不具有较好的性能。在下一节中我们将从自适应滤波的角度来分析问题的实质。3.5 信道辨识与信道均衡相同的算法用于信道辨识与信道均衡时，却有着不同的效果。这使得我们不得不探讨一下，同属于自适应滤波系统的信道辨识和信道均衡到底存在着什麽区别？自适应模拟系统和自适应逆模拟系统 自适应滤波系统分为自适应模拟系统和自适应逆模拟系统。自适应模拟系统与辨识可以用一个自适应系统模拟一个未知的、可以随时间慢变的系统。自适应逆模拟则可消除信号在器件和煤质中传输所受到的影响。信道辨识属于自适应模拟系统，而信道均衡属于自适

29、应逆模拟系统。P(z)P(z)rks rkrkP(z)3.8 自适应模拟系统P(z)rk一个单输入单输出的未知系统( 或称被控系统)r的自适应模拟示于图k3.8，被控系统与自适应滤波器有相同的输入激励。自适应滤波器调整自身以得到一个与未知系统相匹配的输出，通常是得到一个未知系统输出最好的最小均方拟合。这种拟合的程度与自适应系统的可调权系数( 即“自由度”) 有关。如图中所示，假定未知系统的传输函数为P( z)，自适应滤波器达到稳态后的传输函数为H (z)，信号的相关函数和互相关函数定义为*xx(n)Exkxk n式 (3-49)式 (3-50)xd (n) Ex*kdk n上两式的Z变换 xx

30、(z)及 xd (z)分别表示信号的功率谱和互功率谱。若自适应滤波器是一个横向滤波器结构，则由维纳滤波理论，输出均方误差可表示为E| d2 | WH RW 2R ePTW*xx(0)wlwm xx(l m) 2 wl dx( l)式 (3-51)lml其中 W 为滤波器的权向量，R EX*XT式 (3-52)*P Ed X式 (3-53)令上式对权的导数为零，可得最小均方误差时的权向量Wopt ， 对每个权分量，则有2wm xx(k m) 2 xd (k) 0式 (3-54)wk m于是wopt xx(k l) xd (k),k式 (3-55)l上式左边是一个卷积的形式，经Z 变换后成为两部分

31、的乘积，于是，可得到最佳权向量的Z 变换opt式 (3-56)xd(z)xx (z)因此， 最佳权的Z 变换是信号x和 d之间的互功率谱与自适应模拟器的输入x的功率谱之比。由于xd (z) P(z) xx(z) xn(z)式 (3-57)若假定系统噪声nk与输入xk相互独立，则上式右边第二项为零，由此知式 (3-58)Hopt (z) P(z)式 (3-58)表明，当自适应模拟系统具有足够的“自由度”去匹配未知系统的输出或输入时，不可能同时匹配系统噪声nk 。事实上，内部系统噪声表现在系统输出一般可看成是一个加性噪声，并认为该噪声与系统本身的输出是不相关的。若自适应模型为线性组合器，并且它的权

32、值已调到使均方误差达到最小，则其最小均方误差解将主要由被控系统的冲激响应所决定，而不是受噪声存在的影响。但自适应收敛过程还是会受到噪声的影响。另外，系统模拟和辨识还与输入信号谱或统计特性有关，一般要求输入信号谱足够宽或统计相关性足够小，才能得到好的模拟和辨识结果。xknkP(z) s dk rkP(z) xssP(z)3.9 自适应逆模拟系统 rkP(z)P(z)rkrkrk图 3.9是一个自适应系统用作逆模拟的情况。自适应滤波器的输入为未知被控系统的输出，而它的输出将为未知系统输入的最小均方组合。由图可知，逆模型的输入功率谱( z在单位圆上)为xx(z) |P(z)|2ss(z) nn(z)

33、, z ejw式 (3-59)同样假定系统噪声nk 和输入信号sk 独立，由于*xd (z)dx(z)P (z) ss(z)nd (z)式 (3-60)式 (3-61)P* (z) ss(z)zejw因此，逆模拟的最佳传输函数xd (z)P* (z) ss(z)H opt (z)2opt xx(z)| P(z) |2 ss(z) nn (z)若系统内部噪声为零，则上式为H o p (tz)1P(z)式 (3-62)即此时逆模型的传输函数为被控系统传输函数的倒数。事实上， 也可以证明式(3-61)和式 (3-62)对所有的z都成立。从式(3-61)看出，自适应逆模拟和自适应模拟不同，系统噪声nk

34、和被控系统输出信号一起作为自适应滤波器的权输入，因而对权值的最小均方解产生影响，同时也影响自适应收敛的情况。在图 3.9的整个逆模拟系统中，若被控系统为一个全零点 FIR 滤波器，则自适应滤波器模型将是一个全极点的递归滤波器，由此，存在自适应滤波器的稳定性和收敛性问题。其所有极点均应位于z 平面的单位圆之内。若考虑自适应滤波器的因果可实现性，通常可在期待响应dk 通道中插入一定的时间延迟。自适应滤波算法的跟踪性能和收敛性能对于自适应滤波算法而言，跟踪性能和收敛性能是两个很重要的考察指标。跟踪是一个稳态过程。相比之下，收敛是一个暂态过程。因而要训练一个自适应滤波器的跟踪能力，就必须从暂态阶段过渡

35、到稳态阶段，这就需要不断的调整滤波器的参数。而且，收敛的速度与跟踪能力是算法的两个不同的性质。一个具有好的收敛性质的自适应滤波算法不一定具有快速的跟踪能力，相反具有快速跟踪能力的算法也不一定具有好的收敛性能。信道辨识和信道均衡都是自适应滤波系统，但跟踪性能却是一个特殊问题。自适应滤波算法用于系统辨识与用于信道均衡或者加性噪声中的信号恢复却有很大的不同。实际中动态的系统可能来自两种情况。参考信号可能是时变的。譬如， 当自适应横向滤波器用于动态系统的系统辨识时，这种情况就会出现。此时，自适应横向滤波器的输入向量的自相关矩阵保持不变，而输入向量与参考向量的互相关矩阵却是时变形式。自适应滤波器的输入随

36、机过程是动态的。这种情况发生在自适应横向滤波器用于均衡一个时变信道时。此时，自适应横向滤波器的输入向量的自相关矩3.11 信道均衡跟踪性能的数学分析比起系统辨识问题而言要复杂得多。因而，一个时变系统的跟踪性能，不仅取决于所采用的滤波器类型，而且也 是个较特殊问题。 3.5.3 自适应均衡器的均方误差下面让我们来分析一下无限长度线性均衡器的均方误差。我们知道在信号为复值的情况下，均方误差函数被定义为J E |ek |2E | IkI?k |2。J为均衡器系数的二次函数。在第二章中我们曾经给出，基于最小均方误差准则的均衡器，按照正交性原理，应有K*E I kcjrk j rj* l0,l式 (3-

37、63)jK即就是*cj E(rk jrj l) E(Ikrk l),l式 (3-64)j将L*E(rk j rk l ) fn fn l j N0 ljn0式 (3-65)xl j N0 lj ,|l j | L0,其它*E(I krk l)f *l , L l 00,其它式 (3-66)将式 (3-65)和 (3-66)代入(3-64) 并对方程两边取z变换得到式 (3-67)11C(z)F(z)F*(z) N0 F*(z1)MMSE 准则的均衡器传递函数为C(z)F*(z1)*1F(z)F * (z 1) N 0式 (3-68)噪 声 白 化 滤 波 器 被 包 含 在 C(z) 里 时 ， 我 们 可 以 得 到 均 衡 器 的 传 输 函 数C (z)1F(z)F* (z 1) N0式 (3-69)minE(ekIk)E|Ik | j zX (z) N 0 cjE(rk jI*k)j式 (3-69)cjf jj我们注意到(3-69)中的求和为cj和 fj的卷积在0点的值。 因而， 如果bk 定义这两个序列的卷积，那末(3-69)的求和就等于b0。序列bk 的 z变化等于B(z) C(z)F(z)F(z)F (z X(z) dz)F(z)F (z 1) N式 (3-70)X(z)X(z) N0b0 为1 b02j zB(z)dz式 (3-71)通过变量