李亚普诺夫稳定性分析和二次型最佳控制.docx
《李亚普诺夫稳定性分析和二次型最佳控制.docx》由会员分享，可在线阅读，更多相关《李亚普诺夫稳定性分析和二次型最佳控制.docx（16页珍藏版）》请在冰豆网上搜索。
李亚普诺夫稳定性分析和二次型最佳控制
5.6.3二次型最优控制问题
现在我们来研究最优控制问题。
已知系统方程为
（5.20）
确定最优控制向量
（5.21）
的矩阵K，使得性能指标
（5.22）
达到极小。
式中Q是正定（或正半定）Hermite或实对称矩阵，R是正定Hermite或实或实对称矩阵。
注意，式（5.22）右边的第二项是考虑到控制信号的能量损耗而引进的。
矩阵Q和R确定了误差和能量损耗的相对重要性。
在此，假设控制向量是不受约束的。
正如下面讲到的，由式（5.21）给出的线性控制律是最优控制律。
所以，若能确定矩阵K中的未知元素，使得性能指标达极小，则对任意初始状态x（0）而言均是最优的。
图5.6所示为该最优控制系统的结构方块图。
图5.6最优控制系统
现求解最优控制问题。
将式（5.21）代入式（5.20），可得
在以下推导过程中，假设是稳定矩阵，的所有特征值均具有负实部。
将式（5.21）代入（5.22），可得
依照解参数最优化问题时的讨论，取
式中的P是正定的Hermite或实对称矩阵。
于是
比较上式两端，并注意到方程对任意x均应成立，这就要求
（5.23）
根据Lyapunov第二法可知，如果是稳定矩阵，则必存在一个满足式（5.23）的正定矩阵P。
因此，该方法由式（5.23）确定P的各元素，并检验其是否为正定的（注意，这里可能不止一个矩阵P满足该方程。
如果系统是稳定的，则总存在一个正定的矩阵P满足该方程。
这就意味着，如果我们解此方程并能找到一个正定矩阵P，该系统就是稳定的。
满足该方程的其他矩阵P不是正定的，必须丢弃）。
性能指标可计算为
由于假设A-BK的所有特征值均具有负实部，所以。
因此
（5.24）
于是，性能指标J可根据初始条件x（0）和P求得。
为求二次型最优控制问题的解，可按下列步骤操作：
由于所设的A是正定Hermite或实对称矩阵，可将其写为
式中T是非奇异矩阵。
于是，式（5.23）可写为
上式也可写为
求J对K的极小值，即求下式对K的极小值
（见例5.21）。
由于上面的表达式不为负值，所以只有当其为零，即当
时，才存在极小值。
因此
（5.25）
式（5.25）给出了最优矩阵K。
所以，当二次型最优控制问题的性能指标由式（5.22）定义时，其最优控制律是线性的，并由
给出。
式（5.25）中的矩阵P必须满足式（5.23），即满足下列退化方程
（5.26）
式（5.26）称为退化矩阵黎卡提方程，其设计步骤如下：
1、求解退化矩阵黎卡提式（5.26），以求出矩阵P。
如果存在正定矩阵P（某些系统可能没有正定矩阵P），那么系统是稳定的，即矩阵是稳定矩阵。
2、将矩阵P代入式（5.25），求得的矩阵K就是最优矩阵。
例5.9是建立在这种方法基础上的设计例子。
注意。
如果矩阵是稳定的，则此方法总能给出正确的结果。
确定最优反馈增益矩阵K还有另一种方法，其设计步骤如下：
1、由作为K的函数的式（5.23）中确定矩阵P。
2、将矩阵P代入式（5.24）,于是性能指标成为K的一个函数。
3、确定K的各元素，使得性能指标为极小。
这可通过令等于零，并解出的最优值来实现J对K各元素为极小。
这种设计方法的详细说明见例5.11和5.12。
当元素的数目较多时，该方法很不便。
如果性能指标由输出向量的形式给出，而不是由状态向量的形式给出，即
则可用输出方程
来修正性能指标，使得J为
（5.29）
且仍可用本节介绍的设计步骤来求最优矩阵K。
------------------------------------------------------------------
[例5.9]研究如图5.7所示的系统。
假设控制信号为
试确定最优反馈增益矩阵K，使得下列性能指标达到极小
式中
由图5.7可看出，被控对象的状态方程为
式中
图5.7控制系统
以下说明退化矩阵黎卡提代数方程如何应用于最优控制系统的设计。
求解（5.26），将其重写为
注意到A为实矩阵，Q为实对称矩阵，P为实对称矩阵。
因此，上式可写为
该方程可简化为
由上式可得到下面3个方程
将这3个方程联立，解出、、，且要求P为正定的，可得
参照式（5.25）,最优反馈增益矩阵K为
因此，最优控制信号为
（5.28）
注意，由式（5.28）给出的控制律对任意初始状态在给定的性能指标下都能得出最优结果。
图5.8是该系统的方块图。
图5.8图5.7所示对象的最优控制
------------------------------------------------------------------
5.7二次型最优控制问题的MATLAB解法
在MATLAB中，命令
可解连续时间的线性二次型调节器问题，并可解与其有关的黎卡提方程。
该命令可计算最优反馈增益矩阵K，并且产生使性能指标。
在约束方程
条件下达到极小的反馈控制律
另一个命令
也可计算相关的矩阵黎卡提方程
的唯一正定解P。
如果为稳定矩阵，则总存在这样的正定矩阵。
利用这个命令能求闭环极点或的特征值。
对于某些系统，无论选择什么样的K，都不能使为稳定矩阵。
在此情况下。
这个矩阵黎卡提方程不存在正定矩阵。
对此情况，命令
不能求解，详见MATLABPrgram5.1。
MATLABProgram5.1
%——Designofquadraticoptimalregulatorsystem——
%*****DeterminationoffeedbackgainmatrixKforquadratic
%optimalcontrol*****
%*****EnterstatematrixAandcontrolmatrixB*****
A=[-11;02]
B=[1;0];
%*****EntermatricesQandRofthequadraticperformance
%index*****
Q=[10;01];
R=[1];
%*****Toobtainoptimalfeedbackgainmatrix,K,enterthe
%followingcommand*****
K=lqr（A,B,Q,R）
Warning:
Matrixissingulartoworkingprecision.
K=
NaNNaN
%*****lfweenterthecommand[K,P,E]=lqr（A,B,Q,R）.then*****
[K,P,E]=lqr（A,B,Q,R）
Warning;Matrixissingulartoworkingprecision.
K=
NaNNaN
P=
-lnf-lnf
-lnf-lnf
E=
-2.0000
-1.4142
------------------------------------------------------------------
[例5.10]考虑由下式确定的系统
证明：
无论选择什么样矩阵K，该系统都不可能通过状态反馈控制
来稳定（注意，该系统是状态不可控的）。
定义
则
因此特征方程为
闭环极点为
由于极点在s的右半平面，所以无论选择什么样的矩阵K，该系统都是不稳定的。
因此，二次型最优控制方法不能用于该系统。
假设在二次型性能指标中的Q和R为
并且写出MATLABProgam5.1。
所得的MATLAB解为
其中NaN表示“不是一个数”。
每当二次型最优控制问题问题的解不存在时，MATLAB将显示矩阵K由NaN组成。
------------------------------------------------------------------
[例5.11]考虑下式定义的系统
式中
性能指标J为
这里
假设采用下列控制u
确定最优反馈增益矩阵K。
最优反馈增益矩阵K可通过求解下列关于正定矩阵P的黎卡提方程得到
其结果为
将该矩阵P代人下列方程，即可求得最可求得最优矩阵K为
因此，最优控制信号为
利用MATLABProgram5.2也能求解该问题。
MATLABProgram5.2
%——Designofquadraticoptimalregulatorsystem
%*****DeterminationoffeedbackgainmatrixKforquadratic
%optimalcontrol*****
%*****EnterstatematrixAandcontrolmatrixB*****
A=[01;0-1];
B=[0;1];
%*****EntermatricesQandRofthequadraticperformance
%index*****
Q=[10;01];
R=[1];
%TheoptimnalfeedbckgainmatrixK（ifsuchmatrixK
%exists）canbeobtainedbyenteringthefollowingcommand*****
K=lqr（A,B,Q,R）
K=
1.00001.0000
------------------------------------------------------------------
[例5.12]考虑下列系统
式中
性能指标J为
式中
求黎卡提方程的正定矩阵R、最优反馈增益矩阵K和矩阵A-BK的特征值。
利用MATLABProgram5.3，可求解该问题。
MATLABProgram5.3
%--------Designofquadraticoptimalregulatorsystem-------
%*****DeterminationoffeedbackgainmatrixKforquadratic
*****EnterstatematrixAandcontrolmatrixB*****
A=[00;001;-35-27-9];
B-[0;0;1]
%******EntermatricesQandRofthequadraticperformance
%index*****
Q=[100;010;001];
R=[1];
%******TheoptimalfeedbackgainmatrixK,solutionPofRiccati
%equation,andclosed-looppoles（thatis,theeigenvalues
%ofA-BK）canbeobtainedbyenteringthefollowing
%command*****
[K,P,E]=[qr（A,B,Q,R）]
K=
0.01430.11070.0676
P=
4.26252.49570.0143
2.49572.81500.1107
0.01430.11070.0676
E=
-5.0958
-0.9859
升级会员 