第9章 最优控制1

《第9章 最优控制1》由会员分享，可在线阅读，更多相关《第9章 最优控制1（15页珍藏版）》请在装配图网上搜索。

1、9最优控制用经典控制理论设计系统时，是根据给定的频域和时域指标，通过选择适当 结构和参数的校正装置进行调节，本质上是一个试凑的方法，设计质量的好坏很 大程度上依赖于设计人员的实践经验。而最优控制的设计目标是要选择适当的控 制的规律，使控制系统在严格的数学基础上实现系统的性能品质在某种意义下是 最优的。最优控制是现代控制理论的一个重要组成部分。本章主要讨论最优控制问题的基本概念。二次型最优控制问题以及最优观测 器设计。9. 1最优控制问题实例9. 1. 1最优控制问题实例设有一盛放液体的连续搅拌槽，如图9.3 所示。槽内装有搅拌器使槽内液体经常处于完 全混合状态。若槽内液体初始温度为0C,现 需

2、将其温度经1小时后升高到40Co为此在 入口处送进一定量的温度为u(t)的加热液体。(为保持槽内液体恒定，出口处的流出量等于 入口处的流入量)，试寻找加热液体温度u(t) 的变化规律，使槽中的液体温度经1小时后上 升到40C，且散失热量最小。因假定槽内液体处于完全混合状态，由热力学可知，槽内液体温度x(t)的变化率与两种液体的温差(u(t)-x(t)成正比，为简便计，令比例系数为1，于是有 = u (t) - x(t)(9.1-1)dt初始条件x(0)=0(9 .1-2)终值条件x(1)=40(9 .1-3)在1小时内散失掉的热量可用下式表示J= j1 qx2(t)+ru2(t)dt(9.1-

3、4)0其中q和r都是正的常数。该问题的任务是，寻求加热液体温度的最优变化规律u*(t)，使槽内液体由 初态9.1-2转移到终态9.1-3的过程中散失热量9.1-4为最小。9.1.2最优控制问题的一般提法根据上例可概括出，最优控制问题用数学语言来描述时应该包含以下几方面 内容：1. 系统状态方程.X =fx(t),u(t),t(9.1-5)其中，x为n维状态向量，u为r维控制向量，f(.)是 n维函数向量。2. 控制变量的约束条件大多数实际控制系统中的控制变量的取值范围是受限制的，如发动机的推力， 电动机的转矩等都不能超出某一极限，即|u| K在数学上，表示容许控制域为控制空间中的一个集合ueU

4、3. 初始条件和终值条件在最优控制中，初始条件通常是已知的，即X(t0)=X0,而终值条件则要复杂些， 它可以是状态空间中一个确定的点，或状态空间中某一个点集(目标集)中的任 何一个点。表示为：X(t )=X f f如例9.1中，tf规定为1小时。有时终值状态也可能是状态空间中一个运动 的点，如用导弹拦截一个运动的目标，终值是被拦截物运动轨迹上可能的一个点， 此时终值受运动轨迹约束，可表示为。Ci(Xf、tf)=0, C=1, 2,、，m, m0,因此U-0)o式(9.2-2 )中的Q和R被称之为加权矩阵(Weight matrix),它们起着对 式中两个积分项重视程度的调节钮(tuning

5、knob )作用。即，若想重视状态的 控制(如高的响应)，增加Q的各元，若着重考虑控制能量消耗最小，则取较大 的R值。Q和R的选取尚无一般规律可循，一般取决于设计者的经验。常用的是所谓 试行错误法(trial and error method ),即选取不同的Q、日代入计算并通过比 较结果而确定。为计算方便起见，Q、R通常取为对角矩阵。一、问题求解对这个问题的求解可以采用多种方法，这里采用基于Liapunov第二方法的 求解方法。设最优控制律为：u(t)=-kx(t)(9.2-3)代入(9.2-1 )得X =Ax-BKx=(A-BK)x在下面的推导中，设矩阵A-BK是稳定的，即A-BK的特征值

6、具有负实部。将式(9.2-3)代入式(9.2-2)，得J=f ( xtQx+xtKtRKx)dt 0=f xt (Q+KTRK)xdt 0设对任意x都有xt(Q+KtRK)x= - (xtPx) dt式中p为正定实对称矩阵，则可导出xt(Q+KtRK)x=- xtPx-xtP x=-xt (A-BK) tP+P(A-BK)x由Liapunov第二方法可知，对于给定的正定矩阵Q+KtRK，如果矩阵A-BK是 稳定的，则存在正定矩阵P，使得(A-BK)tP+P( A-BK)=-( Q+KtRK)(9.2-4)由此可求出J为：J= f xt (Q+KTRK)xdt=-xTPx |0 0=-xt (

7、3 )p x( 3 )+xt (0)px ( 0 )由于假定A-BK的所有特征值都具有负实部，因此有x(3) T0。这样就得 到J=xt(0)Px(0)(9.2-5)即J可由初始条件x(0)和P得到。由于r为正定实对称矩阵，因此可将r写成R=TtT( 9.2-4 )式中T为一非奇异矩阵。这样式(9.2-4 )可写成(At-KtBt ) P+P ( A-BK ) +Q+KtTtTK=0上式经整理还可写成AtPA+TK- ( T ) -iBtP t TK- (T ) -iBtP-PBR-iBtP+Q=0J对k取极小值就要求下式对k取极小值。XT TK-(T) -iBtP t TK-(T) -iBt

8、Px由于上式是非负的，其最小值为零，即TK=( Tt ) -iBtP因此可求得最优反馈增益矩阵为K=T-1( T ) -iBtP=R-iBtP( 9.2-6 )最优控制u(t)为u(t)=-Kx(t)=-R-1BTPx(t)式(9.2-6 )中的矩阵P必须满足(9.2-4 )或者满足下列方程：AtP+PA-PBR-1BtP+Q=0( 9.2-7 )式(9.2-7 )称为Riccati(黎卡提)方程。二、matlab 实现MATLAB中提供的Lqr和lqry命令可以直接求解二次型最优控制问题(亦称 最优调节器，Optimal regulator )其格式为：K,P,E=Lqr(A,B,Q,R)K

9、,P,E=Lqry(A,B,C,D,Q,R)式中，K为最优反馈增益矩阵；P为对应Riccati方程的正定解，E为A-BK的特 征值，Lqry命令用于求解二次型调节器问题的特例，即目标函数中用输出Y来 代替状态X，也就是说，目标函数为。J=F (yTQy+uTRu)dt0例9. 2给定系统X =Ax+BuY=Cx 式中0-2-3C=11 0 0以高的响应性能为目标，求最优控制律U=-Kx使评价指标3J=j(xtQx+utRu)dt0取最小化解、取100 00Q=010 ;R=0.010 01并代入程序chp9_1.m% Program Chp9_1% Design of quadratic op

10、timal control system% Design of quadratic optimal control system% state matrics A,B,C,D% state matrics A,B,C,DA=0 1 0; 0 0 1; 0 -2 -3;B=0; 0; 1;C=1 0 0;D=0;% Weight matrices Q and RQ=100 0 0; 0 1 0; 0 0 1;R=0.01;% Obtain the optimal state feedback gain matrix KK=lqr(A,B,Q,R)% Optimal contral sysytem

11、 state matrix AA AA=A-B*K;% Unit-step response t=0:0.01:6;% Original systemfigure(1)y,x,t=step(A,B,C,D,1,t);plot(t,y)grid% Optimal contral systempausefigure(2)yo,xo,to=step(AA,B,C,D,1,t);plot(to,yo)grid程序运行结果：图9.2二次型最优控制9.3离散系统的二次型最优控制对于离散系统x(k+1)=Ax(k)+Bu(K)Y(K)=cx(k)+Duc(k)(9.3-1)如果考虑控制步数N-8,则该系统的

12、二次型最优控制问题的评价函数可表示 为J=尹xT(k)Qx(k)+UT(k)Ruc(k)(9.3-2)k=0同样，可证明使(9.3-2)式取极小值的最优控制解可表示为：u(k)=-kx(k)(9.3-3)式中 K=(R+BtPB)-iBtPA-Q=0其中，P是Riccati方程P-AtPA+AtPB( R+BtRB)-1BtPA-Q=0的解MATLAB中的K,P,E=dlqr(A,B,Q,R)K,P,E=dlqr(A,B,C,D,Q,R)等命令支持离散系统二次型最优控制的计算机辅助设计。9.4状态观测器与最优观测器设计状态反馈是现代控制中的一个重要方法，在本章的最优控制、前面讨论的极 点配置以

13、及以后的一些控制方法讨论中，都离不开状态反馈。这意味着、在控制 过程中，需要对各个状态变量变化有很好的了解，而在实际控制系统中，这也意 味着需要对控制变量进行在线检测。然而，实际系统中并非所有的状态变量都很 容易被检测，有的状态变量可能并非物理量，有些状态变量可能很难或根本无法 用物理方法去测量，有时即使状态变量可以被测量，但传感器的过多使用也可能 导致系统性能的下降与成本的上升。因此，要使状态反馈方法能够很好地在工程 中得以实现，需要解决这个问题。本节所讨论的状态观测器便是解决这一问题的 方法之一。所谓状态观测器是一个在物理上可实现的动力学系统，它在待观测系统的输 入和输出的驱动下(这总是可

14、以测量到的)产生一组逼近于待观测系统状态变量 的输出。该动力学系统装量所输出的一组状态变量便可作为该观测系统的状态的 估计值。正因为此，状态观测器也被称为状态估计器。9. 4. 1状态观测器原理对于系统x(t)=Ax(t)+Bu(t)(9 .4 -1)y(t)=Cx(t)考虑一个动力学模型(t ) =A w (t)+Ky(t)+Bu(t)x (t)=Dw (t)+Hy(t)(9.4-2)关联接系统(9.4-1 )与(9.4-2)如图9.3所示图9.3状态观测器结构 如果C1.存在一个矩阵M,并且A M=MA-KC(9.4-3)In=DM+HC( 9.4-4 )和 B =MB(9.4-5)C2.

15、A是一个稳定矩阵它将导致 (t) -x(t) (t-8)模型(9.4-2 )被称为兰恩伯格观测器(Luenberger observer ) 证明：定义一个误差向量e(t)= w (t)-Mx(t)(9.4-6)那么 e (t)= - (t)-M (t)=A w (t)+KCx(t)+ B u(t)-MAx(t)+Bu(t)=A w (t)-(MA-KC)x(t)从而有e (t)= A w (t)- AMx(t)=A e(t)所以 e(t)=exp( A t)e(o)(9.4-7)另外，定义n( t ) = x (t)-x(t)(9.4-8)于是 n (t ) =D w (t)+HCx ( t

16、 ) -x ( t )=Dw (t)-IHCx ( t )=Dw (t)-DMx ( t )=De ( t )所以 n (t ) =Dexp( A t)e(o)因为A是一个稳定矩阵，它将导致n( t )t 0 (t T8)从而有 X (t) Tx(t) (tT8)得证如果让 M=D=In ,H=0于是 A =A-KC, B =BX (t)= w (t)从而 x ( t ) = A X (t)+Ky(t)+ B u(t)即 X (t) = (A-KC) X (t)+Ky(t)+ B u(t)(9.4-9)式(9.4-9 )所示系统即是常见的状态观测器(国外教科书上称为 identify obse

17、rver )9.4.2状态观测器的设计对式(9.4-9 )进行分析可知，状态观测器设计的关健在于矩阵K的求取。 一般可以通过两种方法来求取K。一是给出期望极点矢量求取，这种方法在许多 教科书上均有介绍。另一种是应用Kalman滤波进行最优估计求取，本节中介绍 后一种方法。实际系统在控制过程中，在得到状态反馈的同时，也会不可避免地受到量测 噪声的干绕，所谓最优估计问题，就是指从带有量测噪声的检测状态中获得状态 的最优估计。考虑如下系统x (t)=Ax(t)+Bu(t)+Gv(t)y(t)=cx(t)+w(t)式中，v ( t )为随机噪音干扰输入(System noise ); w(t)为由传感

18、器带来 的随机测量噪音(measurement noise )。进一步，假设v(t)及w(t)为白噪声且 互不相干，于是有以下统计特性：Ev(t)=0;Ew(t)=0 Ev(t)wT (t)=0Ev(t)vT (t)=R kEw(t)wT (t)=Q k这里 Rk0Qk0现在，对于这个包含有噪音v(t)和w(t)的系统，需要做的是根据数据u(t) 和y(t)，找出状态变量x(t)的最优估计值X (t)。参照先前有关兰恩伯格观测器的讨论，可以假设：X (t) = (A-LC)x(t)+Ly(t)+Bu(t)X (o)=xo(9.4-10)式(9.4-10 )意味着，观测器的设计取决于矩阵L的确定

19、，这里，L被称为 增益矩阵。建立一个评价指标：J=Ex(t)- X (t)Tx(t)- X (t)(9.4-11)于是，有关观测器的设计问题可转化为对(9.4-11 )求极小值的问题，这个 过程被称为对系统(9.4-10 )的Kalman滤波。该问题的解被求得为L=PCtR1(9.4-12)式中，L为最优增益矩阵，或称为Kalman增益矩阵，P是下述Riccati方程 的正定解PAT+AP-PCTRk-1CP+Qk=0于是，通过应用最优增益矩阵L,已知输入u (t )和被测得的输出y(t)可构 成一个观测器为X (t) = A-LC X (t)+Ly(t)+Bu(t)J (t) (t)(9.4

20、-13)在MATLAB中，可以用命令Lqe来求解最优滤波器的参数，其格式为L, P，E=Lqe( A，B，C，Q%).式中，L为Kalman增益矩阵，P为对应Riccati方程的解，E为估计器的闭 环特征值。命令estim用来构成连续系统Kalman滤波器，命令格式为Ae,Be,Ce,De=estim(A,B,C,D,L)根据式(9.4-13)可知，上式所求得的滤波器为含有测量噪音的传感器输出Y，而滤波器的作用是根据Y求出估计值J和状态X的估计值，而当应用Kalman 滤波器作状态观测器来使用时，对一些已可以用传感器检测的状态值已不必再去 估计，系统的输入中也还可能有其它确定性的控制输入和外部

21、给定输入，这时， 可以应用更一般的形式的estim命令格式，如下所示。Ae,Be,Ce,De=estim(A,B,C,D,L, sensors, known)式中，Sensors是一个包含传感器输出序号的矢量，Known是包含外部输入矢量 序号的矢量。有关这种情况我们将在下一节讨论，此处，我们先就比较简单的情 况讨论状态估计器的设计。例9.3给定系统 (t)=Ax(t)+Bu(t)Y(t)=Cx(t)式中0-18162 -TA二100，B=0_01-100C=0 18 -162设计系统的状态估计器并进行仿真。解：首先假设系统的输入噪音为零均值，方差为1的高斯白噪声，传感器输 出的测量噪音为零均

22、值，方差为0.01的高斯白噪声，并利用randn命令形成随 机噪音对系统进行仿真，然后根据Rk,Qk和命令Lqe计算最优增益矩阵L，利用 estim命令形成最优估计器Ae,Be,Ce,De。再进行状态值和状态估计值的仿真并绘出图形进行比较。设计程序为chp9-2.m,% Design of the estimera0=0 -18 162; 1 0 0; 0 1 -10;b0=1 0 0;c0=0 18 -162;d0=0;dimf=size(b0);noinp=dimf(2);dimc=size(c0);noout=dimc(1);q0=1;w=sqrt(q0)*randn(100,noinp

23、);r0=0.01;nu=sqrt(r0)*randn(100,noout);t=0:0.1:(100-1)*0.1;y,x=lsim(a0,b0,c0,d0,w,t);y=x*c0，+nu;L,P,E=lqe(a0,b0,c0,q0,r0);L,P,Eae,be,ce,de=estim(a0,b0,c0,d0,L);xh=lsim(ae,be,ce,de,y,t);subplot(2,2,1);plot(t,x(:,1),t,xh(:,2)title(，(1) x1 & xh1 versus t，)subplot(2,2,2);plot(t,x(:,2),t,xh(:,3)title(，(2) x2 & xh2 versus t，)subplot(2,2,3);plot(t,x(:,3),t,xh(:,4)title(，(3) x3 & xh3 versus t，)xerr=x-xh(:,2:4);subplot(2,2,4),plot(t,xerr)title(，(4) xerr versus t，)x2、图9.4为状态与其估计值的对照曲线，图中(1 )、( 2 )、( 3 )分别x1、x1； x2； x3 x3的对照，(4)为估计误差x=x-x的情况。(1) x1 & xh1 versus t图9.4状态与其估计值的对照曲线