CN105446320B - 一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 - Google Patents
一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 Download PDFInfo
- Publication number
- CN105446320B CN105446320B CN201510955103.8A CN201510955103A CN105446320B CN 105446320 B CN105446320 B CN 105446320B CN 201510955103 A CN201510955103 A CN 201510955103A CN 105446320 B CN105446320 B CN 105446320B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- msubsup
- msup
- 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
- 239000000725 suspension Substances 0.000 title claims abstract description 71
- 238000001514 detection method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 78
- 239000013598 vector Substances 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 230000004888 barrier function Effects 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 238000006073 displacement reaction Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
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
- G05B23/0254—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 based on a quantitative model, e.g. mathematical relationships between inputs and outputs; functions: observer, Kalman filter, residual calculation, Neural Networks
-
- 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/26—Pc applications
- G05B2219/2604—Test of external equipment
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Vehicle Body Suspensions (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明一种基于有限频域的高速列车垂向悬挂系统的故障检测方法,包括以下步骤:1)采集高速列车垂向悬挂系统的状态方程和输出方程;2)构造误差方程;3)设计一种高速列车垂向悬挂系统的故障检测观测器;4)为了使得高速列车垂向悬挂系统对干扰具有鲁棒性,并对故障具有敏感性,引进H∞范数和H‑指标,同时给出一种迭代算法,保证得到的故障检测观测器增益矩阵L为最优结果;5)通过观察高速列车垂向悬挂系统实际输出向量y(t)和估计输出向量的残差判断是否发生故障;若t1时刻r(t1)≠0,则说明t1时刻发生故障,能及时检测出故障。本发明能够有效地克服高速列车在运行过程中的各种干扰信号,并准确、快速地解决高速列车垂向悬挂系统的低频微小渐变故障的检测问题。
Description
技术领域
本发明属于高速列车垂向悬挂系统的故障诊断技术领域,涉及一种基于有限频域的高速列车垂向悬挂系统的故障检测方法。
背景技术
基于高速列车高速性、安全性的特点,目前国家正在大力发展高速列车事业,到2015年底,中国高速铁路营业里程达1.8万公里以上,中国已经拥有全世界最大规模以及最高运营速度的高速铁路网,而悬挂系统是高速列车的重要组成部分,按照作用方向可分为高速列车垂向悬挂系统和横向悬挂系统,其支撑着车体和转向架,分别称为二系和一系悬挂系统,还起到缓冲由轨道不平顺所引起的轮轨作用力,控制列车行驶方向,保持舒适性等作用。因此,高速列车悬挂系统是监控列车运行状态的必要组成部分,而其故障检测对保障列车正常运行尤为重要,但由于其悬挂系统建模的复杂性或者模型的高维性使得对悬挂系统的故障检测研究比较少。
由于存在扰动噪声等干扰信号使得故障检测系统会出现误报的情况,而高速列车在运行过程中必然遇到各种干扰信号,所以必须对其进行鲁棒故障检测从而减小误报率。另外,现在大多数高速列车悬挂系统的故障诊断方法是在全频域下进行的,但微小渐变故障由于变化缓慢一般属于低频段,且初始值通常很小,现有的基于全频域的故障诊断方法无法快速地检测出故障,以上局限性也体现在高速列车垂向悬挂系统上。
发明内容
本发明的目的在于克服上述现有技术的不足之处,提出一种基于有限频域的高速列车垂向悬挂系统的故障检测方法,能够有效地克服高速列车在运行过程中的各种干扰信号,并准确、快速地解决高速列车垂向悬挂系统的低频微小渐变故障的检测问题。
本发明具体采用以下技术方案解决上述技术问题:
本发明一种基于有限频域的高速列车垂向悬挂系统的故障检测方法,其特征在于,包括以下步骤:
1)采集高速列车垂向悬挂系统的状态方程和输出方程;
2)构造误差方程和故障检测观测器增益矩阵;
3)根据采集到的高速列车垂向悬挂系统的状态方程和输出方程以及构造的误差方程和故障检测观测器增益矩阵,设计一种高速列车垂向悬挂系统的故障检测观测器;
4)为了使得被检测的高速列车垂向悬挂系统对干扰具有鲁棒性,并对故障具有敏感性,引进H∞范数和H-指标,求得故障检测观测器增益矩阵L;同时为了保证得到的故障检测观测器增益矩阵是最优的结果,给出一种迭代算法;
5)通过观察高速列车垂向悬挂系统实际输出向量y(t)和估计输出向量的残差来判断是否发生了故障;若t1时刻r(t1)≠0,则说明t1发生了故障,能够及时检测出故障。
在所述步骤2)中,所述的构造误差方程的过程是:
通过定义得到所述的构造误差方程为:
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,
在所述步骤3)中,所述的设计一种连续时间高速列车垂向悬挂系统的故障检测观测器具体如下:
其中,A,B,C为高速列车垂向悬挂系统的系统矩阵;u(t)∈Rp是系统输入向量,x(t)∈Rn,y(t)∈Rm分别表示系统的状态向量和输出向量,分别表示状态和输出的估计向量;r(t)是残差信号,L是待设计的观测器增益矩阵。
进一步的,在所述步骤4)中,为了使得被检测的高速列车垂向悬挂系统对干扰具有鲁棒性,对故障具有敏感性,引进H∞范数和H-指标,求得故障检测观测器增益矩阵,其具体过程为:
401)通过定义得到如下误差动态方程:
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,A,B1,B3,C,D1,D3为已知的具有合适维数的常值矩阵,A为高速列车垂向悬挂系统的系统矩阵,C为高速列车垂向悬挂系统的输出矩阵,B1,B3为高速列车垂向悬挂系统的故障分布矩阵,D1,D3为高速列车垂向悬挂系统的干扰分布矩阵,L是待设计的观测器增益矩阵;x(t)表示系统的状态向量,表示状态的估计向量,那么e(t)表示状态向量的估计误差向量,表示误差估计向量的时间导数,f(t)表示故障输入向量,d(t)表示干扰输入向量,r(t)是残差信号。
基于对干扰鲁棒,对故障敏感的设计原则,引入H∞范数和H-指数,使得下式成立:
其中,
σmin表示矩阵Grf(jω)的最小奇异值,σmax表示矩阵Grd(jω)的最大奇异值,采用H-指数描述故障对残差信号的最小影响,H∞范数描述故障对残差信号的最大影响,分别为故障信号和干扰信号的频率。
402)求取故障检测观测器增益矩阵,从而故障检测观测器满足上面的设计原则,给出定理如下:
给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1′,V2′,V3′,V4′使得下面线性矩阵不等式成立:
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
Ψ1=P2+X-XTA-V1′TC+YTC
Ψ2=-XTB3-V1′TD3+YTD3
Ψ4=CTV4′-XTB3-V3′TD3+YTD3
其中,Y=LTX,即故障检测观测器增益L=(YX-1)T,是已知的实数,代表故障和干扰信号的频率范围;在低频故障和干扰信号的情况下,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
进一步的,在所述步骤4)中,所述的为了使得到的故障检测观测器增益是最优的结果,给出一种迭代算法,其过程是:
403)给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1′,V2′,V3′,V4′使得下面线性矩阵不等式成立;
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
Ψ1=P2+X-XTA-V1′TC+YTC
Ψ2=-XTB3-V1′TD3+YTD3
Ψ4=CTV4′-XTB3-V3′TD3+YTD3
其中,Y=LTX,即故障检测观测器增益L=(YX-1)T,是已知的实数,指故障信号的频率范围,指干扰信号的频率范围。在低频故障和干扰信号的情况下,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
404)求故障检测观测器增益时,为了使J=β2/β1尽可能小,按照下面步骤进行:
404-1)求解线性矩阵不等式(1)(2)得到可行性解β1,β2;
404-2)选择合适的参数ε,δ,使得β1=β1+ε,β2=β2-δ,重新代入线性矩阵不等式(1)(2)中看是否有可行解,若有可行解迭代继续进行本步骤;
404-3)若步骤404-2)无可行解,则输出L,那么可以保证求得的故障检测观测器增益L即为使得性能最优的增益。
相比现有技术,本发明含有以下优点和有益效果:
1.由于微小渐变故障变化范围属于低频域,干扰有各种频域范围,目前常用的故障检测观测器是针对全频域范围设计的,没有考虑故障和干扰信号的频率范围,这就会存在很大的保守性问题,针对这个问题,本发明提出一种基于有限频域的高速列车悬挂系统的故障检测方法,尤其给出了低频故障情况下的公式,同时为了增强对故障检测的灵敏性和对干扰信号的鲁棒性,通过最优化方法找到最优指数,更好地满足性能要求,使得对微小渐变故障的检测更加准确和快速。
2.引入两个指标β1,β2,提出一种最优化方法使其满足最优性能,使得对高速列车悬挂系统的微小渐变故障检测更快速、更准确。
附图说明
图1为本发明的一种基于有限频域的高速列车垂向悬挂系统的故障检测方法的流程框图。
图2为本发明所述的一种高速列车垂向悬挂系统模型示意图。其中:变量φb为车体点头角位移,zb,zl,zt分别表示车体和前后转向架质心垂向位移,ul,ut为控制力,其余参数详见表1.
图3为本发明所述的一种高速列车垂向悬挂系统故障检测观测器结构图。
图4为本发明所述的一种高速列车垂向悬挂系统的故障检测报警示意图。
图5为本发明的一种实施例的(A-LC)特征根的分布图。
图6(a)、6(b)为本发明的一种实施例的两个输出向量的故障检测仿真曲线图,其中报警阈值是0.2。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
本发明的原理是:以高速列车垂向悬挂系统为研究对象,针对低频故障,提出了一种基于有限频的故障检测方法,引进两个指标β1,β2,通过优化指标使得J=β2/β1尽可能小,从而使得故障检测更加快速和准确。
本发明一种基于有限频域的高速列车垂向悬挂系统的故障检测方法,其特征在于,包括以下步骤:
1)采集高速列车垂向悬挂系统的状态方程和输出方程;
2)构造误差方程和故障检测观测器增益矩阵;
3)根据采集到的高速列车垂向悬挂系统的状态方程和输出方程以及构造的误差方程和故障检测观测器增益矩阵,设计一种高速列车垂向悬挂系统的故障检测观测器;
4)为了使得被检测的高速列车垂向悬挂系统对干扰具有鲁棒性,并对故障具有敏感性,引进H∞范数和H-指标,求得故障检测观测器增益矩阵L;同时为了保证得到的故障检测观测器增益矩阵是最优的结果,给出一种迭代算法;
5)通过观察高速列车垂向悬挂系统实际输出向量y(t)和估计输出向量的残差来判断是否发生了故障;若t1时刻r(t1)≠0,则说明t1发生了故障,能够及时检测出故障。
在所述步骤2)中,所述的构造误差方程的过程是:
通过定义得到所述的构造误差方程为:
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,
在所述步骤3)中,所述的设计一种连续时间高速列车垂向悬挂系统的故障检测观测器具体如下:
其中,A,B,C为高速列车垂向悬挂系统的系统矩阵;u(t)∈Rp是系统输入向量,x(t)∈Rn,y(t)∈Rm分别表示系统的状态向量和输出向量,分别表示状态和输出的估计向量;r(t)是残差信号,L是待设计的观测器增益矩阵。
进一步的,在所述步骤4)中,为了使得被检测的高速列车垂向悬挂系统对干扰具有鲁棒性,对故障具有敏感性,引进H∞范数和H-指标,求得故障检测观测器增益矩阵,其具体过程为:
401)通过定义得到如下误差动态方程:
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,A,B1,B3,C,D1,D3为已知的具有合适维数的常值矩阵,A为高速列车垂向悬挂系统的系统矩阵,C为高速列车垂向悬挂系统的输出矩阵,B1,B3为高速列车垂向悬挂系统的故障分布矩阵,D1,D3为高速列车垂向悬挂系统的干扰分布矩阵,L是待设计的观测器增益矩阵;x(t)表示系统的状态向量,表示状态的估计向量,那么e(t)表示状态向量的估计误差向量,表示误差估计向量的时间导数,f(t)表示故障输入向量,d(t)表示干扰输入向量,r(t)是残差信号。
基于对干扰鲁棒,对故障敏感的设计原则,引入H∞范数和H_指数,使得下式成立:
其中,σmin表示矩阵Grf(jω)的最小奇异值,σmax表示矩阵Grd(jω)的最大奇异值,采用H_指数描述故障对残差信号的最小影响,H∞范数描述故障对残差信号的最大影响,分别为故障信号和干扰信号的频率。
402)求取故障检测观测器增益矩阵,从而故障检测观测器满足上面的设计原则,给出定理如下:
给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1′,V2′,V3′,V4′使得下面线性矩阵不等式成立:
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
Ψ1=P2+X-XTA-V1′TC+YTC
Ψ2=-XTB3-V1′TD3+YTD3
Ψ4=CTV4′-XTB3-V3′TD3+YTD3
其中,Y=LTX,即故障检测观测器增益L=(YX-1)T,是已知的实数,代表故障和干扰信号的频率范围;在低频故障和干扰信号的情况下,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
进一步的,在所述步骤4)中,所述的为了使得到的故障检测观测器增益是最优的结果,给出一种迭代算法,其过程是:
403)给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1′,V2′,V3′,V4′使得下面线性矩阵不等式成立;
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
Ψ1=P2+X-XTA-V1′TC+YTC
Ψ2=-XTB3-V1′TD3+YTD3
Ψ4=CTV4′-XTB3-V3′TD3+YTD3
其中,Y=LTX,即故障检测观测器增益L=(YX-1)T,是已知的实数,代表故障和干扰信号的频率范围;在低频故障和干扰信号的情况下,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
404)求故障检测观测器增益时,为了使J=β2/β1尽可能小,按照下面步骤进行:
404-1)求解线性矩阵不等式(1)(2)得到可行性解β1,β2;
404-2)选择合适的参数ε,δ,使得β1=β1+ε,β2=β2-δ,重新代入线性矩阵不等式(1)(2)中看是否有可行解,若有可行解迭代继续进行本步骤;
404-3)若步骤404-2)无可行解,则输出L,那么可以保证求得的故障检测观测器增益L即为使得性能最优的增益。
以下为高速列车垂向悬挂系统的故障检测方法仿真实例:
设高速列车垂向悬挂系统的状态变量输入变量u(t)=[ul1 ul2 ut1 ut2]T,输出变量y(t)=[zb φb]T,则上面给出的四个微分方程可用状态空间方程表示为:
其中,变量φb为车体点头角位移,zb,zl,zt分别表示车体和前后转向架质心垂向位移,ul,ut为控制力,式中:
其中,在矩阵A中,有:
仿真模型各参数如下表1:
假设该系统发生传感器故障,故障矩阵及干扰分布矩阵如下:
令
d1(t)=d2(t)=sin(3t)e-0.2t,
建立一般故障模型如下:
为了检测故障,本发明设计了如下故障检测观测器:
其中,分别表示状态和输出的估计向量;r(t)是残差信号,L是待设计的观测器增益。
定义:则误差方程如下:
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,
给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1′,V2′,V3′,V4′使得下面线性矩阵不等式成立:
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
Ψ1=P2+X-XTA-V1′TC+YTC
Ψ2=-XTB3-V1′TD3+YTD3
Ψ4=CTV4′-XTB3-V3′TD3+YTD3
其中,Y=LTX,则故障检测观测器增益矩阵L=(YX-1)T,是已知的实数,代表故障和干扰信号的频率范围,指故障信号的频率范围,指干扰信号的频率范围。本发明针对的是低频故障和干扰信号的情况,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
假设故障为微小渐变故障,即频率为低频,利用MATLAB软件中的线性矩阵不等式工具箱可得:两个性能指标β1=1.005,β2=0.527,所以J=β2/β1=0.525,而基于全频域的故障检测观测器且没有采用优化算法的指标β1=0.843,β2=0.755,J=β2/β1=0.896,明显大于本发明提出算法下的性能指标要求,由此可见基于有限频的高速列车悬挂系统的故障检测更加准确和快速。
为了验证本发明故障检测的效果,进行以下仿真验证。假设传感器发生微小渐变故障,即在第10秒发生故障。
图5为A-LC的特征根分布图,其特征根全位于左半平面,可见实部均小于零,故设计的误差动态系统是稳定的。图6为基于故障检测观测器卡死故障的检测曲线,阈值取为0.2,由图可知,故障是在第10.8秒和第11秒时检测出来,可见能及时检测出微小渐变故障。其中,图6(a)为高速列车垂向悬挂系统车体质心垂向位移输出,图6(b)为高速列车垂向悬挂系统车体点头角位移输出。
本发明可以有效的进行故障检测,不仅解决了基于全频域保守性的问题,而且性能指标的优化使得对微小渐变故障的检测更加准确和快速。
Claims (2)
1.一种基于有限频域的高速列车垂向悬挂系统的故障检测方法,其特征在于,包括以下步骤:
1)采集高速列车垂向悬挂系统的状态方程和输出方程;
2)构造误差方程和故障检测观测器增益矩阵;
所述的构造误差方程的过程是:通过定义得到所述的构造误差方程:
<mrow>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,
3)根据采集到的高速列车垂向悬挂系统的状态方程和输出方程以及构造的误差方程和故障检测观测器增益矩阵,设计一种高速列车垂向悬挂系统的故障检测观测器;
所述的设计一种高速列车垂向悬挂系统的故障检测观测器的过程如下:
<mrow>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,A,B,C为高速列车垂向悬挂系统的系统矩阵;u(t)∈Rp是系统输入向量,x(t)∈Rn,y(t)∈Rm分别表示系统的状态向量和输出向量,分别表示状态和输出的估计向量;r(t)是残差信号,L是待设计的观测器增益矩阵;
4)为了使得被检测的高速列车垂向悬挂系统对干扰具有鲁棒性,并对故障具有敏感性,引进H∞范数和H-指标,求得故障检测观测器增益矩阵L;同时为了保证得到的故障检测观测器增益矩阵是最优的结果,给出一种迭代算法;
5)通过观察高速列车垂向悬挂系统实际输出向量y(t)和估计输出向量的残差来判断是否发生了故障;若t1时刻r(t1)≠0,则说明t1发生了故障,能够及时检测出故障;
在所述步骤4)中,为了使得被检测的高速列车垂向悬挂系统对干扰具有鲁棒性,对故障具有敏感性,引进H∞范数和H-指标,求得故障检测观测器增益矩阵,其具体过程为:
401)通过定义得到如下误差动态方程:
<mrow>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
r(t)=Ce(t)+D1f(t)+D3d(t)
其中,A,B1,B3,C,D1,D3为已知的具有合适维数的常值
矩阵,A为高速列车垂向悬挂系统的系统矩阵,C为高速列车垂向悬挂系统的输出矩阵,B1,B3为高速列车垂向悬挂系统的故障分布矩阵,D1,D3为高速列车垂向悬挂系统的干扰分布矩阵,L是待设计的观测器增益矩阵;x(t)表示系统的状态向量,表示状态的估计向量,那么e(t)表示状态向量的估计误差向量,表示误差估计向量的时间导数,f(t)表示故障输入向量,d(t)表示干扰输入向量,r(t)是残差信号;
基于对干扰鲁棒,对故障敏感的设计原则,引入H∞范数和H-指数,使得下式成立:
其中,
σmin表示矩阵Grf(jω)的最小奇异值,σmax表示矩阵Grd(jω)的最大奇异值,采用H-指数描述故障对残差信号的最小影响,H∞范数描述故障对残差信号的最大影响,指故障信号的频率范围,指干扰信号的频率范围;
402)求取故障检测观测器增益矩阵,从而故障检测观测器满足上面的设计原则,给出定理如下:
给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1',V2',V3',V4'使得下面线性矩阵不等式成立:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>X</mi>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<mi>A</mi>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<mi>A</mi>
<mo>-</mo>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<mi>C</mi>
<mo>+</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mi>Y</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
<mtd>
<msubsup>
<mi>V</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>I</mi>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<mi>C</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>4</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>&beta;</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>D</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
<mrow>
<msub>
<mi>&Phi;</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msubsup>
<mi>V</mi>
<mn>1</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mi>I</mi>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<mi>C</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>4</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>&beta;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>D</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Ψ1=P2+X-XTA-V1'TC+YTC
Ψ2=-XTB3-V1'TD3+YTD3
<mrow>
<msub>
<mi>&Psi;</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<msub>
<mi>B</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
其中,Y=LTX,即故障检测观测器增益L=(YX-1)T,是已知的实数,代表故障和干扰信号的频率范围;在低频故障和干扰信号的情况下,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
2.根据权利要求1所述的一种基于有限频域的高速列车垂向悬挂系统的故障检测方法,其特征在于,在所述步骤4)中,所述的为了使得到的故障检测观测器增益是最优的结果,给出一种迭代算法,其过程是:
403)给定性能指标β1,β2,系统矩阵A,B1,B2,B3,C,D1,D2,D3具有合适维数,如果存在对称正定矩阵P0,Q1,Q2,对称矩阵P1,P2和适合维数的矩阵X,Y,V1,V2,V3,V4,V1',V2',V3',V4'使得下面线性矩阵不等式成立;
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>X</mi>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<mi>A</mi>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<mi>A</mi>
<mo>-</mo>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<mi>C</mi>
<mo>+</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mi>Y</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>X</mi>
<mo>+</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
</mrow>
</mtd>
<mtd>
<msubsup>
<mi>V</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>I</mi>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<mi>C</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Phi;</mi>
<mn>4</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>&beta;</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>D</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Φ1=P1+X-XTA-V1 TC+YTC
Φ2=-XTB1-V1 TD1+YTD1
<mrow>
<msub>
<mi>&Phi;</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mi>T</mi>
</msup>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msubsup>
<mi>V</mi>
<mn>1</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mi>I</mi>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>3</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<mi>C</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>2</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&Psi;</mi>
<mn>4</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>&beta;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>D</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>V</mi>
<mn>4</mn>
<mrow>
<mo>&prime;</mo>
<mi>T</mi>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Ψ1=P2+X-XTA-V1′TC+YTC
Ψ2=-XTB3-V1'TD3+YTD3
其中,Y=LTX,即故障检测观测器增益L=(YX-1)T,是已知的实数,代表故障和干扰信号的频率范围;在低频故障和干扰信号的情况下,则存在故障检测观测器使得误差方程稳定并且满足下面有限频性能指标:
404)求故障检测观测器增益时,为了使J=β2/β1尽可能小,按照下面步骤进行:
404-1)求解线性矩阵不等式(1)(2)得到可行性解β1,β2;
404-2)选择合适的参数ε,δ,使得β1=β1+ε,β2=β2-δ,重新代入线性矩阵不等式(1)(2)中看是否有可行解,若有可行解迭代继续进行本步骤;
404-3)若步骤404-2)无可行解,则输出L,那么可以保证求得的故障检测观测器增益L即为使得性能最优的增益。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510955103.8A CN105446320B (zh) | 2015-12-17 | 2015-12-17 | 一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510955103.8A CN105446320B (zh) | 2015-12-17 | 2015-12-17 | 一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105446320A CN105446320A (zh) | 2016-03-30 |
CN105446320B true CN105446320B (zh) | 2018-02-16 |
Family
ID=55556640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510955103.8A Active CN105446320B (zh) | 2015-12-17 | 2015-12-17 | 一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105446320B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107168279B (zh) * | 2017-05-11 | 2019-03-12 | 浙江工业大学 | 一种基于h∞预演控制的车辆主动悬架系统的控制方法 |
CN113641104A (zh) * | 2021-08-23 | 2021-11-12 | 江南大学 | 动态事件触发下釜式反应器的有限频域故障检测方法 |
CN115128957B (zh) * | 2022-07-20 | 2023-03-14 | 西南交通大学 | 一种基于迭代学习的重载列车运行控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176159A (zh) * | 2011-02-28 | 2011-09-07 | 哈尔滨工业大学 | 一种基于状态观测器和等价空间的卫星姿控系统故障诊断装置及方法 |
CN103529830A (zh) * | 2013-11-05 | 2014-01-22 | 南京航空航天大学 | 一种基于有限频域飞行控制系统渐变故障的诊断设计方法 |
CN104155968A (zh) * | 2014-07-17 | 2014-11-19 | 南京航空航天大学 | 一种针对高速列车悬架系统执行器的微小故障诊断方法 |
CN104458298A (zh) * | 2014-12-09 | 2015-03-25 | 南京航空航天大学 | 基于多模型的高速列车悬架系统多执行器故障检测与隔离方法 |
CN104697804A (zh) * | 2015-03-24 | 2015-06-10 | 清华大学 | 用于检测与分离列车主动悬架系统间歇故障的方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6721683B2 (en) * | 2002-03-08 | 2004-04-13 | Insightek, Llc | Pump motor diagnosis |
-
2015
- 2015-12-17 CN CN201510955103.8A patent/CN105446320B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176159A (zh) * | 2011-02-28 | 2011-09-07 | 哈尔滨工业大学 | 一种基于状态观测器和等价空间的卫星姿控系统故障诊断装置及方法 |
CN103529830A (zh) * | 2013-11-05 | 2014-01-22 | 南京航空航天大学 | 一种基于有限频域飞行控制系统渐变故障的诊断设计方法 |
CN104155968A (zh) * | 2014-07-17 | 2014-11-19 | 南京航空航天大学 | 一种针对高速列车悬架系统执行器的微小故障诊断方法 |
CN104458298A (zh) * | 2014-12-09 | 2015-03-25 | 南京航空航天大学 | 基于多模型的高速列车悬架系统多执行器故障检测与隔离方法 |
CN104697804A (zh) * | 2015-03-24 | 2015-06-10 | 清华大学 | 用于检测与分离列车主动悬架系统间歇故障的方法及系统 |
Non-Patent Citations (1)
Title |
---|
基于观测器的传感器故障重构方法及其应用;何静 等;《兵工学报》;20090630;第30卷(第6期);第672-676页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105446320A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105446320B (zh) | 一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 | |
CN104155968B (zh) | 一种针对高速列车悬架系统执行器的微小故障诊断方法 | |
CN107329083A (zh) | 针对高速列车牵引电机非线性系统传感器的故障诊断方法 | |
Liu et al. | Sideslip angle estimation of ground vehicles: a comparative study | |
CN106202886A (zh) | 基于模糊粗糙集与决策树的轨道电路红光带故障定位方法 | |
CN103529830B (zh) | 一种基于有限频域飞行控制系统渐变故障的诊断设计方法 | |
CN104267716B (zh) | 一种基于多智能体技术的分布式飞行控制系统故障诊断设计方法 | |
CN104458297B (zh) | 基于非线性随机模型的列车悬架系统传感器故障检测方法 | |
CN103901882A (zh) | 一种列车动力系统在线监控故障诊断的系统及其方法 | |
CN104724098B (zh) | 一种城轨列车制动系统故障诊断方法 | |
CN103838959A (zh) | 偏最小二乘回归应用于配电网谐波源定位与检测中的方法 | |
CN103294895B (zh) | 基于证据推理的航迹航路分类方法 | |
CN103900835B (zh) | 基于模型与规则的高速列车转向架蛇形运动实时监测方法 | |
CN103675565B (zh) | 超/特高压交流输电线路工频参数测量方法 | |
CN104835102B (zh) | 列车测速方案评价方法和装置 | |
CN107315421A (zh) | 一种时延无人机编队的分布式速度传感器故障诊断方法 | |
CN102072789A (zh) | 一种地面测试铁道车辆轮轨力的连续化处理方法 | |
CN107453484A (zh) | 一种基于wams信息的scada数据校准方法 | |
CN107885927A (zh) | 一种铁路桥梁运营状态预警方法 | |
CN109766912A (zh) | 一种基于卡尔曼滤波与支持向量机的输电线路覆冰评估方法和系统 | |
CN105203130B (zh) | 一种基于信息融合的船舶组合导航系统故障诊断方法 | |
CN104965506B (zh) | 一种基于可调参数分布式飞行控制系统实时故障诊断方法 | |
CN103278744A (zh) | 一种基于t型线路等效及测后模拟思想的同塔双回输电线路纵联保护的方法 | |
CN104215379B (zh) | 斜坡式码头趸船系泊缆绳载荷监测预警系统及方法 | |
CN104155543B (zh) | 一种基于功率损耗法的输电线路参数辨识方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |