最优控制习题及参考答案

合集下载

最优控制胡寿松版部分习题答案

最优控制胡寿松版部分习题答案

2-5 求通过(0)1x =,(1)2x =,使下列性能泛函为极值的极值曲线*()x t :2(1)ft t J x dt =+⎰解:由题可知,始端和终端均固定被积函数21L x =+,0L x ∂=∂,2L x x ∂=∂, 2d L x dt x∂⋅=∂ 代入欧拉方程0L d L x dt x∂∂-⋅=∂∂,可得20x =,即0x = 故1x c = 其通解为:12x c t c =+代入边界条件(0)1x =,(1)2x =,求出11c =,21c = 极值曲线为*()1x t t =+2-6 已知状态的初值和终值为(1)4x =,()4f x t =式中f t 自由且f t >1,试求使下列性能泛函达到极小值的极值轨线*()x t :211[2()()]2ft J x t x t dt =+⎰ 解:由题可知,2122L x x =+,()4f t ψ=,()14x =,()4f x t = 欧拉方程:L 0d L x dt x∂∂-=∂∂ 横截条件:()00t x =x ,()()f f x t t ψ=,()0fTt L L x x ψ∂⎛⎫+-= ⎪∂⎝⎭易得到2dxdt= 故12x t c =+ 其通解为:()212x t t c t c =++根据横截条件可得:()()()122121114424f f f f f x c c x t t c t c x t t c ⎧=++=⎪⎪=++=⎨⎪=+=⎪⎩解以上方程组得:12569f t c c =⎧⎪=-⎨⎪=⎩将f t ,1c ,2c 代入J 可得5*201500502150233J x x dt =+=-=⎰ 极值轨线为()*269x t t t =-+2-7 设性能泛函为120(1)J x dt =+⎰求在边界条件(0)0x =,(1)x 自由情况下,使性能泛函取极值的极值轨线*()x t 。

解:由题可知,21L x =+,()00x =,()1x 自由欧拉方程:L 0d L x dt x∂∂-=∂∂ 横截条件:()00t x =x ,L 0ft x∂=∂,0fTt L L x x ∂⎛⎫+= ⎪∂⎝⎭易得到()x t a =其通解为:()x t at b =+代入边界条件()f x t a =,()00x =,1f t =,求出0a =,0b = 将f t ,a ,b 代入J 可得()1*211J x dt =+=⎰极值轨线为()*0x t = 2-9 求使泛函22211220(2)J x x x x dt π=++⎰为极值并满足边界条件1(0)0x =,2(0)0x =1()12x π=,2()12x π=- 的极值轨线*1()x t 和*2()x t 。

最优控制习题答案

最优控制习题答案

最优控制习题答案最优控制习题答案最优控制是一门研究如何在给定的约束条件下,使某个系统的性能指标达到最优的学科。

在实际应用中,最优控制被广泛应用于工程、经济、生态学等领域。

然而,最优控制问题通常非常复杂,需要运用数学方法进行求解。

在本文中,我们将探讨一些最优控制习题,并给出相应的答案。

习题一:一辆汽车行驶在一条直线上,其速度v(t)满足以下微分方程:m*dv(t)/dt = F(t) - kv(t),其中m为汽车的质量,F(t)为外部施加的力,k为阻力系数。

求使得汽车行驶时间最短的外部力F(t)。

解答:首先,我们需要确定行驶时间的数学表达式。

设汽车的初始速度为v0,行驶时间为T,根据题意,我们可以得到以下约束条件:v(0) = v0,汽车的初始速度为v0;v(T) = 0,汽车的最终速度为0。

根据最小时间原理,我们可以建立一个最优控制问题的数学模型,即求解以下极值问题:minimize T,使得v(T) = 0;subject to m*dv(t)/dt = F(t) - kv(t),v(0) = v0。

通过拉格朗日乘子法,我们可以得到最优控制问题的解析解。

最终,我们可以得到外部力F(t)的解析表达式。

习题二:一个农民想要将一块矩形土地分成两块,使得两块土地的总面积最大。

该农民只能在土地的一条边上建立一道直线围栏。

求最优划分方案。

解答:设矩形土地的长为a,宽为b。

我们需要确定如何划分土地,使得两块土地的总面积最大。

根据题意,我们可以得到以下约束条件:2a + b = L,L为围栏的长度。

根据最大面积原理,我们可以建立一个最优控制问题的数学模型,即求解以下极值问题:maximize A = ab,使得2a + b = L。

通过拉格朗日乘子法,我们可以得到最优控制问题的解析解。

最终,我们可以得到最优划分方案的解析表达式。

习题三:一架飞机要从A地飞往B地,途中需要经过一个位于C地的雷达站。

飞机的速度恒定为v,雷达站可以通过调整飞机的航向角度来监测飞机的位置。

最优控制第五章习题答案

最优控制第五章习题答案

1. ·2.已知二阶系统的状态方程122()(),()()x t x t x t u t ==性能泛函3222221212120111[(3)2(3)][2()4()2()()()]222J x x x t x t x t x t u t dt =+++++⎰求最优控制。

解:把状态方程和性能指标与标准状态方程和标准性能指标比较,可得0,101,02,11,,,,0,010,21,42A B P Q R ⎡⎤⎡⎤⎡⎤⎡⎤=====⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦考虑到()K t 是对称阵,设11121222,(),k k K t k k ⎡⎤=⎢⎥⎣⎦代入黎卡提方程1()()()()()()()()()()()T T K t K t A t A t K t K t B t R t B t K t Q t -=--+-即1112111211121112111212221222122212221222,,,,,0,10,002,12[0,1],0,01,0,,1,1,4,k k k k k k k k k k k k k k k k k k k k ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤=--+-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦令上式等号左右端的对应元相等,得211121211122222212222221224k k k k k k k k k =-=-+-=-+-这是一组非线性微分方程。

由边界条件(3)K P =即11121222(3),(3)1,0(3),(3)0,2k k k k ⎡⎤⎡⎤=⎢⎥⎢⎥⎣⎦⎣⎦ 最优控制为11112112122212222()()(),()2*[0,1]2()2(),()T u t R B K t X t k k x t k x t k x t k k x t -=-⎡⎤⎡⎤=-=--⎢⎥⎢⎥⎣⎦⎣⎦3. )4.能控的系统状态方程为122()(),()()x t x t x t u t ==这是一种双积分系统,其输出为1()x t ,其输入为()u t ,其传递函数为12()1()()x s G s u s s==其性能泛函为222112201[()2()()()()]2J x t bx t x t ax t u t dt ∞=+++⎰其中220a b ->求最优控制。

最优控制期末考试题精华版

最优控制期末考试题精华版
九、(10分)用 算法求下图中s到t的最短路及长度:
解:
设拉格朗日乘子为 、 ,根据非线性K-T条件定理:
考虑到约束条件:
分三种情况考虑式(1)~(5)的解:
1.若(5)式等号不成立,则由(2)式有 ,再代入(1)式得 ,这和(4)式矛盾。因此,(4)式等号一定成立。
2.若(6)式等号不成立,则由(3)式有 ,代入(1)式得
由 和(7)式以及(5)式可得: .
3.若(6)式等号成立,则由(5)、(6)式可解得:

由(4)式、(1)式可知, 不能取 ,而若取 ,则 就应该为 ,此时(1)式不成立。
综上所述,所求非线性规划有唯一K-T点 ,根据最优条件,最优点必定是K-T点,即 。
七、(10分)用惩罚函数法求下列问题的最优解:
八、(10分)试求 的有效解和弱有效解
一、(10分)用图解法求解:
二、(10分)用单纯形法求解:
三、(10分)求函数 的极值。
四、(15分)用0.618法解下列问题:
,初始区间为 ,
五、(15分)设 ,求在以下各点处的最速下降方向:
, ,
解:函数在某一点处的最速下降方向即为函数在该点方向导数最小的方向

由已知条件可得:
所的表达式求出它们的最优解:

最优控制测试题

最优控制测试题

1. 已知受控系统1)0( ),()(==x t u t x,试求)(t u 和f t ,使系统在f t 时刻转移到坐标原点0)(=f t x ,且使dt t u t J ft f ⎰+=022)(最小。

2.设系统状态方程为),()(11t u t x= 0)0()0(21==x x , (t))()(212u t x t x += , 1)1()1(21==x x 。

其中)(1t u 无约束,41)(2≤t u 。

求)(1t u *,)(2t u *,)(1t x *,)(2t x *,使得系统从的0=t 初始态转移到1=t 的末态,并使性能指标 dt t u t u t x J ⎰++=122211)]()()([为最小值。

3. 在重量为10 kg 的禁止物体上,加一个垂直方向的力)(t F ,物体允许的最大加、减速度为5 m/2s ,欲使物体以最短的时间升高100 m ,并停留在这一位置,求)(t F 的变化规律,并求出最短时间。

取重力加速度为2/10s m g =。

4. 生产库存的状态方程为)()()()1(k s k u k x k x -+=+,其中)(),(),(k s k u k x 分别是库存量,生产速度和销售速度。

设生产费用为)(005.02k u ,库存费用为)(k x ,那么四季度的总费用为)]()(005.0[302k x k u J k +=∑=现设初始库存量0)0(=x ,四季度的订货分别为1200)3(,500)2(,700)1(,600)0(====s s s s 。

求最优生产速度)(k u )3,2,1,0(=k 使0)4(=x (年底无积压),并且使总费用J 最小。

5. 已知系统方程为)()( ),()(221t u t x t x t x== ,其中1)(≤t u ,试求最优控制)(t u *,把系统从初始状态1)0( ,1)0(21==x x 转移到原点,且使性能指标dt t u J ft ⎰+=0])(4[取最小。

最优控制课后习题答案

最优控制课后习题答案

最优控制课后习题答案最优控制课后习题答案最优控制是现代控制理论中的重要分支,它研究如何在给定约束条件下,使系统的性能指标达到最优。

在最优控制的学习过程中,课后习题是巩固理论知识、培养解决问题能力的重要环节。

本文将为大家提供一些最优控制课后习题的答案,希望能对大家的学习有所帮助。

1. 线性二次型最优控制问题考虑一个线性时不变系统,其状态方程和性能指标分别为:$$\begin{align*}\dot{x}(t) &= Ax(t) + Bu(t) \\J(u) &= \int_{0}^{T} (x^T(t)Qx(t) + u^T(t)Ru(t))dt\end{align*}$$其中,$x(t)$为系统的状态向量,$u(t)$为控制输入向量,$A$和$B$为系统矩阵,$Q$和$R$为正定矩阵,$T$为最优控制的时间段。

求解该问题的最优控制输入$u^*(t)$。

答案:根据最优控制的原理,最优控制输入$u^*(t)$满足以下的最优性条件:$$\begin{align*}\frac{\partial J}{\partial u}(u^*(t)) &= 2R u^*(t) + 2B^T P(t)x(t) = 0 \\\dot{P}(t) &= -PA - A^T P - Q + PBR^{-1}B^T P\end{align*}$$其中,$P(t)$为状态向量的共轭变量矩阵。

通过求解上述的代数方程和微分方程,可以得到最优控制输入$u^*(t)$和状态向量的共轭变量矩阵$P(t)$。

2. 非线性最优控制问题考虑一个非线性系统,其状态方程和性能指标分别为:$$\begin{align*}\dot{x}(t) &= f(x(t), u(t)) \\J(u) &= \int_{0}^{T} g(x(t), u(t)) dt\end{align*}$$其中,$f(x(t), u(t))$为非线性函数,$g(x(t), u(t))$为性能指标函数。

最优控制作业

最优控制作业
1 x 2ak 2 x 0 1 其值为 0,最大值为正无穷大;当 a 0 且 k 0 时, y 0 和 y , 2 ak 0 1 k 1 x 2ak 2 x 0 1 是 f ( x, y) 的最大值点, y 0 是 f ( x, y) 的极小值点,其值为 0, y 2ak 0 1 k
解:不显含 x ,由欧拉方程可得: L y
L c1 ,代入整理可得 y 2 (1 y 2 ) c1 ,进 y
一步得 y
c1 积分整理可得 y 2 c1 ( x c2 )2 , 再由横截条件:y(0) y0 , 1 , y2
(( y )
2)
L L) 0 , 即为 ( x 9)2 2 y 2 9 ,联立解得。 y xf
(( y )
L L) 0 , 即为 y x 2 5 ,联立解得。 y xf
第4周 1、求最优轨线及控制规律
11 a、 系统 x x u , J [u ] u 2 dt , x(0) 10 , x(1) 0 ,求最优 x (t ) , u (t ) 0 2 11 x x2 x (0) 0 b、系统 1 ,J [u ] u 2 dt , 1 , 终点约束 x1 (1) x2 (1) 1, 求 x (t ) , 0 2 x2 u x2 (0) 0
z (c1 2c4 )sin x (c2 2c3 ) cos x c3 x sin x c4 x cos x ,其中 c1 , c2 , c3 , c4 为待
定系数;若给出具体的横截条件即可解得待定系数。 4、试写出最优问题模型 球面两点最短路径问题和定周面积最大问题 解:(1)、设过球面两定点 x0 ( x(t0 ), y(t0 ), z (t0 )) 和 x f ( x(t f ), y(t f ), z(t f )) 的参数 曲线为 x x(t ) , y y(t ) , z z (t ) ;球面半径为 R ,则球面两点最短路径问题 可 描 述 为 求 泛 函 J [ x(t ), y(t ), z (t )]

最优控制习题及参考答案

最优控制习题及参考答案

1 2f最优控制习题及参考答案习题 1 求通过 x (0) = 1 , x (1) = 2 ,使下列性能指标为极值的曲线:t f J = ∫ (x2 +1)dt t 0解: 由已知条件知: t 0 = 0 , t f = 1d由欧拉方程得: (2x ) = 0dtx = C 1x = C 1t + C 2将 x (0) = 1,x (1) = 2 代入,有:C 2 = 1,C 1 = 1得极值轨线: x *(t ) = t +1习题 2 求性能指标: J =∫1(x 2 +1)dt在边界条件 x (0) = 0 , x (1) 是自由情况下的极值曲线。

解:由上题得: x *(t ) = C t + Cx * (t )由 x (0) = 0 得: C 2 = 0∂L由∂xt =t f= 2x (t f ) = 2C 1 t =t= 0 t于是: x *(t ) = 0【分析讨论】对于任意的 x (0) = x 0 ,x (1) 自由。

20 1 0 ∫ ⎩ λ= −λ有: C = x , C = 0 ,即: x *(t ) = x 其几何意义: x (1) 自由意味着终点在虚线上任意点。

习题 3 已知系统的状态方程为: x1 (t ) = x2 (t ) , x 2 (t ) = u (t )边界条件为: x 1 (0) = x 2 (0) = 1 , x 1 (3) = x 2(3) = 0 ,31 试求使性能指标 J =u 2(t)dt2取极小值的最优控制 u * (t ) 以及最优轨线 x *(t ) 。

⎡ x ⎤解:由已知条件知: f = ⎢ 2⎥⎢⎣ u ⎥⎦Hamiton 函数: H = L + λT fH = 1u 2 + λ x+ λ u⎧λ = 0 由协态方程: ⎨ 12 1 21 22⎧λ = C① 得: ⎨11⎩λ2 = −C 1t + C 2②∂H由控制方程: ∂u= u + λ2 = 0得: u = −λ2 = C 1t − C 2 ③由状态方程: x 2 = u = C 1t − C 2得: x (t ) = 1C t 2− C t + C ④22由状态方程: x 1 = x 21 2 3得: x (t ) = 1 C t 3− 1C t 2+ C t + C ⑤16 122 3 41⎪⎩=− ∫⎡1⎤ ⎡0⎤将 x (0) = ⎪ ⎪ , x (3) = ⎪0⎪ 代入④,⑤,⎣1⎦ ⎣ ⎦10联立解得: C 1 =由③、④、⑤式得:u * (t ) = 10t− 29, C 2 = 2 , C 3 = C 4 = 1 9x * (t ) = 5 t 3−t 2 + t +1 27 x *(t ) = 5 t 2 − 2t +1 29习题 4 已知系统状态方程及初始条件为x =u , x (0) = 1试确定最优控制使下列性能指标取极小值。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

[最优控制习题及参考答案习题1 求通过x(0) = 1 ,x(1) = 2 ,使下列性能指标为极值的曲线:J = ∫(x+1)dt解:由已知条件知:t= 0 ,t = 1d由欧拉方程得:(2x) = 0dtx= Cx = Ct + C将x(0) = 1,x(1) = 2 代入,有:C= 1,C= 1得极值轨线:x(t) = t +1习题2 求性能指标:J =(x+1)dt∫在边界条件x(0) = 0 ,x(1) 是自由情况下的极值曲线。

x(t)解:由上题得:x(t) = C t+ C Array由x(0) = 0 得:C= 0∂L由= 2x(t ) = 2C= 0 t ∂于是:x(t) = 0【分析讨论】对于任意的x(0) = x,x(1) 自由。

∫ ⎩ λ=有: C = x , C = 0 ,即: x (t ) = x 其几何意义: x (1) 自由意味着终点在虚线上任意点。

习题 3 已知系统的状态方程为: x(t ) = x (t ) , x (t ) = u (t )边界条件为: x (0) = x (0) = 1 , x (3) = x (3) =0 ,1试求使性能指标 J =u (t )d t2取极小值的最优控制 u (t ) 以及最优轨线 x (t ) 。

⎩ x ⎩ 解:由已知条件知: f = ⎩ ⎩⎩⎩ u ⎩⎩Hamiton 函数: H = L + λfH = 1u + λ x+ λ u⎩λ = 0由协态方程: ⎩2⎩λ = C①得: ⎩⎩λ= −Ct + C ② ∂H由控制方程: ∂u= u + λ= 0得: u = −λ= Ct − C ③由状态方程: x = u = Ct − C得: x (t ) = 1C t − C t + C④2由状态方程: x = x得: x (t ) = 1 C t − 1C t + C t + C⑤6 2⎩ ⎩=− =− ∫⎩1⎩ ⎩0⎩将 x (0) = ⎩ ⎩ , x (3) = ⎩0⎩ 代入④,⑤,⎩1⎩ ⎩ ⎩10联立解得: C =由③、④、⑤式得:u (t ) = 10t− 29,C = 2 , C = C = 1 9x (t ) = 5t −t + t +1 27 x (t ) = 5t − 2t +1 9习题 4 已知系统状态方程及初始条件为x =u , x (0) = 1试确定最优控制使下列性能指标取极小值。

J = 解: H = xe + ue + λu⎩x= u 列方程: ⎩λ= −2xe ⎩2eu + λ = 0 (x + u )edt①②③由③得, u代入①得,x 1 e λ④21e λ=−2x 1 e λ e λ =− +2将②,③代入,并考虑到 u = xx 1 e (−2xe ) + e (−2ex ) 2整理可得: x + 2x − x = 0) = u =⎩特征方程: s + 2s −1 = 0s = −1+s = −1−于是得: x (t ) = C e + C eλ(t③ −2e ①−2exλ(t ) = −2e(Cse + C s e )由 x (0) = 1 ,得: C + C = 1⑤由 λ(t ) = λ(1) = 0 得: Cse+ Cse = 0⑥⑤、⑥联立,可得 C 、C代回原方程可得 x → u(略)习题 5 求使系统: x = x , x = u由初始状态 x (0) = x (0) = 0出发,在 t = 1 时转移到目标集1 x (1) + x (1) = 1,并使性能指标 J =∫u (t )dt2为最小值的最优控制 u (t ) 及相应的最优轨线 x (t ) 。

解: 本题 f (i),L (i) 与习题 3 同,故 H (i) 相同→方程同→通解同⎩λ= C ,λ= −Ct + C⎩⎩x = 1 C t − 1 C t + C t + C 有: ⎩ 6 2⎩x = 1 C t −C t + C ⎩2⎩⎩u = Ct −C∫⎩0⎩ x (0) = ⎩ ⎩⎩0⎩由,有: C = C = 0①由 x (1) + x (1) = 1,有:1 C– 1 C+ 1 C − C = 1 6 2 22 C −3 C = 1②3 2∂ϕ ∂ψ由 λ(1) = + ⋅γ = 0 ,ψ = x + x −1∂x ∂x⎩1⎩有: λ(1) = ⎩ ⎩ γ = 0 ⇒ λ (1) = λ (1) ⎩⎩1⎩⎩于是: C = −C + C2C = C ③3 6 ②、③联立,得: C =- 、C = -77于是: u = − 3 t + 67 7x =− 1 t + 3 t 14 7x =− 3 t + 6t 14 7习题 6 已知一阶系统: x (t ) = −x (t ) + u (t ) , x (0) = 3(1)试确定最优控制 u (t ) ,使系统在t = 2 时转移到 x (2) = 0 ,并使性能泛函J = (1+ u )dt = min⎩ ⎩⎩ C C4(2)如果使系统转移到 x (t ) = 0 的终端时间 t 自由,问 u (t ) 应如何确定解:H = 1+ u + λu − λ x⎩x = −x + u 列方程: ⎩λ= λ ⎩2u + λ = 0由协态方程得: λ = C e ①1由控制方程: u=− Ce ②2① t1 代入状态方程: x = −x −Ce2= 2,x (2) = 0 ⇒ x (t ) = Ce – 1 C e 4⎩− 1 C = 3 ⎩ 4 ⎩⎩C e − 1 C e = 0 ⎩⎩ 412解得: C = ,e −13eC =e −1 代入②得: u (t ) =−② x (t ) = 2,t 自由6ee −1⎩− 1 C = 3 ⎩ 4 ⎩ C e – 1C e = 0 ⎩⎩⎩H (t ) = 0 ⎩ ⎩解得:C=40 − 6 =∫u (t ) = −习题 7 设系统状态方程及初始条件为x (t ) = u (t ) , x (0) = 1试确定最优控制 u (t ) ,使性能指标1J = t + 2 u dt 为极小,其中终端时间 t 未定, x (t ) = 0 。

解: H = 1u + λu2由协态方程得: λ = 0→ λ = C ①由控制方程: u + λ = 0→ u = −C ②由状态方程: x = u =−C⇒ x (t ) = −Ct + C ③由始端: x (0) = 1→ C = 1由末端: x (t ) = 0→ −Ct +1 = 0④∂ϕ 考虑到: H (t ) =− t– ∂ψ t ⋅γ = −1∂ ∂1有: u + λu = −12 1C − C = −1 ⇒ C = 22C=⑤当C=时,代入④有:t = 1C⎩ 6当 C = −时,代入④有: t= 1= 1,不合题意,故有2C 最优控制u = −习题 8设系统状态方程及初始条件为x (t ) = x (t ) , x (0) = 2 性能指标为x (t ) =u (t ) ,J =1∫ udtx (0) = 12要求达到 x (t ) = 0 ,试求:(1) t= 5 时的最优控制 u (t ) ;(2) t 自由时的最优控制 u (t ) ; 解:本题 f (i ),L (i ),H (i ) 与前同,故有 ⎩ ⎩λ= C ⎩λ= −Ct + C ⎩x = 1C t − 1 C t + C t + C ⎩ ⎩ 2⎩x = 1C t − C t + C ⎩ 2 ⎩⎩⎩u = Ct − C⎩2⎩⎩0⎩ ⎩C = 2 ⎩C = 1 ⎩125 25 ① 由 x (0) = ⎩ ⎩x (5) = ⎩0⎩ ,得: ⎩ C − C + 5C + C = 0⎩1⎩ ⎩ ⎩ ⎩ 62 ⎩ 25C−5C + C = 0⎩⎩26 ∫ 2联立得: C = ,C = ,⇒ u = −② t 自由⎩⎩C = 1⎩ ⎩C = 2 ⎩ 1C t − 1 C t + C t + C = 0 ⎩ ⎩ 2⎩ 1C t − C t+ C = 0 ⎩ 2 ⎩ ⎩⎩H (t ) = 0联立有: C t − 2C t+ 2 = 0 , 无论C 为何值, t 均无实解。

习题 9 给定二阶系统x (t ) = x (t ) + 1 , x (0) = − 14 4 1x (t ) = u (t ) , 1x (0) = −4控制约束为 u (t ) ≤ ,要求最优控制 u (t ) ,使系统在 t = t 2并使时转移到 x (t ) = 0 ,其中 t 自由。

J = u (t )dt = min解: H = u + λ x + 1λ+ λ u4⎩− 1λ λ ≤ 1 ⎩ 2 ⎩ 本题属最小能量问题,因此: u (t ) = ⎩−1λ > 1 ⎩⎩⎩ 1λ< −1⎩2⎩=− ⎩⎩λ= 0 → λ = C 由协态方程: ⎩⎩⎩λ = −λ → λ = −C t + C λ是 t 的直线函数。

当 u (t ) = − 1 λ = 1C t − 1 C时(试取) 2 2 2x (t ) = 1 C t − 1C t + C42x (t ) =1 C t − 1 C t + 1t + C t + C 124 41由始端条件 → C = C =4由末端条件 → 1 C t −1 C t+ 1 t + 1= 0 12 42 4 1 C t − 1 C t+ 1 = 0 42 4 另: H (t ) = 01联立解得: C =,C = 0,t = 39于是,λ1 t⎩λ= 1时,t < 0=−⎩ 9⎩λ= −1时,t = 9在 t 从 0 → 3 段, λ≤ 1满足条件。

故, u1 λ = 1 t2 180 1 2 3 4 t⎩−习题 10 设二阶系统x (t ) = −x (t ) + u (t ) , x (0) = 1 x (t ) = x (t ) , x (0) = 0控制约束为 u (t ) ≤ 1 ,当系统终端自由时,求最优控制 u (t ) ,使性能指标J = 2x (1) + x (1)取极小值,并求最优轨线 x (t ) 。

解:由题意, f ⎩−x + u ⎩ = , ϕ = x+ x ,L = 0 , ⇒ H = λ u − λ x + λ x ⎩ ⎩⎩ x ⎩由控制方程可得: u =⎩+1 ⎩ λ< 0 λ> 0 ⎩λ = λ − λ ⇒ λ = C e + C 由协态方程可得: ⎩⎩λ= 0 ∂ϕ ⎩2⎩⇒ λ= C由 λ(t ) = = ⎩ ⎩ ⇒ C = 1,C= e ∂x (t )⎩1⎩⎩λ = e +1 → 在t > 0的围λ > 1⇒ ⎩故: u = −1t ∈[0,1]⎩λ= 1若需计算 最优轨线 ,只需把 u = −1 代入状态 方程,可 得:⎩x (t ) = 2e −1 ⎩ ⎩ x (t ) = −2e − t + 2 ⎩⎩习题11 设系统状态方程为x(t) = x(t) ,x(0) = x∫⎩1 0⎩ ⎩0 0⎩ ⎩0⎩ !性能指标为 J = 12x (t ) =u (t ) ,(4x +u )dt x (0) = x试用调节器方法确定最优控制 u (t ) 。

相关文档
最新文档