现代控制理论第六章课件

1、第6章 状态观测与状态最优估计 某些状态量，或者由于不具明确的物理意义，或者由于量测手段的限制，在工程实际中不能直接获取它们。状态观测器可实现对状态的重构。而对于存在随机噪声的系统，则必须利用统计方法对状态量进行最优估计。 1 状态重构与状态观测器 一、状态重构问题 输入量u和输出量y总是可以直接量测的，能否通过输入量u和输出量y间接获取状态量的信息。 为此，对输出方程进行逐次微分运算，并代之以状态方程，可得： 写成矩阵方程形式： 矩阵 满秩 ，x有唯一解。但实际应用中不可取。 启示：如果系统满足一定条件，利用系统的输入量和输出量，得到原系统状态量的间接值 ，它在一定的指标下与x(t)等价。

2、称 为状态量x(t)的重构值，将得到重构状态 的系统称为状态观测器，表示为 。等价性指标一般采用渐近等价，即 如果状态观测器的维数与原系统的维数相同，称为全维状态观测器；如果状态观测器的维数小于原系统的维数，称为降维状态观测器。 二、全维状态观测器1观测器的构成用原系统的结构、输入构造一个模拟系统：有：开环型状态观测器 （1）A包含有不稳定的特征值时，即使很小的 也会使 远离x(t)； （2）观测器参数对原系统参数的任何偏离都会产生不利影响。 所以开环型状态观测器不能实际使用。解决的办法是利用输出偏差 进行反馈，反馈矩阵为M。如图： 观测器的状态方程式为： 有望通过设计合适的偏差反馈矩阵M以调

3、整观测器系统矩阵的特征值（观测器极点），实现渐近等价指标下的状态重构。 所以，一个性能优良的观测器应该是所有极点可以任意配置的。这就是观测器的极点配置问题。2极点任意配置条件结论：系统能采用全维状态观测器重构其状态，并且能通过改变M矩阵任意配置观测器极点的充要条件是原系统完全能观。 其转置 特征值不变，即通过 K 矩阵可任意配置特征值； 取 ，即矩阵(AMC)的特征值可通过M矩阵任意配置；（1）判断 的能观性； 显然原系统能观，它对应的全维状态观测器就能通过改变M矩阵任意配置它的极点。3极点配置算法（1）判定 的能观性； （2）如能观，写出原系统的对偶系统 ; （3）利用状态反馈极点配置算法求

4、出期望极点为 的状态反馈系统 的反馈矩阵 ； （4）取 ; （5）得状态观测器为： 对于单输出系统，除了通过对偶系统求解外，也有类似于单输入系统状态反馈极点配置的二种算法。 方法一（解联立方程）：（5）将m代入方程 ，得出全维状态观测器。 解 （1）系统的能观性矩阵为满秩，系统能观；（4）由 （5）得全维观测器为 ：系统的状态变量图为：方法二（利用能观规范型求）： （1）先判断 的能观性，若能观，则往下进行 ； （2）开环系统的特征多项式：（6）由 求得偏差反馈向量m，并代入观测器方程 。（3）由给定的期望极点求得期望的特征多项式 ： （4）按下式求取具有能观规范型形式的状态空间中的偏差反馈向

5、量： （5）求取将原系统化为能观规范型的变换矩阵P； 对于期望极点的位置，仅从渐近收敛速度看，希望极点尽量远离虚轴。但是极点离虚轴太远，会使观测器频带过宽，不利于扼制观测器输入量的高频干扰。要根据工程实际折衷考虑。 一般，系统中总有一部分状态变量是可以直接量测的。从而，只需构造维数小于n的观测器来得出另一部分状态变量（降维状态观测器） 。如果 ，则有q个输出变量是相互独立的，那么由输出方程就能得出q个状态变量。例如极端情况 ，则后q个输出量就是状态变量，可量测；一般情况下，降维状态观测器的最小维数为 。 三、降维状态观测器 1. 降维状态观测器的构成引入非奇异变换，使新状态空间的状态量为： 使

6、新状态空间的输出矩阵为： q维分状态向量 直接由y得出，而 维分状态向量 需要通过观测器重构。由上面式子可写出： 为了重构(n-q)维状态向量 ，只要构造上述子系统的全维状态观测器即可。 由于原系统能观，非奇异变换后仍然能观，它的部分状态变量构成的子系统当然也能观。所以能对上述子系统构造全维状态观测器。有： 上式含有输出的导数项，这对于观测器抗干扰及观测值的唯一性考虑都是不允许的，为此引入一个新的状态量：于是，降维状态观测器的方程可写为： 或者写为： 而状态量的重构值为： 如将非奇异变换矩阵表示为： 则在原状态空间中状态量的重构值为 ： 此即为(n-q)维降维状态观测器，也称 Luenberg

7、er观测器。降维状态观测器结构图为： 2. 降维状态观测器的设计算法 判别 (A,C) 的能观性，并确定q和 n-q：n(n-q)nq 对原系统实施非奇异变换： 写出降维状态观测器方程：并按观测器极点配置算法求出M； 写出状态量 的估计值 ： 经反变换求出原系统状态 x 的估计值 ：即：降维观测器方程为： 可以画出降维状态观测器如下：2 引入观测器的状态反馈控制系统 一、系统的构成 控制系统由三部分组成：被控对象、状态观测器、状态反馈控制。结构图如下：将三部分合在一起，即得含观测器的状态反馈控制系统：二、系统的特性：1、系统的维数=原系统的维数+观测器的维数。系统的特征值集合=状态反馈系统的特

8、征值集合+观测器的特征值集合。系统矩阵为：引入非奇异变换：有：2、由上式还可以得出结论：通过K 配置系统特征值（闭环极点）和通过M配置观测器特征值（极点）是互相分离的，可以完全独立地进行。这就是分离性原理。 可见，系统的特征值由状态反馈系统的特征值和状态观测器的特征值二部分组成。3、观测器的引入不改变原状态反馈系统的传递函数矩阵。上面的讨论给出了 ，同样可得新状态空间的输入矩阵和输出矩阵： 非奇异变换不改变传递函数矩阵，所以有： 进一步分析可知， 具有按能控性分解的形式，能控子系统为 ，观测器部分是不能控的。所以，观测器的引入使状态反馈控制系统不再保持能控性。 分块矩阵的求逆公式： 必发生了零

9、极点相消现象，相消的n个极点是属于观测器的。由于观测器设计保证了其极点的渐近稳定性，所以零极点相消不影响闭环系统的正常运行。 4、观测器为渐近等价，观测器动态特性将影响闭环系统动态特性，要求观测器的动态过程快于闭环系统的动态过程是合理的。通常把观测器特征值的负实部取为状态反馈系统特征值的负实部的23倍。 解：（1）独立设计降维状态观测器；2）构造非奇异变换矩阵Q,使变换后的对于该系统可以通过重新安排状态变量实现，即输出方程：状态方程为：3）降维状态观测器的方程为观测器的特征多项式为：由 求得：4）降维状态观测器的方程为期望的特征多项式应为： 5）状态量 的重构值为： 6）再顺序安排状态变量，得

10、状态量x的重构值 ： （2）独立设计状态反馈控制； 1）首先判别系统的能控性：系统能控；2）由给定的期望闭环极点求得期望的闭环系统特征多项式为：3）由闭环系统动态方程写出的闭环系统特征多项式为：4）由 求得 ： （3）引入状态观测器的状态反馈控制为：被控对象：观测器：控制作用：可画出引入观测器的状态反馈控制系统的状态变量图 ： 3 状态最优估计 一、状态估计问题的描述w(t)为m维随机干扰(噪声)向量，称为系统（输入）噪声； 状态方程：测量方程： v(t)为q维随机干扰（噪声）向量，称为测量噪声； 所谓状态估计就是根据测量值y(t)及随机干扰的统计特性，对系统的状态量进行估计，得出尽量接近状态

11、量真实值x(t)的估计值 。 希望在一定的准则函数下所作出的估计是最好的，即最优估计问题。最优估计的解通过准则函数极小化（或极大化）得出。 用准则函数（或指标函数）来衡量估计的好坏。 不同的准则函数对应得出不同的估计方法。 二、最小二乘估计以误差平方和最小作为准则函数。 对系统进行k次测量，记第i次测量为 ： k次测量后，可得： 以加权误差平方和（按测量值的精度分配权值）最小作为准则函数： 1、静态最小二乘估计 假设了测量过程中x不变 2、动态最小二乘估计实际的控制系统状态量是变化的，变化规律由系统的状态方程决定： 线性定常离散系统，下标表示时间，仅考虑测量噪声 输入为确定性输入时，可设 考虑

12、 ，有： 得：时为：当同时考虑系统噪声可得到类似的结果。 由测量序列 求得估计值 的基础上，通过新测量值 对 的修正得出新估计值 ，解决存储量和计算量不断增大的问题。 3、最小二乘估计的递推算法 以 的动态估计为例：令：考虑用 替代 ，上式为：再令：则有：还可推导出： 递推计算： 如果估计值 是测量值 的线性函数，则只需事先知道系统噪声和测量噪声的一、二阶矩，即线性最小方差估计。 三、线性最小方差估计估计值的方差最小作为准则函数。 一般需要已知系统噪声、测量噪声的概率密度以及它们的联合概率密度，较难满足。 考虑估计值是测量值的线性函数： 估计值的方差： 令： 则：要求 最小，必须有：由此得：

13、在已知状态量和测量量（或者系统噪声和测量噪声）的一、二阶矩时，就能得到状态量的线性最小方差估计值。 四、卡尔曼滤波基于线性最小方差估计的递推算法。为实际应用提供了可能性。 确定性输入，设 为一步转移矩阵 系统噪声和测量噪声为零均值的白噪声，它们相互独立，并与状态量也不相关。 1一步预测与新息：一步预测 一步预测误差 （新息） 加权修正系数 一步预测 一步预测误差 增益矩阵 估计值： 与 、 不相关， ，则：2估计误差方差阵的求取：估计误差： 估计误差的方差阵： 一步预测误差的方差阵 ： 又 与 、 不相关， ，可得：得： 3增益矩阵 的求取： 估计误差方差最小等价于误差方差阵的迹最小，即：而：和 为对称矩阵，所以有：、令其为零，解得： 可得：总结上面可知，卡尔曼滤波由4个基本方程组成，它们是： （1）一步预测误差方差阵： （2）最优滤波增益阵： （3）状态最优估计方程： （4）估计误差方差阵： 其中一步预测方程 ：卡尔曼滤波示意图： 一步预测最优估计基本方程应用时的几点说明（稳定性、初值、发散）； 带控制作用项系统的卡尔曼滤波； 系统噪声与测量噪声相关的卡尔曼滤波