通信信号自适应滤波处理仿真研究（二）

1 / 17通信信号自适应滤波处理仿真研究（二）图()( =)图()( =)观察两个不同步长情况下的误差曲线不难看出，步长越小，误差越小，但收敛速度越慢，为了好的精度，我们在选择时必然牺牲收敛速度。以上就是围绕对 LMS 算法的分析，着重讨论了算法的实现及算法中重要参数 的选择问题。在实际中，噪声功率大小的也会对系统的收敛程度产生影响，噪声功率越大，即信噪比 SNR 越小，误差曲线就会明显增加，这就是更大噪声功率对算法中随机梯度的影响，可以通过下面两个仿真图看出。分别取信噪比 SNR=5 和 SNR=20。 =图()(SNR=5)图()(SNR=20)递推最小二乘(RLS)算法最小二乘法2 / 17设已知 n 个数据 ， ， ，利用图所示的滤波器结构来估计期望信号 ， ， 。对的估计可表示成 式()估计误差 - 式()根据最小二乘法， (n)的最佳值应该使下列累计平方误差性能函数为最小 式() ， 其中 0 前加窗法最小二乘性能函数为 式()其中 。 引入 m 维矢量： 式()，而 维矩阵： 式()式()的最佳值满足方程 式()从而有 式()最终得到最小二乘算法的最后方程 式()递推最小二乘(RLS)算法由于最小二乘法的运算量较大，一般不适合实时滤波，采用递推算法可以减少运算量。由式()有 式()根据式得 式()对矩阵求逆得 式() 其中 为一纯量。 矩阵 式()N 维矢量 ， 为增益系数 式()3 / 17由式和式逆推式可得式()利用式，就可以用递推的方式求 m m 维矩阵 的逆，使运算量降低。式两端乘以 ，利用式可得式()另外，根据式可得 式()将式，式,式代入式就可以得到式()利用式和式的最后两项可简化为 ，而式的前两项中的 即为 。所以由式可得式()这即为递推最小二乘(RLS)算法的递推公式。下图为 RLS 算法的流程图：RLS 算法的参数分析RLS 算法具体实现需要以下计算式 ； ； = ； ；其中 个参数意义与 LMS 相同，新增个参数意义为 ：反相关矩阵； (n)：增益向量； ：遗忘因子。在 RLS 算法中遗忘因子是一个接近 1 但是小于 1 的正数，一般来说介于到 1 之间。使用遗忘因子的目的在于把接近目前时间点的信息乘上越大的权值，而离目前时间4 / 17点越远的信息乘上越小的权值，也就是说，我们重视较近时间点的信息甚与较远时间点的信息。若等于 1，则表示对所有的信息都一样，其权值都是相同的。RLS 算法 Matlab 仿真分析图上图分别为输入信号，输出信号和误差信号的曲线，可以看出输出信号在经过一段时间的自适应调整后，便能基本达到跟踪，滤波的效果。从误差信号曲线也可以看出这点，误差输出经过一段时间就趋于稳定。图上图为误差平方的均值曲线，大约在 t=300 时，误差趋于收敛，系统完成自适应过程。 以上就是围绕对 RLS 算法的分析，着重讨论了 RLS算法推导，具体实现的相关公式以及运用 matlab 软件对其进行仿真。RLS 算法与 LMS 算法的比较分析RLS 算法能够在很短的时间内就趋于收敛，而 LMS 算法则有一个比较长的渐变过程，所以 RLS 的跟踪性能要优于 LMS，这可以从图和图看出。换句话说，RLS 比 LMS 的收敛速度要快。可以通过下图看出：5 / 17上图蓝色是 LMS 收敛曲线，红色为 RLS 收敛曲线。可以看出明显 RLS 收敛性要优于 LMS 算法。而且 LMS 在收敛后波形还有较大波动，而 RLS 就要小的多，基本没有波动，这说明 RLS 的稳态误差也是小于LMS 的，从图和图可以看出。但是由于 LMS 计算量简单，适合于硬件实现，这是 RLS 无法相比的。所以二者各有优劣。(以上 LMS 和 RLS 算法仿真均采用相同的外界信号及采样时间点) 归一化 LMS 算法(NLMS)NLMS 算法实现NLMS 算法是将 LMS 算法中的 值重新定义，让 值会随输入信号之正规化作改变，能提升收敛的稳定性。下面为 NLMS 算法实现所需的计算式：； ； ； ； 各参数的定义和 LMS 算法定义相同，新增参数的定义为 ：很小的正常数，一般取 =1e-10。LMS 算法的稳定度和收敛速率受到 值和参考信号的影响，由于 值为一固定值，因此 LMS 的整体收敛速率就受它的影响，收敛速率对变化较快的信号反应并不理想。而NLMS 算法能改善输入信号对收敛因子的影响， 值随着时间6 / 17n 变化成为 (n)，使之随时变化，从而调节至最佳值。另外为了避免当输入信号过小时造成收敛因子的发散，还加入 值。下图为 NLMS 算法的流程图：SHAPE * MERGEFORMATNLMS 算法的 Matlab 仿真分析图()图()由图和图可以看出 NLMS 算法的自适应滤波过程及误差收敛情况，而且相比 LMS 算法在相同条件下，NLMS 算法要比 LMS 算法收敛更快一些。这是因为 NLMS 算法的计算量与 LMS 相当，但是由于 NLMS 算法的收敛条件与输入信号的特征值无关，故 NLMS 算法比 LMS 算法的收敛速率快。下图可以看出二者收敛的差别：图 LMS 的误差收敛图 NLMS 的误差收敛LMS 算法中大约在 t=600 时开始收敛，而 NLMS 则是7 / 17在大约 t=400 时开始收敛。明显 NLMS 的收敛速度要快与LMS。3 自适应信号处理的应用及 Matlab 仿真通信中的自适应噪声抵消 在通信和其他许多信号处理应用问题中，接受信号中往往伴随着干扰和噪声，从而显著影响接受信号的可靠性，或者导致误码率上升。一般来说，干扰和噪声的存在总是难免的。信号处理技术的核心问题之一就是从受到干扰污染的信号中估计，检测或者恢复出原始信号。而自适应噪声抵消的基本原理就是将被噪声污染的信号与参考信号进行抵消运算，从而消除带噪信号中的噪声。其关键问题是自适应噪声抵消系统的参考信号一定要与待消除的噪声具有一定的相关性，而与要检测或者提取的信号不相关。自适应噪声抵消系统经过自适应系统的控制和调整，系统能够有效地从噪声中恢复出原始信号。作为自适应信号处理领域的重要分支之一，它已经受到了人们的普遍关注并得到了广泛的应用。自适应噪声抵消系统的基本原理 下图为典型自适应噪声抵消系统的原理框图：图 自适应噪声抵消系统8 / 17在图中，原始输入信号 d(n)为有用信号 s(n)与噪声v(n)之和,参考输入信号 x(n)是与 v(n)相关的噪声 u(n)。假定 s(n)，v(n)，u(n)均为零均值平稳随机过程，且满足s(n)与 v(n)及 u(n)互不相关。由图可见，自适应滤波器的输出 为噪声 u(n)的滤波信号。则整个自适应噪声抵消系统的输出 y(n)为 式() 而 式()对式()两边取期望，由于 s(n)与 v(n)及 u(n)互不相关，且 s(n)与 也不相关，故有 式()信号功率 与自适应滤波器的调节无关，因此，调节自适应滤波器使 最小，等价于使 最小。这样由式()，有 式()由此可见，当 最小时， 也达到最小，即自适应噪声抵消系统的输出信号 y(n)与有用信号 s(n)的均方误差最小。在理想情况下，当 时，有 。这时，自适应滤波器自动地调节其权值，将 u(n)加工成 v(n)，与原始输入信号d(n)中的 v(n)相减，使输出信号 y(n)的噪声完全被抵消，而只保留有用信号 s(n)。但是自适应滤波器能够完成上述9 / 17任务的必要条件为：参考输入信号 必须与被抵消的噪声v(n)相关。自适应噪声抵消系统 Matlab 仿真以下仿真采用图的结构，分别运用 LMS，NLMS 和RLS 循环算法进行噪声消除。图仿真得出三种自适应滤波算法提取正弦信号的曲线图。可以看出系统能基本还原出原始信号，达到噪声抵消的效果。但是用 RLS 算法提取的正弦信号质量要好，其中LMS 算法提取的信号效果最差，存在没有滤除的随机噪声部分较多，而 NLMS 算法要比 LMS 的效果要好，但比起 RLS 算法在估计精度上有些波动，存在一定的残余误差，即有一定失调。自适应陷波滤波器在通信系统和其他电子系统中，经常会受到诸如50Hz 工作频率等单频干扰或者窄带干扰的影响。这种干扰的存在，严重影响了信号的接收或者检测的可靠性和正确性，因此必须加以消除。陷波滤波器是消除这种干扰的有力工具，当自适应噪声抵消系统的参考输入为单一频率正弦信号时，则系统可以构成自适应陷波滤波器。自适应陷波滤波器的原理10 / 17自适应陷波滤波器具有陷波中心频率，且该频率与其参考输入的正弦信号的频率相同。另一方面，自适应陷波滤波器还能够随着干扰频率的变化，自动地修正系统自身参数来跟踪这种变化。典型的单一频率自适应陷波滤波器的原理图如图()所示，图()表示一个具有两个自适应实权的自适应噪声对消器。它等效于有一个复权的噪声对消系统，即用两个实权达到同时调整单一频率正弦波的幅度和相位，以消除干扰的目的。假定原始输入信号的类型是任意的，而参考输入是频率为 f 的纯正弦波，即式()图中第一个权的输入直接由参考输入采样得到，而第二权的输入是将第一个权输入移相 产生。即它们可分别表示为式() SHAPE * MERGEFORMAT 图()其中 ，权的迭代用 LMS 算法，如下式所示，权的修正过程如下： 式() 自适应陷波滤波器 Matlab 仿真分析11 / 17图上图可以看出经过正弦信号干扰的原始信号，在通过自适应陷波滤波器后，基本达到噪声消除的效果。上图中第一个图为原始信号，第二个为经过正弦信号干扰后的信号，第三个为消噪后的信号，第四个为误差信号曲线。自适应预测自适应预测的基本思想要得到预测系数，必须获得输入信号采样值的相关函数矩阵，而实际上它不是一个定值，是时变的，所以就要求必须自适应调整预测系数，以保持最佳的预测增益。求相关函数的简单方法是，先采样并存储一个定长时间间隔的信号值，计算这些采样值的自相关函数，然后确定最佳的预测系数。预测器每隔规定的时间间隔更新依次存储的采样数据，并且每次将计算的预测参数发送到接收端。通过上述方法动态调整预测参数，在存储采样值时间间隔较长或每次存贮采样值个数较大的情况下，可以获得很大的预测增益。这就是自适应预测器的基本思想。忽略量化噪声的影响，预测误差函数式()需要说明，考虑到实际系统的可实现性，可以用误差函数的量化值 。调整预测系数使误差函数向负梯度的12 / 17方向变化，即式()式中，sgn是符号函数， 是预测系数自适应速率，需要根据实验确定其最佳值。也可以考虑用平方差值函数确定预测系数，即式()自适应预测器的实现比较复杂，但是，当信号采样值相关距离大或信号统计特性的平稳性不佳，无法获得确切和恒定的相关系数的情况下，自适应预测是较理想的预测方法。在许多情况下，一个宽带信号既受到周期性干扰的污染，又没有无信号的外部参考输入可以利用。此时，可以直接从原始输入引出，接入一具有固定延迟的延迟线，则可得到类似的参考输入支路。这种结构实际上是一个自适应预测器。下图仿真采用的是线性预测滤波方法抑制窄带干扰的算法。 自适应预测的 Matlab 仿真图由上图可以看出预测信号在经过一段自适应过程后能够很好的跟踪接收信号，达到预测效果。13 / 17自适应均衡自适应均衡的基本原理自适应均衡器的工作过程包含两个阶段，一是训练过程，二是跟踪过程。在训练过程中，发送端向接收机发射一组已知的固定长度训练序列，接收机根据训练序列设定滤波器参数，使检测误码率最小。典型的训练序列是伪随机二进制信号或一个固定的波形信号序列，紧跟在训练序列后面的是用户消息码元序列。接收机的自适应均衡器采用递归算法估计信道特性，调整滤波器参数，补偿信道特性失真，训练序列的选择应满足接收机均衡器在最恶劣的信道条件下也能实现滤波器参数调整，所以，训练序列结束后，均衡器参数基本接近最佳值，以保证用户数据的接收，均衡器的训练过程成功了，称为均衡器的收敛。在接收用户消息数据时，均衡器还需要不断跟踪信道特性的变化并随信道特性的变化连续地改变均衡器参数。 自适应均衡器的实现下面讨论自适应均衡器的具体实现。我们知道信道均衡器均衡器的作用是在信道通带内形成一个信道传输函数的逆，而在通带之外它的增益则很小或者为零。因而，由信道和均衡器级联组成的系统在通带内有基本均匀的振幅特性，而带外基本为零，相位响应在带内是频率的线性14 / 17函数。如果条件满足，联合冲激响应就是辛格函数，符号间干扰可以消除。自适应调整也解决了信道本身未知，时变的特性所带来的困难。下图为自适应均衡器的基本结构。SHAPE * MERGEFORMAT图逆模拟用一个自适应横向滤波器(LMS 滤波器)，由于输入 x(k)的信号带宽受信道带宽的限制，因而，自适应滤波器仅需在信道的通带内去均衡信道的振幅和相位特性。如果能知道信道的输入，并考虑到整个系统的延迟，就可得到期望响应 d(k)，但是一般是难于获得的。周期性地中断信息传输，发射一些已知的码序列，便可以进行自适应调整。贝尔电话实验室的拉克提供了一种得到期望响应 d(k)的方法，这种方法用自适应滤波器自身输出提供 d(k)，因此避免了对发射信号任何先验信息的依赖，拉克称该方法为“判决指向学习” 。更确切地说，期望信号 d(k)=sgn y(k)，如图所示，它是由一个量化滤波器产生的。由于数据是二进制的，若不考虑噪声影响，则经适当均衡了的信道在选通时间内的取样输出为+1 或-1，然后将滤波器输出和经量化后的输出比较，产生误差信号 e(k)。由于均衡器输出应15 / 17该在适当的选通时间内唯一地表示各自的辛格脉冲，因而自适应只许在选通时间内进行，这可用与发射信号同步的闸门脉冲对误差信号 e(k)选通来实现。从平均意义上来说，如果量化后的期望响应是正确的，则自适应将沿着正确的方向进行。自适应均衡器的 Matlab 仿真下面将采用上述的实现方法进行 matlab 仿真，可以看出自适应均衡器对消除信道的干扰的作用。图下图是对自适应均衡器在不同信噪比下误码率的仿真，能进一步说明自适应均衡器的作用。采用 50 万点仿真误码率：图上图中虚线部分是没有经过均衡器的误码