《最新大连海事大学-刘巍-高等运筹第三讲PPT课件.ppt》由会员分享，可在线阅读，更多相关《最新大连海事大学-刘巍-高等运筹第三讲PPT课件.ppt（91页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、大连海事大学大连海事大学-刘巍刘巍-高等运筹第高等运筹第三讲三讲第二篇 运筹学中的数学规划第三章第三章 线性规划线性规划第四章第四章 非线性规划非线性规划第五章第五章 锥规划锥规划第六章第六章 矩阵规划矩阵规划第七章第七章 变分不等式与互补问题变分不等式与互补问题第八章第八章 整数规划整数规划第九章第九章 动态规划动态规划第十章第十章 向量优化（多目标优化）向量优化（多目标优化）第十一章第十一章 全局优化全局优化投资决策问题另外一种提法:某公司在一个时期内可用于投资的总资本为b万元,可供选择的项目有n个。假定对第i个项目的投资总额为ai万元，收益总额为ci万元。问如何确定投资方案，使总的利润率
2、达最高？设投资决策变量为：数学模型关于决策变量xi的非线性约束整数规划问题。收益占总投资的比率1.非线性规划模型：非线性规划模型：数学规划模型的一般形式：数学规划模型的一般形式：其中其中,x=(x1,x2,xn)T，f(x),gi(x),hj(x)为为x的实值函数的实值函数简记为简记为MP(Mathematical Programming)退 出前一页后一页2.1 基本概念基本概念可行域和可行解：可行域和可行解：称称为为MP问题的约束集或可行域。问题的约束集或可行域。若若x在在X内，称内，称x为为MP的可行解或者可行点。的可行解或者可行点。退 出前一页后一页简记形式：简记形式：引入向量函数符号
3、：引入向量函数符号：退 出前一页后一页数学规划问题的分类：数学规划问题的分类：若若f(x),gi(x),hj(x)为线性函数，即为为线性函数，即为线性规划线性规划(LP)；若若f(x),gi(x),hj(x)至少一个为非线性，即为至少一个为非线性，即为非线性规划非线性规划(NLP)；对于非线性规划，对于非线性规划，若没有若没有gi(x),hj(x)即即X=Rn,称为称为无约束无约束非线性规划非线性规划或或无约束最优化问题无约束最优化问题；否则称为；否则称为约束非线性规约束非线性规划或约束最优化问题划或约束最优化问题。退 出前一页后一页最优解和极小点最优解和极小点 对于数学规划（对于数学规划（M
4、P），若），若 ，并且有，并且有如果有如果有定义定义:退 出前一页后一页如果有如果有定义定义退 出前一页后一页例例退 出前一页后一页三角形表示的是可行域。三角形表示的是可行域。同心圆表示的是目标函数的等值同心圆表示的是目标函数的等值线。线。最优解为（最优解为（1/2，1/2）最优值为最优值为1/2问题：问题：(1/2,1/2)是整体的还是局部的？是严格的还是非是整体的还是局部的？是严格的还是非严格的？严格的？1/21/2退 出前一页后一页2.非线性规划方法概述非线性规划方法概述微分学方法的局限性：微分学方法的局限性：实际的问题中，函数可能是不连续或者不可微的。实际的问题中，函数可能是不连续或者
5、不可微的。需要解复杂的方程组，而方程组到目前仍没有有效需要解复杂的方程组，而方程组到目前仍没有有效的算法。的算法。实际的问题可能含有不等式约束，微分学方法不易实际的问题可能含有不等式约束，微分学方法不易处理。处理。退 出前一页后一页数值方法的基本思路：数值方法的基本思路：迭代迭代给定初始点给定初始点x0根据根据x0,依次迭代产生点列依次迭代产生点列xkxk的最后一点为最优解的最后一点为最优解xk有限有限xk无限无限xk收敛于最优解收敛于最优解退 出前一页后一页迭代格式迭代格式xkxk+1pk称称pk为第为第k轮轮搜索方向搜索方向，tk为第为第k轮沿轮沿pk方向的方向的步长步长。产生产生tk和和
6、pk的不同方法，形成了不同的算法。的不同方法，形成了不同的算法。退 出前一页后一页定义：下降方向定义：下降方向退 出前一页后一页定义定义解非线性规划问题，关键在于解非线性规划问题，关键在于找到某个方向，使得在此方向找到某个方向，使得在此方向上，目标函数得到下降，同时上，目标函数得到下降，同时还是可行方向。还是可行方向。这样的方向称为这样的方向称为可行下降方向。可行下降方向。退 出前一页后一页1.凸函数及其性质：凸函数及其性质：定义定义退 出前一页后一页2.2 凸函数和凸规划凸函数和凸规划退 出前一页后一页定理定理:关于凸函数的一些结论关于凸函数的一些结论定理定理:是凸集。是凸集。函数函数f在集
7、合在集合S上关于上关于c的水平集的水平集退 出前一页后一页定理定理?还有什么方法判断一个函数是凸函数呢？还有什么方法判断一个函数是凸函数呢？退 出前一页后一页退 出前一页后一页2.凸规划及其性质：凸规划及其性质：凸规划定义：凸规划定义：退 出前一页后一页凸规划性质：凸规划性质：凸规划的任一局部最优解都是它的整体最优解。凸规划的任一局部最优解都是它的整体最优解。凸规划是以后重点讨论的一类非线性规划凸规划是以后重点讨论的一类非线性规划凸函数凸函数线线性性函函数数退 出前一页后一页解：解：（1）目标函数是不是凸函数？）目标函数是不是凸函数？（2）gi(x)是不是凸函数？是不是凸函数？退 出前一页后一
8、页t为实数为实数一维搜索问题指目标函数为单变量的非线性规划问题。一维搜索问题指目标函数为单变量的非线性规划问题。又称线性搜索问题。其模型为：又称线性搜索问题。其模型为：什么叫一维搜索问题？什么叫一维搜索问题？或或一般一维搜索问题一般一维搜索问题有效一维搜索问题有效一维搜索问题退 出前一页后一页2.3一维搜索方法一维搜索方法一维搜索问题的算法分类：一维搜索问题的算法分类：精确一维搜索（最优一维搜索）精确一维搜索（最优一维搜索）非精确一维搜索（可接受一维搜索）非精确一维搜索（可接受一维搜索）本节内容：本节内容：两种精确一维搜索方法：两种精确一维搜索方法：0.618法，法，Newton法。法。两种非
9、精确一维搜索方法：两种非精确一维搜索方法：Goldstein(戈德斯坦）戈德斯坦）法法,Armijo（阿米霍）阿米霍）法。法。退 出前一页后一页1.0.618法（近似黄金分割法）法（近似黄金分割法）问题：问题：凸函数是不是单谷函数？严格凸函数是不凸函数是不是单谷函数？严格凸函数是不是单谷函数？单谷函数是不是凸函数？是单谷函数？单谷函数是不是凸函数？单谷函数单谷函数退 出前一页后一页搜索法求解：搜索法求解：或或基本过程：基本过程：给出给出a,b,使得使得t*在在a,b中。中。a,b称为称为搜索区间搜索区间。迭代缩短迭代缩短a,b的长度。的长度。当当a,b的长度小于某个预设的值，或者导数的绝对的长
10、度小于某个预设的值，或者导数的绝对值小于某个预设的正数，则迭代终止。值小于某个预设的正数，则迭代终止。退 出前一页后一页假定：已经确定了单谷区间假定：已经确定了单谷区间a,bt1t2ababt1t2新搜索区间为新搜索区间为a,t2新搜索区间为新搜索区间为t1,b退 出前一页后一页区间缩小比例的确定：区间缩小比例的确定：区间缩短比例为区间缩短比例为(t2-a)/(b-a)缩短比例为缩短比例为(b-t1)/(b-a)缩短比例缩短比例 满足：满足：每次插入搜索点使得两个区间每次插入搜索点使得两个区间a,t2和和t1,b相等；相等；每次迭代都以相等的比例缩小区间。每次迭代都以相等的比例缩小区间。0.6
11、18法法t1t2ababt1t2退 出前一页后一页确定确定a,b,计算探索点计算探索点t1=a+0.382(b-a)t2=a+0.618(b-a)0.618法解题步骤：法解题步骤：是是否否是是停止，输出停止，输出t1否否以以a,t2为新的搜索区间为新的搜索区间是是停止，输出停止，输出t2否否以以t1,b为新的搜索区间为新的搜索区间退 出前一页后一页例：例：解：解：t1t230t1.第一轮：第一轮：t1=1.146,t2=1.854t200.5退 出前一页后一页2.第二轮：第二轮：t2=1.146,t1=0.708t20=1.1460.53.第三轮：第三轮：t1=0.438,t2=0.708b-
12、t1=1.146-0.4380.51.8540tt2t11.4160tt2t1退 出前一页后一页4.第四轮：第四轮：t2=0.876,t1=0.708b-t1=1.146-0.7080.5输出：输出：t*=t2=0.876为最优解，最优值为为最优解，最优值为-0.0798课下练习：分析上述迭代过程，体会课下练习：分析上述迭代过程，体会0.618法的实质。法的实质。01.416tt1t2退 出前一页后一页2.Newton法法Newton法基本思想：法基本思想：用探索点用探索点tk处的二阶处的二阶Taylor展开式近似代替目标展开式近似代替目标函数，以展开式的最小点为新的探索点。函数，以展开式的最
13、小点为新的探索点。退 出前一页后一页解题步骤：解题步骤：给定初始点给定初始点t1和精度和精度是是是是停止，输出停止，输出t1是是否否停止，解题失败停止，解题失败否否停止，输出停止，输出t2否否退 出前一页后一页例：例：解：解：取取t1=1,计算：计算：迭代过程如下表：迭代过程如下表：1.1370.11630.11693-0.00106141.3258-0.5178-0.5708220.785411退 出前一页后一页3.非精确一维搜索法非精确一维搜索法数值方法的关键是从一个点迭代到下一个点。数值方法的关键是从一个点迭代到下一个点。确定下一个点的关键是确定搜索方向和步长确定下一个点的关键是确定搜索
14、方向和步长如果已经确定了搜索方向如果已经确定了搜索方向pk,则只要确定一个最佳则只要确定一个最佳的步长即可。的步长即可。所谓的最佳步长即是在所谓的最佳步长即是在pk方向上走一个最好的长度方向上走一个最好的长度使得目标函数下降的最多，即下述的最优化问题：使得目标函数下降的最多，即下述的最优化问题：这样的最优化问题不需要太高的精度，只要满足这样的最优化问题不需要太高的精度，只要满足某些更宽松的精度要求即可。某些更宽松的精度要求即可。这样的搜索方法称之为这样的搜索方法称之为非精确一维搜索方法非精确一维搜索方法退 出前一页后一页Goldstein法原理：法原理：yt0bcdaY=(0)+(0)tY=(
15、0)+m2(0)tY=(0)+m1(0)t退 出前一页后一页是是Goldstein算法算法确定确定m1,m2,t0,a=0,b=+(t0)(0)+m1 (0)t0(t0)(0)+m2(0)t0是是停止停止,输出输出t0否否a=a,b=t0,t1=(a+b)/2否否a=t0,b=b,t1=(a+b)/2(若若b=+,则则t1=a)退 出前一页后一页Armijo法原理：法原理：yt0tkMtk退 出前一页后一页本节课讨论本节课讨论n元函数的无约束非线性规划问题：元函数的无约束非线性规划问题：求解此类模型求解此类模型(UMP)的方法称为的方法称为无约束最优化方法无约束最优化方法。无约束最优化方法通常
16、有两类：无约束最优化方法通常有两类：解析法：解析法：要使用导数的方法；要使用导数的方法；直接法：直接法：无须考虑函数是否可导，直接使用函数值。无须考虑函数是否可导，直接使用函数值。退 出前一页后一页2.4无约束最优化方法无约束最优化方法1.无约束问题的最优性条件无约束问题的最优性条件定理定理1定理定理2梯度为梯度为0的点称为函数的的点称为函数的驻点。驻点。驻点可能是极小点，也可能是极大点，也可能即不是极大驻点可能是极小点，也可能是极大点，也可能即不是极大也不是极小，这时称为函数的也不是极小，这时称为函数的鞍点。鞍点。定理定理2说明：说明：UMP问题的局部最优解必是目标函数的驻点。问题的局部最优
17、解必是目标函数的驻点。注：注：退 出前一页后一页定理定理3定理定理4退 出前一页后一页例例解：解：1.先求出目标函数的全部驻点；先求出目标函数的全部驻点；2.利用充分条件判断驻点是不是最优点。利用充分条件判断驻点是不是最优点。退 出前一页后一页关于梯度的复习：关于梯度的复习：梯度是一个向量。梯度是一个向量。n元函数元函数f(x1,x2,xn)在某点在某点x处的处的梯度为：梯度为：梯度的方向与函数梯度的方向与函数f的等值线的一个法线方向相同，的等值线的一个法线方向相同，从较低的等值线指向较高的等值线。从较低的等值线指向较高的等值线。梯度的方向就是函数梯度的方向就是函数f的值增加最快的方向，其相反
18、的值增加最快的方向，其相反方向就是函数值降低最快的方向。方向就是函数值降低最快的方向。2.最速下降法最速下降法退 出前一页后一页最速下降法又称为最速下降法又称为梯度法梯度法，由，由Cauchy于于1847年给出。年给出。最速下降法解决的是最速下降法解决的是具有连续可微的目标函数具有连续可微的目标函数的的UMP问题。问题。最速下降法的基本思想：从当前点最速下降法的基本思想：从当前点xk出发寻找使得出发寻找使得目标函数下降最快的方向，即目标函数下降最快的方向，即负梯度方向负梯度方向。退 出前一页后一页最速下降法计算步骤：最速下降法计算步骤：选区初始点选区初始点x0和精度和精度计算计算是是否否停止，
19、输出停止，输出x0求求p0=计算计算t0,使使计算计算x1=x0+t0 p0退 出前一页后一页例例解：解：退 出前一页后一页说明：说明：观察观察P119的图，可以发现的图，可以发现x1 x0垂直于目标函数的垂直于目标函数的等值线（图中的虚线）在等值线（图中的虚线）在x0的切线；的切线；最速下降方法相邻的两个搜索方向是相互垂直的，最速下降方法相邻的两个搜索方向是相互垂直的，即即x1 x0垂直垂直x1 x2；最速下降法解决最速下降法解决UMP的缺陷：迭代点越靠近最优的缺陷：迭代点越靠近最优解则目标函数下降的速度越慢；解则目标函数下降的速度越慢；优点：迭代点列总是收敛的，而且计算过程简单。优点：迭代
20、点列总是收敛的，而且计算过程简单。退 出前一页后一页本节课讨论约束非线性规划问题本节课讨论约束非线性规划问题MP其中其中,x=(x1,x2,xn)T，f(x),gi(x),hj(x)为为x的实值函数的实值函数求解此类模型求解此类模型(MP)的方法称为的方法称为约束最优化方法约束最优化方法。退 出前一页后一页2.5约束最优化方法约束最优化方法1.约束最优化问题的最优性条件约束最优化问题的最优性条件对于对于MP问题：问题：退 出前一页后一页若若x*有变化，则约束条件可能没有破坏有变化，则约束条件可能没有破坏若若x*有变化，则约束条件一定被破坏有变化，则约束条件一定被破坏令令J表示表示MP的全部等式
21、约束的下标集合，即的全部等式约束的下标集合，即J=1,2q,I表示表示MP的全部不等式约束的下标集合，即的全部不等式约束的下标集合，即I=1,2px*的积极约束的下标集合的积极约束的下标集合退 出前一页后一页定理定理1对于对于若若x*是局部最优解是局部最优解,则则退 出前一页后一页定理定理1的说明：的说明：2、称下述表达式为、称下述表达式为MP的的Kuhn-Tucker条件，简称条件，简称K-T条件条件满足满足K-T条件的点称为条件的点称为MP的的K-T点，定理点，定理1说明说明MP的局部最的局部最优解一定是优解一定是MP的的K-T点。点。为了求出为了求出MP的最优解，可以先找出的最优解，可以
22、先找出MP的的K-T点，再做进一点，再做进一步的判断。步的判断。退 出前一页后一页3、定理、定理1的实例说明的实例说明定理定理1表明：若表明：若(x1,x2)T是局部最优解，是局部最优解，g1和和g2为积极约束，则：为积极约束，则：退 出前一页后一页4.定理定理1的特例的特例1退 出前一页后一页5.定理定理1的特例的特例2退 出前一页后一页6.定理定理1的改进：的改进：对于对于若若x*是局部最优解是局部最优解,则则互补松紧条件互补松紧条件退 出前一页后一页7.实例说明改实例说明改进后的定理进后的定理1：定理定理1改进后表明：若改进后表明：若(x1,x2)T是局部最优解，则：是局部最优解，则：退
23、 出前一页后一页互补松紧条件互补松紧条件退 出前一页后一页定理定理2对于对于注：定理注：定理2表明，在凸性条件下，表明，在凸性条件下，K-T点是整体最优解。点是整体最优解。退 出前一页后一页例：例：写出写出K-T条件条件;求出相应的求出相应的K-T点点;判断判断K-T点是不是点是不是问题的最优解问题的最优解解：解：由于全部函数都是连续可微的，所以应用以下由于全部函数都是连续可微的，所以应用以下K-T条件条件退 出前一页后一页首先写出原首先写出原MP问题的问题的K-T条件：条件：根据定理根据定理1,K-T点还应该满足原问题的约束条件点还应该满足原问题的约束条件互补松紧条件互补松紧条件退 出前一页
24、后一页利用互补松紧条件，可以求出利用互补松紧条件，可以求出K-T点：点：利用定理利用定理2，由于全部函数都连续可微，并且，由于全部函数都连续可微，并且f和和g都是凸函数，都是凸函数，h是线性函数，所以是线性函数，所以K-T点就是整体最点就是整体最优解。优解。退 出前一页后一页2.惩罚函数法惩罚函数法惩罚函数法的基本思想：利用原问题的中的约惩罚函数法的基本思想：利用原问题的中的约束函数构造适当的惩罚函数，并和原问题的目标束函数构造适当的惩罚函数，并和原问题的目标函数相加，得到带参数的增广目标函数，从而将函数相加，得到带参数的增广目标函数，从而将原问题题转换为一系列无约束非线性规划问题。原问题题转
25、换为一系列无约束非线性规划问题。惩罚函数法的分类：罚函数法（外部惩罚法），惩罚函数法的分类：罚函数法（外部惩罚法），障碍函数法（内部惩罚法）障碍函数法（内部惩罚法）退 出前一页后一页(1)罚函数法罚函数法罚函数法基本原理：罚函数法基本原理：考虑：考虑：构造惩罚函数：构造惩罚函数：很大的正数很大的正数无约束最优化问题无约束最优化问题min F(x)=f(x)+p(x)的最优解必定的最优解必定是原问题的最优解。是原问题的最优解。退 出前一页后一页可选的惩罚函数：可选的惩罚函数：惩罚函数法的经济解释：惩罚函数法的经济解释：f(x)为产品成本，约束条件为产品质量约束；为产品成本，约束条件为产品质量约束
26、；如果违反质量约束，就给予一定的惩罚如果违反质量约束，就给予一定的惩罚p(x)；追求的目标就是成本追求的目标就是成本f(x)和惩罚量和惩罚量p(x)的总和最小的总和最小（即构造的无约束最优化问题）；（即构造的无约束最优化问题）；如果惩罚条件很苛刻，最好的结果就是不违反质量如果惩罚条件很苛刻，最好的结果就是不违反质量约束（无约束最优化问题的最优解为约束（无约束最优化问题的最优解为MP的最优解）的最优解）退 出前一页后一页(2)障碍函数法障碍函数法障碍函数法基本原理：障碍函数法基本原理：构造一个新的目标函数，它在可行区域的边界筑构造一个新的目标函数，它在可行区域的边界筑起一道墙；起一道墙；当迭代点
27、靠近边界时，新的目标函数迅速增加；当迭代点靠近边界时，新的目标函数迅速增加；迭代点被档在可行区域的内部；迭代点被档在可行区域的内部；迭代得到的点列就只可能在可行区域的内部。迭代得到的点列就只可能在可行区域的内部。退 出前一页后一页可选的惩罚函数：可选的惩罚函数：考虑：考虑：构造最优化问题：构造最优化问题：或：或：当当x靠近边界时，至少有一个靠近边界时，至少有一个gi(x)趋近于零，则趋近于零，则F(x)将将无限增大，从而使得迭代点保持在可行区域的内部。无限增大，从而使得迭代点保持在可行区域的内部。退 出前一页后一页应用实例：应用实例：供应与选址供应与选址 某公司有6个建筑工地要开工，每个工地的
28、位置（用平面坐标系a，b表示，距离单位：千米）及水泥日用量d(吨)由下表给出。目前有两个临时料场位于A(5,1)，B(2,7)，日储量各有20吨。假设从料场到工地之间均有直线道路相连。（1）试制定每天的供应计划，即从A，B两料场分别向各工地运送多少吨水泥，使总的吨千米数最小。（2）为了进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量各为20吨，问应建在何处，节省的吨千米数有多大？（一）、建立模型（一）、建立模型 记工地的位置为记工地的位置为(ai，bi)，水泥日用量为，水泥日用量为di，i=1,6;料场位置为料场位置为(xj，yj)，日储量为，日储量为ej，j=1,2；从料场；从
29、料场j向工地向工地i的运送量为的运送量为Xij。当用临时料场时决策变量为：Xij，当不用临时料场时决策变量为：Xij，xj，yj。（二）使用临时料场的情形（二）使用临时料场的情形 使用两个临时料场A(5,1)，B(2,7).求从料场j向工地i的运送量为Xij，在各工地用量必须满足和各料场运送量不超过日储量的条件下，使总的吨千米数最小，这是线性规划问题.线性规划模型为：设X11=X1,X21=X 2,X31=X 3,X41=X 4,X51=X 5,X61=X 6X12=X 7,X22=X 8,X32=X 9,X42=X 10,X52=X 11,X62=X 12 编写程序gying1.mMATLA
30、B（gying1）cleara=1.25 8.75 0.5 5.75 3 7.25;b=1.25 0.75 4.75 5 6.5 7.75;d=3 5 4 7 6 11;x=5 2;y=1 7;e=20 20;for i=1:6 for j=1:2 aa(i,j)=sqrt(x(j)-a(i)2+(y(j)-b(i)2);endendCC=aa(:,1);aa(:,2);A=1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1;B=20;20;Aeq=1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0
31、 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1;beq=d(1);d(2);d(3);d(4);d(5);d(6);VLB=0 0 0 0 0 0 0 0 0 0 0 0;VUB=;x0=1 2 3 0 1 0 0 1 0 1 0 1;xx,fval=linprog(CC,A,B,Aeq,beq,VLB,VUB,x0)计算结果为：计算结果为：x=3.0000 5.0000 0.0000 7.0000 0.0000 1.0000 0.0000 0.0000 4
32、.0000 0.0000 6.0000 10.0000fval=136.2275（三）改建两个新料场的情形（三）改建两个新料场的情形 改建两个新料场，要同时确定料场的位置(xj,yj)和运送量Xij，在同样条件下使总吨千米数最小。这是非线性规划问题。非线性规划模型为：设 X11=X1,X21=X 2,X31=X 3,X41=X 4,X51=X 5,X61=X 6 X12=X 7,X22=X 8,X32=X 9,X42=X 10,X52=X 11,X62=X 12 x1=X13,y1=X14,x2=X15,y2=X16 （1）先编写M文件liaoch.m定义目标函数。MATLAB（liaoch）
33、function f=liaoch(x)a=1.25 8.75 0.5 5.75 3 7.25;b=1.25 0.75 4.75 5 6.5 7.75;d=3 5 4 7 6 11;e=20 20;f1=0;for i=1:6 s(i)=sqrt(x(13)-a(i)2+(x(14)-b(i)2);f1=s(i)*x(i)+f1;endf2=0;for i=7:12 s(i)=sqrt(x(15)-a(i-6)2+(x(16)-b(i-6)2);f2=s(i)*x(i)+f2;endf=f1+f2;(2)取初值为线性规划的计算结果及临时料场的坐标:x0=3 5 0 7 0 1 0 0 4 0
34、6 10 5 1 2 7;编写主程序gying2.m.MATLAB（gying2）x0=3 5 0 7 0 1 0 0 4 0 6 10 5 1 2 7;A=1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0;B=20;20;Aeq=1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0
35、 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0;beq=3 5 4 7 6 11;vlb=zeros(12,1);-inf;-inf;-inf;-inf;vub=;x,fval,exitflag=fmincon(liaoch,x0,A,B,Aeq,beq,vlb,vub)(3)计算结果为：x=3,4.9994,4,7,1.0006,0,0,0.0006,0,0,4.9994,11,5.6774,4.9055 7.2499,7.7500fval=89.8851exitflag=1(4)若修改主程序gying2.m,取初值为上面的计算结果:x0=3,4
36、.9994,4,7,1.0006,0,0,0.0006,0,0,4.9994,11,5.6774,4.9055,7.2499,7.7500得结果为:x=3.0000 5.0000 4.0000 7.0000 1.0000 0 0 0 0 0 5.0000 11.0000 5.6958 4.9283 7.2500 7.7500fval=89.8835exitflag=1总的吨千米数比上面结果略优.MATLAB（gying2）(5)若取初值为上面的计算结果:x0=3.0000 5.0000 4.0000 7.0000 1.0000 0 0 0 0 0 5.0000 11.0000 5.6958 4.9283 7.2500 7.7500 则计算结果为:x=3.0000 5.0000 4.0000 7.0000 1.0000 0 0 0 0 0 5.0000 11.0000 5.6958 4.9283 7.2500 7.7500fval=89.8835exitflag=0总的吨千米数89.8835与上面结果一样.第三讲结束第三讲结束