优化理论和最优控制

1、分 数: _ 任课教师签字：_ 华北电力大学研究生结课作业学 年 学 期：2013-2014第二学期课 程 名 称：优化理论和最优控制 学 生 姓 名： 学 号： 提 交 时 间：2014年4月26日 优化理论和最优控制结课总结摘 要： 最优控制理论是现代控制理论的核心，控制理论的发展来源于控制对象的要求。尽50年来，科学技术的迅速发展，对许多被控对象，如宇宙飞船、导弹、卫星、和现代工业设备的生产过程等的性能提出了更高的要求，在许多情况下要求系统的某种性能指标为最优。这就要求人们对控制问题都必须从最优控制的角度去进行研究分析和设计。最优控制理论研究的主要问题是：根据已建立的被控对象的时域数学模

2、型或频域数学模型，选择一个容许的控制律，使得被控对象按预定要求运行，并使某一性能指标达到最优值1。关键字：最优控制理论，现代控制理论，时域数学模型，频域数学模型，控制率Abstract： The Optimal Control Theory is the core of the Modern Control Theory，the development of control theory comes from the requires of the controlled objects.During the 50 years， the rapid development of the scie

3、ntific technology puts more stricter requires forward to mang controlled objects，such as the spacecraft，the guide missile，the satellite，the productive process of modern industrial facilities，and so on，and requests some performance indexes that will be best in mang cases.To the control problem，it req

4、uests people to research ,analyse，and devise from the point of view of the Optimal Control Theory. There are mang major problems of the Optimal Control Theory studying,such as the building the time domains model or the frenquency domains model according to the controlled objects,controlling a contro

5、l law with admitting, making the controlled objects to work according to the scheduled requires, and making the performance index to reseach to a best optimal value. Keywords: The Optimal Control Theroy， The Modern Control Theroy， The Time Domaints Model， The Frequency domains Model，The Control Law0

6、 引言最优控制理论的形成和发展和整个现代自动控制理论的形成和发展十分不开的。在20世纪50年代初期，就有人开始发表从工程观点研究最短时间控制问题的文章，尽管其最优性的证明多半借助于几何图形，仅带有启发性质，但毕竟为发展现代控制理论提供了第一批实际模型。由于最优控制问题引人注目的严格表述形式，特别是空间技术的迫切需求，从而吸引了大批科学家的密切注意。经典变分理论只能解决一类简单的最优控制问题，因为它只对无约束或开集性约束是有效的。而实际上碰到的更多的是容许控制属于闭集的一类最优控制问题，这就要求人们去探索、求解最优控制问题的新途径。下面介绍本课程的主要内容，线性规划：单纯形法和对偶规划；非线性规

7、划：共轭梯度法、最速下降法和牛顿法，还有最优控制问题。1 优化理论的数学模型1.1 基本数学概述线性规划的标准形式：方程解的情况: 有规划数学的基本知识可以知道：二维线性规划问题若有最优解，则最优解一定可在可行域的某个顶点上达到。1.2 一维搜索1.2.1 进退法进退法特点&适用条件：-可以用相同的试点数计算出最精确的解的估计区间.-所用函数为下单峰函数基本算法：-确定试点个数-根据相对精度，得出Fibonacci数Fn-使n是满足的最小数。 -对于初始区间令计算函数值，比较其大小 若，则令，并令否则，令，并令如第3步继续迭代，通式为令，其中是充分小的数在两点中以函数值较小的为近似极小点，相应

8、的函数值为近似极小值，并得最终区间或1.2.2 黄金分割法黄金分割法实际上是试探法的一种，它根据单峰函数构造。设F(x)是搜索区间a，b上的单峰函数。为了进行一维搜索，求一维目标函数的极小点，我们可以采用试探方法来进行。为了逐步缩小单峰区间在区间内任取两点，算函数值为和)，比较这两个函数值的大小，将出现以下三种情况。(1)当时，由于函数单峰性极小点必于区间内，这时可丢掉，把搜索区间缩小为。(2)当时，同理极小点必在区间内把搜索区间缩小为 。(3)当时，这时极小点应在区间内，缩小了区间。若计算出搜索区间内两个点函数值，能把搜索区间缩短，这样不断的重复，就可越来越精确的估出区间的位置，这就是试探法

9、的基本思想。若第一次选取的试点为，则下一步保留区间为a,x2或x1,b，两者的机会是均等的，因此选取试点时希望x2-a=b-x1，实际计算取近似值：黄金分割法是Fibonacci法的极限形式。每次缩小区间的比例是一致的，每次将区间长度缩小到原来的0.618倍。2 线性规划2.1 单纯形法线性规划问题的可行域是 n维向量空间Rn中的多面凸集，其最优值如果存在必在该凸集的某顶点处达到。顶点所对应的可行解称为基本可行解。单纯形法的基本思想是：先找出一个基本可行解，对它进行鉴别，看是否是最优解；若不是，则按照一定法则转换到另一改进的基本可行解，再鉴别；若仍不是，则再转换，按此重复进行。因基本可行解的个

10、数有限，故经有限次转换必能得出问题的最优解。如果问题无最优解也可用此法判别。根据单纯形法的原理，在线性规划问题中，决策变量（控制变量）x1，x2，x n的值称为一个解,满足所有的约束条件的解称为可行解。使目标函数达到最大值（或最小值）的可行解称为最优解。这样，一个或多个最优解能在整个由约束条件所确定的可行区域内使目标函数达到最大值（或最小值）。求解线性规划问题的目的就是要找出最优解。最优解可能出现下列情况之一：存在着一个最优解；存在着无穷多个最优解；不存在最优解，这只在三种情况下发生，即没有可行解或各项约束条件不阻止目标函数的值无限增大（或向负的方向无限增大）。单纯形法的一般解题步骤可归纳如下

11、：把线性规划问题的约束方程组表达成典范型方程组，找出基本可行解作为初始基本可行解。若基本可行解不存在，即约束条件有矛盾，则问题无解。若基本可行解存在，从初始基本可行解作为起点，根据最优性条件和可行性条件，引入非基变量取代某一基变量，找出目标函数值更优的另一基本可行解。按步骤3进行迭代,直到对应检验数满足最优性条件（这时目标函数值不能再改善），即得到问题的最优解。若迭代过程中发现问题的目标函数值无界，则终止迭代。用单纯形法求解线性规划问题所需的迭代次数主要取决于约束条件的个数。现在一般的线性规划问题都是应用单纯形法标准软件在计算机上求解，对于具有106个决策变量和104个约束条件的线性规划问题已

12、能在计算机上解得。2.2 对偶规划原始规划与对偶规划是同一组数据参数，只是位置有所不同，所描述的问题实际上是同一个问题从另一种角度去描述。推论若是原始线性规划的可行解，是对偶线性规划的可行解，则与分别是原始线性规划问题与对偶线性规划问题的最优解。对偶的线性规划都有最优解的充要条件是两者都有可行解。若原始线性规划问题与对偶线性规划问题之一具有无界的目标函数值，则另一个无可行解。若原始线性规划问题与对偶线性规划问题之一有最优解,则另一个也有最优解,并且它们目标函数的最优值相等.3 非线性规划3.1 最速下降法最速下降法又称为梯度法，是1847年由著名数学家Cauchy给出的，它是解析法中最古老的一

13、种，其他解析方法或是它的变形，或是受它的启发而得到的，因此它是最优化方法的基础。作为一种基本的算法，他在最优化方法中占有重要地位。其优点是工作量少，存储变量较少，初始点要求不高；缺点是收敛慢，效率不高，有时达不到最优解。最速下降法迭代公式是迭代步骤如下：（1）给定初点，允许误差0，令k=0。（2）计算搜索方向（3）若，则 ，停止；否则令 ，由一维搜索步长，使得（4）令 ，k=k+1，转步骤（2）。3.2 共轭梯度法共轭梯度法（Conjugate Gradient）是介于最速下降法与牛顿法之间的一个方法，它仅需利用一阶导数信息，但克服了最速下降法收敛慢的缺点，又避免了牛顿法需要存储和计算Hess

14、e矩阵并求逆的缺点，共轭梯度法不仅是解决大型线性方程组最有用的方法之一，也是解大型非线性最优化最有效的算法之一。 在各种优化算法中，共轭梯度法是非常重要的一种。其优点是所需存储量小，具有步收敛性，稳定性高，而且不需要任何外来参数。二次函数的共轭方向法的迭代步骤：已知具有正定矩阵G的二次目标函数 和终止限 。（1）给定初始点下降方向 ，置k=0。（2）作精确一维搜索 ，求步长。（3）令。（4）若 ，则，停；否则，转步骤（5）。（5）取共轭方向使得 =0，i=0，1，k（6）令k=k+1，转步骤（2）。3.3 牛顿法（1）基本思想：用二次函数逼近目标函数，用二次函数的极小值点逼近目标函数的极小值点

15、。（）计算方法将在点展成二阶泰勒级数，即 令，即 若正定，由上式解出，并把它记作得 以此作为迭代公式就是牛顿法。（3）广义牛顿法牛顿法中： 此方法对二次严格凸函数是非常有效的，迭代一步即可求出最优解。一般不能保证点列收敛。广义牛顿法基本思想：，按最佳步长确定和，即 3.4 单纯形法1维单纯形：不在同一超平面上的个点生成的凸多面体。维、维、维单纯形例子。2基本思想：比较目标函数在单纯形的个顶点处的函数值，去掉其中最差点，代之以新点构成新的单纯形。重复上述过程，使单纯形逐步逼近于极小值点。（1）反射设为单纯形的个顶点，记 求反射点， 其中 是去掉后所有顶点的形心。（2）扩张1若，2若，以代替构成新

16、的单纯形；否则，用代替构成新的单纯形，并返回（1）。3若，则以代替构成新的单纯形，并返回（1）。（）收缩1如果，令：，然后压缩求点，2若，将点压缩在与之间，仍上式确定。3若，则以代替得新的单纯形；否则，令 得新的单纯形，返回（1）。如此继续计算，直至满足某个收敛指标为止。3.5 DFPDFP校正是第一个拟牛顿校正是1959年由Davidon提出的后经Fletcher和Powell改进故名之为DFP算法它也是求解无约束优化问题最有效的算法之一.DFP算法基本原理考虑如下形式的校正公式 (7)其中是特定维向量，是待定常数.这时，校正矩阵是.现在来确定.根据拟Newton条件，必须满足（6），于是有

17、或.满足这个方程的待定向量和有无穷多种取法，下面是其中的一种：，注意到和都是数量，不妨取，同时定出，.将这两式代回（5.32）得. （8）这就是DFP校正公式.3.6 罚函数法罚函数法是利用原问题的目标函数和约束条件构造新的目标函数-罚函数, 把约束最优化问题转化为相应的罚函数的无约束最优化问题来求解。罚函数分为内罚函数法、外罚函数法、广义乘子法法。罚函数根据约束条件的不同构造的辅助函数也不相同。不等式约束问题的辅助函数与等式约束的辅助函数情形不同，但构造辅助函数的基本思想是一致的，这就是：在可行点辅助函数等于原来的目标函数值，在不可行点，辅助函数值等于原来的目标函数值加上一个很大的正数。无约

18、束优化问题的最优解趋于一个极限点，这个极限点正是原来的约束问题的最优解。此外，无约束问题的最优解往往不满足原来问题的约束条件，它是从可行域外部趋向原问题的最优点。因此 也称为外罚函数，相应的最优化方法称为外点法或外罚函数法。内点法在迭代总是从内点出发，并保持在可行域内部进行搜索。因此，这种方法适用于不等式约束的问题。4 最优控制最优控制，就是将通常的最优控制问题抽象成一个数学问题，并且用数学语言严格的表示出来，最优控制可分为静态最有和动态最有两类。静态最优是指在稳定情况下实现最优，它反映系统达到稳态后的静态关系。系统中的各变量不随时间变化，而只表示对象在稳定情况下各参数之间的关系，其特性用代数方程来描述。大多数的生产过程受控对象可以用静态最优控制来处理，并且具有足够的精度。静态最有一般可用一个目标函数J=f（x）和若干个等式约束条件或不等式约束条件来描述。要求在满足