基于非光滑观测器的含有死区和迟滞的复合三明治系统状态估计方法
【技术领域】
[0001] 本发明属于非线性系统的状态估计领域，特别涉及一种基于非光滑观测器的含有 死区和迟滞的复合三明治系统状态估计方法。
【背景技术】
[0002] 对于任何控制系统来说，无论是实现系统的闭环极点配置，还是实现系统的解耦 都离不开状态反馈。同时，准确估计系统的状态也是实现系统最优控制和故障诊断的重要 依据。然而，系统的状态并不都是易于直接能检测到的，有些状态变量甚至根本无法检测 至IJ。因此，研究者常常通过构造状态观测器来解决以上问题。但是，对于特定系统，如果不 根据系统实际情况构造观测器进行状态估计，必将导致状态估计的不准确。依据不准确的 估计状态无法实现对系统的准确控制，在故障诊断中，不准确的状态估计也会导致故障的 误报和漏报。因此，构建能准确估计特定系统状态的状态估计观测器和能够实现准确故障 预报的故障预报观测器显得尤为重要。
[0003] 死区、间隙和迟滞是最常见的非光滑非线性特性。死区广泛存在于直流电机，机械 传动系统、液压传动系统、机电一体化系统和功放电路系统中。间隙广泛存在于齿轮机械 传动系统、电动阀、数字电路、传感器和液压系统中。近年来，采用记忆合金、压电陶瓷的智 能执行装置由于具有定位精度高、驱动力大和响应快等优点，被广泛应用在精密加工机床、 航天飞机的柔性机械手臂和天文望远镜等精密设备的定位系统中，但是记忆合金和压电陶 瓷中存在着迟滞这种非光滑非线性，因此，迟滞广泛存在于精密加工机床、航天飞机的柔性 机械手臂和天文望远镜等精密设备的定位系统中。另外，通信系统中常用的射频电路和锁 相环路也包含有迟滞三明治系统。而在工程实际中，死区、间隙和迟滞非线性特性往往不是 孤立的，而是与其他传统环节相联接，夹在两个线性动态环节之间，被称为非光滑三明治系 统。带死区、间隙和迟滞的非光滑三明治系统可统一用图1所示的结构框图来描述。
[0004] 但是，在实际的工程中，很多系统并非是如图1所示的典型非光滑三明治系统，例 如，一个柔性机械手臂由放大电路、直流电机、压电陶瓷和机器手臂的执行机构组成，同样 放大电路可以看作一个线性环节L1，直流电机存在死区看作一个死区环节DZ，压电陶瓷存 在迟滞特性可以看作一个迟滞环节HS，最后，机器手臂的执行机构可以看作是一个线性环 节L2。因此，整个系统可以表示为图2所示的结构，它的非光滑环节由一个死区串联一个迟 滞组成。
[0005] 为此，本发明定义了一类新的非光滑三明治系统：如果非光滑三明治系统的中间 环节不止一个，而是多个非光滑环节串联在一起，这样的非光滑三明治系统称作复合非光 滑三明治系统，图2所示的系统就是含有死区和迟滞的复合三明治系统。
[0006] 针对特定系统构造相应的观测器一直都是控制工程领域的一个研究热点。自 D.J.Luenberger(1971)在本世纪70年代提出著名的Luenberger观测器以来，对于定常线 性系统观测器的设计理论和方法就已成熟。但对于非线性系统却不同，首先，非线性系统的 能观性是一个局部特性；第二，线性系统的能观性与系统输入无关，只取决于系统本身的结 构，而非线性系统的能观性不但与系统结构有关，同时还与系统输入有关。正是由于非线性 系统的复杂性，对于非线性系统我们很难找到一个统一的观测器构造方法，往往是针对某 一类非线性系统构造一类特定的观测器。
[0007] 在特殊系统观测器设计方面，马克茂，马萍（2003)对满足Lipschitz条件的非线 性系统开展了观测器设计研究。朱芳来（2004)对类Lyapunov非线性系统设计了降维观测 器进行状态估计。周少武，章兢（2001)提出一种改进的Lipschitz非线性系统观测器的设 计方法。栾小丽等（2008)构造了非线性不确定时滞系统的观测器，利用该观测器对该类系 统进行鲁棒控制。韩春艳等（2009)针对随机时滞系统构建特定观测器进行状态估计。吴敏 等（2008)完成了含不可观变迀事件图的状态估计器设计。A.Lj.Juloski.等（2002,2003) 针对分段线性系统设^ -种Luenberger型切换观测器进行状态估计。A.Alessandri等 (2001)利用共用Lyapunov函数方法求解分段线性切换系统的观测器增益矩阵。周祖鹏，谭 永红（2011，2012)对单一的带死区、间隙和迟滞三明治系统设计了非光滑观测器进行了状 态估计。
[0008] 但是，以上状态估计方法都是针对特定系统开展观测器设计，在这些方法中都假 设系统是完全能观的或是给定只有一个非光滑环节，而复合非光滑三明治系统在某工作区 间不完全能观且具有两个串联的非光滑环节。因此，在对复合非光滑三明治系统进行状态 估计观测器设计时，必须同时要考虑系统不完全能观和复合后的新问题。因此，本发明将拓 展复杂系统状态观测器理论，解决含有死区和迟滞的复合三明治系统状态变量的准确状态 估计问题，本发明具有一定的科学理论意义与实际应用前景。
【发明内容】
[0009] 针对上述技术的不足，本发明公开了一种基于非光滑观测器的针对含有死区和迟 滞的复合三明治系统的状态估计方法，与带死区，带间隙和带迟滞的单一非光滑三明治系 统比较，采用该方法能解决单一非光滑观测器无法解决的复合问题，与传统的观测器比较， 采用该方法的观测器能更准确地估计系统的状态值。
[0010] 实现本发明目的的技术方案是：
[0011] 基于非光滑观测器的含有死区和迟滞的复合三明治系统状态估计方法，包含如下 步骤：
[0012] 步骤1 :利用关键项分离原则和切换函数，由简单到复杂，借鉴已构建的带死区、 带间隙和带迟滞三明治系统非光滑状态空间方程，构建能准确描述含有死区和迟滞的复合 三明治系统的非光滑状态空间方程；
[0013] 步骤2 :根据步骤1构建的复合非光滑三明治状态空间方程，当系统满足观测器的 存在性条件时，构造能随复合非光滑三明治系统工作区间变化而自动切换的非光滑状态估 计观测器，并给出相应非光滑状态估计观测器的存在条件和收敛性定理。
[0014] 所述步骤1包括如下步骤：
[0015] (1)含有死区和迟滞的复合三明治系统的前段线性环节Li(?)，如式（1)所示：
[0016]
(1i
[0017] 含有死区和迟滞的复合三明治系统的后端线性环节L2(.)，如式（2)所示：
[0018]
(2)
[0019] (2)死区和迟滞的建模
[0020] 死区环节DZ的输入输出之间的数学模型，如式（3)所示：
[0021]
(3)[0022] 迟滞的输入输出关系的数学模型，如式（4)所示：
[0023] 4)
[0024]
[0025] 其中，，尽e/Txl，yiGRlxl，C；.ei?1XBSuGRlxl，ViGRlxl， i= 1，2,Xll和x21分别代表前段线性环节和后端线性环节的第i个状态变量；4 是 状态转移矩阵，在e/r4是输入矩阵，yieR1X1是输出变量，ni表示第i个线性子系统的维 数，ueR1X1是输入变量，VleR1X1是迟滞的输入变量同时也是死区的输出变量，v2eR1X1 是迟滞的输出变量。Zl(k)是第i个间隙的输出，Wl是第i个间隙在构成迟滞中所占的权 重，□)是单个间隙的输入输出函数，n用来构建迟滞的间隙的个数;mlb是间隙的斜率， mld是死区的斜率；Dlb是间隙的宽度，Dld是死区的宽度。
[0026] 所述步骤2包括如下步骤：
[0027](1)构造该系统的非光滑观测器，根据式（1)，（2)，（3)和（4)，构造含有死区和迟 滞的复合三明治系统的状态观测器如下：
[0028]
(5)
[0032] 其中，M e /r'1是反馈矩阵；
[0033] (2)给定式（5)观测器的存在条件和收敛性定理：
[0034] 设定系统满足如下条件：
[0035] 条件1状态变量x是有界的，S卩Vife., ||x(k) ||xb，xb> 0;其中||IIm表示m范 数，11)表示状态变量的最大边界值；
[0036] 条件2观测器的初始误差是有界的，即||e(l) ||eb，eb>0;eb表示观测器初 始误差的最大边界值；
[0037] 条件3前段线性子系统的转移矩阵&的特征值都在单位圆内；
[0038] 定理：对于满足以上三个假设条件的含有死区和间隙的复合三明治系统，可以 构造如式（5)所示的观测器对该系统进行状态估计，如果选择观测器的反馈矩阵1( 2，使得 (A2_K2C22)的特征值均在单位圆内，那么式（5)所示的观测器的估计误差最终会收敛到零， 其中，A2表示后端线性系统的转移矩阵，K2表示后端线性系统的反馈矩阵，C22表示后端线 性系统的输出矩阵。
[0039] 本发明方法的优点在于：
[0040] 1.充分考虑了死区和迟滞的非线性和非光滑特性和两种的串联复合特性，通过在 复合三明治系统中，引入切换函数更准确地描述该类系统，通过复合函数公式（4)建立了 其准确的数学模型以描述该类系统，公式（4)把死区和迟滞非光滑特性连接在一起，解决 了如何准确描述复合三明治的难题。
[0041] 2.在构建的系统的准确数学模型的基础上，构造复合非光滑状态观测器估计其状 态值，其中复合非光滑观测器能够随系统工作区间的切换而自动切换，特别是公式（5)中 包含了连接两个非光滑环节的复合表达关系，成功地解决了复合非光滑