拓端tecdat|matlab实现扩展卡尔曼滤波(EKF)进行故障检测

本文主要是介绍拓端tecdat|matlab实现扩展卡尔曼滤波(EKF)进行故障检测,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

最近我们被客户要求撰写关于扩展卡尔曼滤波的研究报告,包括一些图形和统计输出。

本文展示了如何使用扩展卡尔曼滤波器进行故障检测。本文使用扩展的卡尔曼滤波器对一个简单的直流电机的摩擦力进行在线估计。估计的摩擦力的重大变化被检测出来,并表明存在故障。

电机模型

电机被模拟成具有阻尼系数c,转动惯量J,由一个扭矩u驱动。电机的角速度w和加速度˙w是测量输出。

为了使用扩展卡尔曼滤波器估计阻尼系数c,为阻尼系数引入一个辅助状态,并将其导数设为零。

因此,模型状态,x = [w;c],和测量,y,方程式为:

使用近似值˙x=xn+1-xnTs将连续时间方程转换为离散时间,其中Ts为离散采样周期。这就得到了离散时间模型方程。

指定电机参数

J  = 10;    % 惯性
Ts = 0.01;  % 采样时间

指定初始状态。

x = [...0; ... % 角速度 1]; % 摩擦力

% 以摩擦为状态的电机的状态更新方程
% 状态更新方程式
x1 = [...x0(1)+Ts/J*(u-x0(1)*x0(2)); ...x0(2)];

% 以摩擦为状态的电机的测量方程% 输出。
% y - 电机测量元素[角速度;角加速度]。y = [...x(1); ...(u-x(1)*x(2))/J];

电机经历状态(过程)噪声干扰,q,和测量噪声干扰,r。噪声项是相加的。

过程和测量噪声的平均值为零,E[q]=E[r]=0,协方差Q=E[qq']和R=E[rr']。摩擦状态有很高的过程噪声干扰。这反映了我们希望摩擦力在电机正常运行期间会有变化,并且滤波器能跟踪这种变化。加速和速度状态的噪声很低,但速度和加速度测量的噪声相对较大。

指定过程噪声协方差。

[...1e-6 0; ...   % 角速度0 1e-2];      % 摩擦力

指定测量噪声协方差。

[...1e-4 0; ...  % 速度测量0 1e-4];     % 加速度测量

创建一个扩展的卡尔曼过滤器

创建一个扩展的卡尔曼滤波器来估计模型的状态。我们特别关注阻尼状态,因为这个状态值的急剧变化表明存在故障事件。
创建一个扩展卡尔曼滤波器对象,并指定状态转换和测量函数的雅各布系数。

扩展卡尔曼滤波器的输入参数是之前定义的状态转换和测量函数。初始状态值x0、初始状态协方差、过程和测量噪声协方差也是扩展卡尔曼滤波器的输入。在这个例子中,精确的雅各布函数可以从状态转换函数f和测量函数h中得到。


% 输出
Jac - 在x处计算出的状态雅各布系数% 雅各布系数
Jac = [1-Ts/J*x(2) -Ts/J*x(1); ...0 1];% 电机模型测量方程的雅各布系数% 输出
Jac - 在 x 处计算的测量雅各布系数% 雅各布系数
J = [ ...1 0;-x(2)/J -x(1)/J];

Simulation仿真

为了模拟工厂,创建一个环路,在电机中引入一个故障(虚构的电机剧烈变化)。在模拟回路中,使用扩展的卡尔曼滤波器来估计电机状态,并特别跟踪摩擦状态,检测摩擦力何时发生统计意义上的变化。

电机被模拟成一个脉冲序列,反复加速和减速。这种类型的电机操作对于生产线上的采摘机器人来说是典型的。

在模拟电机时,加入与构建扩展卡尔曼滤波器时使用的Q和R噪声协方差值相似的过程和测量噪声。对于摩擦,使用一个小得多的噪声值,因为除了故障发生时,摩擦大多是恒定的。在模拟过程中人为地诱发故障。

Qv = chol(Q); % 过程噪声的标准偏差
Qv(end) = 1e-2; % 较小的摩擦噪声
Rv = chol(R); % 测量噪声的标准偏差

使用状态更新方程对模型进行仿真,并在模型状态中加入过程噪声。仿真十秒钟后,强制改变电机的摩擦力。使用模型测量功能来模拟电机传感器,并在模型输出中加入测量噪声。

for ct = 1:numel(t)  % 模型输出更新   y = y+Rv*randn(2,1); % 添加测量噪声% 模型状态更新 xSig(:,ct) = x0;% 诱发摩擦力的变化if t(ct) == 10x1(2) = 10; % 步骤变化x1n = x1+Qv*randn(nx  % 加入过程噪声

Significant friction change at 10.450000

为了从电机测量值中估计电机状态,使用扩展卡尔曼滤波器的预测和纠正命令。

% 使用扩展卡尔曼滤波器进行状态估计x_corr = correct(ekf,y,u(ct),J,Ts); % 根据当前测量结果修正状态估计。predict(ekf,u(ct),J,Ts); % 根据当前状态和输入预测下一个状态。

为了检测摩擦力的变化,使用4秒的移动窗口来计算估计的摩擦力平均值和标准偏差。在最初的7秒之后,锁定计算的平均值和标准差。这个最初计算出的平均值是摩擦力的预期无故障平均值。7秒后,如果估计的摩擦力与预期的无故障平均值相差超过3个标准差,这就意味着摩擦力有了明显的变化。为了减少噪音和估计摩擦力的影响,在与3个标准差的界限比较时,使用估计摩擦力的平均值。

% 计算估计平均值和标准偏差。else% 存储计算出的平均数和标准差,不需要%重新计算。fMean(ct) = fMean(ct-1) % 使用预期的摩擦力平均值和标准偏差来检测%摩擦力变化。estFriction = mean(xSigEst(2, if fChanged(ct) && ~fChanged(ct-1) % 检测摩擦变化信号的上升沿|fChanged|

使用估计的状态来计算估计的输出。计算测量输出和估计输出之间的误差,并计算出误差统计。误差统计可用于检测摩擦力的变化。这一点将在后面详细讨论。

 kurtosis(ySigEst(1,idx)-ySig(1,idx));  kurtosis(ySigEst(2,idx)-ySig(2,idx))];

扩展的卡尔曼滤波器性能

请注意,在10.45秒时检测到了一个摩擦变化。我们现在描述一下这个故障检测规则是如何得出的。首先检查仿真结果和过滤器的性能。

figure, plot(t,  Sig(1,:)  Sig(2,:));

模型的输入输出响应表明,很难直接从测量信号中检测出摩擦力的变化。扩展的卡尔曼滤能够估计状态,特别是摩擦状态。比较真实的模型状态和估计状态。估计的状态显示了对应于3个标准差的置信区间。

plot(t, True(1,:), t,  Est(1,:), ...

请注意,滤波器的估计值跟踪了真实值,而且置信区间仍然有界。检查估计误差可以更深入地了解滤波器。

plot(t,True(1,:)-Est(1,:)

误差图显示,滤波器在10秒的摩擦力变化后进行了调整,并将估计误差降低到了零。然而,误差图不能用于故障检测,因为它们依赖于对真实状态的了解。将测量的状态值与加速度和速度的估计状态值进行比较,可以提供一种检测机制。

plot(t,Sig(1,:-Est(1,:)

加速度误差图显示,在引入故障的10秒左右,平均误差有微小的差异。查看误差统计,看看是否可以从计算的误差中检测出故障。加速度和速度误差预计是正态分布的(噪声模型都是高斯的)。因此,加速度误差的峰度可能有助于识别由于摩擦力的变化和由此产生的误差分布从对称变为不对称的变化情况。

plot(t,Kur(1,:)

忽略估计器仍在收敛和收集数据的前4秒,误差的峰度相对稳定,在3(高斯分布的预期峰度值)附近有微小的变化。因此,在这个应用中,误差统计不能被用来自动检测摩擦力的变化。在这个应用中,使用误差的峰度也是很困难的,因为过滤器正在适应并不断地将误差推向零,只给出了一个误差分布与零不同的短暂时间窗口。

因此在这个应用中,使用估计的摩擦力的变化提供了自动检测电机故障的最好方法。来自已知无故障数据的摩擦力估计值(平均值和标准偏差)提供了摩擦力的预期界限,当这些界限被超过时,很容易检测出来。下面的图强调了这种故障检测方法。

plot(t,x,[nan t],[Mean+3*STD,Mean-3*STD]

摘要

这个例子展示了如何使用扩展的卡尔曼滤波器来估计一个简单的直流电动机的摩擦力,并使用摩擦力估计值进行故障检测。


这篇关于拓端tecdat|matlab实现扩展卡尔曼滤波(EKF)进行故障检测的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/154536

相关文章

python实现svg图片转换为png和gif

《python实现svg图片转换为png和gif》这篇文章主要为大家详细介绍了python如何实现将svg图片格式转换为png和gif,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录python实现svg图片转换为png和gifpython实现图片格式之间的相互转换延展:基于Py

Python利用ElementTree实现快速解析XML文件

《Python利用ElementTree实现快速解析XML文件》ElementTree是Python标准库的一部分,而且是Python标准库中用于解析和操作XML数据的模块,下面小编就来和大家详细讲讲... 目录一、XML文件解析到底有多重要二、ElementTree快速入门1. 加载XML的两种方式2.

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

Python实现图片分割的多种方法总结

《Python实现图片分割的多种方法总结》图片分割是图像处理中的一个重要任务,它的目标是将图像划分为多个区域或者对象,本文为大家整理了一些常用的分割方法,大家可以根据需求自行选择... 目录1. 基于传统图像处理的分割方法(1) 使用固定阈值分割图片(2) 自适应阈值分割(3) 使用图像边缘检测分割(4)

Android实现在线预览office文档的示例详解

《Android实现在线预览office文档的示例详解》在移动端展示在线Office文档(如Word、Excel、PPT)是一项常见需求,这篇文章为大家重点介绍了两种方案的实现方法,希望对大家有一定的... 目录一、项目概述二、相关技术知识三、实现思路3.1 方案一:WebView + Office Onl

C# foreach 循环中获取索引的实现方式

《C#foreach循环中获取索引的实现方式》:本文主要介绍C#foreach循环中获取索引的实现方式,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录一、手动维护索引变量二、LINQ Select + 元组解构三、扩展方法封装索引四、使用 for 循环替代

Spring Security+JWT如何实现前后端分离权限控制

《SpringSecurity+JWT如何实现前后端分离权限控制》本篇将手把手教你用SpringSecurity+JWT搭建一套完整的登录认证与权限控制体系,具有很好的参考价值,希望对大家... 目录Spring Security+JWT实现前后端分离权限控制实战一、为什么要用 JWT?二、JWT 基本结构

Java实现优雅日期处理的方案详解

《Java实现优雅日期处理的方案详解》在我们的日常工作中,需要经常处理各种格式,各种类似的的日期或者时间,下面我们就来看看如何使用java处理这样的日期问题吧,感兴趣的小伙伴可以跟随小编一起学习一下... 目录前言一、日期的坑1.1 日期格式化陷阱1.2 时区转换二、优雅方案的进阶之路2.1 线程安全重构2

Android实现两台手机屏幕共享和远程控制功能

《Android实现两台手机屏幕共享和远程控制功能》在远程协助、在线教学、技术支持等多种场景下,实时获得另一部移动设备的屏幕画面,并对其进行操作,具有极高的应用价值,本项目旨在实现两台Android手... 目录一、项目概述二、相关知识2.1 MediaProjection API2.2 Socket 网络