自适应滤波器翻译作业

1、第八章快速横向LMS滤波算法8.1简介在大量的算法解决最小二乘问题递归形式的方法中快速横向递 归最小二乘（FTRLS）算法是非常具有吸引力，因为其能减少计算复杂 度。FTRLS算法可以通过求解同时向前和向后的线性预测问题，连同其他两个横向过滤器：过程估计量和一个辅助滤波器的期望信号向量 有一个作为其第一和唯一的非零元素 （例如,d（0）= 1）。与格型算法相 比,FTRLS算法只需要时间递归方程。然而，需要得到一些FTRLS算 法的关系，可参考前面一章 LRLS算法。FTRLS算法考虑快速的横 向滤波器RLS的算法更新的解决方法。因为顺序固定，更新横向自 适应滤波器系数向量在每个计算中都迭代。

2、格型算法的向后和向前的派生关系可以用于预测所派生的FTRLS算法。由此产生的算法计算复杂度在实际中实现N使它们特别具有吸引力。相比格型算法，FTRLS算法的计算复杂度较低，由 于没有权向量更新方程。特别是，FTRLS算法通常需要7 n到11 n 每输出样本，乘法和除法则需要LRLS 14n到29 n计算。因此，FTRLS 算法被认为是最快的解决方案的实现 RLS的问题1-7。在工程实践领域相继提出几种不同的 FTRLS算法，所谓的快速 卡尔曼算法1,这的确是一个早期的快速横向 RLS算法,计算11n次乘 法和除法的复杂运算在每次输出示例。在后面的研究阶段开发领域的 快速横向算法，快速后验误差序

3、列的技术 (fa)2,快速横向滤波器 (FTF)3算法提出了要求，同样需要7n乘法和每次除法的输出样本。 FTF算法是具有最低的复杂性的 RLS算法，不幸的是，这些算法对量 子化效应非常敏感，如果有一些步骤没被采取将会变得不稳定。在这一章,FTRLS算法的一种特殊形式将被提到，基于那些被提 的网格算法所派生出来的。众所周知，量子化错误在FTRLS算法 中是指数发散1-7。自从FTRLS算法不稳定的行为用有限精度算 法实现的时候，我们讨论实现FTRLS数值稳定的算法，并提供一个特 定算法的描述8,10。8.2递归最小二乘预测快速算法探索一些结构性的信息数据以达到低计算的复杂性。 在特定情况下的快

4、速 RLS算法本文中讨论达到减少计算复杂度的情 况下，由输入信号连续推迟样本中相同的信号。在本例中，模式的快速算法是相似的，向前和向后预测这些过滤器是必不可少的部分算法。 建模的预测执行任务的输入信号，因此允许替换矩阵方程的矢量和标 量关系。派生的FTRLS算法，解决方案的RLS向前和向后的预测问题需要 权向量递归方程。在本节中，这些解决方案进行了综述强调 FTRLS算 法相关的结果。如前所述，我们将借一些派生的前一章对点阵算法。 是值得的提及,FTRLS可以被介绍通过一个独立的推导，基于格型的推 导在这点可能更加深刻的当然更直截了当的。821向前预测关系瞬时向前后验Nth-order预测作为

5、预测误差打(札 TV) = x(k)-wJ(k.N)x(k-l,N)U&N + l)后验和先验的向前预测误差之间的关系，首次提出了方程(7-49) 和为了方便在这里重复诃N)丁伙一 1,N)一个简单的处理方程(7.73)，导致以下的最小加权最小二乘误差时间的更新，这种方法将用于FTRLS算法:歆应临N)=入氛皿伙一 1,N)十引(层同样的从等式(7.73)我们可以获得，需要的等式在FTRLS算法中对(k - lV)mm应N + 1)=欝仏百於-AMJmin '1*可以通过执行前一章的方程(7.40)提出更新方程预测抽头系数矢 量w/fk, N) = w/(fc 1, TV)十

6、0(上一1, N)ef(k1 TV)在这里c>(k lj TV) = S)(k 一 1, A)x(k 1. A )将会看到，向量的更新© (k-1,N) © (k,N + 1)是需要更新落后的预测系数向量。同时，最后一个元素的© K N + 1)是用于更新反向预测先验误差和获得丫 K N)。向量© K N + 1)可以通过自右乘方程(7.56),双方在即时k和系数N通过x(k, N + 1)=x(k) xT(k-1,N)。结果可以表示为协,N十1)=刑_ 的十号.N) I恥“)然而，不方便使用FTRLS算法因为上面的方程产生反向预测部分，它将导致额

7、外的计算。解决方案是使用另一种递归涉及亠(k N +1)7)k，N 1r代替.k，N 1 (具体参照问题氛k、TV + 1)=后产生的递归可以派生一些代数运算方程(8.6)和(8.3)(8.5),得到- U +-W/(A - lV)可(站N)正向预测抽头系数向量应该被更新使用，这样Wy(A；, N) = Wy(A： 1,N)十 4>(k h N)Ef(k N)8.2反向预测关系在本节中，关系涉及用于FTRLS反向预测问题算法。后验概率预测与先验概率预测误差之间的关系可以表示为£t(fcT N) eb(k, N)j(kN)我们也知道对于不同转换因素的比率表示为丁仇N + l)=兀

8、壯k-l,N)见前一章的方程(7.79)我们为了方便重写了最后的平等方程(7.70)，得到* (k, N) = >4t (fc - 1, N)这个等式也可以这样写肝(屁M怨心依1，N)一现在我们回想一下，反向预测滤波器的更新的时间可以写成wb(fc,AT) = wb(fc- 1,AT) + (k,N)eb(k,N)=轨(k-LTV) + 驱 N)為仏,N)以下类似的方法，得到方程(8.7)，首先两边的方程(7.59)，在即后乘时k和N,通过x(k,N+1)=(k,N) xT x(k-N)，并使用关系(8.10),(8.11),(8.13)我们有褓N)0二0(札冷+1) -皿;1肿卜劇注意

9、,在这个等式的最后一个元素k, N 1已经在方程(8.7)计A算。在任何情况下，值得一提的是，最后一个元素的k, N 1或者可以表达祗+i（匕 N 4- 1）=通过方程(8.9),(8.15),在方程(8.12)和(8.10),我们可以得到+ 兀+i（EN + l血（札N）二了（札 N）丁（机 N + 1 I将方程(8.9)代入上面的方程,我们可以归纳出更新方程，并用于FTRLS算法N) = 了 j(底 N + 1 _ 瓦+息 N + l)eb(k.N)有关后验与先验的预测问题和转换因子 丫 (k,N)的更新方程现在 可用。我们可以通过期望信号d(k)进行派生解决估计的更一般的问题 相关的过程

10、，称为过程评估。8.3过程评估对于所有先前提出了自适应滤波器算法，得到FTRLS算法是很有 用的，可以匹配一个期望信号d(k)的最小化加权方差。从先验误差e(k,N) = d(k) - wr(fc - 1, N)x(k, N)我们可以计算后验误差E(k,N) = e(k,N)y(k,N)在传统的RLS算法，更新的时间输出联合过程的抽头系数估计量可以执行w(A：, N) = w(A： 1, N) + (爲N)e(為 N 'w(A： 1： N) + </>(A 7V)e伙，A I现在所有的更新方程可用来描述快速横向RLS算法。的FRLS 算法由方程(8.1)-(8.3),(8.

11、7)-(8.8)和(8.4)提出相关预测；方程 (8.15),(8.17),(8.9),(8.11),(8.14)和(8.13)相关的预测和落后的东西转换 因子;(8.18)-(8.20)与过程估计量有关。FTRLS算法在逐步形成算法8.1。FTRLS算法的计算复杂度7(N)+ 14乘法/输出示例。FTRLS算 法的关键特性是它不需要矩阵乘法。正因为如此 ,FTRLS算法的实现 每输出样本顺序相乘N的复杂性。初始化过程包括设置反向预测的抽头系数，前进预测和过程评估过滤器为零，即wj( 1?7V) = Wb(1, N) = w( 1, N) = 0向量 -1, n设置0假设的输入和期望信号零 k

12、 < 0,即prewindowed数 据。转换因子应该初始化丁(一1)N) = 1算法8.1快速横向RLS算法科了=叫（1, TV）=叫一1.= 0机一 1,W) =0, 7(-UJ = 1(1. .V) = e (a smaJI positive constant)Prediction PurlDo fort'3th k >0,£/(fe, JV) =- 1, N)银切,(比N) = A£：协” (k-1,N) +引(鼠N冷快N) N) = W/(fc -1,7V) + 0(it - lt N)今住,N)磁 JV + 1)= 0tN)j minj_!P

13、min 住SI 仆n + i)=:肝二加 iz eb N)=人bmin (血1、川)电!+1 (职 * 十】) 邙一】闪 W)=邙 j(匕 N+ 1) 心"+ 1)殆(k,N)(k,N) = eb(k, AT)7(tb jV) f 敗 N) = Ain (fc-l：JV) + Eb(k, N)eb(k,N) '如a0叫(比 N) =-1,N) + 冰时)5(釈 N)，址nin=$(肌 N + 1+1 (札 N + 1)(S.2)(S3)<8.8)(S.7)<8.4)(8.15)(3.17)(8J1)(8.14)(S.13)Joint-Process Estimat

14、ion(S.1S)(8.19)(5.20)e(h,N) = d(k) - wr(fc - m&N)EktN) = e(Jt, Nf(k, N)Ar) = w(fc 1, JV) + 讽息 AL)(t, A End因为在初始化期间先验和后验误差之间没有区别。加权最小二乘 误差应该初始化与一个正的常数。"孟-1川）之爲（-1川）为了避免除零在第一次迭代。引入这个初始化参数的原因表明 ，它 应该是一个小的价值。原因，然而，对于稳定的价值不应小（见本章末尾 的例子）。应该提到，有确切的初始化程序的快速横向 RLS过滤器，目的是最小化目标函数的瞬间在初始化期间3这些程序在初始化期间探索

15、事实数据样本的数量在d(k)和x(k)小于N + 1。因此，目标函数可以是零，因为比需要更多的参数。3的 确切的初始化过程取代了计算密集型回来时相当简单替换算法和自 适应滤波器系数和零初始化。这个过程也可以广义的情况下一些非零 抽头系数的初始值是可用的。正如前面提到的，一些快速RLS算法基于横向实现存在，这里介绍 的一个对应于所谓的在3提出了 FFT。大量的替代算法引入的问题。8.4稳定快速横向 RLS算法尽管速度横向算法提出了文学提供一个不错的解决方案固有的计 算复杂度负担传统的RLS算法，这些算法用有限精度算法实现时不稳 定。增加字并不能解决不稳定的问题。唯一的采用更长的字的效果是, 该算

16、法将不再有分歧。解决这个问题早些时候由重新启动算法选择的 累积误差变量时达到规定的阈值3。虽然过去再启动过程将使用信 息，由此产生的表现不佳是由于不连续的信息在相应的确定性的相关 矩阵。不稳定行为的原因快速横向算法固有的正反馈机制。这个解释了这个想法，如果一些特定的测量数值错误，他们可以方便地反馈为了使 负面反馈误差传播动力学中占主导地位。幸运的是，一些测量的数值错误可以通过引入快速算法计算冗余。 这种计算冗余可能涉及使用两 个不同的公式计算一个给定的数量。在有限精度实现中，结果的数量通过这些公式计算值不相等和他们的区别是一个很好的测量数量的累积误差。这个错误可以反馈为了稳定算法。 关键问题是

17、确定的数量 应该引入计算冗余的误差传播动力学可以稳定。早期提出的解决方案 67,只有一个数量选择引入冗余。之后，这是表明，至少有两个量要求 为了保证稳定的FTRLS算法9。另一个相关的问题是，这个错误应该 内反馈算法。注意，任何时候可以选择在不影响算法的行为实现无限 精度时，自反馈误差为零。自然选择是错误反馈回相关的物理量的表 达式。这意味着对于每个数量，介绍了冗余,其最终价值是计算的两种 形式的组合。FTRLS算法可以看作是一个离散时间非线性动态系统9:有限精度实现中使用时，量化误差将会上升。在这种情况下，内部的数量将 摄动与无限精确数量相比。非线性系统建模误差传播时，可以被描述,如果适当的

18、线性化，允许误差传播机制的研究。使用一个平均分析，这是有意义的固定的输入信号，可以得到一组系统的特点是它的特征值 的动态行为类似于k时的误差传播行为th和(1-入予0。通过这些特 征值,可以确定反馈参数以及数量选择引入冗余。这里的目标是修改 不稳定模式通过错误的反馈以让他们稳定9。幸运的是，它被发现在9,可以修改和稳定不稳定模式引入错误的反馈。不稳定模式可以修 改通过引入冗余丫 KN)和eb(k,N)。这些数量可以计算使用不同的关 系，以便区分它们包含在一个额外的索引他们的描述。先验向后误差可以被描述的替代形式1) = Amin(Ai TV + 1)為(层-uj (k 1, N) 1龙(札 J

19、V + 1)3) = eh(k, N, 2)股 + 為(息 1)1 叫=efr(fc)N, 1) + Kf 瞅爲 JV, 2) et(fc, TV, 1)第一种形式是受雇于FTRLS算法和第二种形式对应的内积实现 先验向后误差。第三形式对应于一个线性组合的两种形式，这些形式反馈确定数值差别的最终值，w(k,N,3),它将使用在不同的地方稳定算 法。对于每个Ki，i=1,2,3中，我们选择一个不同的值，以保证相关特征 值小于1。转换因子丫 KN)可能是第一个参数显示算法变得不稳定的迹象。 这个参数也可以通过不同的计算关系。必须保证所有这些替代关系模式的误差传播系统变得稳定。第一个方程给出 1,

20、N,3) + 札 N + l)e/(fct2V)在 k,N 1里的第一个元素是0 k,N 1。上述等式源于(8.4),(8.3),(8.2)和(8.7)等式以及迭代。第二个表达式的转换因子来源于方程(8.14)和给出的= 7 N + 1丁 1) (札N + 1)寂爲3(船M 3)第三个表达式是b j(眾皿3) = 1 +$(底N)双机TV)在方程(8.27),转换因子是用不同的方式表达，这是第一个提出FTRLS算法9。第二种形式已经使用一个先验落后的错误和冗余。 第三种形式可由方程(7.48)晶格RLS算法(参见问题10)。另一种关系稳定快速横向算法利用涉及到最低最小二乘误差。从方程(8.3

21、)和(8.7),我们可以写魏昭札小尸=H21 mi £ef(k, N)£f(kN) 胶爲/斤一1川)珠J肌N) $0(乩”)今(息¥)从(8-6)我们可以推断出，”北=%典N)=沐匕心匕N + 1,1)带有这个关系，我们可以获得所需的方程=入一惟盒金 _ i,n)t格 n+1,1血 n+1)选择的丫(k,N + 1,1)是用来保持系统错误的工作状态的稳定9。使用方程转换因子和冗余的先验向后误差，我们可以获得稳定的快速横向RLS算法(SFTRLS)逐步实现给定的算法8.2。参数心，i=1,2,3 确定通过计算机模拟搜索9的最佳值发现Kj =1.5, K2 =2.5,

22、 K3 =1。在9 还发现，数值表现对于Ki的最优值毫无反应，选择最佳值对于一个给定 的情况对各种环境和工作算法设置情况(例如，对于不同的遗忘因子的 选择)。SFTRLS算法相关的另一个问题涉及的范围值 入的稳定保证。1?< A < 12(N + 1) 一 忑大量仿真实验结果9表明,该范围这里的N是滤波器的阶数，实验验证最优数值选择当 入的值选择04?vn算法8.2稳定快速横向RLS算法lnitializaliunWy ( 1T 1V) =1, Ar) = lt N) = 00(L/V) = fl. 7(-1,3)=1Kl = 1.5, H2 = 2.5, 3 = 1£

23、la small positive consiant)Prvdiciion Ptirtlk)for each k > 0,与(札 N) = Ej(札 Ny)(k13)*ui陋N + l=冷化_屮)J + 7用7N) L(k,/V + 1,1) = yl(k- ltJV, 3) +X(>+f(ktN) 吩丽临呵L F险曲佻- IN- WN + l孟临N + l) 巧(此 N) = 1 N) + 松1, N打(爲 N)跖他/VJ) =入住(fc -1, JV)jv+i(*j JV +1) 环(爲Nd)=嚣征一1 xxV + 1) 甄此M3)=电&化M浚)怙 +jV, 1)1 -

24、 kJ furi= 1T2,3们啥阴2)=下啥N 41,1) 亦i(爲N 4-1)%忍爲M町 阪M，N、3)=叫血N讹,Nj二1北褊他旳=轧#-1川)+臨冷如決凤3) 呱評卜敞E)-砥曲"+ 1) F叮脚 w&(fc, N) = wb(fc 1, N) + 0( A')Eb,L(A?tA 3) 们(札N=1*&他眄紅人尙(8.2)(8.7)(«.27)(8.30)(88)(815)(8.25)(8.28)(8 11)(8.14(8.138.29cint-Process lisLiraatione(t, N) = d(k) - wT (k - 1, N

25、)xk, N)£(Jt,jV) = e( N)7札M3) fl(fctJV)=稱(k +&&7V)e(&JVtS.18)(8J0)End入值的范围以及它可以非常接近最优值高阶滤波器。这可能是个潜在的限制SFTRLS算法的使用，特别是在不稳定环境中较小值入是必需的。SFTRLS算法的计算复杂度是订单9 n乘法/输出示例。有另一种 算法计算复杂度的订8n(参见问题9)。在离开之前这一节中，值得一提的是快速横向RLS算法的一个很 好的解释。FTRLS算法可以看作四个横向滤波器并行工作，互相交换 数量,如在图8.1。第一个过滤器是远期预测滤波器 利用x(k- 1,N)

26、作 为输入信号矢量，wf(k,N)系数向量，并提供数量£ f(k,N),英孚(k,N), E dfmin(k,N)作为输出。第二个过滤器反向预测滤波器，利用x(k,N)作为输入信号矢量，世行KN)系数向 量，并提供大量£ b(k,N),eb(k,N), Ed bmin(k,N)作为输出。第三个过滤 器是一个辅助滤波器的系数由-? © KN),其输入信号向量x(k,N),和 其输出参数丫 - 1(k,N)。对于这个过滤器，所需的信号矢量常数和等于 (1 0 0。0T。第四个和最后一个过滤器的过程估计量输入信号向量 x(k,N)的系数向量w(k,N),并提供大量的(k

27、,N)和e(k,N)作为输出。x(k)分比川）® (A. W eb (k N)£ (ktN)e(t N)图8.1快速横向RLS算法框图例8.1小节中描述的系统辨识问题顺序稳定快速横向算法来求解这一章。主要目的是检查算法在有限精度实现时的稳定性。解答根据方程（8.31），下界为入在这种情况下是0.9375。选择一个值入=0.99。稳定快速横向算法应用于解决识别问题和测量MSE是 0.0432。使用二2,我们运行的算法与有限的精度和结果总结在表8.1。没有发现不稳定的迹象为入二0.99。这些结果产生的系综平均 200 实验。比较的结果表8.1与表5.2和7.2所示,SFTRLS

28、算法也有类似 的性能比传统和格型RLS算法,量化误差积累。问题是该算法在大多 数情况下保持稳定。对于SFTRLS,对于大订单的过滤器我们剩下的有 限选择入的值的范围。同时，发现在我们的实验中，初始化参数的选择 起着重要的作用在这个算法的性能实现有限的精度。在某些情况下，甚至当入的值是在推荐范围内，该算法不收敛很小。通过增加的价值， 我们通常的收敛时间增加同时该算法稳定。例8.2小节中描述的通道均衡的例子（363）也用于模拟测试SFTRLS算 法。我们使用25阶均衡器和遗忘因子入二0.99。解答为了解决平衡问题稳定快速横向 RLS算法初始化二0.5。这里给出的 结果所产生的系综平均200实验。生

29、成的MSE的学习曲线显示在图 8.2中,测量和MSE是0.2973。的整体性能SFTRLS算法对于这个特 定的例子是一样好RLS算法，如lattice-based算法。200300400herjtionsh k5006008.5结束语在这一章里，我们已经提出了一些快速横向 RLS算法。这类算法 比传统和格型RLS算法计算效率更高。一些模拟例子包括 SFTRLS 算法使用的地方。有限字长模拟是读者的特殊兴趣。大量的替代FTRLS算法以及理论结果可以在3中找到。归一化版本 的FTRLS算法的推导也可能并没有解决 。目前的一章，这一结果参 考4。已知最计算有效FTRLS算法不稳定。允许的错误反馈方法

30、简 要介绍稳定的FTRLS算法。完整的推导并给出理由错误反馈的方法9。在不稳定的环境中，它可能是有用的采用时变遗忘因子。因此希 望得到FTRLS算法允许使用变量入。这个问题第一次被在11。然而 计算提出了更有效的解决方案在8介绍了数据加权的概念来取代误差加权的概念FTRLS算法具有潜在的应用。特别是，信号的问题可以从环境嘈 杂的版本的传输信号和噪音过滤版本相同的传输信号是一个有趣的 应用程序。在这个问题，延迟和未知的滤波器系数估计。加权平方误 差最小化，同时考虑延迟和未知系统参数。这个问题的联合估计可以 解决采用优化FTRLS算法12。8.6引用1. D. D. Falc oner and L

31、. Ljun g, “Applicatio n of fast Kalma n estimati on to adaptive equalization” IEEE Trans. on Communications, vol. COM-26, pp. 1439-1446, Oct. 1978.2. G. Carayannis, D. G. Manolakis, and N. Kalouptsidis, “ A fast sequential algorithm for leastsquaresfiltering and prediction,” IEEE Trans. on Acoust.,

32、Speech, and Signal Processing, vol.ASSP-31, pp. 1394-1402, Dec. 1983.3. J. M. Cioffi and T. Kailath, “ Fast, recursive-least-squares tran sversal filters for adaptive filters, ” IEEE Trans. on Acoust., Speech, and Signal Processi ng, vol. ASSP-32, pp. 304-337, April 1984.4. J. M. Cioffi and T. Kaila

33、th, “Windowed fast transversal filters adaptive algorithms with normalization,” IEEE Trans. on Acoust., Speech, and Sig nal Processi ng, vol. ASSP-33, pp. 607-627,June 1985.5. S. Ljung and L. Ljung, “Error propagation properties of recursive least-squaresadaptation algorithms,” Automatica, vol. 21,

34、pp. 157-167, 1985.6. J.-L. Botto and G. V. Moustakides, “Stabilizing the fast Kalman algorithms,” IEEE Trans. On Acoust., Speech, and Sig nal Process ing, vol. 37, pp. 1342-1348, Sept. 1989.7. M. Bellanger,“ Engineering aspects of fast least squares algorithms intran sversal adaptivefilters, ” Proc.

35、 IEEE In tern. Conf. on Acoust., Speech, Sig nal Process ing, pp. 49.14.1-49.14.4,8. D. T. M. Slock and T. Kailath, “Fasttransversal filters with data sequenee weighting, ” IEEE Trans. on Acoust., Speech, and Signal Processi ng, vol. 37, pp. 346-359, March 1989.9. D. T. M. Slock and T. Kailath, “Num

36、erically stable fast transversal filters for recursive least squares adaptive filtering, ” IEEE Trans. on Sig nal Processi ng, vol. 39, pp. 92-113, Jan. 1991.10. J. G Proakis, C. M. Rader, F. Li ng, and C. L. Nikias, Adva need Digital Sig nal Process in g,MacMilla n, New York, NY, 1992.11. B. Toplis

37、 and S. Pasupathy, “Track ing improveme nts in fast RLS algorithms using a variable forgetting factor,” IEEETrans. on Acoust., Speech, and Sig nal Process in g, vol. 36, pp. 206-227,Feb. 1988.12. D. Boudreau and P. Kabal, “ Joint-time delay estimation and adaptive recursive least squares filteri ng”

38、 IEEE Trans. on Signal Process ing, vol.41, pp. 592-601, Feb. 1993.8.7 问题1. 如下给出0(职 N) = A' )x(fe, N)Sd(晁1,2V)x(fct TV)A + T(k, N)S)(k 1, jV)x(fc. N)提示:使用矩阵求逆引理SD(k,N)。2. 如下给出出汕1需fjv仗)可(札N) _ -祇N) _出“Mil)如)加啊-盂更丽-2、+)wf,N(k)代表了 wf的最后一个元素(k,N)。3. 使用一个适当的混合关系的晶格 RLS算法基于后验和FTRLS算法,得到一个快速横向滤波器系数精确的初

39、始化过程。4. 显示下面的关系是有效的，假设输入信号prewindoweddetS£)(A N + 1)1detSD(-l,.V) =0min(AV)detS°（息 N + l）1detSD（N）=級.（kN）5. 如下给出7-i（a.心一ER川I）7 l'丿 AdctRD（A - 1,A7）提示：（let I + AB = det I + BA6使用的问题4和5的结果，证明入化爲(小)7推出方程(8.7)和(8.14)。还要表明，使用© KN)会增加FTRLS算法的 计算复杂度。8如果有一个避免使用转换因子丫 KN),有必要使用内部产品获得后验误差的快速

40、算法。推出一个没有转换因子的快速算法。9在问题6关系换为丫(k,N,3)SFTRLS关系派生算法，入叽()絆“(心)了(人 N )=描述结果相乘算法和显示它需要顺序8 N /输出示例。10. 推出等式(8.29)。11. FTRLS算法应用于预测信号x(k)= sin辿 匸。入=0.98计算误差143丿和前十的抽头系数迭代。12.SFTRLS算法应用于预测信号x(k)=sin岂匸。入二0.98计算误差143丿和前十的抽头系数迭代。13应用FTRLS算法于确定一个7阶未知系统的系数WT =0.0272 -0.0221 - 0.06210.11910.6116 - 0.3332 - 0.0190

41、- 0.0572输入高斯白噪声信号方差 衬二1,测量噪声也独立的高斯白噪声输入信号方差二2 = 0.01。模拟上面描述过程和测量超额 MSE当，=0.97和=0.98 时。14. 重复问题13,输入信号是一个一阶马尔可夫过程，二=0.98。15. 重做问题13使用FTRLS和SFTRLS定点实现算法。使用12位小数部分的信号和参数表示。16. 假设一个15阶FIR数字滤波器与乘数系数确定下面通过一个自适 应数字滤波器的使用FTRLS算法相同的顺序。假设定点运算，模拟识 别问题的以下规格描述：额外的噪音：白噪声方差(T2n = 0.0015系数字长:bc= 16位信号字长:bd = 16位输入信

42、号：高斯白噪声方差(T2x = 0.7入二 0.98wT = 0.0219360 0.0015786 - 0.0602449 - 0.0118907 0.13753790.0574545- 0.3216703- 0.5287203- 0.29577970.00020430.290670-0.0353349 - 0.0068210 0.0026067 0.0010333 - 0.0143593 画出有限和无限精确的实现的学习曲线。17. 重复上述的问题用SFTRLS算法。同样减少之前使用的字长直到明显增加(10%)多余的MSE观察输出。18. 重复问题16 SFTRLS算法，利用入二0.999和

43、入二0.960并评论结果。19. 将SFTRLS算法用于执行远期预测所产生的信号x(k)应用零均值高斯白噪声与输入单元方差的线性滤波器的传递函数0.5(1 - 0.827-2)(1 - l.Sz-1 - 0.87-2)计算得到的预测误差传递函数的零和比较线性滤波器的极点20. 计算均衡信道的脉冲响应hk) = 0.96fe + (-0.9)fc其中k=0,1,2,3,4,5,传输信号的零均值高斯白噪声与单位方差和自适应滤器的输入信噪比为 30 dB。使用SFTRLS算法顺序计算到100。第四章、基于LMS的改进型算法4.1 简介有许多自适应滤波器算法是来自前一章讨论传统的LMS算法。改进型算法

44、的目标是减少计算复杂度或收敛时间。在这一章主要介绍几个基于LMS算法的提出和分析，即量化误差算法1-11频域或变 换域 LMS 算法12卜14,归一化 LMS 算法15丄MS-Newton 算法 16-17,和仿射投影算法19-25。这一章还简要地讨论了几个相关 的算法的主要算法。量化误差算法减少计算复杂性代表了LMS算法的误差信号与短字或一个简单的2的幕数。LMS-Newton算法的收敛速度是独立于输入信号相关矩阵的特征 值扩散。这改进是通过使用逆的估计输入信号相关矩阵，计算复杂度大幅增加。归一化LMS算法利用变量收敛因子，最小化了瞬时错误。这样一 个收敛因子通常降低了收敛时间，但增加了失调

45、。在频域算法，变换应用到输入信号，以允许减少转换后的信号相关 矩阵的特征值分布与输入信号相关矩阵的特征值扩散。LMS算法应用到更好的条件转换信号实现更快的收敛。快速仿射投影算法重用旧的数据导致收敛，当输入信号是高度相关，导致一个家庭可以权衡计算复杂度的算法收敛速度。4.2量化误差LMS算法LMS算法的计算复杂度主要是因为乘法的执行自适应滤波器的系数更新和计算输出。在应用程序自适应过滤需要在高速、 回波消 除、通道等均衡，重要的是要尽量减少硬件的复杂性。简化LMS算法的第一步是应用量化误差信号，生成根据量化误差 算法更新滤波器系数w(A + 1) = w(fe) + 2/_tQe(k)x(k)Q

46、( )代表一个量化操作。量化函数离散的值，有界，不减少的。量 化的类型代表了量化误差算法。如果收敛系数卩是2的幕，号码，系数更新可以实现简单的乘法、基本组成的变化和补充。在许多应用程序中，如在全双工数据传输回波消除2和均衡的渠道与二进制数据3,输入 信号x(k)是一个二进制信号，即+ 1和-1假设值。在这种情况下，自适 应滤波器可以实现没有任何复杂的乘法。误差的量化实际上意味着最 小化目标函数，即通过Fe(k).般渐变类型算法系数执行更新通过w(t + 1) = w(fe)-对于一个线性组合器上面的方程可以改写为仏 +1)二-W(血)+ “(Q因此，目标函数最小化量化误差的算法是这样Oe(k)=2Qe 伙 I在这里Fe(k)是2 Qe(k)对e(k)的2倍积分。注意，链式法则应用于方程(4.3)是无效的不连续点Q( )在这里Fe(k)不是可微。量化误差算法和LMS算法的表现显然是不同的。一些广泛使用的分析量化误差算法将在下面写出表4.1符号误差算法Algoritlim 4.1Sign-Error AlgorithmInitializatio