第三章 极小值原理及应用 经典变分法缺陷: 1、应用前提:a 、控制量 u(t)的取值不受任何限制,没有任何 不等式约束。 b 、 f、L、

Slides:



Advertisements
Similar presentations
旋毛虫病  旋毛虫为毛首目毛形科的线虫,是一种人 畜共患病。幼虫寄生于肌肉中称肌旋毛虫, 成虫寄生于小肠称肠旋毛虫。它是多宿主 寄生虫。除猪、人以外,鼠类、狗、猫、 熊、狼等均可感染,目前已有 65 种哺乳动 物可感染此病。
Advertisements

稳恒磁场习题课. 类比总结 1. 产生 静止电荷运动电荷 2. 被作用 电荷 与电荷运动 状态无关 只对运动电荷作用 3. 表观 性质 力力  作功 力力 4. 基本物 理量 5. 基本 定理 基本 性质 表一 场的产生与性质静电场稳恒磁场.
1 主題三 網路常見衝突事件 的防範 3-1 認識網路兩性交往常見的衝突事件 3-2 瞭解處理兩性網路交往衝突之注意事項 3-3 認識處理兩性網路交往常見的衝突事件的 有效方法 有效方法.
XX啤酒营销及广告策略.
德 国 鼓 励 生 育 的 宣 传 画.
泛黄的春联还残留在墙上 依稀可见几个字岁岁平安 在我没回去过的老家米缸 爷爷用楷书写一个满 黄金葛爬满了雕花的门窗 夕阳斜斜映在斑驳的砖墙 铺着榉木板的屋内还弥漫 姥姥当年酿的豆瓣酱 我对着黑白照片开始想像 爸和妈当年的模样 说着一口吴侬软语的姑娘缓缓走过外滩 消失的旧时光一九四三 在回忆的路上时间变好慢.
中考热点 空气中氧气含量的测定习题课.
2013届高考复习方案(第一轮) 专题课件.
专题19 自然灾害与防治.
第2节《动量守恒定律》 张映平.
第三組做的報告 . 組長:紀美朱 組員:謝寶岳 羅芳婷 李依芳 李銘賢 還有內容.
未成年少女墮胎的法律問題.
人體的奧妙 姓名:林佩儀 班級:六年4班 指導老師:趙美華老師.
第十六章 狭义相对论基础 经典力学对于解决宏观物体的低速运动是卓有成效的,我们在日常生产、生活和科学实验中所遇到的物体的运动速度,大都远比光速为小,所以经典力学的适用范围很广泛。但是自从十九世纪中期以来,电子和放射性被发现以后,人们对于物质世界的认识开始进入微观领域,分子、原子和电子等微观粒子成为我们的研究对象。一个宏观物体,例如一块石头,质量以kg计;而微观粒子,例如一个电子,质量以10-31kg计。宏观物体的运动速度以m·s-1计,而微观粒子的运动速度却可以达到接近光速(3×108.
第二章 复式记账原理*** 主要内容、重点难点: 1.会计要素与会计等式*** 2.会计科目与账户*** 3. 借贷记账法***
舌尖上的昭通.
第四章 种群和群落 第33讲 种群的特征及其数量变化.
簡報大綱 前言 為何會有異質採購最低標 異質採購最低標法令規定 各種決標方式之履約成果分析.
氧气的制法 装置 原理 练习 随堂检测.
恒泰期货研究所2016年 期债暴跌告一段落,短期波动降低 国债期货周报
文明史观 文明史观,通常被称为文明史研究范式,是研究历史的一种理论模式。人类社会发展史,从本质上说就是人类文明演进的历史。
1、分别用双手在本上写下自己的名字 2、双手交叉
南美洲 吉林省延吉一高中 韩贵新.
税收新政 -2008年度.
將軍澳循道衛理小學 申請中一自行分配學位 家長晚會
2007年11月考试相关工作安排 各考试点、培训中心和广大应考人员:
死與生的自我掌握.
欢迎大家来到生命科学课堂.
专题4 地表变化及影响.
主题一 主题二 模块小结与测评 主题三 考点一 主题四 考点二 主题五 考点三 主题六 考点四 命题热点聚焦 考点五 模块综合检测 考点六.
分式的乘除(1) 周良中学 贾文荣.
第四章 制造业企业 主要经济业务核算.
高三生物备考复习指导 稳 准 恒 温州中学 张微华 2009年3月29日.
第五章 血液生理 第一节 概述 第二节 血液的组成和理化性质 第三节 血细胞生理 第四节 血液凝固 第五节 免疫系统 第六节 血型与输血原则.
《思想品德》七年级下册 教材、教法与评价的交流 金 利 2006年1月10日.
动画分镜头技巧 梁思平.
实验设计中的因变量检测 乐清中学 霍晓珍.
食用受污染三鹿牌婴幼儿配方奶粉相关的 婴幼儿泌尿系统结石的超声诊断.
北京中日友好医院骨科 刘朝晖 程立明 郭万首 王冉东 邹海波 张铁立 韦宏宇 孙伟
我国三大自然区.
“这是一道选择题,请看题板:由于他( )成一个商人,日本鬼子没有认出他来。
专题二 识图题增分技巧.
B F C D G E B E A 下图是沿20°经线所作的地形剖面示意图
物 资 供 应 简 报 第三期 2014年3月 中铁二局物资重庆分公司项目物资简报.
我是情緒管理小高手 黃玲蘭老師.
第五讲 地球公转及其地理意义 一、地球公转的 1.公转特征 (1)轨道:近似正圆的 轨道,太阳位于其中一个 上。
北师大版七年级数学 5.5 应用一元一次方程 ——“希望工程”义演 枣庄市第三十四中学 曹馨.
麻疹的院内感染控制 敦煌市医院院感科 梁荣.
动物激素的调节及其在农业生产中的应用(B级)
海洋存亡 匹夫有责 ——让我们都来做环保小卫士 XX小学三(3)班.
深圳市晨兴餐饮投资管理有限公司 招商手册.
機車第六篇 事故預防 單元二 行駛中注意事項.
《美国的两党制》选考复习 温州第二高级中学 俞优红 2018年6月14日 1.
实验一 双绞线的制作与应用.
一、认真审题,明确作图目的。 二、作图按投影规律准确无误。 三、图线粗细分明。 四、需要保留作图线的一定保留。
力 学 第二章 杨维纮 中国科学技术大学 近代物理系.
04 第四章 應用幾何 4-1 概說 4-2 認識尺度符號 4-3 等分線段、圓弧與角 4-4 垂直線與平行線 4-5 多邊形
性騷擾之調查與防治 主講人:龜山分局 家防官 劉淑卿.
教学建议 学习目标 § 6.1 矩阵的概念 § 6.2 矩阵运算 § 6.3 矩阵的初等行变换与矩阵的秩 § 6.4 线性方程组的消元解法
圓心角 A 劣弧 優弧 C O B D 對 的圓心角 AOB 顧震宇老師 台灣數位學習科技股份有限公司.
3、运动快慢的描述—— 速度.
铰链四杆机构 5 机械2组 郭广新.
第六章 静电场 第2课时 电场力的性质.
义务教育课程标准实验教科书 八年级 上册 第三章 3 生活中的旋转.
用待定系数法求二次函数的解析式 乌鲁木齐市第六十七中学:王磊.
3.2 生长素的生理作用 孔楠 南京师范大学
第2节 大气的热力状况 基础知识回顾 重点难点诠释 经典例题赏析.
自动控制原理.
● o. ● o 涂色部分是扇形 o ● 与“扇形”相关的概念. 概念1:“弧” 概念2:“圆心角”
第十二章 机械波 4、波的反射和折射.
Presentation transcript:

第三章 极小值原理及应用 经典变分法缺陷: 1、应用前提:a 、控制量 u(t)的取值不受任何限制,没有任何 不等式约束。 b 、 f、L、 第三章 极小值原理及应用 经典变分法缺陷: 1、应用前提:a 、控制量 u(t)的取值不受任何限制,没有任何 不等式约束。 b 、 f、L、 等函数对其自变量有充分可微性。 2、实际控制要求: a 、控制量u受不等式约束,如: ,i=1,2,3…… b 、性能指标有时并不完全可微

如:燃料最优控制: 若采用经典变分:

若采用经典变分法: 不再适用,求不出解来 实际应为 极小值原理 若在容许控制范围内,J或H有极值且唯一,用极小值 原理与经典变分法,所得 结论一致。

一、<定理>极小值原理:[时变系统] 时变受控系统 为容许控制 ,其中控制向量 , 域, U(t)是在 内取值的任何分段连续函数,为使状态向量由初始 转移到末端 , 满足约束: , 未定, 并使性能指标达 到极小值。 设 和 是如上J为最小的最优解, 为最优状态轨 线,则必存在不 为0的n维向量 ,满足: 2、边界条件: 1、规范方程:

二、极小值原理的意义: 3、与 对应的哈密顿函数H取极小值。 即:设 为满足 状态方程和协状态方程的最优解。 在 中。把H仅看作U的函数,若J为最小,必要条 件为 使得 仅看作U的函数时也取最小值。 极小值原理的证明:应用数学基础较多,有些书中用很大篇幅进行 证明,省略。 二、极小值原理的意义: 1 、容许控制条件放宽 变分法:在整个控制域,对U没有约束 且即使U不受限制, 有时 计算不易。 极小值原理:H在U的约束闭集中取极小值。 变分法仅为极小值原理的一个特例。

2、最优控制 使哈密顿函数H取极小值,极小值原理由此得名。 这一原理是苏联学者 “庞特里亚金”等人首先提出,而后加以证明得。 在证明过程中: 与H得符号与这里所定义的相反。 ∴所以有的文献中也称为“极大值原理”。 3、H对u没有可微要求,因此应用拓宽。 4、 极小值原来是求取最优控制的必要条件,非充分条件。 即:满足极小值原理不一定J取极小值,需进一步判断。 一般:对于实际系统 有最优解 有唯一解 最优解

三、几种边界条件得讨论: 上面所讨论的是 和 已知。 受约束, 自由的最一般 情况。若 和末端状态不同,只需改变极小值原理的边界条件即可。 1) 已知, 边界条件为: 2) 给定, 自由, 未给定, 边界条件: 确定 3) 已知, 给定,末端受约束 边界条件为: 若 自由:外加:

四、例题分析 :设一阶系统状态方程: 试求使性能指标: 为极小值的最优控制 及最优性能指标 解:定常系统, 固定,末端自由问题 x(0)=5 控制约束: 试求使性能指标: 为极小值的最优控制 及最优性能指标 解:定常系统, 固定,末端自由问题 根据极小值原理,使H绝对极小相当于使J为极小 所以 由协状态方程:

由横截条件: 显然:当 时, 产生切换 所以

由x(0)=5代入,得 所以 令t=0.307可得0.307≤t≤1时x(t)的初始条件: 解得 所以 将 代入J可得:

例2: 求 a)对U没有约束 b) |u| 解:a)

解得: b) |u| 由极小值原理: 当t=1时 在[0,1]区间 所以

五、极小值原理中哈密顿函数H的性质讨论 用途:对于所求解的最优控制的验证,或帮助求解最优控制及 1、线性定常系统: 固定, 包括 则: 常数 。 { H中不显函t} (与末端状态无关) 自由, 沿最优控制轨线: (与末端状态无关) 因为 中不显函t所以 常数 又因为 自由,

2、对于时变系统: 固定: 证明:见胡寿松P91页 自由: ,末端 若末端自由:

第四节最小值原理在实际中的应用 几个典型例子: 1.时间最优控制问题 2.最小燃料消耗问题 3.最小能量控制问题 4.线性调节问题 介绍重点: 时间最优控制问题(其他求解思想与此类似)

一、时间最优控制问题 所谓时间最优控制,就是把系统从初始状态转移到目标状态的时间作为性能指标,即使转移时间为最短。 这也是发展得最早的最优控制问题之一。

1、问题提出(时变系统) 并设 f 和 B对X(t)和t 连续可微。 已知受控系统 u: r×1 控制向量 f :n×1 函数向量 X:n×1 状态向量 u: r×1 控制向量 f :n×1 函数向量 B:n×r 函数值矩阵 控制向量约束条件: 末端状态: g:p ×1维函数向量 目标函数: : 自由 问题:寻求最优控制u*(t),使系统由初态到终态, 目标函数J 为最小

应用最小值原理进行问题的求解 步骤: ⑴列写哈密顿函数 ⑵由控制方程求u*(t) ∵u有约束, ∴H在u*上取得极小值,即: 令 q:r ×1维向量函数 [注: ]

则有: j =1,2…r 最优控制u*(t)是使 为极小,则: +1 -1 不定 可见:当 时, 有确定值,正常情况 +1 -1 奇异 t -1 不定 可见:当 时, 有确定值,正常情况 当 时, 不定, 奇异情况

我们仅研究正常情况 u*(t)写成符号函数sgn{ }形式 则 j =1,2…r 向量形式:u*(t)=-sgn{q*(t)} ⑶根据规范方程: 及初始条件和横截条件: 可求得x*(t)及

⑷求最优控制u*(t) →砰一砰控制 2、砰一砰控制定理: 要求控制量始终为最大或最小 设u*(t)是上述问题提出的解,x*(t), 是相应的状态轨线和协状态轨线。若问题正常(非奇异),则 这是一个继电器控制方式,称为砰一砰控制

3、线性定常系统的最小时间控制问题的解法: ⑴如何确定最优控制u*(t) 设线性定常系统的状态方程为: 其中,X:n ×1维状态向量 u:控制变量 A,B分别为n ×n,n ×1矩阵 约束条件:         末端条件: 求   ,使系统状态从    转移到       所用时间最短,即使      为最小

⑵问题的求解 ①首先列写哈密顿函数: ②根据极小值原理分析可得: 注: 为标量函数,题意要求 ③有规范方程:

代入   得:   可见,  的值完全由 的符号决定   但是, 的确定是不容易的。因为它还和系统的状态变量有关系。通常采用的方法是: 先设一个 ,求出  ,求出  ,判定 若为0,则  即为所求;否则修正  重复上述过程

⑶开关次数定理: 设线性系统 是正常的(不存在非奇异问题),若矩阵A的特征值均为实数,假定时间最优控制存在,并令其为 则u*(t)的切换次数最多不超过(n-1)次,n为系统的维数。 以下将根据极小值定理,开关次数定理及相平于状态空间分析,求u* 例题分析1: 时间最优控制问题 求u*(t)

解:对象为二阶线性系统[双积分模型的时间最优控制](应用最普通最广泛的一种) 由规范方程: 则

由开关次数定理知:切换一次,设切换时间为ts,则令 C1,C2的取值要求:保证 由开关次数定理知:切换一次,设切换时间为ts,则令 为了求出ts,必须首先找出状态在 平面上的转移轨线。 ts tf

由 设u=1,则 X2 则: s t p 其中 如图(a)所示,为一组抛物线, 当K=0时经过原点[pos]

为一组抛物线,如图(b),当K1=0时过原点[NOT] 若u=-1,则 X2 N o X1 T u=-1 为一组抛物线,如图(b),当K1=0时过原点[NOT]

显然:若 初始状态在NO或在PO上,可进一步转移到目标原点,称NOP为开关曲线 由题意假设 它落在u=-1相应抛物线组中的一条上,即AQB,这时在u=-1的作用下, 状态由 沿AQB 转移到B,进行切换,B位于PO上,一步可到达原点。 X2 N u=-1 A[1,1] o X1 u=+1 B p

从初始状态到达末端状态的轨迹为AQBO, 即 u*= 进而,可求出转移时间ts及最优时间 因此,问题的解为: ①先以u=-1控制到达Po曲线上的B点 ②以u=+1沿开关曲线Po到达原点 从初始状态到达末端状态的轨迹为AQBO, 即 u*= 进而,可求出转移时间ts及最优时间 把状态轨线控制序列分成若干段,逐步算出所需时间,最后相加。 求 及ts 在AQB段,u=-1, 切换次数为1 -1,+1 到达B点:t=ts,

BO段:u=+1, 当 时, ,则 在B点应有: 联立求解: 即:

例题分析2:二阶积分系统的最小时间控制系统 最小时间控制问题:求u*(t),使系统由初态 转移到末端状态 的时间为最小,且满足 解:⑴列写哈密顿函数: ⑵求解协状态方程 设 ,则:

显然,由⑵知, 为一条直线,其形式有可能为4种 ⑶确定控制序列: 显然,由⑵知, 为一条直线,其形式有可能为4种 +1 -1 -1 u u u u 因此,u相应的控制序列为:{-1},{+1,-1}, {-1,+1} {+1}

⑷状态轨线: 由⑶知,u有4种可能的取值,其值为±1,代入状态方程: 注:

利用上式,消去中间变量t,可导出x1和x2的关系为: 如图:u=+1为实 u=-1为虚 X2 A X1 u=-1 u=+1 B

⑸确定开关曲线:使系统状态直接回到末端状态的曲线AO和BO 总的开关曲线:AOB 显然: X2 A O X1 u=-1 AOB将状态平面分为两部分 和 B 显然:

⑹确定最优控制作用u* u*与初始状态 有关 分析: ①若 位于BO上,则u*= +1; ②若 位于AO上,则u*= -1; ③若 位于 内,则u*=[ -1,+1]; ④若 位于 内,则u*=[ +1,-1]; ③④在开关曲线上为转折点

设有一升降机W,它的质量为1,升降机一方面 受重力g的作用,另一方面受控制器的作用力u(t)的作用,且 例3:升降机的快速降落问题: 设有一升降机W,它的质量为1,升降机一方面 受重力g的作用,另一方面受控制器的作用力u(t)的作用,且 (M>g是常数) 设x(t)为升降机离开地面的距离, 当t= 时, [离地面距离] [垂直运动速度 ] 问题:求u*(t),使升降机最快的到达地面,并且到达地面时的速度为零。 即: u 最小, 自由 W X(t) g

解:建立升降机系统的数学模型,F=ma 即: 令: 即: 哈密顿函数: 显然,为了使H为最小,则 即: 不确定

协状态方程: 即: 常数 相应于 的4种可能,u*的取值有4种可能 {+M},{-M},{+M,-M},{-M,+M} 因此,下面只研究u=±M时升降机的状态轨线 设u=M,则状态方程为: …① …② ①/②: 是一组抛物线, 图中实线箭头表示状态运动的方向 在此族曲线中,只有 到达原点,

设u=-M,同理可得: 如图虚线所示 只有 到达原点, 开关曲线 r将相平面分为两部分,在r下半部的记为 ,包括 在r上半部的记位 ,包括 ∵u*只取+M或-M,切换最多一次,因此可得到结论: 〈ⅰ〉初始状态 在 上, 状态沿 回原点

〈ⅱ〉当 在曲线 上时 , 状态沿 回原点 〈ⅲ〉当 时, 沿相应的虚线抛物线运动到 时, 沿 回到原点。 马上切换 〈ⅳ〉当 时, ,沿相应的实抛物线运动到 时, 马上切换 ,沿 回到原点。 总之:

它在地面之上,∴ ,处于相平面的右半部分,且设 对于实际问题升降机的分析: 它在地面之上,∴ ,处于相平面的右半部分,且设 a〉若 ,而 时状态沿实抛物线运动与 轴交于N ,这意味着升降机到达地面时,速度不为0,不合要求。 <b> 当 即开始以最大推力向下最用, 使升降机尽快下降。当其状态检测到达 时,马上改变控制,使它以 的最大推力向上作用,这样升降机将以速度0到达地面。 N

从上例可以看出:快速最优控制有如下特点: <ⅰ > u*要么最大,要么最小。 < ⅱ> u*的取值经过有限的(n-1)次(可为最多次)数切换可到达平衡点。 < ⅲ > u*的取值仅在开管曲线上切换。 注意:时间最优控制的应用中,有些实际问题并不要求将相点控制到状态空间原点,而是到某一集合,其分析方法与上类似 (若二阶系统为一般的二阶系统,特征值为实数时,分析方法类似;为复数或纯虚数时,开关次数定理不成立,问题较为复杂,如无阻尼振荡二阶系统。

二、燃料最优控制问题 节约能源,减少燃料消耗在国民经济各部门中都是一项重要的技术经济课题。在航空和宇航中使用的原料是由地面起飞时带到空间去的。在空中携带的燃料是有限的,要保证长时间的飞行计划,就希望空中的控制系统消耗的燃料最小,而燃料的消耗一般是和控制力u的大小成正比的。U有正有负。因袭燃料消耗的性能指标 : 也可以以升降机系统分析,只是 相应于时间最优控制,要求到达地于所用时间最小, 相应于燃料最优控制,要求达到目的地时所用燃料最小

1、数学描述[以二阶级分模型的燃料最优控制为例] 系统: 约束: 要求:系统从初始状态 转移到(0,0) 使 最小, 给定。 解:应用极小值原理

正常:仅在有限个点上 奇异:至少在一段时间[t1,t2]间隔内 正常:u*可取+1,-1,0随着t增大,u*在三个值上 切换,是一种三位控制{开关控制}。 奇异:不能用极小值,死区函数。 为使H为最小,则使    为最小 分析:①若u=+1,则         若使H最小,则    ②若u=-1,则    ③若

由: 和相应的最优控制 之间的关系: 显然,燃料最优控制也是 开关式控制,控制器应为 一个具有死区的继电器。 +1 -1 +1 -1 +1 和相应的最优控制   之间的关系: -1 +1 -1 +1 显然,燃料最优控制也是 开关式控制,控制器应为 一个具有死区的继电器。 -1 +1 ta tb tf -1

和 的计标 当   时,                           相平面上一组抛物线[实线] 当   时,                           相平面上一组抛物线[虚线] 当   时,                          

以下两个图形画出了不同初始状态转移轨线 仅为<1>进行分析: 在t=ta处应满足: 相对于X2而言,a点 相对于X1而言,a点 [1,1] a =0 b =-1 =1 =-1 - =1 =1 =-1 a b a b =0 =0 仅为<1>进行分析:   在t=ta处应满足: 相对于X2而言,a点 相对于X1而言,a点

在t=tb处应满足: 解方程可得ta,tb的值

习题1:设系统为 求最短时间控制 及最短时间 提示: 开关曲线: 对于AB段, 对于BO段, 切换点为B A[10,0) =1 =-1 ts 求最短时间控制   及最短时间 提示: 开关曲线: A[10,0) 对于AB段, =1 =-1 ts 对于BO段, B 切换点为B

当t=ts时 BO段:u*=+1 当 时,X1=X2=0 ,则:

在B点应有: 联立求解:

习题2分析:设线性状态方程为: 边界条件: 容许控制为: 求最短时间控制u*(t)及开关曲线(做出大致图形) 分析: 根据最小值原理: 则: 切换周期为

当u*=+1时, 是一组同心圆,圆心为(0,1) 同理,当u*=-1时,可得: 是一组同心圆,圆心为(0,-1) 只有NO右半圆及MO坐半圆弧能够到达原点,,u*的切换周期为 ,曲线如图。

当相点运动到 或 上的任意一点时,均可在相应的控制律u=+1或u=-1作用下,沿 当X2>1时, ∴ X1 ↑,X2 ↓ 当X2<1时,X1 ↓,X2 ↓ 所以箭头如图 X0 当相点运动到 或 上的任意一点时,均可在相应的控制律u=+1或u=-1作用下,沿 或 最快地到达原点。 现在改查最优轨线的倒数第二段。设u*(t)的最后一次切换发生在 上的A点,则倒数第二段的控制必有:u=-1,其最优轨线必为(0,-1)为圆心的圆弧。

由于时间持续不超过 ,故改圆弧的长度最多等于半圆,到达A’点,发上第二段转换进而进入倒数第三段 。 由于A点为 上的任一点,因此A’点形成以(-3,0)为圆心,1为半径的半圆 。显然: 是u=-1到u=+1的开关曲线,而 则为u=+1到u=-1的开关曲线。同理可取: ,一次类推,可得一系列圆弧,可谓开关曲线。

极小值原理的证明: 、基础证明: 针对定常系统 末端自由, 得出的极小值原理的结论,<定理> 二、对于时变系统 及 等情况,可通过 引入新状态变量的方法,将时变系统化为定常系统,利用定常系统极小值原理定理的结论进行证明。

极小值原理的应用(时间最优) 例: 已知无阻尼振荡二阶系统的状态方程为: 其中 试求最优控制 使系统由任意初态 以最短时间转移到状态空间原点。 解:由极小值原理,可求取最优控制的必要条件为: 正则方程:

边界条件: 极小值条件: 解协状态方程为: 所以 最优控制特点: a、 只在某些孤立时刻为0,不存在奇异段,故 为砰-砰控制 。 b、 的切换次数与系统阶数无关。 c、除首尾两端外,最优控制每隔π时间切换一次。

下面分析开关曲线: 首先考虑相平等方程: 若 则: 是一组(1,0)为圆心的同心圆。 若 则: 是一组(-1,0)为圆心的同心圆。 若 则: 是一组(1,0)为圆心的同心圆。 若 则: 是一组(-1,0)为圆心的同心圆。 方向如图: -1 1

显然,只有c=1及 两条曲线可到达末端而考虑到最优控制最优一段的时间 间隔≤π,则最优轨线最后一段必位于下列两条半圆形开关线上。 当相点运动到 上的任一点时,均可在相应的控制律U=+1或U=-1作用下,沿 很快地到达原点。 现在考查最优轨线的倒数第二段。 设 的最优一次切换发生在 的A点,则倒数第二段的控制必为: 轨迹为(-1,0)为圆心的圆弧。考虑到第二段在时间上不大于π。故设圆弧最多等于 半圆,到达 发生倒数第二段转换,进入倒数第三段。 最优控制在某曲线上进行切换的曲线称为开关曲线。

由于A点可为 上的任一点,所以 点形成(-3,0)为圆心,1为半径的半圆。 显然 : 到 的开关曲线 : 到 的开关曲线 同理:对亢于 : 到 的开关曲线 : 到 的开关曲线 同理:对亢于 可得: 依上述过程类推可得一系列圆弧: :U=+1到U=-1的开关曲线:U=+1到U=-1的开关曲线 66

这些圆弧的全体构成了所求问题的开关曲线: 开关曲线r将相平等分为两部分 所以 起点 的最优轨线

E A B C D AB弧: ,(-1,0)为圆心, 为半径的圆弧交于开关曲线B BC弧: ,(-1,0)为圆心, 为半径圆弧,交开关曲线于C CD弧: ,(+1,0)为圆心, 为半径,交开关曲线于D DE弧: ,(+1,0)为圆心, 为半径,交开关曲线于E。 EO弧: 回到终点。 所以总的控制作用: 共转换四次。

习题:已知线性定常系状态方程: 其中 , 求 使系统由任意初态 以最短时间转移到目标集: 习题:已知受控系统: ,目标集: 求满足约束条件 的时间最优控制函数,求开关曲线

注:在时间最优控制中,我们知道: 可知: 之间的关系 由前分析知: 时,可由极值条件确定 ,正常情况; 时,可为满足约束条件的任意值,为不定状态,异步情况。 但是,奇异状态并不表示时间最优控制不存在,只表明用极小值原理 不能确定最优解,需采用奇异最优控制方法,以下介绍:

<定义1>若在区间[ ]内,存在时间的可数集合: 即: 使得对所有的 均有: 则称时间最优控制是正常的。 , <定义2>若在区间 内,存在一个(或多个)子区间 使得对所有 ,有 则称时间最优 控制异步。 奇异区间。

如何判定系统是正常的,还是奇异的。 设计时间最优控制之前总希望知道问题是否有解?是否有唯一解?问题是正常 的,还是奇异的。初次之外,我们还希望了解时间最优控制的共同特点和性质。 这种一般规律的认识和了解会有助于具体系统的设计计算。 然而:对任意的非线性系统和任意的目标集,没有明确结论。 对于线性定常系统,可以回答上述问题,(目标集假设为坐标原点) 至于线性时变系统及一般性目标集问题,只有其中的部分结论适用。

<问题1>:已知线性时不变系统, 时完全能控的 求满足下列不等式或约束的r维容许控制向量U(t), 使系统 由已知初态 转移到状态空间原点的时间最短 ,根据极小值原理, <问题>最优控制的必要条件如下: 或 为B的第j列向量

从上述必要条件出发,可得一些有用的结论: <定理1>当且仅当 个矩阵 中至少有一个奇异矩阵时<问题1>是奇异的。 证明:由已知条件: 由6式知, 否则1=0错 若问题正常,则对于给定的初协态 ,可唯一确定砰-砰控制 怎样知道是正常还是奇异呢?推证定理。 假定<问题1>是奇异的,至少存在一段时间 使某 对所有 均成立:

由此: 考虑到A与 可前后交换顺序,则有: 令: 则关于n维待定向量 的代数方程组可写成: 所有 由于 为奇异矩阵,为使 ,则 必为奇异矩阵, 即: 奇异控制问题的必要条件。可以证明其为充分条件 <定理1>得证: 由设定理可进一步得出<问题>为正常得充分必要条件

<定理2>:当且仅当 全部为非 奇异矩阵,则时间最优控制是正常得。 <定理1>和<定理2>得推证过程都没有设计到目标集,因此,不论目标集如何,只 要受控系统是线性时不变得,因此两个定理可用。 将满足<定理2>得系统叫做正常系统。正常受控系统,其时间最优控制问题也是正 常得,对于正常问题,由下列存在性与唯一性定理。 <定理3>若受控系统 是正常的,且时间最优控制存在,则最优控制 必定唯一。 证明:见“百年学书”p176页。

另外,我们知道,一个完全能控的线性定常系统: 必需满足 n:系统维数 若把系统表征为: 其中 控制分量 正常问题要求 都是完全能控。 即: 说明:每一个控制分量 均能单独使受控系统由任意初态在有限时间内转 移到坐标原点。 据此,常可很容易地判断问题的时间最优控制是否属于正常情况。 显然:一个输入完全能控的线性不变系统,其时间最优控制问题也一定是正常的。 燃料最优控制的一般情况,接<之二本>

<问题>已知线性定常系统: 求最优控制 ,使系统由任意初态 转移到目标集: 且使性能指标: 为最小,T未知。 分析: 若记: 为B的第j列向量, 则H种与U(t)由关的部分R(u)为:

根据极小值原理, 应使H或R(u)取极小,则: 求出: 这就是燃料最优控制。 如何判定燃料最优控制是正常还是奇异? <定理><问题>为正常得充分条件为,对所有j=1,2,3,……r,均有 其中 <问题>为奇异得必要条件为 :对于某个或某些 有: 证明从略

注意:在燃料最优控制中,区分正常情况与时间最优控制不同。首先: 对时间最优:系统正常时,最优控制问题一定是正常的。 2.对燃料最优:即使系统正常( ),如果系统矩阵A是奇异得(A有零 特征值,即系统中含有积分环节),问题仍可能属于奇异情况。 只有当系统是正常得,且A有事非奇异矩阵,才能保证燃料最优控制有正常解。 3、另外(1)式为系统正常得充分条件,次条件不满足时, 系统仍可能有正常解(有可能正常或有可能奇异)视初始状态而定。

习题:设二阶系统 1) 试证明系统由初态: 所消耗燃料为最小得最优控制 为: 2)欲求系统由初态X(0)最快地转移到终态

2、二阶空间控制系统的状态方程为: 不等式控制约束 ,试求使系统由初态 达到平衡状态 的最短时间最优控制。 关于“二次积分模型”的燃料最优控制问题的进一步讨论: 系统: 求 ,使系统由任意初态( )转移到状态空间原点,且使性能指标: 为最小值,T自由。

解 :求解最优解的必要条件: 1)正则方程: 则: 2)边界条件: 3)极小值条件: 4)H函数变化率: 则: 仅在有限个点上为1,则正常; 在一段区间上为1,则奇异。

最优控制必为三位式控制,且至多有两次切换,候选解为:{0},{+1},{-1}, {+1,0,-1} 具体分析: 解协状态方程: 常数, 的不同,系统有可能为正常或奇异。 <1> 奇异情况: 成立,必有: 奇异。 若 ,使系H的变化规律 无法用极小值原理求解。 时, 是时间的线性函数,这时至多有两个 <2> 当 点满足 正常情况 最优控制必为三位式控制,且至多有两次切换,候选解为:{0},{+1},{-1}, {+1,0,-1} ,{-1,0,+1},{+1,0},{-1,0},{0,-1},{0,+1}由于 结尾的三种控制序列不可能为最优 控制。

因为有状态方程知:是一组不通过原点的平行线或轴上的孤立点 所以可能的最优控制序列为:六种可能:{+1},{-1},{0,-1},{0,+1}{+1,0,-1} {-1,0,+1} 为了进一步分析燃料最优控制解的性质,转向相平等分析: 当u=+1,u=-1时,有初态转移到原点的两条轨线为: 及 轴将相平等分为四部分:

当系统初始状态位于不同区域时,解大不相同。 1) 位于 上, 是唯一的燃料最优控制,且 2) 位于 时, 时燃料最优控制,且时间最短 分析: 的可能选择:{0,+1},{-1,0,+1} 分别计算两种轨线的燃料最优值,取其 最小者,即为所求。 同理: 位于 时, 为转移时间 最短的燃料最优控制

分析:设 求解积分并代入 初、末条件可得: 对系统 进而 : 说明:燃料消耗的下限为: 。因此,只要能找到 ,使系统由初态 转移到(0,0),并且所消耗燃料为 ,则此控制必为最优控制 设 下计算 两种控制所消耗燃料 是燃料最优 所以不是燃料最优

若 采用的控制方式为: 可以证明,其仍然三位燃料最优控制 只是:并不是时间也最优。 A B C D E F 仅有ABO为燃料最优,且所用时间为最短。 位于 内时, 设 若燃料最优控制存在,必有: 则在控制候选函数中,仅有u={-1,0,+1} A C D 可将系统转移到坐标原点。 但是,此时 ,不是燃料最优控制 所以无解,但存在 燃料最优

若选用 ={-1,0}则系统在u=-1作用下,将状态转移到 轴上时所消耗燃料正好 为,相应状态记为(a,0);然后取 ,系统不在消耗燃料,但状态会保持在(a,0)不动, 不会到达坐标原点。所以 ={-1,0}不满足要求 但可得到启示,进而有如下结论: 燃料最优控制解。 <命题>:若初态 ,则对任给 ,使系统 转移到(0,0),且燃料消耗量为: 的燃料最优控制为 ={-1,0,+1} 证明方法见书:对于u={-1,0,+1}分别确定各段的燃料最优值J,最后求和。 因此 足够小时, 燃料最优接近于燃料最优。 存在问题: 燃料最优存在,但由C点D点的状态转移时间为: 是很长的。

初态 时,有相同的结论。 综上所述,可得燃料最优控制为: 注意: 上述控制律同时也保证了转移时间最小。

若不考虑转移时间,则 有唯一解. 有无穷多解 无解 燃料最优控制问题虽存在最优解,但因状态转移时间 过长,从而使 系统总转移时间过长。 时间-燃料最优控制.自学,分析方法类似。

作业1:设二阶系统方程: 其中   。试确定将系统由已知初态最快地转移到     的最优控制函数    作业2:已知系统的状态方程 试求最优控制   使系统曲任意初态最快的转移到         的状态。写出开关曲线方程,并绘出开关曲线的图形。

作业3:给定二阶系统 控制约束:     试确定使系统由已知初态最快地转移到坐标原点的最优控制 作业4:设二阶系统: 控制约束     试确定使系统由已知初态最快地转移到坐标原点的最优控制和开关曲线。

作业6:系统状态方程 求时间最优控制函数,使系统由       转移到终端状态        ,并求开关曲线,绘出图形。

沿最优轨线,哈密顿函数: