CN110543162B - 一种强噪声下的运动控制系统多重故障辨识方法 - Google Patents
一种强噪声下的运动控制系统多重故障辨识方法 Download PDFInfo
- Publication number
- CN110543162B CN110543162B CN201910669435.8A CN201910669435A CN110543162B CN 110543162 B CN110543162 B CN 110543162B CN 201910669435 A CN201910669435 A CN 201910669435A CN 110543162 B CN110543162 B CN 110543162B
- Authority
- CN
- China
- Prior art keywords
- matrix
- attack
- observer
- representing
- motion control
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/24—Pc safety
- G05B2219/24065—Real time diagnostics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
一种强噪声下的运动控制系统多重故障辨识方法,对运动控制系统进行建模;得到系统状态空间方程并离散化;利用小波变换去噪;构建中间观测器;通过矩阵不等式求解中间观测器增益;通过中间观测器估计出攻击。构建观测器,并通过矩阵不等式求解中间观测器增益。通过与标称的中间观测器、H∞观测器分别进行2组对比,验证本方法的有效性。本发明的方法考虑了运动控制系统发生网络攻击的现象,并且本发明不只局限于此实例,同时其估计结果可以满足实际应用的精度与实时性要求。
Description
技术领域
本发明属于网络安全技术领域,具体针对强高斯白噪声情况下存在的运动控制系统多重故障辨识问题,提出了一种小波变换结合中间观测器的解决方法。
背景技术
随着“工业4.0”的提出以及工业化进程的加速,运动控制系统发挥的作用日益显著。在系统运行的过程中,不可避免地会受到故障的影响,此时若能对故障作出辨识,无疑有着重要的作用。在实际的工业生产中,由于受到影响较强的高斯白噪声的干扰,微弱的输出信号完全“淹没”在噪声中,从而致使观测器难以对故障作出准确估计。
在强噪声下对运动控制系统进行故障辨识,传统方法有着一定的局限性。在现有技术中,有鲁棒观测器、滑模观测器、中间观测器。这些观测器均具有一定的鲁棒性。然而标称的中间观测器,虽然具有一定的鲁棒性,但是仅仅通过调节参数,无法得到理想的效果。H∞按照比例调整,如果噪声太强,同样无法取得理想的效果。滑模观测器需要故障的先验信息,例如故障导数的上界、故障本身的上界,但是这些信息在实际情况中无法获取。
发明内容
基于上述问题,本发明提出了一种强噪声下的运动控制系统多重故障辨识方法,它偏向工程,更适用于实际工业情况。具体地说,它首先利用小波变换对带有强噪声的输出信号进行去噪处理;其次引入一个中间变量,构造一个中间观测器来同时估计状态和故障。同时与标称的中间观测器、H∞观测器分别进行2组对比实验,验证该方法的有效性和优越性。
本发明为解决上述技术问题提供了如下解决方案:
一种强噪声下的运动控制系统多重故障辨识方法,包括以下步骤:
步骤1),确定运动控制系统传递函数;
通过系统辨识,确定运动控制系统传递函数如式(1)所示:
其中G(s)为运动控制系统的传递函数,K、Ts为辨识出来的参数;
步骤2),建立运动控制系统状态空间方程并离散化,过程如下:
2.1)将上述传递函数转换为状态空间方程并对其离散化,并考虑系统中存在执行器攻击、传感器攻击以及白噪声的情况:
其中A为系统的状态矩阵,B为输入矩阵,x表示系统状态量,u为系统输入,wi表示白噪声,fa表示执行器攻击,fs表示传感器攻击,E和D表示攻击增益矩阵;
步骤3),利用小波变换对含有攻击及白噪声的输出信号进行去噪;
小波变换方法是依据数据自身的时间尺度特征,将一个复杂信号分解为若干频率不同的简单信号;利用小波变换对含有攻击及白噪声的输出信号进行分解,并剔除高频分量,将剩下的分量合并,即得到去噪后的输出信号;
步骤4),构建中间观测器并通过矩阵不等式求解观测器增益,过程如下:
4.1)定义中间变量如所示:
4.2)基于中间变量,设计中间观测器如(5)所示:
4.3)构建矩阵如式所示:
其中*表示对称元素,P1表示待求解的正定矩阵,H表示待求解的矩阵,P2表示带求解标量,w为调优参数,ε为给定标量,I表示单位阵,Π11,Π12,Π22表示中间矩阵:
4.4)求解矩阵不等式Π<0,得到P1、H,中间观测器增益L如式所示:
L=P-1H (7)
其中上标“-1”表示矩阵的逆,从而由中间观测器(5)实现对执行器攻击和传感器攻击的准确估计;
步骤5),构建H∞观测器并通过矩阵不等式求解观测器增益,过程如下:
5.1)设计鲁棒观测器为:
5.2)构建矩阵如式所示:
本发明一种强噪声下的运动控制系统多重故障辨识方法,通过对含有攻击及白噪声的输出信号进行小波变换去噪,并构建中间观测器,来对执行器攻击和传感器攻击进行估计。同时与标称的中间观测器、H∞观测器分别进行2组对比实验,验证该方法的有效性和优越性。
与现有鲁棒观测器技术相比,本发明的有益效果是:在强高斯白噪声影响下,其攻击辨识效果较鲁棒观测器精度更高,并且可以通过调节特定的调节参数来继续改善性能。
附图说明
图1a是利用小波变换结合中间观测器对执行器攻击fa的估计效果;
图1b是利用标称的中间观测器对执行器攻击fa的估计效果;
图2a是利用小波变换结合中间观测器对传感器攻击fs的估计效果;
图2b是利用标称的中间观测器对传感器攻击fs的估计效果;
图3a是利用小波变换结合中间观测器对执行器攻击fa的估计效果;
图3b是利用鲁棒观测器对执行器攻击fa的估计效果;
图4a是利用小波变换结合中间观测器对传感器攻击fs的估计效果;
图4b是利用鲁棒观测器对传感器攻击fs的估计效果。
具体实施方式
为使本发明的目的、技术方案和优点更加清晰,下面结合附图和实际实验对本发明的技术方案作进一步描述。
参照图1a~图4b,一种强噪声下的运动控制系统多重故障辨识方法,包括以下步骤:
1)确定运动控制系统传递函数;
2)建立运动控制系统状态空间方程并离散化;
3)利用小波变换对含有攻击及白噪声的输出信号进行去噪;
4)构建中间观测器并通过矩阵不等式求解观测器增益;
5)构建H∞观测器并通过矩阵不等式求解观测器增益,以作为效果的对比。
进一步,所述步骤1)中,确定运动控制系统传递函数:
通过系统辨识,确定运动控制系统传递函数如式(1)所示:
其中G(s)为运动控制系统的传递函数,K=0.08373,Ts=0.02433为辨识出来的参数。
再进一步,所述步骤2)中,建立运动控制系统状态空间方程并离散化,过程如下:
2.1)将上述传递函数转换为状态空间方程并对其离散化,并考虑系统中存在执行器攻击、传感器攻击以及白噪声的情况:
更进一步,所述步骤3)中,利用小波变换对含有攻击及白噪声的输出信号进行去噪:
对含有攻击及白噪声的输出信号进行小波变换分解,剔除高频噪声分量,将剩余分量重新整合,从而将部分白噪声从输出信号中剔除。
进一步,所述步骤4)中,构建中间观测器并通过矩阵不等式求解观测器增益,过程如下:
4.1)定义中间变量如所示:
4.2)基于中间变量,设计中间观测器如所示:
4.3)构建矩阵如式所示:
其中*表示对称元素,P1表示待求解的正定矩阵,H表示待求解的矩阵,P2表示带求解标量,调优参数ω=80,给定标量ε=1,I表示单位阵,Π11,Π12,Π22表示中间矩阵:
4.4)求解矩阵不等式Π<0,得到
中间观测器增益L如式所示:
L=P-1H (7)
所述步骤5)中,构建H∞观测器并通过矩阵不等式求解观测器增益的过程如下
5.1)设计鲁棒观测器为:
5.2)构建矩阵如式所示:
求解不等式(9),取r=33,得到
F=[133.4388 0.0643]
从实验结果可以看出,本发明中提出的方法能在强噪声的干扰下较为准确地对攻击进行实时估计,其运行结果可以满足实际应用的精度与实时性要求。
以上结合附图详细说明和陈述了本发明的实施方式,但并不局限于上述方式。在本领域的技术人员所具备的知识范围内,只要以本发明的构思为基础,还可以做出多种变化和改进。
Claims (1)
1.一种强噪声下的运动控制系统多重故障辨识方法,其特征在于,所述方法包括以下步骤:
步骤1),确定运动控制系统传递函数;
通过系统辨识,确定运动控制系统传递函数如式(1)所示:
其中G(s)为运动控制系统的传递函数,K、Ts为辨识出来的参数;
步骤2),建立运动控制系统状态空间方程并离散化,过程如下:
2.1)将上述传递函数转换为状态空间方程并对状态空间方程离散化,并考虑系统中存在执行器攻击、传感器攻击以及白噪声的情况:
其中A为系统的状态矩阵,B为输入矩阵,x表示系统状态量,u为系统输入,wi表示白噪声,fa表示执行器攻击,fs表示传感器攻击,E和D表示攻击增益矩阵;
步骤3),利用小波变换对含有攻击及白噪声的输出信号进行去噪;
小波变换方法是依据数据自身的时间尺度特征,将一个复杂信号分解为若干频率不同的简单信号,利用小波变换对含有攻击及白噪声的输出信号进行分解,并剔除高频分量,将剩下的分量合并,即可得到去噪后的输出信号;
步骤4),构建中间观测器并通过矩阵不等式求解观测器增益,过程如下:
4.1)定义中间变量如所示:
4.2)基于中间变量,设计中间观测器如(5)所示:
4.3)构建矩阵如式所示:
其中*表示对称元素,P1表示待求解的正定矩阵,H表示待求解的矩阵,P2表示带求解标量,w为调优参数,ε为给定标量,I表示单位阵,Π11,Π12,Π22表示中间矩阵:
4.4)求解矩阵不等式Π<0,得到P1、H,中间观测器增益L如式所示:
L=P-1H (7)
其中上标“-1”表示矩阵的逆,从而由中间观测器(5)实现对执行器攻击和传感器攻击的准确估计;
步骤5),构建H∞观测器并通过矩阵不等式求解观测器增益,过程如下:
5.1)设计鲁棒观测器为:
5.2)构建矩阵如式所示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910669435.8A CN110543162B (zh) | 2019-07-24 | 2019-07-24 | 一种强噪声下的运动控制系统多重故障辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910669435.8A CN110543162B (zh) | 2019-07-24 | 2019-07-24 | 一种强噪声下的运动控制系统多重故障辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110543162A CN110543162A (zh) | 2019-12-06 |
CN110543162B true CN110543162B (zh) | 2021-02-26 |
Family
ID=68709821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910669435.8A Active CN110543162B (zh) | 2019-07-24 | 2019-07-24 | 一种强噪声下的运动控制系统多重故障辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110543162B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112398840B (zh) * | 2020-11-06 | 2022-03-08 | 国网山东省电力公司 | 基于自适应滑模观测器的电力信息物理系统攻击定位方法 |
CN113093540B (zh) * | 2021-03-31 | 2022-06-28 | 中国科学院光电技术研究所 | 一种基于小波阈值去噪的滑模扰动观测器设计方法 |
CN113359438A (zh) * | 2021-05-18 | 2021-09-07 | 浙江工业大学 | 一种基于二维增益调节机制的两轴雕刻机故障估计方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540028A (zh) * | 2012-01-13 | 2012-07-04 | 广东电网公司电力科学研究院 | 发电机定子绝缘局部放电缺陷的自动模式识别方法 |
CN103529830A (zh) * | 2013-11-05 | 2014-01-22 | 南京航空航天大学 | 一种基于有限频域飞行控制系统渐变故障的诊断设计方法 |
CN104238357A (zh) * | 2014-08-21 | 2014-12-24 | 南京航空航天大学 | 一种近空间飞行器的容错滑模控制方法 |
CN105066996A (zh) * | 2015-07-20 | 2015-11-18 | 东南大学 | 自适应矩阵卡尔曼滤波姿态估计方法 |
CN108345711A (zh) * | 2018-01-09 | 2018-07-31 | 山东科技大学 | 基于事件驱动的动车组鲁棒传感器间歇故障检测方法 |
CN108629132A (zh) * | 2018-05-10 | 2018-10-09 | 南京邮电大学 | DoS攻击下的故障检测滤波器和控制器的协同设计方法 |
CN108628331A (zh) * | 2018-05-09 | 2018-10-09 | 北京航空航天大学 | 一种空间对抗环境下传感器受攻击的航天器姿控方法 |
CN108733031A (zh) * | 2018-06-05 | 2018-11-02 | 长春工业大学 | 一种基于中间估计器的网络控制系统故障估计方法 |
CN108803325A (zh) * | 2018-06-06 | 2018-11-13 | 黄山学院 | 永磁同步电机伺服系统鲁棒有限时间控制方法 |
CN109150639A (zh) * | 2018-11-05 | 2019-01-04 | 江南大学 | 一种高速率通信网络影响下时变系统的有限时域h∞控制方法 |
CN109557818A (zh) * | 2019-01-10 | 2019-04-02 | 南京航空航天大学 | 具有执行器和传感器故障的多智能体跟踪系统的滑模容错控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7480377B2 (en) * | 2003-12-31 | 2009-01-20 | Intel Corporation | Dual adaptive filter apparatus and method |
-
2019
- 2019-07-24 CN CN201910669435.8A patent/CN110543162B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540028A (zh) * | 2012-01-13 | 2012-07-04 | 广东电网公司电力科学研究院 | 发电机定子绝缘局部放电缺陷的自动模式识别方法 |
CN103529830A (zh) * | 2013-11-05 | 2014-01-22 | 南京航空航天大学 | 一种基于有限频域飞行控制系统渐变故障的诊断设计方法 |
CN104238357A (zh) * | 2014-08-21 | 2014-12-24 | 南京航空航天大学 | 一种近空间飞行器的容错滑模控制方法 |
CN105066996A (zh) * | 2015-07-20 | 2015-11-18 | 东南大学 | 自适应矩阵卡尔曼滤波姿态估计方法 |
CN108345711A (zh) * | 2018-01-09 | 2018-07-31 | 山东科技大学 | 基于事件驱动的动车组鲁棒传感器间歇故障检测方法 |
CN108628331A (zh) * | 2018-05-09 | 2018-10-09 | 北京航空航天大学 | 一种空间对抗环境下传感器受攻击的航天器姿控方法 |
CN108629132A (zh) * | 2018-05-10 | 2018-10-09 | 南京邮电大学 | DoS攻击下的故障检测滤波器和控制器的协同设计方法 |
CN108733031A (zh) * | 2018-06-05 | 2018-11-02 | 长春工业大学 | 一种基于中间估计器的网络控制系统故障估计方法 |
CN108803325A (zh) * | 2018-06-06 | 2018-11-13 | 黄山学院 | 永磁同步电机伺服系统鲁棒有限时间控制方法 |
CN109150639A (zh) * | 2018-11-05 | 2019-01-04 | 江南大学 | 一种高速率通信网络影响下时变系统的有限时域h∞控制方法 |
CN109557818A (zh) * | 2019-01-10 | 2019-04-02 | 南京航空航天大学 | 具有执行器和传感器故障的多智能体跟踪系统的滑模容错控制方法 |
Non-Patent Citations (2)
Title |
---|
Distributed Attack Detection and Secure Estimation of Networked Cyber-Physical Systems Against False Data Injection Attacks and Jamming Attacks;Yanpeng Guan;《IEEE TRANSACTIONS ON SIGNAL AND INFORMATION PROCESSING OVER NETWORKS》;20180331;第48-59页 * |
基于扩张状态观测器的轮式移动机器人全阶滑模控制;郭一军等;《重庆邮电大学学报(自然科学版)》;20170615(第03期);第382-388页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110543162A (zh) | 2019-12-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110543162B (zh) | 一种强噪声下的运动控制系统多重故障辨识方法 | |
Wu et al. | Indirect adaptive robust control of nonlinear systems with time‐varying parameters in a strict feedback form | |
CN110531616B (zh) | 一种有色噪声下的网络化运动控制系统攻击辨识方法 | |
Xiao et al. | Stability of stochastic delay switched neural networks with all unstable subsystems: A multiple discretized Lyapunov-Krasovskii functionals method | |
CN109033021B (zh) | 一种基于变参收敛神经网络的线性方程求解器设计方法 | |
Zheng et al. | Asynchronous nonfragile H∞ filtering for discrete-time nonlinear switched systems with quantization | |
CN116992927A (zh) | 一种欺骗攻击下的无线传感器网络分布式保概率估计方法 | |
Wang et al. | Adaptive filtering-based recursive identification for time-varying Wiener output-error systems with unknown noise statistics | |
Song et al. | A method for simulation model validation based on Theil’s inequality coefficient and principal component analysis | |
Wang et al. | A variable step size for maximum correntropy criterion algorithm with improved variable kernel width | |
CN110311898B (zh) | 基于高斯径向基函数分类器的网络化数控系统中间人攻击检测方法 | |
Ding et al. | Adaptive modified input and state estimation for linear discrete-time system with unknown inputs | |
Lima et al. | On the steady-state MSE performance of the set-membership NLMS algorithm | |
Xiong et al. | Improved stabilization for continuous dynamical systems with two additive time‐varying delays | |
Chen et al. | Refined instrumental variable identification of continuous-time OE and BJ models from irregularly sampled data | |
Yu et al. | Observability forms for switched systems with zeno phenomenon or high switching frequency | |
Lyu et al. | Variance-constrained resilient H∞ filtering for time-varying nonlinear networked systems subject to quantization effects | |
CN114326382A (zh) | 一种具有欺骗攻击的自适应弹性跟踪控制方法及系统 | |
Zhao et al. | Robust transcale state estimation for multiresolution discrete‐time systems based on wavelet transform | |
Zhang et al. | Consistent parameter estimation and convergence properties analysis of Hammerstein output-error models | |
Xia et al. | Application of Kalman filter in microseismic data denoising based on identified signal model | |
CN114019944B (zh) | 一种fdi攻击下网络化控制系统状态与故障的联合区间估计方法 | |
JP2020027245A5 (ja) | 情報処理方法、情報処理装置およびプログラム | |
CN117033910B (zh) | 一种海面高精度信号的处理方法及系统 | |
Qin | Estimation and identification of nonlinear parameter of motion index based on least squares algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |