线性规划及非线性规划.ppt

由此得到规划模型: 注意到, 该规划是个非线性规划! 分别编写函数文件和脚本文件: 运行后得到问题的局部最优解: 相应的目标函数值为 与原问题比较, 节省路程 图形为 解的分析: 该解似乎不太合理! 经过多次试验, 问题的最优解应该为 最优值为 * s.t. (subject to),称为约束条件. * 注意：这样的表示方式是MatLab固有的形式。 可以缺省，但是只能从右边开始缺省，如中间需要缺省，则可输入空矩阵 * * * 程序Ch6_1.m * 程序Ch6_2.m * 程序Ch6_3.m * 程序cx1.lg * 程序cx2.lg * 程序cx2.lg * 程序cx3.lg * 程序cx4.lg * 程序cx6.lg * 程序cx7.lg H是二次项的系数矩阵。如ax1^2+bx2^2-cx1x2; H=[2a –c;-c 2b]. * * 程序Ch6_4.m * 程序Ch6_5.m * 程序Ch6_6.m * 程序Ch6_6.m * 程序Ch6_7.m * 程序Ch6_8.m 例 求解最小指派问题 相应程序如下: 例 求下面指派问题的最小值解: 例 求下面指派问题中的最小值解 其中某个人完成两项工作. 例 求解运输问题 有一连锁超市系统, 系统中4个存储站和10个门市部, 仓储站的库存量为196, 187, 179, 176; 10个门市部当前 的需求量分别为75, 69, 72, 83, 66, 65, 74, 62, 81, 56, 从仓储站到门市部的运输成本如下表: 4 4 7 5 5 6 8 7 6 5 4 7 8 9 8 5 7 5 4 6 6 3 8 7 5 3 4 7 6 9 4 5 2 8 5 6 7 8 9 6 7 8 7 1 10 9 8 7 6 5 4 3 2 1 求相应的运输方案, 使总运输成本为最小. 模型为 且为整数. 三、非线性规划 若在一个优化模型中, 目标函数为二次函数, 而约束条 件为线性函数, 则该模型称为二次规划. 1.二次规划 一个常见的二次规划标准形式为 求解二次规划的MatLab函数为: 格式 例 求解二次规划 解 输入下面语句 结果为 例 求解二次规划 解 输入下面语句 结果为 ⑵非线性规划 非线性规划的一般形式 求解函数为 格式 例 求下面问题的最优解, 取初始点为 将约束条件改写为标准形式: 首先编写函数文件 再编写脚本文件 运行后得到问题的解: 例 求解非线性规划 初始点分别取 首先编写函数文件 再编写脚本文件 运行后得到问题的解: 四、选址问题 问题 某公司有6个建筑工地要开工, 每个工地的位置 用平面坐标 来表示及水泥日用量 给出. 目前有 两个临时料场位于 日储量为 假 设从料场到工地之间均有直线相连, 试制定每天的供应 计划, 使总的吨公里数最小. 11 6 7 4 5 3 7.75 6.5 5 4.75 0.75 1.25 7.25 3 5.75 0.5 8.75 1.25 6 5 4 3 2 1 工地 工地的位置及水泥日用量 模型建立 记工地的位置为 水泥日用量为 料场 的日储量为 从料场 到工地 的运送量 为 则问题的目标函数为 相应的约束条件为 由此得到规划模型: 先画出相应的散点图: 编写脚本文件, 求最小值解. 运行后得到最优解为: 及相应的最优值: 进一步的问题是: 若选择新建两个料场, 则这两个料场应建在何处, 能 路程为最短, 又节省了多少路程? 模型建立 在上面模型中, 只需将原料场的坐标取为未知变量即 可, 因而有 解 输入下面语句 运行结果为 5.整数规划与 规划 在线性规划中, 若变量的取值均为整数时, 相应的规 划称为整数规划; 若仅要求其中部分为整数时, 相应的 规划为混合型规划. 在MatLab中, 无相应的求解程序! 二、应用Lingo软件求解线性规划 ⑴程序的安装和运行 在Lingo官方网站上, 有Lingo软件的学生版. 安装完毕 后, 在Windows主菜单中会产生Lingo程序相应的启动图 标: 启动程序后, 进入Lingo主窗口. 程序输入处 程序运行 例 求解线性规划 解 启动Lingo, 在主