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

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

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

3、则 I?(t) 变为I?(t)N1c j (t1) y(tj )j 0C N (t1)YN (t )式(3-2)其中 CN (t1) 和 YN (t ) 分别为均衡器系数c j (t1)， j0,1, N1 和输入信号y(t j ) ， j0,1, , N1的列向量。类似的，在 DFE 均衡器结构中， 均衡器系数 c j(t ) ， j0,1, N1的前 K1 1个系数为前向滤波器系数，剩下的K 2NK 11为反馈滤波器系数。用来预测I?(t) 的数据为 rt K1 , , rtj ,1jK 2 为判决器先前作出判决的数1 , I t1 , I tK 2 ，其中 I t据。这里，我们忽略判决器

4、判错的情况，因而I t j ,1jK 2 。同时为方便起I tj见定义第三章递归最小二乘算法19vt K1j (0jK 1 )y(t j )式 (3-3)I t K1j (K1jN 1)因此YN (t) y(t), y(t1), , y(tN 1) rt K 1 , rt 1 , rt , I t 1 , I t K 2 式 (3-4) 3.2 RLS 自适应算法RLS 算法对于 I?(t) 的估计可以从下面的式子得到。假定我们的观测向量为YN (n) ， n 0,1, t ，我们期望得到均衡器的系数向量C N (t) 使得均方误差的加权平方和tw t n | eN (n, t) |2(n)式

5、 (3-5)n0最小。其中误差定义为eN (n,t )I (n) C N (t)YN (n)式(3-6)w 代表遗忘因子， 0 w 1。这样我们对过去的数据引入了一个指数权， 这对于信道特性为时变的情况非常合适。关于权向量 CN (t ) 的(n) 最小化便得到下面的线性方程RN (t )CN (t ) DN (t)式 (3-7)其中 RN (t) 为信号的自相关矩阵，定义为tw t nY * N (n)Yn (n)RN (t )式 (3-8)n0D N (t ) 为互相关向量tD N (t )wt n I ( n)Y * N (n)式 (3-9)n 0式(3-7)的解为著名的 WienerH

6、opf 方程CN (t ) RN1 (t )DN (t )式 (3-10)为了避免复杂的求逆运算，引入一 NN 矩阵20P(t)CN1 (t)式(3-11)N由式 (3-8)有tw t n YN* (n)YN ( n)RN (t)n 0wRN (t1)YN* (t )YN (t)式 (3-12)又由矩阵求逆引理有：111RN1 (t 1)YN* (t )YN (t )RN1 (t 1)式 (3-13)RN (t )wRN(t 1)wYN (t )RN1 (t1)YN*(t)在上式中定义 PN (t)RN1 (t ) ，令N(t)Y(t )P(t1)Y*(n)式 (3-14)NNNK N (t

7、)PN (t1)YN* (t )式 (3-15)wN (t )N (t ) 为一标量， K N (t ) 为一 N 维矢量，称为 Kalman 增益向量。则PN (t )1PN (t1)K N (t)YN (t) PN (t1)式 (3-16)w假定我们在式 (3-16)两边右乘以 YN* (t ) ，PN (t)YN* (t )1 PN (t1)YN* (t)K N (t )YN (t) PN (t1)YN* (t )w1N (t )K N (t )K N (t )N (t ) wwK N (t )式 (3-17)因此， Kalman 增益向量可以被定义为PN() *() 。t YNt由于C

8、N (t ) PN (t )D N (t )D N (t )wDN (t1)I (t)YN* (t)式 (3-18)我们得到C N (t )1 PN (t 1) K N (t )YN (t) PN (t 1) wD N (t 1) I (t )YN* (t )w1 I (t )PN (t 1)YN* (t)PN (t 1)D N (t 1)wK N (t)YN (t) PN (t1)D N (t1)第三章递归最小二乘算法211 I (t )K N (t )YN (t )PN (t 1)YN* (t)wC N (t 1) K N (t) I (t) YN (t)CN (t 1)式 (3-19)Y

9、N (t)CN (t 1) 为均衡器在 t 时刻的输出，也就是I?(t )YN (t)CN (t1)式 (3-20)而eN (t ,t1) I (t )I?(t) eN (t )式 (3-21)为期望信号与估计信号之间的误差。因此，CN (t ) 可以根据下式来递推更新CN (t)CN (t1)K N (t )eN (t )式 (3-22)式( 3-22) 表明：t 时刻最佳的 CN (t) 值可由 t 1时刻的最佳 CN (t1) 值加一修正量得到。这就是递推最小二乘算法或 Kalman 算法。将上述在推导过程中出现的各式予以整理， 可得到正规 RLS 算法的计算步骤。由于此算法为迭代型，故

10、应在已得迭代式组外，还注意在计算的初始部分设置合理的初始值组。根据经验设定则一般可得到较快的收敛效果。由于矩阵RN(t )类似于统计自相关矩阵，而向量DN (t ) 近似于互相关向量。应该注意到RN (t) 不是一个Toeplitz矩阵，对于较小的t ， RN (t) 可能处于病态条件；因而通常初始时需要在RN(t )上加上一个I N ，为一个1 的正常数，I N 为单位阵。由于对于过去的信号引入了指数权，加上I N 的作用将随着时间增加而减弱。正规 RLS 算法的计算步骤如下：步骤 1：初始化：令 CN (0)YN (0)0，式 (3-23)PN (t)I N (一般取1 的正数 )，n 0

11、式(3-24)步骤 2：更新 t t1eN (t )I (t)YN (t)CN (t 1)式 (3-25)N (t )YN (t )PN (t1)YN* (t)式 (3-26)K N (t)PN (t1)YN (n)式 (3-27)wN (t )PN (t)1 PN (t 1)K N (t)YN (t )PN (t1)式 (3-28)wCN (t)CN (t1)K N (t )eN (t )式 (3-29)在一次迭代当中，正规 RLS 算法所需的计算量为乘法3N 25N 次，除法 1 次，22加(减)法 2.5N 21.5N 次。我们看到均衡器系数随时间的变化量等于预测误差乘以Kalman 增

12、益向量。由于 K N (t ) 为 N 维，K N (t) 的每一个元素有效地控制着均衡器每一个系数，因而能够得到快速的收敛性质。相反，最陡梯度算法数的更新可表示为( steepest-descent algorithm) 均衡器系C N(t )C N (t1)YN*(t)eN(t )式 (3-30)唯一变化的参数为步长。图 3.1 给出了这两个算法初始收敛速度的比较，信道选自3, 具有固定参数f 00.26 ，f10.93 ，f 20.26 。信道的特征值为max/min11 。均衡器的所有系数在初始迭代时置为0。最陡梯度算法的步长选为0.020 。与最陡梯度算法相比 RLS 算法具有较快的

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

14、上信道参数的变化。为了解决这一问题，在实践当中，人们提出了许多改进的RLS 算法。例如指数遗忘的加窗RLS 算法，避免了增益向量变成0。这一算法的优点是它对于信道参数的变化总是能够起到预防的作用；然而也因为非0 的增益向量使得该算法对第三章递归最小二乘算法23信道的扰动和噪声都非常敏感。另外一个方法是一旦检测到信道的变化，就重新初始化迭代协方差矩阵 P(t ) ，如何检测信道参数的变化就成为该算法的关键。在实际操作中，我们可以通过设置适当的门限来检测信道的突跳。一旦迭代误差超过该门限， RLS 算法便被重新初始化。我们称此方法为复位 RLS( ResetRLS) 算法。3.3.2 SPRLS

15、算法在文献 4 中 Park 和 Jun 提出了 SPRLS 算法 ( SelfPerturbing Least SquaresAlgorithm ) 。它是一种基于前向预测误差来调整RLS 算法的迭代矩阵，前向预测误差越大，意味着信道发生变化的可能性越大。我们再次给出正规RLS 算法：?(t )?(t1) k(t )e(t )式 (3-31)e(t)y(t )?T (t1)(t )式 (3-32)k(t)P(t 1)(t )式 (3-33)1T (t) P(t1)(t)P(t)P(t1) k (t ) T (t )P(t 1)式 (3-34)P(0)I /式 (3-35)?(0)0式 (3-

16、36)其中， t 为离散采样时刻， ?(t) 要预测的滤波器权向量，y(t) 是期望的输出信号，()* T()() ， * 为信道的真值向量，n(t)为量测噪声，e(t)为前向预y ttn t测误差， k (t ) 为自适应增益向量。SPRLS 算法在 (3-34) P(t ) 的迭代式中加入一项与预测误差有关的项，即P(t)P(t1) k(t ) T (t) P(t 1)N I NTe2 (t1) I式 (3-37)其中为常数， e(t ) 为前向预测误差，定义为e(t )y(t)?T (t) (t)式 (3-38)NINT 函数为四舍五入取整函数， 为灵敏增益系数，根据系统量测噪声的大小调

17、整其大小， I 为单位阵。在迭代矩阵中加上一项与前向预测误差有关的项，我们称之为自扰动项。当信道发生变化时 , 迭代矩阵受到前向预测误差的扰动，变成非 0 从而跟踪信道的变化。该误差平方在 SPRLS 算法中起着关键的作用，当 e2 (t ) 0.5 时，这一扰动项为 0，因而 控制着启动该扰动项的最小误差限。24当信道参数发生突跳时，该算法比其它普通的RLS 算法性能要好。这一算法在信道参数变化很大而噪声很小时确实很有效，然而当噪声电平很高时上述算法几乎不能工作，原因在于信道参数的变化以及噪声的影响都可能导致较大的前向预测误差。针对上述问题， J.Jiang 和 R.Cook 在文献 5中提

18、出了一个新的方案 , 也就是我们下面将要介绍的 MRLS 算法。 3.3.3 MRLS 算法在 MRLS 算法当中，迭代矩阵的更新基于一个观测向量与预测向量相关矩阵的函数。由于相关性，这种算法不仅在噪声存在的情况下非常稳健，而且能够快速的跟踪信道参数的变化。正规 RLS 算法的迭代公式已经给出。我们可以看出当迭代时间很长时，增益向量 k(t) 已经趋于0，即使信道参数已经改变，(t ) 和 (t1) 也非常接近，也就是说 RLS 算法失去了跟踪能力。因此J.Jiang和 R.Cook 提出只要信道参数变化就在 P(t ) 的迭代式中加上一项正的对角阵，即P(t)P(t1)k (t ) (t )

19、P(t1)NINT | f y, y, n , M , (t1)| I式 (3-39)?1t11t 1?n , M ,(t2(i)?2式 (3-40)f y, y,1)yM i t 1y(i ) y(i )nM i t 1 MM、 定义如同 SPRLS 算法。 M 代表了运行窗口的大小，它取决于噪声的强度。一旦该扰动项被激活，它便一直起作用直到f 函数的返回值接近于 0。从式 (3-40)可以看出， f 函数是包括y(t ) 的自相关函数，?y(t) 与 y(t) 的互相关函数，以及噪声方差的代数和。由于n(t)是高斯白噪声，与(t) 不相关，很容易看出，不管噪声强度多大， 当预测参数已经收敛

20、到真值时， 一个相当大的 M 就使得 f 函数的返回值为 0。当然， M 的选择应当正比于噪声的方差。另一方面，在任何噪声电平下，当?不等于* 时，f 函数的值非 0。也就是说 f 函数对于信道参数的变化很敏感而对于噪声不敏感。由于采用了运行窗口M ，该算法对于信道参数变化的检测有一个时延。M 越大，收敛的速度越慢。事实上这一算法也是在检测时延和抗噪声之间的一种折衷。结果表明，在噪声很大的情况下，这种改进的RLS 算法可以很好的跟踪时变信道。因为该算法要调整M 的大小，为方便起见我们称这种改进算法为MRLS 算法。 3.3.4 ISPRLS 算法虽然 MRLS 算法的性能要优于SPRLS 算法

21、，但其灵敏参数以及 M 窗口的大小要根据噪声的强度进行选择，仍然很不方便。于是KwangSeop Eom等在文献第三章递归最小二乘算法256中提出另一种基于Kalman 滤波的改进型 SPRLS 算法。我们称之为ISPRLS算法 ( Improved Self-Perturbing Recursive Least Squares Algorithm) 。它不需要在不同的噪声条件下改变灵敏参数及窗口大小，能够很好的克服噪声的影响，快速的跟踪信道。我们同样假定P(t) 为 NN 协方差矩阵，(t) 为输入向量， y(t ) 为系统输出向量， ?(t) 为滤波器的权向量， e(t ) 为前向预测误差

22、， Q(t ) 为一自扰动项。 ISPRLS 算法如下：k(t)P(t1)(t )式 (3-41)1T (t) P(t1)(t)e(t )y(t)?(t1)(t)式 (3-42)?(t)?(t1) k(t )e(t )式 (3-43)P(t ) Ik(t )T (t ) P(t1) Q(t)式 (3-44)Q(t) 对于决定该算法的性能是否良好是一个关键的因素。因此Q(t) 的设计是主要要解决问题。为了避免根据噪声的大小来确定参数，该算法基于Kalman 滤波的原理设计出一个新的自扰动项。假定在某一时刻 tm ，系统的真值参数从(0)0 变到 (t)m ，时变系统可以表示为(t )(t 1)m

23、(t ) (tm)式 (3-45)y(t )T (t )(t )n(t )式 (3-46)其中 (nm) 为 Kronec kerdelta 函数。 n(t) 为方差为2 的量测噪声。由 Kalman滤波原理可以得到式 (3-44)中的 Q(t ) 为Q (t )NI NTE e2(t )2I|2|式 (3-47)其中 2 () 可以由2(t)代替，E etee2 (t )e2 (t1)(1)e2 (t)式 (3-48)为实常数， 01。在式 (3-47)和 (3-48)中，和 不需要根据噪声的大小进行选择。 因为分母上的2 会自动根据噪声的大小来调整 Q(t ) 。前面介绍的四种改进型的RL

24、S 算法应用的背景均为信道估值( 辨识 ) ，鉴于自适应辨识与自适应均衡都属于自适应滤波系统，我们也期望这些改进型的RLS 算26法应用于均衡器时可以获得较好的性能。在下面一节，我们将给出各种算法的计算机仿真以及结果分析。3.4 计算机模拟及结果分析本节通过计算机模拟来比较各种RLS 算法的跟踪能力和收敛性能。在上述SPRLS、 MRLS 、 ISPRLS 三种具有相同思想的改进型 RLS 算法中，由于 SPRLS 对于噪声过于敏感，相比 MRLS 、ISPRLS 算法性能较差，故我们选取正规 RLS，Reset-RLS，MRLS ，ISPRLS 这几种算法分别用在信道估值和信道均衡中，并进行

25、了比较。 3.4.1 计算机仿真仿真 1，我们选择文献 4中的例子，在不同的信噪比下比较 Reset-RLS 和 ISPRLS 的均方误差性能。假定未知信道的参数在第 70 次迭代时由 0.1 0.2 0.3 0.4 0.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 -0.2 -0.3，Reset-RLS在信道变化时重新启动迭代，而ISPRLS 在不同信噪比情况下参数不用改变。图 3.2均方误差曲线图 3.3均方误差曲线第三章递归最小二乘算法27可以看出， Reset-RLS与 ISPRLS 对于突变信道的跟踪能力大体相当。

26、但这时Reset-RLS需要对于信道参数的变化实时地检测， 门限的选取便成为该算法的关键因素。而 ISPRLS 算法只是在迭代矩阵中加进一项扰动项， 该扰动项会对信道参数的变化作出反应，而无需另外进行专门的检测。仿真 2，我们选择正规RLS， MRLS ，ISPRLS 三种算法进行比较，信道条件同仿真 1，仍在第 70 次迭代时发生突变。图 3.4 参数误差曲线图 3.5 参数误差曲线28图中给出的是信道参数估计的误差曲线，参数预测误差被定义为*? ，我们取 2范数，即 xx12x22xn2。可以看出，在第70 次迭代信道发生变化时，MRLS 算法和 ISPRLS 算法可以迅速的跟上信道的变化

27、， 具有较快的收敛能力。 MRLS 和 ISPRLS 算法相比正规 RLS 对于信道的突跳具有快速的跟踪能力和收敛性能。 MRLS 和 ISPRLS 相比，由于 MRLS 的运行窗口 M 的限制，在刚开始迭代时，它对于信道的跟踪有一段时延。而且在信噪比较低时( 10dB、0dB) ，由于在 MRLS 和 ISPRLS 算法中多加的扰动项受到噪声的影响，在信道突跳前的一段时间，两者收敛时的参数误差要大于正规RLS 算法。因而较好的选择应该将正规RLS 与 MRLS 或 ISPRLS 结合起来，在信道不变时，用正规RLS 就能达到较好的收敛效果，一旦信道发生变化，启动MRLS 或 ISPRLS 快

28、速的跟踪上信道的变化，可见 MRLS 和 ISPRLS 算法对于时变信道快速的跟踪能力。我们将它们用在均衡系统中，下面的仿真给出了算法的模拟结果。仿真 3，将 MRLS ，ISPRLS 算法应用于信道均衡，采用线性均衡器，信号为2PSK，信道采用复值信道，在第 70次 迭 代 时 信 道 参 数 从 1+0j0+0j0.5*(0.996+0.087j)0+0j 0.3*(0.985+0.174j)变 为 1+0j 0.3*(0.996+0.087j)0+0j0.2*(0.985+0.174j)0+0j 。图 3.6 均方误差曲线第三章递归最小二乘算法29图 3.7 均方误差曲线从曲线可以看出，

29、在高信噪比情况下( 30dB、 20dB) ，ISPRLS 这种改进型算法，在信道发生突跳时具有较快的跟踪能力和较好的收敛性能，MRLS 算法略优于 RLS 算法，但二者的跟踪能力明显不如ISPRLS 算法。在低信噪比情况下 ( 10dB、0dB), MRLS 和 ISPRLS 算法却显示不出它们对于正规 RLS 算法性能的改善。信道的变化早已淹没在很强的噪声之中， 区分不开。这与在仿真 2 当中 MRLS 和 ISPRLS算法的优良性能似乎还有一定的差距。我们期望的改进型RLS算法用与均衡系统时，在低信噪比下仍不具有较好的性能。在下一节中我们将从自适应滤波的角度来分析问题的实质。3.5 信道

30、辨识与信道均衡相同的算法用于信道辨识与信道均衡时，却有着不同的效果。这使得我们不得不探讨一下，同属于自适应滤波系统的信道辨识和信道均衡到底存在着什麽区别？ 3.5.1 自适应模拟系统和自适应逆模拟系统自适应滤波系统分为自适应模拟系统和自适应逆模拟系统。自适应模拟系统与辨识可以用一个自适应系统模拟一个未知的、可以随时间慢变的系统。自适应逆模拟则可消除信号在器件和煤质中传输所受到的影响。信道辨识属于自适应模30拟系统，而信道均衡属于自适应逆模拟系统。P( z)噪声nkd k未知系统r rkek输入xxkH ( z)ykssx自适应模拟sP( z)xP( z)P( z)sP( z) rks rk r

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

32、x ( z) 及 xd (z) 分别表示信号的功率谱和互功率谱。若自适应滤波器是一个横向滤波器结构， 则由维纳滤波理论， 输出均方误差可表示为E| d 2 |W H RW2R ePT W xx (0)wl* wm xx (lm) 2 wl dx ( l )式 (3-51)lml其中 W 为滤波器的权向量，REX*XT式 (3-52)P E d * X 式 (3-53)令上式对权的导数为零， 可得最小均方误差时的权向量Wopt ，对每个权分量， 则有第三章递归最小二乘算法312wm xx (k m) 2 xd (k) 0式 (3-54)wkm于是woptxx (kl )xd (k ),k式 (3

33、-55)l上式左边是一个卷积的形式，经Z 变换后成为两部分的乘积，于是，可得到最佳权向量的 Z 变换H opt ( z)xd ( z)式 (3-56)xx ( z)因此，最佳权的 Z 变换是信号 x 和 d 之间的互功率谱与自适应模拟器的输入x 的功率谱之比。由于xd ( z)P(z) xx ( z)xn ( z)式 (3-57)若假定系统噪声 nk 与输入 xk 相互独立，则上式右边第二项为零，由此知H opt (z)P( z)式 (3-58)即自适应模拟器的传输函数就是未知系统的传输函数。式 (3-58)表明，当自适应模拟系统具有足够的“自由度”去匹配未知系统的输出或输入时，不可能同时匹配

34、系统噪声 nk 。事实上，内部系统噪声表现在系统输出一般可看成是一个加性噪声，并认为该噪声与系统本身的输出是不相关的。若自适应模型为线性组合器，并且它的权值已调到使均方误差达到最小，则其最小均方误差解将主要由被控系统的冲激响应所决定，而不是受噪声存在的影响。但自适应收敛过程还是会受到噪声的影响。另外，系统模拟和辨识还与输入信号谱或统计特性有关，一般要求输入信号谱足够宽或统计相关性足够小，才能得到好的模拟和辨识结果。P( z)噪声H (z)nk输 入逆模型被控系统P(z) rksskxk rkP( z)ekP( z)sykdk rk rkxP( z) xx 图 3.9 rks自适应逆模拟系统ss

35、P(z)P( z)P( z) rk rk rk32图 3.9 是一个自适应系统用作逆模拟的情况。 自适应滤波器的输入为未知被控系统的输出，而它的输出将为未知系统输入的最小均方组合。由图可知，逆模型的输入功率谱 ( z 在单位圆上 ) 为xx ( z)| P( z) |2ss ( z)nn ( z),ze jw同样假定系统噪声 nk 和输入信号 sk独立，由于xd (z)dx* ( z)P* (z)ss( z)nd ( z)P* ( z)ss (z)zejw因此，逆模拟的最佳传输函数H opt ( z)xd (z)P* (z)ss ( z)xx (z)| P( z) |2ss( z)nn ( z

36、)式(3-59)式 (3-60)式 (3-61)若系统内部噪声为零，则上式为H o p (t z)1式 (3-62)P(z)即此时逆模型的传输函数为被控系统传输函数的倒数。 事实上，也可以证明式 (3-61) 和式 (3-62)对所有的 z 都成立。从式 (3-61)看出，自适应逆模拟和自适应模拟不同，系统噪声 nk 和被控系统输出信号一起作为自适应滤波器的权输入，因而对权值的最小均方解产生影响，同时也影响自适应收敛的情况。 在图 3.9 的整个逆模拟系统中， 若被控系统为一个全零点 FIR 滤波器，则自适应滤波器模型将是一个全极点的递归滤波器，由此，存在自适应滤波器的稳定性和收敛性问题。其所

37、有极点均应位于z 平面的单位圆之内。若考虑自适应滤波器的因果可实现性，通常可在期待响应 dk 通道中插入一定的时间延迟。 3.5.2 自适应滤波算法的跟踪性能和收敛性能对于自适应滤波算法而言，跟踪性能和收敛性能是两个很重要的考察指标。跟踪是一个稳态过程。相比之下，收敛是一个暂态过程。因而要训练一个自适应滤波器的跟踪能力，就必须从暂态阶段过渡到稳态阶段，这就需要不断的调整滤波器的参数。而且，收敛的速度与跟踪能力是算法的两个不同的性质。一个具有好的收敛性质的自适应滤波算法不一定具有快速的跟踪能力，相反具有快速跟踪能力的算法也不一定具有好的收敛性能。信道辨识和信道均衡都是自适应滤波系统，但跟踪性能却

38、是一个特殊问题。自适应滤波算法用于系统辨识与用于信道均衡或者加性噪声中的信号恢复却有很大的不同。实际中动态的系统可能来自两种情况。第三章递归最小二乘算法331参考信号可能是时变的。 譬如，当自适应横向滤波器用于动态系统的系统辨识时，这种情况就会出现。此时，自适应横向滤波器的输入向量的自相关矩阵保持不变，而输入向量与参考向量的互相关矩阵却是时变形式。2自适应滤波器的输入随机过程是动态的。这种情况发生在自适应横向滤波器用于均衡一个时变信道时。此时，自适应横向滤波器的输入向量的自相关矩阵以及输入向量与参考向量的互相关矩阵都是时变形式。因而，时变信道均衡的噪声时变信道输入自适应滤波器图 3.10 信道

39、辨识噪声输入时变信道自适应滤波器图 3.11 信道均衡跟踪性能的数学分析比起系统辨识问题而言要复杂得多。因而，一个时变系统的跟踪性能，不仅取决于所采用的滤波器类型，而且也是个较特殊问题。 3.5.3 自适应均衡器的均方误差下面让我们来分析一下无限长度线性均衡器的均方误差。我们知道在信号为复值的情况下，均方误差函数被定义为JE | ek|2E | I kI?k|2 。 J 为均衡器系数的二次函数。在第二章中我们曾经给出，基于最小均方误差准则的均衡器，按照正交性原34理，应有KE I kc j rk j r j*l 0,ljK即就是c j E(rk j r *j l )E( I k rk* l )

40、,lj将LE(rk j rk*l )f n* f n l jN 0 ljn0xl jN0lj ,| l j | L0,其它E(I k rk*l )f *l , Ll00,其它式 (3-63)式(3-64)式 (3-65)式 (3-66)将式 (3-65)和(3-66)代入 (3-64) 并对方程两边取 z 变换得到C(z) F (z) F * ( z 1 )N0 F * (z 1 )式 (3-67)因而基于 MMSE 准则的均衡器传递函数为F * ( z 1 )式 (3-68)C ( z)1) N0F ( z) F * ( z当噪 声白化滤波器 被包 含在 C (z) 里时 ，我 们可 以得到均 衡器 的传 输函数1C ( z)F ( z) F * (z 1 ) N 01式 (3-69)X ( z) N 0我们注意到 C ( z) 与基于最大失真准则的均衡器传递函数式(2-7) 仅在于 (3-69)表达式中分母上的噪声功率谱密度 N 0 。当 N 0 与信号相比非常小的情况下， 使得峰值失真 D (c) 最小的均衡器系数与使得均方误差J 最小的均衡器系数相当。 也就是，当 N00 时，最大失真准则与最小均方误差准则可以得到相