matlab肌电信号平滑滤波_最速离散跟踪-微分器:信号滤波、提取微分信号

本文介绍了使用最速离散跟踪-微分器对肌电信号进行平滑滤波的方法。通过Matlab仿真,展示了该微分器如何有效抑制噪声,平滑信号并提取微分信号,适用于处理包含随机噪声和单帧跳变的信号。调整参数r和h可以平衡跟踪效果和微分信号的平滑度。
摘要由CSDN通过智能技术生成

一、背景介绍

在实际工程问题中, 测量信号经常有不连续或带有随机噪声的问题, 但经常需要提取连续信号以及对应的微分信号。

比如用编码器做速度检测时, 需要由数字的绝对位移信号提取速度信号,并进一步提取加速度信号等。

95907d14328e1b48b3be6fcfb0ef35cc.png

20世纪80年代末期,系统与控制专家韩京清在吸收并丰富PID控制精髓(基于误差来消除误差)的基础上,针对PID控制的缺点(对增益变化太敏感、快速性与超调之间的矛盾、没有合适的微分器、对变化的扰动抑制能力不显著等),相继开发出非线性跟踪微分器、扩张状态观测器、非线性误差反馈机制等对信号处理及控制具有独特功能和效率的新型非线性动态结构,统称为自抗扰控制技术。

二、最速离散跟踪-微分器

关于跟踪微分器的理论推导见韩京清的《自抗扰控制技术》这本书,这里对最速最速离散跟踪-微分器的结构及公式做一个简单介绍,并通过仿真观察最速离散跟踪-微分器的信号滤波和提取微分信号的效果。

53e3a89f81f85b3a6d071a9ff5455d97.png

三、Matlab代码

仿真测试步骤:

  1. 创建测试信号
  2. 增加随机噪声和跳变的噪声
  3. 调用最速离散跟踪-微分器处理信号
  4. 观察信号滤波和提取微分信号的效果

matlab代码如下:

33014b2dc1b552358e2d0edbafc17fab.png
0958ae3446db2af025b7af0654acb779.png
60521e645daef8b795ff315e1b39b850.png

四、仿真结果

增加噪声后的信号如下图所示,可以看到信号不仅有随机噪声,而且存在单帧跳变。

cefd4db174faae408631f932defca9df.png

带噪声的信号

调用最速离散跟踪-微分器处理信号的结果如下:

5b6765f9d4a55d3e330984c29343c6bc.png

最速离散跟踪-微分器的处理结果

其中,蓝色点线为最速离散跟踪-微分器的处理结果。可以看到,

  • 跟踪信号可平滑的模拟原始信号,几乎不受单帧跳变的影响,有一定滞后;
  • 而微分信号,相比差分的方法(Diff Method)可有效抑制噪声的干扰。

在实际应用过程中,最速离散跟踪-微分器只需调整r和h两个参数。而h通常设为信号的采样周期,通过调整r值的大小来调整跟踪和微分信号的效果。

  • r越大,跟踪效果越好,但微分信号会增加高频噪声;
  • r越小,跟踪和微分信号越平滑,但滞后现象越严重。
3f89cbb7e720af7e7344a37ee4f492bb.png

r较小的效果

五、结论

跟踪微分器可用于信号的滤波,并有效平滑信号的单帧突变;同时可用于提取微分信号,相比用差分法来计算信号微分的方法,跟踪-微分器能有效抑制测量信号的噪声干扰,且不易受信号单帧跳变的影响。

matlab中的微分方程-matlab中的微分方程.doc 1510 matlab中的微分方程 第1节  Matlab能够处理什么样的微分方程? Matlab提供了解决包括解微分方程在内的各种类型问题的函数: 1. 常规微分方程(ODEs)的初始值问题 初值问题是用MATLAB ODE求解解决的最普遍的问题。初始值问题最典型的是对非刚性度(?nonstiff)问题应用ODE45,对刚性度(?stiff)问题采用ODE15S。(对于stiffness的解释,请参照“什么是Stiffness”一节。) 2. 微分-代数方程(DAEs)的初值问题 在那些守恒定律规定一些变量之间满足常数关系领域经常遇到这类问题。Matlab 可以用ODE15S 或者 ODE23T解决索引(index)为1的DAEs。(对于索引的解释,请参阅“DAEs与他们的索引”一章。) 3. 边界值问题(BVPs) 这种通常要求微分方程在两边都具有特殊的条件组成。尽管他们通常不象IVPs那样经常遇到,但是他们也是工程应用中比较常见的问题。可以利用函数BVP4C来解决这类问题。 4. 时延微分方程(DDEs) 这类微分方程包含了独立变量的延迟。他们在生物与化学模型这类大量的应用中遇到,可以通过DDE23来解决这类问题。 5. 偏微分方程(PDEs) 采用PDEPE可以解决一维时空的抛物面与椭圆方程的初值、边界值的问题。而那些对更加多的一般的偏微分方程感兴趣的可以利用PDE工具箱。 更多的matlab的综合应用技术的信息请参阅Solution8314。 更多的有关matlab采用的各种求解的算法的信息请查看下面的URLs: ● ODE 函数 ● BVP 函数 ● DDE 函数 ● PDE 函数 第2节 可以从什么地方获得更多的指导与附加信息?    可以从MATLAB Center、网站的新闻组、文件交换点可以获得一系列资料,可以进一步解释MATLAB解决各种方程(ODE,DAE,BVP,DDE)的求解的算法和使用。你可以下载各种方程的文章与手册,他们通常带有大量的实例。   可以从 matlab自带的帮助文件的 Mathematics|Differential Equations下找到使用指导。   Cleve Moler的《Numerical Computing with MATLAB》的第七章详细讨论了OEDs的解法,并附带有大量的实例与简单的问题练习。    第3节 对ODE求解的语法存在有些什么变化? 在MATLAB6.5(R13)中应用ODE求解求解的首选语法是: [t,y]=odesolver(odefun,tspan,y0,options,parameter1,parameter2,…,parameterN); odesolver 是你采用的求解,例如ODE45或者ODE15S。odefun是微分方程的定义函数,所以odefun定义独立参数(典型的是时间t)的导数y‘ 以及y和其他的参数。在MATLAB6.5(R13)中,推荐使用函数句柄作为odefun。 例如,ode45(@xdot,tspan,y0),而不是用 ode45('xdot',tspan,y0)。 请看采用函数句柄的好处的文档: 采用函数句柄传递你定义MATLAB求解计算的量、例如大规模矩阵或者Jacobian模式的函数。 如果你喜好采用字符串儿传递你的函数,matlab求解将回溯匹配。 在老的matlab版本里,通过传递标志来规定求解的状态和恰当的计算。在MATALB6.0以及其后的版本中,这就没有必要了,可以从matlab自带的文档中发现这个差别。 如果里采用的matlab的ODE求解的老的语法,你可以看看我们FTP站点上的各种求解的老的实例: ftp://ftp.mathworks.com/pub/doc/papers/ 前面的站点包含了BVP,DAE与DDE这三个方向的采用老的语法的实例。你可以在下面的站点中找到应用ODE45与ODE23的实例: ftp://ftp.mathworks.com/pub.mathworks/toolbox/matlab/funfun 你可以在MATLAB Center的文件交换站点查看这些例子的更新版本。 第4节  如何减小ODE的阶次? 求解一阶ODE的代码是很直接的。然而,二阶或者三阶的ODE不能够直接应用求解。你必须先将高阶的ODE改写成一阶的ODEs系统,使得它可以采用MATLAB ODE求解。 这是一个如何将二阶微分方程改写成两个一阶微分
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值