第二篇最优控制理论习题答案
最优控制胡寿松版部分习题答案
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 。
最优控制期末
. . . . . . . .
. . . . . . .
. .
. .
. .
. .
. .
.
课后习题解答
最优控制理论与系统
December 28, 2014
2 / 46
概念及定理
概念
概念
1
求使泛函 J 取极值的极值轨线 x∗ (t),其中 x(0) = 1, ∫ 1 ˙ 2 ) dt J= (1 + x
0
解:本题为定常系统,tf 固定,末端自由,积分型性能指标的泛函 极值问题。
. . .
. .
.
. . . . . . . .
. . . . . . . .
. . . . . . .
. .
. .
. .
. .
. .
.
课后习题解答
最优控制理论与系统
December 28, 2014
2 / 46
概念及定理
0
解:本题为末端时刻自由,初态和末态两端固定、积分型性能泛函的最 优控制问题。令哈密顿函数为: H = 1 + u2 + λ(−x + u), 由协态方程 ˙ (t) = − ∂ H = λ(t) ⇒ λ(t) = c1 et λ ∂x
. . . . . .
. . . .
. . . . . . . .
. . .
. .
.
. . . . . . . .
. . . . . . . .
. . . . . . .
. .
. .
. .
. .
. .
.
课后习题解答
最优控制理论与系统
December 28, 2014
最优控制习题及参考问题详解
标准文档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 + C由 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) 自由。
2 0 1∫⎩ λ = −λ有: 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 )dt 2取极小值的最优控制 u *(t ) 以及最优轨线 x *(t ) 。
⎡ x ⎤解:由已知条件知: f = ⎢ 2⎥⎢⎣ u ⎥⎦Hamiton 函数: H = L + λT f H = 1u 2 + λ x + λ u⎧λ = 0由协态方程: ⎨ 12 121 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 ) = 1C 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 +127 x *(t ) = 5 t 2 − 2t +1 29习题 4 已知系统状态方程及初始条件为x =u , x (0) = 1试确定最优控制使下列性能指标取极小值。
最优控制习题参考解答
而
再由条件
由
所以, ,
3-6 已知线性系统的状态方程 ,其中
, , ,
给定 , ,求 ,使泛函 最小。
解:
系统状态方程: ,
代入初始条件得:
①
②
由①②得, , ,所以
3-7
3-8 设 为平面中一个质点,它自始点 向原点 的转移过程中,其瞬时速度为该质点所在位置的函数 ,质点的运动方程为
§2.6 习题
2.2
解:
已知 ,
当 ,
当 ,
2.4
解:
设 ,
,最短曲线为
2.5
解:
, , ,
,
又由横截条件得:
, 极值轨线为
2.6
解:
, ,
, 极值轨线为
第三章 作业题
1、解:由题意知末端状态
将其带入J得
则
欧拉方程
得
从而得到
因为
再由无约束泛函极小值的充分条件
故所求得 可使性能泛函取极小值
2、解:
构造函数H:
本题为末端时刻固定,末端状态自由,控制受约束的最优化问题。
式中 为轨线方向和 轴之间的夹角,证明:质点在沿最短时间轨线运动时控制变量 须满足
解:由题意得:
两边同时对 求导可得:
得证
4-1解:
本题为末端时刻固定,末端状态自由,控制受约束的最优化问题。
4-4解:
本题为末端固定,控制有约束的最优化问题。
4-6解:
本题为末端时刻固定,末端状态受约束的最优化问题。
4-7解:
正则方程:
控制方程:
由:
横截条件:
程:
控制方程:
状态方程:
第二部分最优控制理论习题答案
x c1 , x c1t c2 (通解形式)
由边界条件
x(0) c2 1 c1 1 , 解之,得 x(1) c1 c2 2 c2 1
*
故最优轨线为 x (t ) t 1
2-2、求一阶系统 x(t ) u(t ), x(0) 1,当性能指标为 J
构造哈密顿函数: H x1 u12 u2 2 1 u1 2 ( x1 u2 ) 最优轨线 x1* (t ) t , x2* (t ) 0.5t 2 0.5t 最优控制 u1* 1 , u2* 0.5
2-8、 设二阶系统状态方程为 x1 x1 u,
x2 (1) ,
1
故
H H 1 2 , 2 0, x1 x2
1 1 2 11 c1et c2 Nhomakorabea c2
由横截条件
1 (1)
0, 2 (1) 1 x1 (1) x2 (1)
那么 所以
代入边界条件
x(0) 1 (c1 c2 1)
,
(1) 0
(c1e c2e 0)
1
, [终端横截条件 t f
] x(t f )
得 c1 0.12, c2 0.88
最优轨线 x* (t ) 0.12et 0.88et 最优控制 u* (t ) 0.12et 0.88et
-1
最优轨线是齐次方程 x x 的解
x(t ) [ A BR1BT K ]x(t )
由 x(0) 1 ,解得: x (t ) e
* t
所以: u (t ) e
最优控制与状态估计2
tf
, t ) 及 x (t ) 在 [t0 , t f ] 上连续可微, t0 和 t f 给定, 其中, L( x , x
x (t f ) x f ,x (t ) R n ,则极值轨线 x * (t ) 满足如下欧 已知 x (t0 ) x0 , 拉方程
L d L 0 x d t x
J [ x (t )] J [ x0 (t )] 就称泛函J [ x (t )]在
x (t ) x0 (t ) 处是连续的。
2、泛函的变分 所谓泛函 J [ x (t )]的宗量 x (t ) 的变分是指两个函数间的差。
δ x x (t ) x0 (t )
n
x (t ), x0 (t ) R n
x (t ) x * (t ) x (t )
于是泛函J 的增量J
J
(t ) x * (t ) x (t ) x
可计算如下(以下将*号省去)
tf t0
tf
x , t L x, x , t dt L x x , x
华东理工大学
ECUST
泛函的变分等于
J x(t ) x 0
3、泛函变分的规则 1) 2) 3) 4)
δ( L1 L2 ) δ L1 δ L2
δ( L1 L2 ) L1 δ L2 L2 δ L1
δ
b a
, t]d t L[ x , x
b a
f ( x, u, t ) x
初始状态
(6) (7)
x (t )
t t 0
x (t0 )
其中,x 为n 维状态向量; u 为r 维控制向量; f 为n 维向量函数。 要求在控制空间中寻求一个最优控制向量 u(t ),使以下性能指标
《控制理论CAI教程》答案
(s)
=
1 C2s
I
(s)
Uo (s) = R2I (s) + Uc2 (s)
1 Uc2(s) Uo(s)
C2 s
R2
第二章 控制系统的数学模型
Uc1(s) R1
R1C1s +1
-
1 I(s)
-
Ls + R2
1 Uc2(s)
C2 s
R2
Ui(s) 1
− R1 R1C1s +1 1 Ls + R2
R2
第三章 控制系统的时域分析法
3-3 s(τs+1)(2s+1)+K(s+1)=0
2τs3+(2+τ)s2+(1+K)s+K=0
s3 2τ
1+K
s2 2+τ K
s1 (2+τ)(1+K) -2τK
s0 K
2τ>0
K>0 (2+τ)(1+K) -2τK>0
第三章 控制系统的时域分析法
τ
τ>0
(1-K) τ+2(K+1)>0
特征多项式:s3 +as2 +bs+d =0
设另一极点为e(e<5),则特征多项式:(s+e)(s2+2s+2)
第三章 控制系统的时域分析法
二阶系统,R(s)=1/s,σp=0.2,ts=1.8s 时,试确定K1、τ值。 当输入信号分别为: r(t)=1(t),r(t)=t,r(t)=1/2t2时,
Uo (s) = P1Δ1 + P2Δ2
Ui (s)
最优控制试题答案(在职研究生班)2
广西工学院在职研究生班课程《最优控制》参考答案一、简答题1、系统数学模型、边界条件与目标集、容许控制、性能指标。
2、积分型性能指标,末值型性能指标,综合型性能指标3、控制向量不受约束,且是时间的连续函数。
4、控制向量受到约束,哈密顿函数对控制向量的偏导不存在时。
5、状态调节器问题;输出调节器问题;跟踪问题。
6、不论初始状态和初始决策如何,当把其中的任何一级和状态再作为初始级和初始状态时,其余的决策对此必定也是一个最优控制。
二、计算题(70分)1、解 本题 t f 固定,末态自由。
由题意 ∙+=21x L欧拉方程2=-=∂∂-∂∂∙∙∙x L dtd xL x解得 ()21c t c t x += 由边界条件及横截条件021==∂∂∙=∙x xLf t解得 c1=0 ,c2=0 故所求极值曲线为 ()0=*t x2、解 本题是求解最短曲线问题,可以将性能指标设定为曲线长度函数的积分,当该指标为最小时,所得的曲线即为最短曲线。
根据几何知识,在直角坐标系中弧线元的长度表示为dtdx dt ds x21)()(22∙+=+=设性能指标为 dt J tftox⎰∙+=21由题意可知,tf 固定,末态固定,21xL ∙+=,由欧拉方程0=∂∂-∂∂∙xL dtd xL ,22c x=∙(常量)解得 x(t)=ct+d根据边界条件,可得c=1,d=0,故所求曲线为:()t t x =*3、解 本题为定常系统,tf 固定,末端自由,末值型指标,控制受约束的最后控制问题,可采用极小值原理求解。
由题意知,性能指标为末值型的,即 [])1()1(2)(22x x tf x +=ϕ 令哈密顿函数 H=1211)(x u x λλ++-协态方程022=∂∂-=∙x H λ,2c =λ2112111,c e c x H t+=-=∂∂-=∙λλλλ,横截条件()()1,1211=+=-t e t t λλ 求出 c1=e 1-t +1,c2=1,则有()()1,1211=+=-t e t t λλ极值条件 u ⎩⎨⎧<>-==*0,10,1)sgn()(111λλλt因为()111+=-t e t λ>0,t []1,0∈,故可确定 10,1)(<≤-=*t t u4、解 根据性能指标的形式,可知本题是线性二次型问题,且是有限时间状态调节器问题。
最优控制理论与系统胡寿松版课后习题答案
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 Lx 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 xx ψ∂⎛⎫+-= ⎪∂⎝⎭易得到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 =⎧⎪=-⎨⎪=⎩ 还有一组解⎪⎩⎪⎨⎧===12121c c t f (舍去,不符合题意f t >1)将f t ,1c ,2c 代入J 可得3140)3(4)212(5025.2*=-=+=⎰⎰•t dt x x J . 极值轨线为()*269x t t t =-+2-7 设性能泛函为120(1)J x dt =+⎰求在边界条件(0)0x =,(1)x 自由情况下,使性能泛函取极值的极值轨线*()x t 。
最优控制课后习题答案
最优控制课后习题答案最优控制课后习题答案最优控制是现代控制理论中的重要分支,它研究如何在给定约束条件下,使系统的性能指标达到最优。
在最优控制的学习过程中,课后习题是巩固理论知识、培养解决问题能力的重要环节。
本文将为大家提供一些最优控制课后习题的答案,希望能对大家的学习有所帮助。
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))$为性能指标函数。
最优控制(2)
则满足末态要求的最优轨线方程可表示为
取u*= -1,也可得到满足末态要求的最优轨线方程 曲线 , 组成曲线 ,称为开关曲线,表示为
开关曲线将相平面分成两部分R+和R-
则时间最优控制为
4.2.4 最小能量控制
设线性定常系统
求满足下列不等式约束的容许控制:
使系统从初始状态x0转移到x(tf)=xf,并使性能指标
由横截条件 解出
由极小值条件
由于
可得到
t=1时,u*(t)应该为 零,即不存在最优 控制
定理 对于如下时变系统、末值型性能指标、末端自由、 控制受约束的最优控制问题
式中末端时刻固定或自由,假设同前,则对于最优解 u*,x*,tf*,必存在非零的 (t ) ,使如下必要条件成立: 1) 正则方程
则对于最优解u*,x*,tf*,必存在非零的 (t ),使如下必要条 件成立: 1) 正则方程 其中
2) 边界条件与横截条件
3) 极小值条件
4) 沿最优轨线哈密尔顿函数变化率(tf自由时用)
例子:
解:已知
由协态方程 可得到
2 (t ) c2 , 1 (t ) c1e c2
t
其中
2) 边界条件与横截条件
3) 极小值条件
4) 沿最优轨线哈密尔顿函数变化率(tf自由时用)
于是该问题就变成了如下定常问题:
(16)
利用定常系统的结论,可知协态方程为
即 (17)
横截条件为
即 极小值条件为 (18)
将式(16)代入可得
即得结论3)。沿最优轨线哈密尔顿函数变化率
将(18)代入可得到本定理的结论4)。
极小,其中 tf 固定。
构造
定义开关向量函数
最优控制习题及参考答案
最优控制习题及参考答案6212最优控制习题及参考答案习题 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 )63x f由 x (0) = 0 得: C 2= 0∂L由 ∂xt =tf= 2x (t f ) = 2C 1 t =t = 0t0 1于是: x *(t ) = 0【分析讨论】对于任意的 x (0) = x,x (1)自由。
6421∫ ⎩λ =有: C = x , C = 0 ,即: x *(t ) = x其几何意义: x (1) 自由意味着终点在虚线上任意点。
习题 3已知系统的状态方程为:x 1 (t ) = x 2 (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 )dt 2取极小值的最优控制 u *(t ) 以及最优轨线 x *(t ) 。
⎡ x ⎤解: 由已知条件知: f = ⎢ 2⎥⎢⎣ u ⎥⎦Hamiton 函数: H = L + λTfH = 1u 2+ λ x + λ u⎧λ = 0由协态方程: ⎨12121 2 2⎧λ = C① 得: ⎨1 1⎩λ2 = −C 1t + C2 ② ∂H由控制方程:∂u= u + λ2 = 0 得: u = −λ2= C 1t − C 2③由状态方程:x2 = u = C1t −C2得:x (t) = 1 C t2 −C t + C ④2 2由状态方程:x1 = x21 2 3得:x (t) = 1 C t3 −1 C t 2 + C t + C ⑤1 6 12 23 465661⎪⎩=− ∫⎡1⎤ ⎡0⎤将x (0) = ⎢ ⎢,x (3) = ⎢0⎢代入④,⑤, ⎣1⎦⎣ ⎦ 10联立解得: C 1 =由③、④、⑤式得:u * (t ) = 10t − 29,C 2 = 2 , C3=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. 求使得2()ln f x x x =-最小的x 值。
解:'1()20f x x x =-=求得可能的极值点是x = "21()2f x x =--恒小于0. 所以使得2()ln f x x x =-最小的x2. 求使221122()10124f X x x x x =---为极值的极值点x 。
解:12'12'12201201280x x f x x f x x =--==--=由上述两个方程得出的可能极值点为[]***12,0,0T T X x x ⎡⎤==⎣⎦二阶导数矩阵为*"20,1212,8X f --⎡⎤=⎢⎥--⎣⎦用塞尔维斯特判据来检验,有200-<, 20,12det 16012,8--⎡⎤=>⎢⎥--⎣⎦故*"X f 为负定,在[]*0,0T X =处,()f X 为极大。
3求.使222123121323()55484f X x x x x x x x x x =+++--为极值的极值点x 。
解:123'123'213'31210480244010840x x x f x x x f x x x f x x x =+-==+-==--=由上述三个方程得出的可能极值点为 []****123,,0,0,0T TX x x x ⎡⎤==⎣⎦ 二阶导数矩阵为*"10,4,84,2,48,4,10X f -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦用塞尔维斯特判据来检验,有100> 10,4det 04,2⎡⎤>⎢⎥⎣⎦10,4,8det 4,2,4808,4,10-⎡⎤⎢⎥-=>⎢⎥⎢⎥--⎣⎦故*"X f 为正定,在[]*0,0,0TX =处,()f X 为极小。
4.求使2212min ()45,f X x x =+且12()2360g X x x =+-=。
解:作拉格朗日函数12221212(,)()()45(236)L x x f X g X x x x x λλ=+=+++-极值的必要条件 11221282010302360L x x L x x L x x λλλ∂=+=∂∂=+=∂∂=+-=∂联立求解上面三个方程可得307λ=- 可能的极值点坐标为11514x =,297x = 根据问题的性质可以判断极小值存在且是唯一的。
最优控制理论_第二章
为正定矩阵 为正定矩阵。
二:有约束条件的函数极值问题 设 元函数f(x1,x2),x1和x2必须满足下列方程: 设二元函数 必须满足 列方程 g(x1,x2)=0 为求函数f(x1,x2)的极值,并找出其极值点(x1*,x2*),作一辅助函 数 拉格朗日函数: 数-拉格朗日函数:
L( x1 , x 2 , ) f ( x1 , x 2 ) g ( x1 , x 2 )
0
t0
不同函数F的欧拉方程为:
F [ x (t ), t ]
F 0 x
2F 0 x 2 x
(t )] F[ x
(t ), t ] F[ x
2F 2F 0 x 2 xt x
2F 2 F F x 0 2 x t x x
* * ( x1 , x2 )
正定
其中
f
* x1 x1
2 f ( x1 , x2 )
2 x1
f
f
* x1 x 2
* x2 x2
2 f ( x1 , x 2 ) * * x1x 2 ( x1 , x2 )
2 f ( x1 , x 2 )
2 x 2 * * ( x1 , x2 )
X(t)
(t )
tf t0
0
X1(t) X2(t)
故边界条件为: x(t0)=x0, x(tf)=xf
t0
X3(t) tf t
(2)自由始端和自由终端
X(t)
F t0 x
0
F tf x
0
t0 tf t
(3)自由始端和固定终端
F t0 x
X(t)
0
x(tf)=xf
最优控制理论 第二章
x(t ) f ( x(t ), u (t ), t ) 1
x(t 0 ) x0 2
t0 , t f
终端状态满足:
目标函数:
已知。
N [ x(t f ), t f ] 0 3
J [ x(t ), u (t ), t ]dt 4
J
tf t0
[
x(t )
x(t ) x x
*
x (t )
x (t )
x x*
]dt
泛函取极值的必要条件:
J 0
即:
tf
t0
[ x(t ) x(t )]dt 0 8 x(t ) x(t )
tf
t0
X
x(t )
tf t0
0 横截条件,又称为边界条件
3.横截条件的分析 <1> x(t0 ), x(t f )都固定,图a
x(t0 )
x(t )
许多状态轨线
求出最优
即 x* (t ) x(t ) 0 0
x (t f ) x(t f )
*
即 x(t 0 ) 0 x(t f ) 0 <2> x(t0 ) 固定, x(t f )自由 图 b
t1
F
证明见书。
< 定理2 > :若可微泛函 J [ y( x)] 在 y0 ( x) 上达到极值, 则在 y y0 ( x)上的变分等于0,即 J 0 证明较简单,见书。
变分规则:<1 > ( F1 F2 ) F1 F2
< 2> ( F1 F2 ) F1F2 F2F1
最优控制习题及参考问题详解
标准文档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 + C由 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) 自由。
2 0 1∫⎩ λ = −λ有: 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 )dt 2取极小值的最优控制 u *(t ) 以及最优轨线 x *(t ) 。
⎡ x ⎤解:由已知条件知: f = ⎢ 2⎥⎢⎣ u ⎥⎦Hamiton 函数: H = L + λT f H = 1u 2 + λ x + λ u⎧λ = 0由协态方程: ⎨ 12 121 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 ) = 1C 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 +127 x *(t ) = 5 t 2 − 2t +1 29习题 4 已知系统状态方程及初始条件为x =u , x (0) = 1试确定最优控制使下列性能指标取极小值。
第二篇最优控制理论习题答案
第二篇最优控制理论习题答案:2-1、求通过x(0)=1,x(1)=2,并使性能指标120(1)J xdt =+∫&为最小的曲线x(t)。
解:本题属于无约束(无状态方程约束),始端和终端均固定的泛函极值问题,可用变分法求解。
被积函数 21,0,2,2L L d LL xx x x x dt x∂∂∂=+==⋅=∂∂∂&&&&&& 代入欧拉方程 0L d Lx dt x ∂∂−⋅=∂∂&, 得20x =&&, 即0x =&& 1xc =&, 12x c t c =+ (通解形式) 由边界条件 212(0)1(1)2x c x c c ==⎧⎨=+=⎩, 解之,得1211c c =⎧⎨=⎩ 故最优轨线为 *()1x t t =+2-2、求一阶系统()(),(0)1xt u t x ==&,当性能指标为12201()2J x u dt =+∫取最小值时的最优控制与最优轨线。
解:本题属于有约束,始端固定;终端时间f t 固定,()f x t 自由,控制u 无限制的泛函极值问题,可用变分法求解。
构造哈密顿函数 222211() ()22H x u u L x u λ=++=+注: 协态方程 H x x λ∂=−=−∂&, 即x λ=−& ① 极值条件/控制方程 0Hu uλ∂=+=∂, 即u λ=− ②由系统的状态方程 xu =&及②式,,x x λλ=−=−&&&& ③ 由①式及③式,得 xx =&& 故12()ttx t c e c e −=+ 12()()t t t xt c e c e λ−=−=−+& 代入边界条件11212(0)1(1)0,(1)(0)x c c c e c e λ−==+=−+=, [终端横截条件()f t f x t φλ∂=∂] 得 120.12,0.88c c ==最优轨线 *()0.120.88t tx t e e −=+ 最优控制 *()0.120.88t tu t e e −=−2-5、有一开环系统,包含放大倍数为4的放大器和一个积分环节。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
所以: u*(t) = −e−t
代入 Riccati 代数方程
−KA − AT K + KBR−1BT K − Q = 0 −K ⋅0 − 0⋅ K + K ⋅1⋅1⋅1⋅ K −1 = 0
u (t )
1
x (t )
s
‐1
故 K = ±1 取 K = 1 最优控制 u*(t) = −R−1BT Kx(t) = −x(t) 最优轨线是齐次方程 x& = −x 的解
故最优轨线为 x*(t) = t +1
∫ 2-2、求一阶系统 x&(t) = u(t), x(0) = 1 ,当性能指标为 J = 1 1 (x2 + u2 )dt 取最小值时的最优控
20
制与最优轨线。
解:本题属于有约束,始端固定;终端时间 t f 固定,x(t f ) 自由,控制 u 无限制的泛函极值问题,
可用变分法求解。
构造哈密顿函数 H = 1 (x2 + u2 ) + λu 2
注:L = 1 (x2 + u2 ) 2
协态方程 λ& = − ∂H = −x , 即 x = −λ& ① ∂x
极值条件/控制方程 ∂H = u + λ = 0 , 即 u = −λ ② ∂u
由系统的状态方程 x& = u 及②式, x& = −λ, &x& = −λ& ③
故最优控制 u * (t) = −1,∀t ∈[0,1)
∞
∫ 2-9、设线性系统为 x&(t) = u(t), x(0) = 1性能指标为 J = (x2 + u2 )dt 试求最优控制 u(t) ,使性 0
能指标 J 取最小值。
解:本题属于线性定常系统状态调节器问题,
由系统状态方程及性能指标 A = 0, B = 1,Q = 1, R = 1
2x&,
d dt
⋅
∂L ∂x&
=
2&x&
代入欧拉方程
∂L ∂x
−d dt
⋅
∂L ∂x&
=0,
得 2&x& = 0 ,
即 &x& = 0
x& = c1 , x = c1t + c2 (通解形式)
由边界条件
⎧ ⎨ ⎩
x(0) = x(1) = c1
c2 = + c2
1 =
2
,
解之,得
⎧⎨⎩ cc12
=1 =1
− xT e2T
+ x0e2T − e−2T
x* (t )
=
xT − x0e−2T e2T − e−2T
⋅ e2t
+
−xT + x0e2T e2T − e−2T
⋅ e−2t
u* (t )
=
1 4
x&* (t )
=
xT − 2(e2T
x0e−2T − e−2T
)
⋅ e2t
+
xT − x0e2T 2(e2T − e−2T
达到最小值,试求输入的控制 u(t)。
解:本题属于有约束,始端和终端均固定的情况。
由题意: x& = 4u x(0) = x0 x(T ) = xT
u
4
x
s
构造哈密顿函数: H = x2 + 4u2 + λ ⋅ 4u
λ& = − ∂H = −2x 即 ∂x
∂H = 8u + 4λ = 0 即 ∂u
λ& = −2x λ = −2u
(1)
=
∂φ ∂x2 (1)
=
1
那么 c1e + c2 = 0 , c2 = 1, ∴c1 = −e−1
所以 λ1(t) = 1 − et−1
由极值条件
H*
=
min
u∈Ω
H
得
u
*
(t)
=
− sgn {λ1 (t )}
不难发现,
λ1(0) = 1 − e−1 > 0 , λ1(1) = 0
即 λ1(t) = 1 − et−1 > 0,t ∈[0,1)
由①式及③式,得 &x& = x 故 x(t) = c1et + c2e−t
λ(t) = −x&(t) = −c1et + c2e−t
代入边界条件
x(0) = 1 ,
(c1 + c2 = 1)
λ(1) = 0 (−c1e + c2e−1 = 0) ,φ ] ∂x(t f )
原理求解。
令 H = λ1(−x1 + u) + λ2x1,
φ = x2 (1) ,
λ&1
=
− ∂H ∂x1
= λ1
− λ2,
λ&2
= − ∂H ∂x2
=0,
故 λ&&1 = λ&1 − λ&2 = λ&1
λ1 = c1et + c2
λ2 = c2
由横截条件
λ1(1)
=
∂φ ∂x1(1)
=
0, λ2
第二篇最优控制理论习题答案:
1
∫ 2-1、求通过 x(0)=1,x(1)=2,并使性能指标 J = (x&2 +1)dt 为最小的曲线 x(t)。 0
解:本题属于无约束(无状态方程约束),始端和终端均固定的泛函极值问题,可用变分法求解。
被积函数
L = x&2 +1,
∂L = 0, ∂x
∂L ∂x&
=
由 x& = 4u 及(1)(2)两式,得: &x& = 4x
(1) (2)
x(t) = c1e2t + c2e−2t
代入边界条件
⎧
⎨ ⎩
x(T
x(0) = c1 + c2 = x0 ) = c1e2T + c2e−2T =
xT
解方程组得: c1
=
xT − x0e−2T e2T − e−2T
,
c2
=
)
⋅ e−2t
2-8、设二阶系统状态方程为 x&1 = −x1 + u, x&2 = x1 , x1(0)=1,x2(0)=0, u ≤ 1,终端 x(tf)自由。
试确定最优控制 u(t) ,使下列性能指标 J= x2(1) 取最小值。
解:本题为控制受限制; t f 给定, x(t f ) 自由;末值型性能指标的最优控制问题,可用最小值
得 c1 = 0.12, c2 = 0.88
最优轨线 x*(t) = 0.12et + 0.88e−t 最优控制 u*(t) = 0.12et − 0.88e−t
2-5、有一开环系统,包含放大倍数为 4 的放大器和一个积分环节。现加入输入 u(t),
T
∫ 要将系统从 t=0 时的 x0 转移 到 t=T 时的 xT,并使性能泛函 J = (x2 + 4u2 )dt 0