优化建模与LINDOLINGO优化软件.ppt

MATLAB优化工具箱能求解的优化模型 LINDO 公司软件产品简要介绍 LINDO和LINGO软件能求解的优化模型 LINDO/LINGO软件的求解过程 建模时需要注意的几个基本问题 需要掌握的几个重要方面 使用LINDO的一些注意事项 “ ”（或“ ”）号与“ ”（或“ ”）功能相同 变量与系数间可有空格 甚至回车 , 但无运算符 变量名以字母开头，不能超过8个字符 变量名不区分大小写（包括LINDO中的关键字） 目标函数所在行是第一行，第二行起为约束条件 行号 行名 自动产生或人为定义。行名以“）”结束 行中注有“!”符号的后面部分为注释。如: ! It’s Comment. 在模型的任何地方都可以用“TITLE” 对模型命名（最多72个字符），如： TITLE This Model is only an Example 使用LINDO的一些注意事项 变量不能出现在一个约束条件的右端 表达式中不接受括号“ ”和逗号“,”等任何符号, 例: 400 X1+X2 需写为400X1+400X2 表达式应化简，如2X1+3X2- 4X1应写成 -2X1+3X2 缺省假定所有变量非负；可在模型的“END”语句后用“FREE name”将变量name的非负假定取消 可在 “END”后用“SUB” 或“SLB” 设定变量上下界 例如： “sub x1 10”的作用等价于“x1 10” 但用“SUB”和“SLB”表示的上下界约束不计入模型的约束，也不能给出其松紧判断和敏感性分析。 14. “END”后对0-1变量说明：INT n 或 INT name 15. “END”后对整数变量说明：GIN n 或 GIN name 二次规划 QP 问题 LINDO可求解二次规划 QP 问题，但输入方式较复杂，因为在LINDO中不许出现非线性表达式 需要为每一个实际约束增加一个对偶变量（LAGRANGE乘子），在实际约束前增加有关变量的一阶最优条件，转化为互补问题 “END”后面使用QCP命令指明实际约束开始的行号，然后才能求解 建议总是用LINGO解QP [注意]对QP和IP: 敏感性分析意义不大 状态窗口（LINDO Solver Status） 当前状态：已达最优解 迭代次数：18次 约束不满足的“量” 不是“约束个数” ：0 当前的目标值：94 最好的整数解：94 整数规划的界：93.5 分枝数：1 所用时间：0.00秒（太快了，还不到0.005秒） 刷新本界面的间隔:1 秒 选项设置 Report/Statistics LINDO行命令、命令脚本文件 LINGO模型 — 例：选址问题 选址问题：NLP 边界 状态窗口 平面示意图 模型假设 卡车在一个班次中不应发生等待或熄火后再启动的情况； 在铲位或卸点处由两条路线以上造成的冲突问题面前，我们认为只要平均时间能完成任务，就认为不冲突。我们不排时地进行讨论； 空载与重载的速度都是28km/h，耗油相差很大； 卡车可提前退出系统，等等。 符号 xij ：从i铲位到j号卸点的石料运量 （车） 单位： 吨； cij ：从i号铲位到j号卸点的距离 公里； Tij :从i号铲位到号j卸点路线上运行一个周期平均时间 分； Aij ：从号铲位到号卸点最多能同时运行的卡车数 辆； Bij ：从号铲位到号卸点路线上一辆车最多可运行的次数 次； pi：i号铲位的矿石铁含量 p 30,28,29,32,31,33,32,31,33,31 % qj : j号卸点任务需求，q 1.2,1.3,1.3,1.9,1.3 *10000 吨 cki ：i号铲位的铁矿石储量 万吨 cyi ：i号铲位的岩石储量 万吨 fi :描述第i号铲位是否使用的0-1变量，取1为使用；0为关闭。 优化模型 计算结果（LINGO软件） 计算结果（派车） 最大化产量 CUMCM其他优化赛题 主要参考文献 飞机精确定位模型 第2类模型: 考虑误差因素 作为硬约束 Min x; Min y; Max x; Max y. 以距离为约束，优化角度误差之和（或平方和）； 或以角度为约束，优化距离误差. 非线性规划 ？ ？ 仅部分考虑误差! 角度与距离的“地位”不应不同！ 有人也可能会采用其他目标，如： 误差非均匀分布！ 飞机精确定位模型 误差一般服从什么分布？ 正态分布！ 不同的量纲如何处理？ 无约束非线性最小二乘模型 归一化处理！ shili0702.m 飞机坐标 978.31，723.98 , 误差平方和0.6685 4 角度需要进行预处理，如利用 Matlab的atan2函数, 值域 -pi, pi 第3类模型: 考虑误差因素 作为软约束 ; 且归一化 飞机精确定位模型 小技巧: