线性连续系统最优控制补充讲义(2004-03-09)3.6可化为规范形式的LQ问题3.6.1具有规定衰减速度（稳定度）的调节器系统稳定条件，极点在左半复平面。+ Im衰减速度：系统离虚轴最近的闭环极点与虚轴间的距离，。越大，系统的非零初态响应的衰减速度愈快。+Re若闭环系统的极点都在距离虚轴为的直线左边，则称闭环系统有至少不低于。的衰减速度。最优化问题I：er00(0+ 0r0)lamin J =-x()= Ax(0)+ Bu(t)x(0)= xos.t式中，R正定，O半正定，α为正常量，（A,B）为能控对，（A,D）为能观对，Q= DDTu(l)=eu(0)令,x(0)=ex(0)则有=e"x+aex= e" (Ax+ Bu)+ ae" x=(A+aαl)e"x+Be"u转化为最优化问题IⅡI：
线性连续系统最优控制补充讲义 （2004-03-09） 3.6 可化为规范形式的 LQ 问题 3.6.1 具有规定衰减速度（稳定度）的调节器 系统稳定条件，极点在左半复平面。 衰减速度：系统离虚轴最近的闭环极点与 虚轴间的距离，。 越大，系统的非零初态 响应的衰减速度愈快。 若闭环系统的极点都在距离虚轴为的直 线左边，则称闭环系统有至少不低于 的衰减 速度。 最优化问题 I： ( ) ( ) ( ) ( ) = + 0 2 2 1 min J e x t Qx t u t Ru t dt t T T s.t x (t) = Ax(t)+ Bu(t) ( ) 0 0 x = x 式中，R 正定，Q 半正定，为正常量，（A,B）为能控对，（A,D）为能观对， T Q = DD 令， x(t) e x(t) t = u(t) e u(t) t = 则有 ( ) (A I)e x Be u e Ax Bu e x x e x e x t t t t t t = + + = + + = + 转化为最优化问题 II： Re Im
[(00x(0)+7(0Ra(0)min J=!x=(A+αl)x+ Bus.tx(0)= x(0)= xo可以证明：[A,B]完全能控→[A+αl,B]完全能控[A,D]完全能观→[A+αl,D]完全能观。则，最优化问题ⅡI有唯一解：u()=-Kax(0)其中Ka=R-B'PaPa(A+αl)+(A +αl)Pa-P.BR-'B'Pa+Q=0又由于u=e-"u=-Ke-"x()=-K,x(l)所以，两个最优控制问题的反馈增益阵是相同的，又由于最优控制问题II的闭环系统文=(A+αl-BK。)x渐近稳定，即lim x(t)= lim e"x(t) = 0所以，x()也是渐近稳定的，且其衰减速度不低于e-"。即闭环系统x=(A-BK。)x的极点均落在复平面-α直线的左边。注意，这时并不能保证衰减速度恰好是α，且大于α多少也不能保证。可以由小到大改变α，通过试验的方法找到合适的α值，保证有规定的衰减率3.6.2具有非零设定点的调节器线性定常系统
( ) ( ) ( ) ( ) = + 2 0 1 min J x t Qx t u t Ru t dt T T u s.t ( ) ( ) ( ) 0 0 0 x x x x A I x Bu = = = + + 可以证明： [A,B]完全能控 A+I,B 完全能控。 [A,D]完全能观 A+I,D 完全能观。 则，最优化问题 II 有唯一解： u(t) K x(t) = − 其中 K R B P −1 T = ( ) ( ) 0 1 + + + − + = − P A I A I P P BR B P Q T T 又由于 u e u K e x(t) K x(t) t t = = − = − − − 所以，两个最优控制问题的反馈增益阵是相同的，又由于最优控制问题 II 的闭环系统 x (A I BK )x = + − 渐近稳定，即 lim ( ) = lim ( ) = 0 → → x t e x t t t t 所以， x(t) 也是渐近稳定的，且其衰减速度不低于 t e − 。即闭环系统 x (A BK )x = − 的极点均落在复平面-直线的左边。 注意，这时并不能保证衰减速度恰好是，且大于多少也不能保证。可以 由小到大改变，通过试验的方法找到合适的值，保证有规定的衰减率。 3.6.2 具有非零设定点的调节器 线性定常系统
x= Ax+ Bu=Cx过去讨论的都是稳态（平衡态）为零状态的问题。设给定点为ysp，稳态时各量应满足Jsp=CxdO = Axa + Bud其中，ud,xd为希望的稳态时控制和希望的状态。取指标函数min J=[6(0)-y,)o(()-y)+(u(0)-ua) R(u(0)-ua)t令[x(0)= x(0)- xa(0)= y(0)- yp[a(0)= u()- ud则[X = Ax + Bu=CX最优控制问题转化为[sr(0)o()+ar(0)Ra(0)tmin2[X=AX+Bix(0)=xo - xds.t=CX这是一个规范化输出调节器问题。解得：i()=-Kx()K= R-'B'PPA+ ATP-PBR-'BIP+C'QC = 0则
= = + y Cx x Ax Bu 过去讨论的都是稳态（平衡态）为零状态的问题。 设给定点为 ysp，稳态时各量应满足 d d sp d Ax Bu y Cx = + = 0 其中，ud, xd 为希望的稳态时控制和希望的状态。 取指标函数 J (y(t) y ) Q(y(t) y ) (u(t) u ) R(u(t) u )dt d T s p d T s p u = − − + − − 2 0 1 min 令 ( ) ( ) ( ) ( ) ( ) ( ) = − = − = − d sp d u t u t u y t y t y x t x t x ~ ~ ~ 则 = = + y Cx x Ax Bu ~ ~ ~ ~ ~ 最优控制问题转化为 J y (t)Qy(t) u (t)Ru(t)dt T T u + 0 ~ ~ ~ ~ ~ 2 1 min s.t = = + y Cx x Ax Bu ~ ~ ~ ~ ~ ( ) d x = x − x 0 0 ~ 这是一个规范化输出调节器问题。 解得： u(t) Kx(t) ~ ~ = − K R B P −1 T = 0 1 + − + = − PA A P PBR B P C QC T T T 则
u(t)= -Kx(0)+ ua + Kx =-Kx(0)+ ul其中第一项由最优化问题解得：[min J =(Oy+u' Ruat2J0(不考虑非零设定值)x=Ax+Bus.t.ly= cx确定ua:X=(A-BK)x+ Bua代入状态方程lim x(0)= 0由于系统渐进稳定O = (A- BK)xa + Bua由于(A-BK)的特征值都在左半复平面，(A-BK)非奇异。则xa = -(A- BK)" Bulayp=-C(A-BK)'Bua（这是一个线性方程组，需讨论其有解、无解、有无穷解的条件。）当r=m，控制维数=输出维数时，有如下定理。定理：已知线性定常系统x=Ax+Buly=Cx其中y和u具有相同的维数，对于任何渐近稳定的定常控制律：u(0)=-Kx(0)+ ud令W。(s)是开环传递函数阵W(s)=C(sI-A)"BW,(s)是闭环传递函数阵W.(s)=C(sl-A+BK)"B当且仅当，W。(s)的“分子多项式”没有s=0的零点，则W.(s)是非奇异矩
( ) ( ) ( ) d d ud u t = −Kx t + u + Kx = −Kx t + 其中第一项由最优化问题解得： ( ) = = + = + y cx x Ax Bu st J y Qy u Ru dt T T . . 2 1 min 0 （不考虑非零设定值） 确定 d u ： 代入状态方程 ( ) Bud x = A− BK x + 由于系统渐进稳定 lim ( ) = 0 → x t t ( ) d Bud 0 = A− BK x + 由于(A-BK)的特征值都在左半复平面，(A-BK)非奇异。 则 ( ) d A BK Bud x = − − −1 ( ) sp C A BK Bud y = − − −1 (这是一个线性方程组，需讨论其有解、无解、有无穷解的条件。) 当 r = m ，控制维数=输出维数时，有如下定理。 定理：已知线性定常系统 其中 y 和 u 具有相同的维数，对于任何渐近稳定的定常控制律： ( ) ( ) ud u t = −Kx t + 令 W (s) 0 是开环传递函数阵 W (s) C(sI A) B 1 0 − = − W (s) c 是闭环传递函数阵 Wc (s) C(sI A BK) B −1 = − + 当且仅当， W (s) 0 的“分子多项式”没有 s = 0 的零点，则 W (s) c 是非奇异矩 = = + y Cx x Ax Bu
阵，若选u=W(0)ys可保证稳态时y(o)=ysp，从而实现非零给定点调节。yspw(0)cx=Ax+BuH例：（解书p389）受控系统：x= Ax+Buy=Cx其中0001000100A=B:0010000-343 4.6332540]00C=[1取01000000Q=C'C:R=r=10500000000解最优控制LQR得K =[30.245.7360.31600.01473]得u=-Kx对零设定点对非零设定点：u=-Kx+W-(0). ysp而
阵，若选 ( ) d c sp u =W 0 y 可保证稳态时 ( ) sp y = y ，从而实现非零给定点调节。 例：（解书 p389）受控系统： y Cx x Ax Bu = = + 其中 − − = 0 0 343 4.63 0 0 0 1 0 0 1 0 0 1 0 0 A = 3254 0 0 0 B C = 1 0 0 0 取 = = 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 Q C C T 5 R = r = 10 解最优控制 LQR 得 K = 30.24 5.736 0.3160 0.01473 得 u = −Kx——对零设定点 对非零设定点： ( ) c sp u = −Kx +W y − 0 1 而 (0) −1 Wc a x = Ax + Bu C + - K ysp u x y