自适应滤波器的算法研究及DSP仿真实现

滤波是电子信息处理领域的一种最基本而又极其重要的技术。在有用信号的传输过程中,通常会受到噪声或干扰的污染。利用滤波技术可以从复杂的信号中提取所需要的信号,同时抑制噪声或干扰信号,以便更有效地利用原始信号。滤波器实际上是一种选频系统,它对某些频率的信号予以很小的衰减,让该部分信号顺利通过;而对其他不需要的频率信号则予以很大的衰减,尽可能阻止这些信号通过。在电子系统中滤波器是一种基本的单元电路,使用很多,技术也较为复杂,有时滤波器的优劣直接决定产品的性能。   本文将介绍自适应滤波器,它是数字信号的基础,大神就看看可以了。

1 自适应滤波器简介
    自适应滤波器属于现代滤波器的范畴,自适应滤波器是相对固定滤波器而言的,固定滤波器属于经典滤波器,它滤波的频率是固定的,自适应滤波器滤波的频率则是自动适应输入信号而变化的,所以其适用范围更广。在没有任何关于信号和噪声的先验知识的条件下,自适应滤波器利用前一时刻已获得的滤波器参数来自动调节现时刻的滤波器参数,以适应信号和噪声未知或随机变化的统计特性,从而实现最优滤波。所谓自适应滤波,就是利用前一时刻已获得的滤波器参数等结果,自动地调节现时刻的滤波器参数,以适应信号和噪声未知的或随时间变化的统计特性,从而实现最优滤波。自适应滤波器实质上就是一种能调节其自身传输特性以达到最优化的维纳滤波器。

2 自适应滤波原理
    自适应滤波器的特性变化是由自适应算法通过调整滤波器系数来实现的。一般而言,自适应滤波器由两部分组成,一是滤波器结构,二是调整滤波器系数的自适应算法。自适应滤波器的结构采用FIR或IIR结构均可,由于IIR滤波器存在稳定性问题,因此一般采用FIR滤波器作为自适应滤波器的结构。图1给出了自适应滤波器的一般结构。
    图1为自适应滤波器结构的一般形式,图中x(n)为输入信号,通过参数可调的数字滤波器后产生输出信号y(n),将输出信号y(n)与标准信号(或者为期望信号)d(n)进行比较,得到误差信号e(n)。e(n)和x(n)通过自适应算法对滤波器的参数进行调整,调整的目的使得误差信号e(n)最小。

自适应滤波器设计中最常用的是FIR横向型结构。图2是横向型滤波器的结构示意图。

其中:x(n)为自适应滤波器的输入;w(n)为自适应滤波器的冲激响应:w(n)={w(O),w(1),…,w(N-1)};y(n)为自适应

3 自适应滤波算法
    自适应滤波器除了包括一个按照某种结构设计的滤波器,还有一套自适应的算法。自适应算法是根据某种判断来设计的。自适应滤波器的算法主要是以各种判据条件作为推算基础的。通常有两种判据条件:最小均方误差判据和最小二乘法判据。LMS算法是以最小均方误差为判据的最典型的算法,也是应用最广泛的一种算法。
    最小均方误差(Least Mean Square,LMS)算法是一种易于实现、性能稳健、应用广泛的算法。所有的滤波器系数调整算法都是设法使y(n)接近d(n),所不同的只是对于这种接近的评价标准不同。LMS算法的目标是通过调整系数,使输出误差序列e(n)=d(n)-y(n)的均方值最小化,并且根据这个判据来修改权系数,该算法因此而得名。误差序列的均方值又叫“均方误差”(Mean Sqluare Error,MSE)。
    理想信号d(n)与滤波器输出y(n)之差e(n)的期望值最小,并且根据这个判据来修改权系数wi(n)。由此产生的算法称为LMS。均方误差ε表示为:

对于横向结构的滤波器,代入y(n)的表达式:

其中:R=E[X(n)XT(n)]为N×N的自相关矩阵,它是输入信号采样值间的相关性矩阵。P=E[d(n)X(n)]为N×1互相关矢量,代表理想信号d(n)与输入矢量的相关性。在均方误差ε达到最小时,得到最佳权系数

它应满足下式:

这是一个线形方程组,如果R矩阵为满秩的,R-1存在,可得到权系数的最佳值满足:W*=R-1p。用完整的矩阵表示为:

显然φx(m)=E[x(n)x(n-m)]为x(n)的自相关值,φxd(R)=E[x(n)d(n一k)]为x(n)与d(n)互相关值。在有些应用中,把输入信号的采样值分成相同的一段(每段称为一帧),再求出R,P的估计值得到每帧的最佳权系数。这种方法称为块对块自适应算法。如语音信号的线性预测编码LPC就是把语音信号分成帧进行处理的。R,P的计算,要求出期望值E,在现实运算中不容易实现,为此可通过下式进行估计:

用以上方法获得最佳W*的运算量很大,对于一些在线或实时应用的场合,无法满足其时间要求。大多数场合使用迭代算法,对每次采样值就求出较佳权系数,称为采样值对采样值迭代算法。迭代算法可以避免复杂的R-1和P的运算,又能实时求得近似解,因而切实可行。LMS算法是以最快下降法为原则的迭代算法,即W(n+1)矢量是W(n)矢量按均方误差性能平面的负斜率大小调节相应一个增量:W(n+1)=W(n)-μ▽(n),这个“是由系统稳定性和迭代运算收敛速度决定的自适应步长。▽(n)为n次迭代的梯度。对于LMS算法▽(n)为下式E[e2(n)]的斜率:

由上式产生了求解最佳权系数W*的两种方法,一种是最陡梯度法。其思路为:设计初始权系数W(0),用W(n+1)=W(n)一μ▽(n)迭代公式计算,到W(n+1)与W(n)误差小于规定范围。其中▽(n)计算可用估计值表达式:

上式K取值应足够大。如果用瞬时一2e(n)X(n)来代替上面对-2E[e(n)X(n)]的估计运算,就产生了另一种算法——随机梯度法,即Widrow-Hoft的LMS算法。此时迭代公式为:

W(n+1)=W(n)+2ue(n)X(n)
    以后讨论的LMS算法都是基于WidrOW-Hoff的LMS算法。上式的迭代公式假定滤波器结构为横向结构。对于对称横向型结构也可推出类似的迭代公式:
    W(n+1)=W(n)+2ue(n)[X(n)+X(n一N+1)]
4 自适应滤波算法的理论仿真
    使用Matlab编程,采用自适应滤波器技术实现信噪分离,也就是去噪。程序如下:

程序运行的结果如图3所示。

通过CCS软件环境,把滤波程序烧录到DSP芯片中,在CCS DSK C5000环境下输出仿真结果:输入信号为余弦信号和随机噪声的叠加。程序正确运行后,观察运行结果,得出如图4,图5所示的仿真图。

说明:本文来源网络;文中观点仅供分享交流,不代表本公众号立场,转载请注明出处,如涉及版权等问题,请您告知,我们将及时处理。

Please clik the advertisement and exit

重点

如何下载 《华为软件编程规范总则(C语言,C++,JAVA)》高清PDF电子书

(0)

相关推荐

  • WEBRTC AEC3算法原理

    延迟估计算法 原理 AEC3的延迟估计算法与AEC的非线性处理的延迟估计算法思想一致,因为回声能量是呈指数衰减,所以计算滤波器能量最大块作为延迟估计值,但是比AEC的延迟估计算法复杂的多 AEC3延迟 ...

  • ANC主动降噪理论及Matlab代码实现

    根据系统是否有参考信号传感器可将ANC系统大致的分为前馈型和反馈型. 前馈控制是产生次级噪声之前就通过传感器测量初级噪声的频率以获取参考信号. 反馈控制不需要测得参考信号就产生次级噪声进行相消干涉 反 ...

  • 【信号去噪】基于低通和自适应滤波LMS去噪matlab源码

    一.简介 自适应滤波是近30年以来发展起来的关于信号处理技术的方法.它是在维纳滤波.Kalman滤波等线性滤波基础上发展起来的一种最佳滤波方法.由于它具有更强的适应性和更优的滤波性能,从而在工程实际中 ...

  • 压缩感知Compressed Sensing基本原理

      压缩感知技术(CS,Compressed Sensing)因其能够大幅度提升MRI采集效率.降低采集时间而被大家所知,但它究竟是如何实现在尽量保证图像真实信息的情况下实现成像速度提升的呢?带着好奇 ...

  • 【学术论文】子阵级LCMV循环优化自适应波束形成算法研究

    摘要: 在传统LCMV波束形成器以及子阵空间部分自适应阵的基础上,提出了一种新颖的降维方法.首先将大规模阵列按照子阵划分的某种规则划分为若干组子阵列,每一组子阵列使用相同的权值.在权值优化过程中,每一 ...

  • 离散状态事件驱动仿真方法及自适应预估校正算法

    清华大学电机系电力系统及发电设备安全控制和仿真国家重点实验室的研究人员杨祎.赵争鸣.檀添.李帛洋.袁立强,在2017年第12期<电工技术学报>上撰文指出,在电力电子系统分析中,需要对带间断 ...

  • 低信噪比下的红外弱小目标检测算法研究综述

    低信噪比下的红外弱小目标检测算法研究综述

  • 基于数据关联的多雷达点迹融合算法研究

    0 引言 现代战争已成为高科技的信息战争,随着反辐射导弹.低空突防.隐身技术和电子对抗技术的发展,仅仅依赖单部的雷达已经无法满足作战需求.目前装备的雷达广泛应用航迹融合技术,但是多雷达组网的点迹融合技 ...

  • 空中侦察图像连续拼接中的平滑算法研究(中国知网 空中侦察)

    63961部队 北京100012 摘    要: 侦察图像在现代战争中具有重要意义, 利用图像处理技术对序列空中侦察图像进行拼接所生成的全场景图能为战场毁伤评估提供重要信息.对配准后的图像重叠部分进行 ...

  • 基于EC集合预报动态百分位的降水订正算法研究

    贵州特殊的地理位置及复杂的地形特征,形成了独有的天气气候特征,因此降水的时空分布比较复杂,呈现出更加明显的不均匀性.贵州的降水预报准确率一直较低,尤其是分级降水预报暴雨以上量级的预报准确率很低,202 ...

  • 基于单幅图像的2D转3D算法研究

    最近,3D影片盛行,3D电视技术也层出不穷,3D技术在带给大家非凡的视觉冲击同时,也在告诉大家这背后隐藏了太多的商机. 目前的3D技术大体分为两种:软件技术和硬件技术.所谓硬件技术就是在拍摄影片时采用 ...

  • 影像仪自动寻边算法研究

    作者:牧象仪器 来源:www.szmxyq.cn 时间:2017-07-27 10:02        影像仪的应用主要在于对机床加工出来的零件进行测量,确保加工过程的精确性.在对一些尺寸进行检测的过 ...

  • 中医状态辨识算法研究

    <中医状态学>是中医治未病理论的奠基之作,为治未病.健康管理提供了理论指导和顶层设计.状态与健康状态辨识是中医状态学的重要组成内容,算法模型的科学性.准确性直接影响着健康状态辨识的发展.建 ...