第五章-章--线性系统二次型指标的最优控制

第五章线性系统二次型指标的最优控制——线性二次型问题返回主目录5.1引言5.2线性二次型问题的提法5.3终端时间有限时连续系统的状态调节器问题5.4稳态时连续系统的状态调节器问题5.5离散系统的线性二次型问题5.6伺服跟踪问题5.7设计线性二次型最优控制的若干问题5.8小结5.1引言用极小值原理解非线性系统的最优控制将导致非线性两点边值问题，这类问题求解是很困难的。即使系统是线性的，但当指标函数是最短时间、最少燃料这种形式，要求得到最优控制的解析表达式，并构成反馈控制（即把表示为的函数）也是非常困难的。)(tX)(tU返回子目录的确定归结为求解一个非线性矩阵黎卡提（Riccati）微分方程或代数方程。而黎卡提方程的求解已研究得很透彻，有标准的计算机程序可应用，因此，求解既规范又方便。这种问题简称为线性二次型(LinearQuadratic简称LQ)问题，目前应用得十分广泛，是现代控制理论最重要的结果之一。)(tG下面我们将看到，若系统是线性的，指标函数是二次型的（指标函数是和的二次函数），则可以求得线性最优反馈控制律。)(tX)(tU)()()(tXtGtU线性二次型问题的实用意义还在于：例如，在飞行器的轨迹优化问题中，根据飞行器的状态方程（一般是非线性的）用极小值原理计算出名义的最优控制和最优状态轨迹，设分别用和表示。)(0tX)(0tU把它所得到的最优反馈控制与非线性系统的开环最优控制结合起来，可减小开环控制的误差，达到更精确的控制的目的。因为状态方程只能是对飞行器实际动力学特性的近似描绘，这里存在着模型误差，把加到飞行器上去，所产生的实际状态将不同于（这里我们还未考虑作用在飞行器上的其它扰动作用）。（这里我们还未考虑作用在飞行器上的其它扰动作用）。)(0tX)(tX)(0tU令状态误差为，我们要使愈小愈好，为此，可根据构成一个最优反馈控制，作为校正信号加到上去，得到的实际控制信号将使飞行器尽可能沿着飞行。)(0tX)()()(0tUtUtU)(0tU)(tU)(tX)()()(0tXtXtX)(tX由于、应该比较小，它们将满足线性的状态方程，所以可用线性二次型问题设计出反馈控制。我们可用图5-1表示上面的思想。)()()(tXtGtU)(tU)(tX)(0tX)(tX)(tX)(tU)(tU)(0tU))00tXtU（名义状态（名义控制计算机存储（飞行器）实际系统线性二次最优反馈控制图5-1线性二次最优反馈控制的应用5.2线性二次型问题的提法一般情况的线性二次型问题可表示如下：其中，为维状态向量，为维控制向量，为维输出向量。设不受约束。)(tU)(tY)(tX)(tUmln设线性时变系统的方程为)()()()()(tUtBtXtAtX（5-1）)()()(tXtCtY（5-2）返回子目录其中，为维理想输出向量。寻找最优控制，使下面的性能指标最小l)(tZfttTTffTdttUtRtUtetQtetPeteuJ0)()()()()()(21)()(21)(（5-4）令误差向量为)(te)()()(tYtZte（5-3）其中，是对称半正定常数阵，是对称半正定阵,是对称正定阵。一般将、、取成对角阵。)(tR)(tQPll)(tR)(tQllPmm下面对性能指标中的每一项作一说明。因为正定阵，则当，就有。例如0)()()(tUtRtUT0)(tU)(tRJ)(00)()(21trtrtR)()()(21tututUdttutrtutrdttUtRtUffttttT00)]()()()([21)()()(21222211设，，则为正定阵，于是)(tR0)(2tr0)(1tr它与消耗的控制能量成正比，消耗得越多，则性能指标值越大。故性能指标中这一项表示了对消耗控制能量的惩罚。、可看作加权系数，如认为的重要性大于，则可加大。将选成时间函数，是为了对不同时刻的加权不一样。实际上，为了简单起见常选用常数阵。R)(tU)(tR)(1tr)(2tu)(1tr)(2tr)(1tuJ为半正定阵，则当，就有，表示误差平方和积分，故这项表示对系统误差的惩罚。表示对终端误差的惩罚，当对终端误差要求较严时，可将这项加到性能指标中。)()(21ffTtPetedttetQtefttT)()()(2100)()()(tetQteT0)(te)(tQ总之，性能指标最小表示了要用不大的控制量来保持较小的误差，以达到能量和误差的综合最优。)(uJ这时（单位阵），理想输出，则，这时，问题归结为用不大的控制量使保持在零值附近。因而称为状态调节器问题。)()()(tetXtY0)(tZItC)()(tX下面讨论几种特殊情况：1）调节器问题。例如电机转速调节系统中，由于外加电压波动使转速偏离要求值，通过施加控制使转速偏差趋于零。这时，，这时要用不大的控制量使跟踪，因而称为跟踪问题。例如，用雷达跟踪飞行器的运动，通过控制使跟踪误差趋于零。)(tZ)(tY)()()(tYtZte0)(tZ2）伺服机问题。5.3终端时间有限时连续系统的状态调节器问题要求寻找最优控制，使最小。这里无约束。、为对称半正定阵，为对称正定阵。终端时间为有限值。ft)(tR)(tQP)(tUJ)(tU)()()()()(tUtBtXtAtX00)(XtX(5-5))()(21ffTtPXtXJfttTTdttUtRtUtXtQtX0)()()()()()(21(5-6)J考虑下面的系统状态方程和性能指标返回子目录5.3.1用极小值原理求解上面的问题因无约束，故等同于用经典变分法求解。取哈密顿函数为)(tU协态方程为最优解的必要条件如下：)()()()()(tUtBtXtAtT)()()()()()(21tUtRtUtXtQtXHTT（5-7）)()()()(ttAtXtQXHT（5-8）)()()(1ttBtRUT因正定，故存在，由上式可确定最优控制。为寻求最优反馈控制律还需把与状态联系起来。)(tX)(tU)(1tR)(tR)(tU0)()()()(ttBtUtRUHT（5-9）控制方程为我们再一次遇到了两点边值问题（已知和），如前所述，一般要试凑再积分协态方程使满足要求。)(0tX)(ft)(0t)(ft但这里处理的是线性微分方程，可找到更简单的解法。从(5-10)可见，协态和状态在终端时刻成线性关系。ft)(tX)(t)()()(21)()()(fffTffftPXtPXtXtXtXt(5-10)横截条件为然后再来求出（这种方法称为扫描法）。将（5-11）代入（5-9），再代入（5-5），得)(tK)()()(tXtKt（5-11）)()()()()()()()(1tXtKtBtRtBtXtAtXT(5-12))()()()()()()()()()(tXtKtAtXtQtXtKtXtKtT(5-13)由（5-11）和（5-8）可得这启发我们假定：上式对任意都应成立，故方括号内的项应为零，这就得出)(tX0)()]()()(tXtQtKtAT)()()()()()()()([1tKtBtRtBtKtAtKtKT)()()()()()(1tQtKtBtRtBtKT()()()()()TKtKtAtAtKt(5-14)将(5-12)代入(5-13)可得上式是的非线性矩阵微分方程，称为黎卡提（Riccati）矩阵微分方程。一般来说得不出的解析表达式，但可用计算机程序算出的数值解。为了求解，要知道它的边界条件。比较(5-11)和(5-10)可知)(tK)(tK)(tK)(tK因此可从到逆时间积分黎卡提微分方程，求出。由(5-9)和(5-11)就可构成最优反馈控制ft0t)(tKPtKf)((5-15)又称为最优反馈增益矩阵。最优反馈系统的结构图如图5-2所示。)()()()(1tKtBtRtGT)()()()()()()(1tXtGtXtKtBtRtUT(5-16))(tU)(tX)(tK)()(1tBtRT)()()()()(tUtBtXtAtX图5-2最优反馈系统的结构图注意到与状态无关，故可在系统未运行前，将先计算出来（称为离线计算），把它存储在计算机中。)(tK)(tX)()()(1tKtBtRT在系统运行时，将从计算机存储元件中取出，与同一时刻测量到的相乘，就可构成最优控制。)()()(1tKtBtRT)(tX)(tU由此可见，系统运行时的计算量（称为在线计算量）只是一个乘法计算，故可用简单的微计算机来完成。5.3.2矩阵黎卡提微分方程的求解及的性质1、ttKttKdttdKt)()(lim)(0于是可用下面的差分方程来近似黎卡提微分方程)()()()({)()(tKtAtAtKttKttKT)}()()()()()(1tQtKtBtRtBtKT（5-17）矩阵黎卡提微分方程是非线性的，一般不能求得闭合形式的解。在数字机上求解时，可用一阶差分代替微分2、求解上式时，以为初始条件，取为负的小量，从到逆时间递推计算，即可出。PtKf)(tft0t)(tK是对称矩阵，即，表示转置。这可证明如下：因为、、都是对称的，将(5-14)式转置一下，可得)(tK)()(tKtKTTP)(tQ)(tR)()()()()()(tKtAtAtKtKtKTTTTT)()()()()()(1tQtKtBtRtBtKTTT因此和一样满足同一黎卡提方程，并且边界条件一样，即。于是，由微分方程解的唯一性可知)(tKT)(tK)()(ffTtKPtK)()(tKtKT利用这个对称性，求维的元时，只需积分个方程即可。)(tK2)1(nnnn3、即使系统是定常的，即系统矩阵A，输入矩阵B为常数阵，加权阵和也是常数阵，但仍为时变阵。QR)(tK这从是黎卡提微分方程的解可看出。时变时，反馈控制增益也时变，在实现时总是不太方便。)(tK)(tK下一段将看到，对线性定常系统，若终端时间，且系统满足一些附加条件时，将变为常数阵。ftK)(tK例5-1设系统状态方程为0)0(22xux1)0(121xxx（5-18）ftdttutxJ0221)()(21（5-19）寻找最优控制使下面的性能指标为最小。()ut解考虑到是对称阵，设)(tK为简单起见，上式右端省略了自变量。把上面的、、、和代入黎卡提方程（5-14）式，可得ABQR)(tKt1,0001,0,10,0010RQPBA（5-20）22121211)(kkkktK（5-21）把状态方程（5-18）和（5-5）式相比较，把性能指标（5-19）和（5-6）式相比较，可得22121211221212112212121101000010kkkkkkkkkkkk212111kk22121112kkkk令上式等号左右端的对应元相等，得000110102212121122121211kkkkkkkk22212222kkk（5-23）222121222112212112122,,1kkkkkkkkk（5-22）由到逆时间积分上面的非线性微分方程组，即可求得。于是最优控制为ft0t)(),(),(221211tktktk)()(101)()()(21221212111txtxkkkktXtKBRtuT)()()()(222112txt