CN110850844A - 一种多维执行器广义最小可检测故障的计算方法 - Google Patents

一种多维执行器广义最小可检测故障的计算方法 Download PDF

Info

Publication number
CN110850844A
CN110850844A CN201911096771.4A CN201911096771A CN110850844A CN 110850844 A CN110850844 A CN 110850844A CN 201911096771 A CN201911096771 A CN 201911096771A CN 110850844 A CN110850844 A CN 110850844A
Authority
CN
China
Prior art keywords
discrete
fault
estimation error
actuator
lpv system
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.)
Granted
Application number
CN201911096771.4A
Other languages
English (en)
Other versions
CN110850844B (zh
Inventor
王学谦
闵波
徐峰
梁斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen International Graduate School of Tsinghua University
Original Assignee
Shenzhen International Graduate School of Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen International Graduate School of Tsinghua University filed Critical Shenzhen International Graduate School of Tsinghua University
Priority to CN201911096771.4A priority Critical patent/CN110850844B/zh
Publication of CN110850844A publication Critical patent/CN110850844A/zh
Application granted granted Critical
Publication of CN110850844B publication Critical patent/CN110850844B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/0208—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterized by the configuration of the monitoring system
    • G05B23/0213—Modular or universal configuration of the monitoring system, e.g. monitoring system having modules that may be combined to build monitoring program; monitoring system that can be applied to legacy systems; adaptable monitoring system; using different communication protocols
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明提供一种多维执行器广义最小可检测故障的计算方法，包括如下步骤：S1：利用龙伯格观测器建立离散LPV系统健康状态估计误差的动态方程；S2：利用所述离散LPV系统健康状态估计误差的动态方程构造健康状态估计误差的mRPI集合；S3：建立多维执行器故障同时发生时所述离散LPV系统状态估计误差的动态方程与集合描述；S4：构建所述离散LPV系统的健康残差集和故障残差集；S5：计算多维执行器故障同时发生时的广义最小可检测故障。给出多维执行器故障同时发生时系统状态估计误差的动态方程与集合描述的建立方法；多维执行器故障广义MDF的定义和计算方法。提供了单个执行器故障MDF计算的扩展，适用范围更广。

Description

一种多维执行器广义最小可检测故障的计算方法
技术领域
本发明涉及故障检测技术领域，尤其涉及一种多维执行器广义最小可检测故障的计算方法。
背景技术
由于现代工业中系统安全性和可靠性的要求日益增加，故障检测(FD)越来越受到学者们的广泛关注。故障的发生影响实际系统的正常运行，FD就是确定系统中是否发生了故障。不变集理论是一种重要的基于模型的FD方法，它主要关注检验系统实时测量的残差信号与通过模型计算出的健康残差集之间的一致性，它适用于模型已知的实际系统。在健康情况下，系统达到稳定状态时，测量的残差轨迹最终会收敛到健康残差集。一旦测量的残差信号超出了健康残差集，表明系统中发生了故障，此时，测量的残差轨迹最终会收敛到故障残差集。因此，只要健康残差集与故障残差集之间是分离的，系统在稳态阶段时就一定可以实现FD。
基于不变集的FD的核心在于构造健康和故障不变残差集，RPI(robustpositively invariant)集合是一种广泛使用的不变集，用来限定受有界不确定性干扰影响的系统状态；mRPI(minimal robust positively invariant)集合是包含在所有RPI集合内的最小的RPI集合。针对线性时不变(LTI:linear time-invariant)系统，已经有一些成熟的RPI集合构造方法；针对线性变参数(LPV:linear parameter varying)系统，有学者提出了一种新型的mRPI集合构造方法，并基于健康残差集与故障残差集分离的条件研究了单个加性执行器故障的最小可检测故障(MDF)计算问题。MDF是判断基于不变集的FD方法性能的指标，其计算问题具有重要意义。
在现有技术中，只研究了单个加性执行器故障MDF的计算方法。比如，采用多面体(polytope)集描述带有有界扰动的LPV系统的状态估计误差不变集，得到健康和故障不变残差集，通过分离健康和故障残差集计算出MDF，则基于不变集的FD方法可以检测出任何大于此MDF的加性执行器故障。基于不变集的LPV系统的状态估计误差mRPI集合构造方法只能描述单个执行器故障发生时状态估计误差的动态方程，无法描述多维执行器故障同时发生时状态估计误差的动态方程，适用范围有限；基于健康与故障残差集分离的单个加性执行器故障MDF计算方法只考虑了单个执行器故障的情况，无法解决多维执行器故障MDF的计算问题，适用范围有限。
发明内容
本发明为了解决现有的问题，提供一种多维执行器广义最小可检测故障的计算方法。
为了解决上述问题，本发明采用的技术方案如下所述：
一种多维执行器广义最小可检测故障的计算方法，包括如下步骤：S1：利用龙伯格观测器建立离散LPV系统健康状态估计误差的动态方程；S2：利用所述离散LPV系统健康状态估计误差的动态方程构造健康状态估计误差的mRPI集合；S3：建立多维执行器故障同时发生时所述离散LPV系统状态估计误差的动态方程与集合描述；S4：构建所述离散LPV系统的健康残差集和故障残差集；S5：计算多维执行器故障同时发生时的广义最小可检测故障。
优选地，建立所述LPV系统健康状态估计误差的动态方程包括：受多维乘性执行器故障影响的所述离散LPV系统为：
xk+1＝A(θk)xk+B(θk)Guk+Ewk
yk＝C(θk)xk+Fηk
其中，k表示离散时间，
Figure BDA0002268581440000021
和
Figure BDA0002268581440000022
分别是所述离散LPV系统矩阵、输入矩阵和输出矩阵，
Figure BDA0002268581440000023
是调度变量，包含在顶点θi生成的凸集Θ＝Conv{θ1,θ2,...,θN}中；
Figure BDA0002268581440000024
分别是所述离散LPV系统状态向量和输出向量；
Figure BDA0002268581440000025
是输入向量，和
Figure BDA0002268581440000027
是所述离散LPV系统的扰动和噪声，分别包含在已知凸集
Figure BDA0002268581440000028
Figure BDA0002268581440000029
中；
Figure BDA00022685814400000210
和是常参数矩阵；对角矩阵表征多维乘性执行器故障，
Figure BDA00022685814400000213
当第i个执行器发生故障时有0≤Gi<1，当所有执行器均健康时G是单位矩阵记为I；
利用如下所述龙伯格观测器：
Figure BDA0002268581440000031
Figure BDA0002268581440000032
其中，分别是估计状态向量和估计输出向量，
Figure BDA0002268581440000034
是观测器增益矩阵，在健康情况下，即G＝I时，所述离散LPV系统的状态估计误差为
Figure BDA0002268581440000035
其动态方程如下：
ek+1＝(A(θk)-LC(θk))ek+Ewk-LFηk。
优选地，构造健康状态估计误差的mRPI集合包括如下步骤：S21：构造初始凸集
Figure BDA0002268581440000036
假设所述离散LPV系统健康状态估计误差的动态方程是稳定的，则所述初始凸集
Figure BDA0002268581440000037
用如下公式获得：其中，存在ξ∈(0,1)，p*∈N，以及集合包含集合S，即
Figure BDA00022685814400000310
使得对于任意k≥p*，有
Figure BDA00022685814400000311
成立；S22：构造初始RPI集合
Figure BDA00022685814400000312
根据所述初始凸集
Figure BDA00022685814400000313
通过计算如下迭代等式：
Figure BDA00022685814400000314
其中，表示闵可夫斯基和，A(·)是一个集合映射函数，
Figure BDA00022685814400000316
Conv(·)表示集合的凸包，
Figure BDA00022685814400000317
若存在一个有限的k*∈N使得
Figure BDA00022685814400000318
则即是所述离散LPV系统健康状态估计误差的动态方程的一个初始RPI集合；S23：构造mRPI集合Ω∞：得到所述初始PRI集合
Figure BDA00022685814400000320
后，针对所述离散LPV系统健康状态估计误差的动态方程的集合序列Ωk：其中，
Figure BDA00022685814400000322
则在每一步迭代，Ωk都是RPI集合并且有：
Figure BDA00022685814400000323
此外，有：
Figure BDA00022685814400000324
Ω∞即是所述离散LPV系统健康状态估计误差的动态方程的mRPI集合。
优选地，建立多维执行器故障同时发生时所述离散LPV系统健康状态估计误差的动态方程包括：在多维执行器故障同时发生时，首先考虑不含扰动wk和噪声ηk的所述离散LPV系统：
Figure BDA0002268581440000041
其中，
Figure BDA0002268581440000042
Figure BDA0002268581440000043
并且在故障情况下G不是单位矩阵I；此时，相应的状态估计误差为
Figure BDA0002268581440000045
的动态方程为：
Figure BDA0002268581440000046
，因为：
Figure BDA0002268581440000047
其中，
Figure BDA0002268581440000048
ui是输入向量uk的第i个分量，又
Figure BDA0002268581440000049
所以有：
Figure BDA00022685814400000410
其中，Bi(θk)是矩阵B(θk)的第i列，则所述动态方程可转化为：
Figure BDA00022685814400000411
优选地，建立多维执行器故障同时发生时所述离散LPV系统健康状态估计误差的集合描述包括：由多维执行器故障同时发生时离散LPV系统状态估计误差的动态方程，根据叠加原理得到单个执行器故障下状态估计误差的动态方程为：
Figure BDA00022685814400000412
其中，假定第i个执行器发生了故障：
Figure BDA00022685814400000413
构造所述单个执行器故障下状态估计误差的动态方程的mRPI集合，假定分量ui包含在凸集
Figure BDA00022685814400000414
中，列向量Bi(θk)也包含在凸集中，则有：
Figure BDA00022685814400000416
假设所述单个执行器故障下状态估计误差的动态方程是稳定的，mRPI集合为
Figure BDA00022685814400000417
其中
Figure BDA00022685814400000418
根据不变集的性质，所述离散LPV系统达到稳定状态时有
Figure BDA00022685814400000419
叠加得：
Figure BDA00022685814400000420
结合式所述离散LPV系统健康状态估计误差的动态方程和所述多维执行器故障同时发生时离散LPV系统状态估计误差的动态方程，进一步得到带有扰动wk和噪声ηk的受多维乘性执行器故障影响的状态估计误差的动态方程：
Figure BDA0002268581440000051
其中，
Figure BDA0002268581440000052
相应的集合描述为：
Figure BDA0002268581440000053
其中，ε＝Ω∞表示式所述离散LPV系统健康状态估计误差的动态方程ek的mRPI集合，所述集合描述即为多维执行器故障同时发生时系统状态估计误差的集合描述。
优选地，所述离散LPV系统健康情况下的残差向量为：
Figure BDA0002268581440000054
所述健康残差集为：
Figure BDA0002268581440000055
其中，
Figure BDA0002268581440000056
优选地，多维执行器故障同时发生时所述离散LPV系统的残差向量为：
Figure BDA0002268581440000057
所述故障残差集为：
Figure BDA0002268581440000058
其中，Rf是多维执行器故障同时发生时的故障残差集，它依赖于所有的fi，i＝1,2,…,nu，定义多维执行器故障同时发生时故障向量
优选地，所述健康残差集R与所述故障残差集Rf之间不相交，即
Figure BDA00022685814400000510
的情况下，定义广义最小可检测故障为：
min||f||1
Figure BDA00022685814400000511
Figure BDA00022685814400000512
其中，||f||1表示故障向量f的1范数；上式等价于：
Figure BDA0002268581440000061
Figure BDA0002268581440000062
Figure BDA0002268581440000063
约束
Figure BDA0002268581440000064
描述了n维欧氏空间中的一个超立方体，优化问题的最优解描述了一个相对于各分量fi的超平面，记为其中λ即为多维执行器故障同时发生时的广义最小可检测故障；
Figure BDA0002268581440000066
的幅值大于所述广义最小可检测故障的幅值即认为所述离散LPV系统发生了故障。
优选地，所述广义最小可检测故障的优化问题转换为线性规划问题求解：
Figure BDA0002268581440000067
Figure BDA0002268581440000068
其中，
Figure BDA0002268581440000069
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上任一所述方法的步骤。
本发明的有益效果为：提供一种多维执行器广义最小可检测故障的计算方法，通过龙伯格观测器建立LPV系统健康状态估计误差的动态方程，建立mRPI集合，然后进行多维乘性执行器故障集理论分析，给出多维执行器故障同时发生时系统状态估计误差的动态方程与集合描述的建立方法；多维执行器故障广义MDF的定义和计算方法。提供了单个执行器故障MDF计算的扩展，适用范围更广。
附图说明
图1是本发明实施例中一种多维执行器广义最小可检测故障的计算方法的示意图。
图2是本发明实施例中建立状态估计误差ek的mRPI集合Ω∞的方法示意图。
图3是本发明实施例中nu＝3时的广义MDF的示意图。
图4是本发明实施例中健康残差集R与故障残差集Rf分离的示意图。
图5(a)和图5(b)是本发明实施例中在线FD的示意图。
图6是本发明实施例中健康残差集和故障残差集的示意图。
具体实施方式
为了使本发明实施例所要解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
需要说明的是，当元件被称为“固定于”或“设置于”另一个元件，它可以直接在另一个元件上或者间接在该另一个元件上。当一个元件被称为是“连接于”另一个元件，它可以是直接连接到另一个元件或间接连接至该另一个元件上。另外，连接既可以是用于固定作用也可以是用于电路连通作用。
需要理解的是，术语“长度”、“宽度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明实施例和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多该特征。在本发明实施例的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
如图1所示，本发明提供一种多维执行器广义最小可检测故障的计算方法，包括如下步骤：
S1：利用龙伯格观测器建立离散LPV系统健康状态估计误差的动态方程；
S2：利用所述离散LPV系统健康状态估计误差的动态方程构造健康状态估计误差的mRPI集合；
S3：建立多维执行器故障同时发生时所述离散LPV系统状态估计误差的动态方程与集合描述；
S4：构建所述离散LPV系统的健康残差集和故障残差集；
S5：计算多维执行器故障同时发生时的广义最小可检测故障。
在步骤S1中，考虑如下一般的受多维乘性执行器故障影响的离散LPV系统：
Figure BDA0002268581440000081
其中，k表示离散时间，
Figure BDA0002268581440000082
和
Figure BDA0002268581440000083
分别是离散LPV系统矩阵、输入矩阵和输出矩阵，是调度变量，包含在顶点θi生成的凸集Θ＝Conv{θ1,θ2,…,θN}中；
Figure BDA0002268581440000085
分别是离散LPV系统状态向量和输出向量；
Figure BDA0002268581440000086
是输入向量，
Figure BDA0002268581440000087
和
Figure BDA0002268581440000088
是离散LPV系统的扰动和噪声，分别包含在已知凸集
Figure BDA0002268581440000089
中；和
Figure BDA00022685814400000811
是常参数矩阵；对角矩阵
Figure BDA00022685814400000812
表征多维乘性执行器故障，当第i个执行器发生故障时有0≤Gi<1，当所有执行器均健康时G是单位矩阵记为I。
利用如下龙伯格FD观测器：
Figure BDA00022685814400000814
其中，
Figure BDA00022685814400000815
分别是估计状态向量和估计输出向量，
Figure BDA00022685814400000816
是观测器增益矩阵。在健康情况下(即G＝I)，离散LPV系统(1)的状态估计误差为
Figure BDA00022685814400000817
其动态方程如下：
ek+1＝(A(θk)-LC(θk))ek+Ewk-LFηk (3)
在步骤S2中，为了降低计算复杂度，以下所有的计算都是基于处理动态方程(3)的mRPI集合的凸包Ω∞进行的，为了方便，直接称Ω∞为mRPI集合。
如图2所示，针对动态方程(3)，采用如下步骤建立状态估计误差ek的mRPI集合Ω∞：
S21：构造初始凸集
Figure BDA00022685814400000818
假设动态方程(3)是稳定的，则初始凸集
Figure BDA00022685814400000819
可以用如下公式获得：
Figure BDA0002268581440000091
其中，存在ξ∈(0,1)，p*∈N，以及集合
Figure BDA0002268581440000092
包含集合S，即
Figure BDA0002268581440000093
使得对于任意k≥p*，有
Figure BDA0002268581440000094
成立。
S22：构造初始RPI集合
Figure BDA00022685814400000917
得到初始凸集
Figure BDA0002268581440000095
后，通过计算如下迭代等式：
Figure BDA0002268581440000096
其中
Figure BDA0002268581440000097
表示闵可夫斯基和，A(·)是一个集合映射函数，
Figure BDA0002268581440000098
Conv(·)表示集合的凸包，
Figure BDA0002268581440000099
若存在一个有限的k*∈N使得
Figure BDA00022685814400000910
则即是动态方程(3)的一个初始RPI集合。
S23：构造mRPI集合Ω∞；
得到初始PRI集合
Figure BDA00022685814400000912
后，针对动态方程(3)，集合序列Ωk：
其中,
Figure BDA00022685814400000914
则在每一步迭代，Ωk都是RPI集合并且有：
Figure BDA00022685814400000915
此外，有：
Figure BDA00022685814400000916
Ω∞即是动态方程(3)的mRPI集合，通过以上步骤可以得到。
在步骤S3中，建立多维执行器故障同时发生时所述离散LPV系统状态估计误差的动态方程包括：
在故障情况下，首先考虑如下不含扰动wk和噪声ηk的系统：
其中，
Figure BDA0002268581440000102
并且在故障情况下G不是单位矩阵I。此时，相应的状态估计误差为
Figure BDA0002268581440000103
Figure BDA0002268581440000104
的动态方程为：
Figure BDA0002268581440000105
为了将(10)进行转化，注意：
Figure BDA0002268581440000106
其中，
Figure BDA0002268581440000107
ui是输入向量uk的第i个分量。
又
Figure BDA0002268581440000108
所以有：
Figure BDA0002268581440000109
其中，Bi(θk)是矩阵B(θk)的第i列，则动态方程(10)可转化为：
Figure BDA00022685814400001010
建立多维执行器故障同时发生时所述离散LPV系统状态估计误差的集合描述包括：
事实上，式(13)表示多维执行器故障同时发生时状态估计误差的动态方程，根据叠加原理，单个执行器故障下状态估计误差的动态方程为：
Figure BDA00022685814400001011
其中，假定第i个执行器发生了故障，
Figure BDA00022685814400001012
我们首先考虑构造动态方程(14)的mRPI集合。假定分量ui包含在凸集
Figure BDA00022685814400001013
中，列向量Bi(θk)也包含在凸集
Figure BDA00022685814400001014
中，则有：
Figure BDA00022685814400001015
假设动态方程(14)是稳定的，根据步骤2的方法，(14)的mRPI集合为
Figure BDA0002268581440000111
其中
Figure BDA0002268581440000112
根据不变集的性质，系统达到稳定状态时有
Figure BDA0002268581440000113
叠加得：
Figure BDA0002268581440000114
结合式(3)和(13)，我们可以进一步得到带有扰动wk和噪声ηk的受多维乘性执行器故障影响的状态估计误差的动态方程：
Figure BDA0002268581440000115
其中，
Figure BDA0002268581440000116
相应的集合描述为：
Figure BDA0002268581440000117
其中，ε＝Ω∞表示式(3)ek的mRPI集合，式(18)即为多维执行器故障同时发生时系统状态估计误差的集合描述。
步骤S4中，建立健康残差集和故障残差集包括:
结合(1)和(2)，我们可以得到健康情况下的残差向量：
Figure BDA0002268581440000118
则健康残差集为：
Figure BDA0002268581440000119
其中,同样地，我们可以得到多维执行器故障同时发生时的残差向量：
相应的集合描述为：
Figure BDA00022685814400001112
其中,Rf是多维执行器故障同时发生时的故障残差集，它依赖于所有的fi，i＝1,2,…,nu，因此我们定义多维执行器故障同时发生时的故障向量
Figure BDA0002268581440000121
显然，单个执行器故障fi下的故障残差集为
Figure BDA0002268581440000122
它仅依赖于fi。根据基于不变集的FD的检测标准，我们需要实时检验rk∈R是否成立，如果某个时刻开始不成立(即
Figure BDA0002268581440000123
)，那么表明系统此时发生了故障；否则，我们相信系统仍然工作在健康状态下。一旦系统中发生了故障，那么残差向量rk会收敛到故障残差集Rf。因此，只要健康残差集R与故障残差集Rf之间不相交，即
Figure BDA0002268581440000124
就一定能对发生的故障实现准确的检测。
步骤S5中，考虑约束MDF就是满足该约束的最小故障，因此，给出广义MDF的定义如下：
Figure BDA0002268581440000126
其中，||f||1表示故障向量f的1范数，优化问题(23)等价于：
Figure BDA0002268581440000127
事实上，约束
Figure BDA0002268581440000128
描述了n维欧氏空间中的一个超立方体，优化问题(24)的最优解描述了一个相对于各分量fi的超平面，记为其中λ即为多维执行器故障同时发生时的广义MDF。
因为
Figure BDA00022685814400001210
所以有
Figure BDA00022685814400001211
其中
Figure BDA00022685814400001212
和
Figure BDA00022685814400001213
分别描述了上述超立方体的一个顶点。所有满足
Figure BDA00022685814400001214
的值表示超立方体内一系列平行的超平面。点
Figure BDA00022685814400001215
与点
Figure BDA00022685814400001216
之间的连线垂直于所有的超平面。对于广义MDF的含义，只要
Figure BDA0002268581440000131
的幅值大于广义MDF的幅值，就认为系统中发生了故障，它描述的是健康残差集R与故障残差集Rf恰好分离的边界。
如图3所示，nu＝3时的广义MDF的示意图。
优化问题(24)可以转化成一个简单的线性规划问题求解：
Figure BDA0002268581440000132
其中，
Figure BDA0002268581440000133
总的来说，优化问题(23)、(24)、(25)给出了多维执行器故障广义MDF的定义和计算方法，它包含了单个执行器故障MDF的计算方法。对于广义MDF，只要所有发生的故障幅值之和大于广义MDF的幅值，就一定能利用基于不变集的方法实现FD。
针对公式(1)的离散LPV系统，设计如下参数：
θk＝[θk(1)；θk(2)]，
A(θk)＝[0.85，0.2θk(2)；0.1θk(1)，0.75]
B(θk)＝[0.5678θk(1)，0.4265θk(1)，0.2464θk(1)；0.3492θk(2)，0.7347θk(2)，0.6356θk(2)]，
C(θk)＝[0.04θk(2)，0.3；0.05，0.01θk(1)]，
E＝[0.6324，0.2785；0.0975，0.5469]
F＝[0.8147，0.1270；0.9058.0.9134]
W＝{w∈R2|||w||∞≤0.02}
V＝{η∈R2|||η||∞≤0.02}
U1＝{μ∈R|1.5≤μ≤1.8}
U2＝{μ∈R|0.8≤μ≤1.2}
U3＝{μ∈R|1.0≤μ≤1.6}
Θ＝Conv{[0.5，0.5]T，[1，0.5]T，[1，1]T，[0.5，1]T}
L＝[0.5，0；0，1]
通过本发明的方法计算得到广义MDF为λ＝0.5426(f1＝0.5426,f2＝0,f3＝0)，相应的健康残差集R与故障残差集Rf分离结果如图4所示，图中左边是故障残差集Rf，右边是健康残差集R。
假定如下多维执行器故障情形：从k＝0到k＝40时刻，系统处于健康状态，在k＝41时刻给系统注入三个执行器故障f1＝0.2426,f2＝0.1,f3＝0.2，则在线FD结果如图5(a)和图5(b)所示，图5(a)和5(b)分别表示残差向量rk的两个分量rk(1)和rk(2)，对应健康残差集R的两个分量R(1)和R(2)，其中直线表示健康残差集R的上下界，曲线表示实时测量的残差信号rk，从图5(a)和图5(b)可以看出，k＝44时刻测量的残差信号rk超出了健康残差集R，表明系统此时发生了故障。
此外，可以得到相应的多维执行器故障同时发生时的故障残差集Rf以及单个执行器故障fi下的故障残差集为Ri，i＝1,2,3，如图6所示。从图6中可以看出，由于单个执行器故障f1，f2，f3都较小，健康集R(右边)和故障集Ri,i＝1,2,3(中间三个)都有交集，不满足FD的条件；而健康集R和故障集Rf(左边)没有交集，满足FD的条件，体现了广义MDF的优势，它可以检测出多个微小故障的叠加，只要所有发生的故障幅值之和大于广义MDF的幅值，就一定能利用基于不变集的方法实现FD。
本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干等同替代或明显变型，而且性能或用途相同，都应当视为属于本发明的保护范围。

Claims (10)

1.一种多维执行器广义最小可检测故障的计算方法，其特征在于，包括如下步骤：
S1：利用龙伯格观测器建立离散LPV系统健康状态估计误差的动态方程；
S2：利用所述离散LPV系统健康状态估计误差的动态方程构造健康状态估计误差的mRPI集合；
S3：建立多维执行器故障同时发生时所述离散LPV系统健康状态估计误差的动态方程与集合描述；
S4：构建所述离散LPV系统的健康残差集和故障残差集；
S5：计算多维执行器故障同时发生时的广义最小可检测故障。
2.如权利要求1所述的多维执行器广义最小可检测故障的计算方法，其特征在于，建立所述LPV系统健康状态估计误差的动态方程包括：
受多维乘性执行器故障影响的所述离散LPV系统为：
xk+1＝A(θk)xk+B(θk)Guk+Ewk
yk＝C(θk)xk+Fηk
其中，
Figure FDA0002268581430000011
和分别是所述离散LPV系统矩阵、输入矩阵和输出矩阵，
Figure FDA0002268581430000013
是调度变量，包含在顶点θi生成的凸集Θ＝Conv{θ1,θ2,...,θN}中；
Figure FDA0002268581430000014
分别是所述离散LPV系统的状态向量和输出向量；
Figure FDA0002268581430000015
是输入向量，和
Figure FDA0002268581430000017
是所述离散LPV系统的扰动和噪声，分别包含在已知凸集
Figure FDA0002268581430000018
Figure FDA0002268581430000019
中；
Figure FDA00022685814300000110
和是常参数矩阵；对角矩阵表征多维乘性执行器故障，
Figure FDA00022685814300000113
当第i个执行器发生故障时有0≤Gi<1，当所有执行器均健康时G是单位矩阵记为I；
利用如下所述龙伯格观测器：
Figure FDA00022685814300000115
其中，k表示离散时间，
Figure FDA00022685814300000116
分别是估计状态向量和估计输出向量，
Figure FDA0002268581430000021
是观测器增益矩阵，在健康情况下，即G＝I时，所述离散LPV系统的状态估计误差为
Figure FDA0002268581430000022
其动态方程如下：
ek+1＝(A(θk)-LC(θk))ek+Ewk-LFηk。
3.如权利要求2所述的多维执行器广义最小可检测故障的计算方法，其特征在于，构造健康状态估计误差的mRPI集合包括如下步骤：
S21：构造初始凸集
Figure FDA0002268581430000023
假设所述离散LPV系统健康状态估计误差的动态方程是稳定的，则所述初始凸集
Figure FDA0002268581430000024
用如下公式获得：
Figure FDA0002268581430000025
其中，存在ξ∈(0,1)，p*∈N，以及集合包含集合S，即
Figure FDA0002268581430000027
使得对于任意k≥p*，有
Figure FDA0002268581430000028
成立；
S22：构造初始RPI集合
Figure FDA0002268581430000029
根据所述初始凸集
Figure FDA00022685814300000210
通过计算如下迭代等式：
Figure FDA00022685814300000211
Figure FDA00022685814300000212
其中，
Figure FDA00022685814300000213
表示闵可夫斯基和，A(·)是一个集合映射函数，
Figure FDA00022685814300000214
Conv(·)表示集合的凸包，
Figure FDA00022685814300000215
若存在一个有限的k*∈N使得
Figure FDA00022685814300000216
则
Figure FDA00022685814300000217
即是所述离散LPV系统健康状态估计误差的动态方程的一个初始RPI集合；
S23：构造mRPI集合Ω∞：得到所述初始PRI集合
Figure FDA00022685814300000218
后，针对所述离散LPV系统健康状态估计误差的动态方程的集合序列Ωk：
Figure FDA00022685814300000219
其中，
Figure FDA00022685814300000220
则在每一步迭代，Ωk都是RPI集合并且有：
Figure FDA00022685814300000221
此外，有：
Figure FDA0002268581430000031
Ω∞即是所述离散LPV系统健康状态估计误差的动态方程的mRPI集合。
4.如权利要求3所述的多维执行器广义最小可检测故障的计算方法，其特征在于，建立多维执行器故障同时发生时所述离散LPV系统健康状态估计误差的动态方程包括：
在多维执行器故障同时发生时，首先考虑不含扰动wk和噪声ηk的所述离散LPV系统：
Figure FDA0002268581430000032
Figure FDA0002268581430000033
其中，
Figure FDA0002268581430000034
并且在故障情况下G不是单位矩阵I；此时，相应的状态估计误差为 的动态方程为：
Figure FDA0002268581430000037
因为：
Figure FDA0002268581430000038
其中，
Figure FDA0002268581430000039
ui是输入向量uk的第i个分量，又
Figure FDA00022685814300000310
所以有：
Figure FDA00022685814300000311
其中，Bi(θk)是矩阵B(θk)的第i列，则所述动态方程可转化为：
Figure FDA00022685814300000312
5.如权利要求4所述的多维执行器广义最小可检测故障的计算方法，其特征在于，建立多维执行器故障同时发生时所述离散LPV系统健康状态估计误差的集合描述包括：
由多维执行器故障同时发生时离散LPV系统状态估计误差的动态方程，根据叠加原理得到单个执行器故障下状态估计误差的动态方程为：
Figure FDA0002268581430000041
其中，假定第i个执行器发生了故障：
Figure FDA0002268581430000042
构造所述单个执行器故障下状态估计误差的动态方程的mRPI集合，假定分量ui包含在凸集
Figure FDA0002268581430000043
中，列向量Bi(θk)也包含在凸集中，则有：
Figure FDA0002268581430000045
假设所述单个执行器故障下状态估计误差的动态方程是稳定的，mRPI集合为其中
Figure FDA0002268581430000047
根据不变集的性质，所述离散LPV系统达到稳定状态时有
Figure FDA0002268581430000048
叠加得：
Figure FDA0002268581430000049
结合式所述离散LPV系统健康状态估计误差的动态方程和所述多维执行器故障同时发生时离散LPV系统状态估计误差的动态方程，进一步得到带有扰动wk和噪声ηk的受多维乘性执行器故障影响的状态估计误差的动态方程：
Figure FDA00022685814300000410
其中，
Figure FDA00022685814300000411
相应的集合描述为：
Figure FDA00022685814300000412
其中，ε＝Ω∞表示式所述离散LPV系统健康状态估计误差的动态方程ek的mRPI集合，所述集合的描述即为多维执行器故障同时发生时系统状态估计误差的集合描述。
6.如权利要求5所述的多维执行器广义最小可检测故障的计算方法，其特征在于，所述离散LPV系统健康情况下的残差向量为：
Figure FDA0002268581430000051
所述健康残差集为：
Figure FDA0002268581430000052
其中，
Figure FDA0002268581430000053
7.如权利要求6所述的多维执行器广义最小可检测故障的计算方法，其特征在于，多维执行器故障同时发生时所述离散LPV系统的残差向量为：
Figure FDA0002268581430000054
所述故障残差集为：
其中，Rf是多维执行器故障同时发生时的故障残差集，它依赖于所有的fi，i＝1,2,…,nu，定义多维执行器故障同时发生时故障向量
Figure FDA0002268581430000056
8.如权利要求7所述的多维执行器广义最小可检测故障的计算方法，其特征在于，所述健康残差集R与所述故障残差集Rf之间不相交，即
Figure FDA0002268581430000057
的情况下，定义广义最小可检测故障为：
min ||f||1
Figure FDA0002268581430000058
Figure FDA0002268581430000059
其中，||f||1表示故障向量f的1范数；
上式等价于：
Figure FDA00022685814300000510
Figure FDA00022685814300000511
Figure FDA00022685814300000512
约束
Figure FDA0002268581430000061
描述了n维欧氏空间中的一个超立方体，优化问题的最优解描述了一个相对于各分量fi的超平面，记为
Figure FDA0002268581430000062
其中λ即为多维执行器故障同时发生时的广义最小可检测故障；
的幅值大于所述广义最小可检测故障的幅值即认为所述离散LPV系统发生了故障。
9.如权利要求8所述的多维执行器广义最小可检测故障的计算方法，其特征在于，所述广义最小可检测故障的优化问题转换为线性规划问题求解：
Figure FDA0002268581430000064
其中，
Figure FDA0002268581430000066
10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-9任一所述方法的步骤。
CN201911096771.4A 2019-11-11 2019-11-11 一种多维执行器广义最小可检测故障的计算方法 Active CN110850844B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911096771.4A CN110850844B (zh) 2019-11-11 2019-11-11 一种多维执行器广义最小可检测故障的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911096771.4A CN110850844B (zh) 2019-11-11 2019-11-11 一种多维执行器广义最小可检测故障的计算方法

Publications (2)

Publication Number Publication Date
CN110850844A true CN110850844A (zh) 2020-02-28
CN110850844B CN110850844B (zh) 2020-09-22

Family

ID=69601329

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911096771.4A Active CN110850844B (zh) 2019-11-11 2019-11-11 一种多维执行器广义最小可检测故障的计算方法

Country Status (1)

Country Link
CN (1) CN110850844B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102980580A (zh) * 2012-11-16 2013-03-20 北京理工大学 基于张量积多胞鲁棒h2滤波的无陀螺卫星姿态确定方法
CN103149046A (zh) * 2013-03-07 2013-06-12 北京化工大学 一种基于专家思维的多维度故障诊断方法
CN105093933A (zh) * 2015-06-23 2015-11-25 北京航天自动控制研究所 一种确定lpv变增益控制器的方法
CN105094114A (zh) * 2015-06-23 2015-11-25 北京航天自动控制研究所 一种确定lpv控制器的稳定性的方法
CN107329083A (zh) * 2017-07-28 2017-11-07 南京航空航天大学 针对高速列车牵引电机非线性系统传感器的故障诊断方法
CN108803316A (zh) * 2018-03-09 2018-11-13 南京航空航天大学 针对多智能体系统执行器故障的主动容错控制方法
CN109308518A (zh) * 2018-09-13 2019-02-05 北京理工大学 一种基于概率神经网络的监测系统及其平滑参数优化方法
CN109508475A (zh) * 2018-10-17 2019-03-22 北京航空航天大学 一种基于多维卡尔曼滤波修正的故障主动预测方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102980580A (zh) * 2012-11-16 2013-03-20 北京理工大学 基于张量积多胞鲁棒h2滤波的无陀螺卫星姿态确定方法
CN103149046A (zh) * 2013-03-07 2013-06-12 北京化工大学 一种基于专家思维的多维度故障诊断方法
CN105093933A (zh) * 2015-06-23 2015-11-25 北京航天自动控制研究所 一种确定lpv变增益控制器的方法
CN105094114A (zh) * 2015-06-23 2015-11-25 北京航天自动控制研究所 一种确定lpv控制器的稳定性的方法
CN107329083A (zh) * 2017-07-28 2017-11-07 南京航空航天大学 针对高速列车牵引电机非线性系统传感器的故障诊断方法
CN108803316A (zh) * 2018-03-09 2018-11-13 南京航空航天大学 针对多智能体系统执行器故障的主动容错控制方法
CN109308518A (zh) * 2018-09-13 2019-02-05 北京理工大学 一种基于概率神经网络的监测系统及其平滑参数优化方法
CN109508475A (zh) * 2018-10-17 2019-03-22 北京航空航天大学 一种基于多维卡尔曼滤波修正的故障主动预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JUNBO TAN 等: "A Novel RPI Set Computation Method for Discrete-time LPV Systems with Bounded Uncertainties", 《THE 15TH IEEE INTERNATIONAL CONFERENCE ON CONTROL AND AUTOMATION (IEEE ICCA 2019)》 *
周萌: "基于多目标观测器设计的故障检测与分离", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Also Published As

Publication number Publication date
CN110850844B (zh) 2020-09-22

Similar Documents

Publication Publication Date Title
Toker et al. Tracking performance limitations in LTI multivariable discrete-time systems
US10611025B2 (en) Method an device for evaluating performance of industrial control loops based on full loop reconstruction simulations
Polycarpou et al. Learning methodology for failure detection and accommodation
Zhang et al. Integrated active fault-tolerant control using IMM approach
Eriksson et al. A method for quantitative fault diagnosability analysis of stochastic linear descriptor models
Rahnavard et al. Robust actuator and sensor fault reconstruction of wind turbine using modified sliding mode observer
CN110850844B (zh) 一种多维执行器广义最小可检测故障的计算方法
Tayari et al. Active fault tolerant control design for lpv systems with simultaneous actuator and sensor faults
de Oca et al. Fault estimation and virtual sensor FTC approach for LPV systems
Cristofaro et al. Fault‐tolerant control allocation for overactuated nonlinear systems
CN113820954B (zh) 一种广义噪声下复杂非线性系统的容错控制方法
Henry et al. A new multi-objective filter design for guaranteed robust FDI performance
Pereira et al. Bank of H∞ filters for sensor fault isolation in active controlled flexible structures
CN114237038A (zh) 基于复系数增益消除周期性扰动的故障检测系统和方法
Zhang et al. An anti-interference dynamic integral neural network for solving the time-varying linear matrix equation with periodic noises
Rodrigues et al. Design of an active fault tolerant control for nonlinear systems described by a multi-model representation
Zhao et al. An identification approach to nonlinear state space model for industrial multivariable model predictive control
Bakhshande Observer-based robust nonlinear control design
Ghanipoor et al. Robust fault estimators for nonlinear systems: An ultra-local model design
de Oca et al. Reliable Fault-Tolerant Control Design for LPV Systems using Admissible Model Matching
Tan et al. Computation of Minimal Detectable Fault under Hybrid Stochastic and Deterministic Framework
Chowdhury et al. Scalable Coherence in Large Scale Second-Order Networks Using High-Gain Observer
Zhang et al. Sensor Threat Isolation for Cyber-Physical Systems
de Oca et al. Fault-tolerant control design using LPV admissible model matching: application to a two-degree of freedom helicopter
Tang Linear Cyber-Physical System Security–Detection and Correction of Adversarial Attacks

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