范婷婷,牟宗磊,何 童
(1.山东科技大学 电气与自动化工程学院,山东 青岛 266590;2.中央财经大学 统计与数学学院,北京 102206)
安全性和可靠性是衡量现代工业系统性能指标的重要因素,近些年来随着现代工业系统安全生产要求不断提高,故障诊断与容错控制技术受到各界的广泛关注[1-2]。故障类型不同,运用的故障诊断与容错控制方法也不尽相同。根据故障发生位置,可将系统故障分为执行器故障、传感器故障和元部件故障等。对于含有执行器或传感器故障的系统,可以通过增加辅助变量来构造观测器[3],或将系统扩维后设计观测器[4-5],在此基础上利用反步法或重构控制率来实现故障容错[6-9]。根据故障的持续时间不同,系统故障又可分为永久故障、瞬时故障和间歇故障,其中间歇故障是工业系统中发生概率较高的一类,但是间歇故障容错控制方法的研究相对较少[10]。文献[11]利用增广卡尔曼滤波器来同时估计线性离散随机系统状态和间歇故障;文献[12]运用等价空间方法将残差和扰动解耦,通过假设检验来提高间歇故障的检测精度;文献[13]利用双层卡尔曼滤波来联合估计系统状态、无故障约束状态和间歇故障的故障诊断等。但上述工作均未对发生故障的系统作进一步的控制设计,应针对间歇故障进行容错控制研究。
此外,工业过程大多采用周期性采样和传输技术来控制系统,这种数据采集和传递方式既增加通讯消耗又导致能量和计算资源的浪费,因此事件触发控制机制成为传统周期性采样技术的替代和补充[14]。相关研究工作包括在模型预测控制的基础上添加事件触发机制,减少控制序列的更新次数,从而减轻网络传输负担,提高计算优越性[15-17];在含有执行器故障情况下,通过事件触发机制来降低通讯频率,并重构控制率使系统满足一致最终有界[18];针对含有外部扰动、参数不确定性和未知执行器故障的卫星姿态控制系统,通过事件触发机制来降低信息传输频率,减轻网络传输负担,并通过反步法来补偿故障和外部扰动带来的系统影响[19]。但上述工作均未考虑传感器故障情况。
本研究以含有间歇传感器故障的线性系统为研究对象,通过事件触发主动容错控制技术优化系统的故障容错性能。主要创新点在于对含有不同间歇传感器故障模态的系统,通过求解多个线性矩阵不等式,使系统在不同故障模态下具有公共的观测器、控制器增益以及事件触发采样参数,降低了含有事件触发机制的容错控制设计难度,事件触发的容错控制策略可以针对间歇故障频繁发生和消失的特点,在故障发生和消失时增加控制率调整的频率,改善系统控制效果,并有效缓解系统网络通讯压力。
符号说明:Rn表示n维欧式范数,diag(p1,p1,…,pm)表示对角矩阵,AT表示矩阵A的转置,A-1表示矩阵A的逆,I和0表示适当维度的单位矩阵和零矩阵。
考虑如下含有间歇传感器故障的线性离散系统:
(1)
式中:x(k)∈Rn,u(k)∈Rq和y(k)∈Rm分别为系统状态、输入向量和输出向量;f(k)为传感器测量处发生的常值间歇故障;Fi为不同故障模态下的方向系数矩阵;A、B、C为已知合适维度的参数矩阵;i∈[1,2,…,M]为故障模态,共有M种类型的故障模态,对αi给定以下假设:
假设1αi是只能取0或1的已知标量,且每个时刻最多只有αi=1,其余αj=0(j≠i)。
为减少通讯负担,在输出传感器端引入事件触发机制,事件触发的条件为:
‖y(k*)-y(k)‖>δ‖y(k)‖,
(2)
式中,y(k)是当前时刻的传感器测量值,y(k*)是上一次触发时刻的传输值,δ(0<δ<1)是待设计的事件触发采样参数。当式(2)成立时,将当前时刻的传感器测量值y(k)传递至观测器;当式(2)不成立时,则将上一次触发时刻输出值传递至观测器。
假设2传感器故障幅值有界且已知,即f(k)≤θ。
首先构造故障诊断观测器
(3)
(4)
式中,K为待设计的控制率参数矩阵。
定义
则根据式(1)和式(3)可以得到误差系统如下:
(5)
考虑到间歇故障在每个故障间隔期间的幅值不变,则只有在间歇故障发生和消失的瞬间,间歇故障增量Δf(k)=f(k+1)-f(k)不为零,其余时间Δf(k)均恒等于零。将控制率式(4)代入系统状态方程,并将状态方程和误差系统式(5)合并,扩展成增广系统:
(6)
定理1如果存在对称矩阵P>0,矩阵N1、N2和标量β>0满足以下矩阵不等式:
(7)
证明:设计李雅普诺夫函数
V(k)=ηT(k)Pη(k),
(8)
则
(9)
将状态方程代入式(2)并结合范数不等式得:
(10)
又因为
(11)
(12)
将式(12)代入式(9),并根据假设2可将式(9)写成:
(13)
即
(14)
其中,
(15)
如果Φi<0,通过求解M种线性矩阵不等式能够得到公共解P,则可以证明增广系统在无故障时渐近稳定,在发生间歇故障时一致最终有界。式(15)等价于:
(16)
式(16)由Schur补[20]可得:
(17)
利用无人艇线性系统来验证上述理论推导的有效性,采样时间T=0.2 s,仿真时间40 s,给定系统各参数矩阵如下[21]:
仿真中,把不同幅值的间歇传感器故障归为同一模态类型,没有发生任何故障的正常系统状态看作另一种模态类型,即假设仿真中仅有发生传感器故障和未发生故障两种模态类型,给出两种模态下的方向矩阵分别为:
故障发生和消失的间隔以及不同时刻的故障幅值描述如式(18)所示。
(18)
通过仿真平台LMI工具箱求得的事件触发采样参数为δ=0.506,各增益矩阵为:
仿真结果如图1、图2所示,其中图1是真实故障幅值与观测器的故障估计值,图2是故障的估计误差。
图1 传感器故障与故障估计Fig. 1 The estimation of sensor fault
图2 真实故障与估计误差Fig. 2 The error of fault estimation
对比图1、图2可以看出,估计值在故障幅值发生改变的短时间内误差较大,但依然能够快速跟上,因此本研究设计的观测器有良好的估计效果。
图3、图4分别是加入设计的控制率后的系统状态及系统状态误差。
图3 容错控制状态Fig. 3 The system states with FTC
图4 真实状态与估计误差Fig. 4 The error of states estimation
分析图3、图4可以看出,系统稳定之后,状态只在每次间歇故障发生和消失的瞬间即系统模态发生改变的瞬间有所波动,在无故障和故障间隔内,状态都能够很好地稳定在零点。
由图5分析可知:系统开始工作的0~4 s内,振荡较大,触发次数较多;4~8 s内发生故障,触发次数依然较多;8~14 s、18~24 s、28~34 s以及38~40 s内系统稳定且没有故障发生,触发次数明显减少;14~18 s、24~28 s、34~38 s内发生故障,触发次数相对较多。可以看出,加入事件触发机制后系统信息的传输频率明显减少,网络的通讯压力有效降低。
图5 事件触发采样图Fig. 5 The event-triggered sampling
综上,本控制方法对故障有更加快速和优越的镇定效果,且加入的事件触发机制有效缓解了系统网络的通讯压力。
本研究提出一种基于事件触发线性离散系统间歇传感器故障的估计和主动容错控制方法。通过将状态、状态及故障估计误差进行扩维,利用线性矩阵不等式同时解得公共且合适的观测器、控制器增益以及事件触发采样参数,使得系统在发生不同模态的间歇故障时,灵活改变控制率重构频率,减少控制器反复调整的次数,缓解控制器性能损耗速度,有效降低了传感器到观测器的传输频率,减轻了通讯压力。仿真实例也证明该方法在估计和控制性能上的效果。下一步将对故障特性和处理方法做更合理和宽泛的优化,并在此基础上考虑非线性系统的间歇故障诊断与容错控制问题,研究加入多个触发机制的事件触发控制问题。