《应用lingo、matlab软件求解线性规划》由会员分享，可在线阅读，更多相关《应用lingo、matlab软件求解线性规划（37页珍藏版）》请在金锄头文库上搜索。
1、1.5 应用LINGO、MATLAB软件求解线性规划1.5.1 应用LINGO软件求解线性规划 一、LINGO使用简介LINGO软件是美国的LINDO系统公司（Lindo System Inc）开发的一套用于求解最优化问题的软件包。LINGO除了能用于求解线性规划和二次规划外，还可以用于非线性规划 求解以及一些线性和非线性方程（组）的求解等。LINGO软件的最大特色在于它允许优化模型中的决策变量为整数，而 且执行速度快。LINGO内置了一种建立最优化模型的语言， 可以简便地表达大规模问题，利用LINGO高效的求解器可快 速求解并分析结果，这里简单介绍LINGO的使用方法。LINGO可以求解线性
2、规划、二次规划、非线性规划、整数规划、图论及网络优化和排队论模型中的最优化问题等。一个LINGO程序一般会包含集合段、数据输入段、优化目标和约束段、初始段和数据预处理段等部分，每一部分有其独特的 作用和语法规则，读者可以通过查阅相关的参考书或者LINGO 的HELP文件详细了解，这里就不展开介绍了。LINGO的主要功能特色为：既能求解线性规划问题，也有较强的求解非线性规划问题的能力；输入模型简练直观；运算速度 快、计算能力强；内置建模语言，提供几十个内部函数，从而 能以较少语句，较直观的方式描述大规模的优化模型；将集合 的概念引入编程语言，很容易将实际问题转换为LINGO模型； 并且能方便地与
3、Excel、数据库等其他软件交换数据。 LINGO的语法规定： （1）求目标函数的最大值或最小值分别用MAX=或 MIN=来表示； （2）每个语句必须以分号“；”结束，每行可以有许多语 句，语句可以跨行； （3）变量名称必须以字母(AZ)开头，由字母、数字 (09)和下划线所组成，长度不超过32个字符，不区分大 小写； （4）可以给语句加上标号，例如 OBJMAX=200*X1+300*X2； （5）以惊叹号“！”开头，以分号“；”结束的语句是注释 语句; （6）如果对变量的取值范围没有作特殊说明，则默认所 有决策变量都非负； （7）LINGO模型以语句“MODEL：”开头，以“END”结 束
4、，对于比较简单的模型，这两个语句可以省略。在LINGO的MODEL窗口内输入如下模型： model: max=2*x1+3*x2; x1+2*x260; 0.1*x1+0.05*x2+0.02*x3+0.2*x4+0.05*x53; 0.05*x1+0.1*x2+0.02*x3+0.2*x4+0.08*x58; X1+x2+x3+x4+x552; 例1.5.2 用LINGO求解例1.1.3食谱问题。 解 例1.1.3食谱问题的数学模型为（1.1.6）求解输出结果如下: Global optimal solution found at iteration: 4 Objective value:
5、22.40000 Variable Value Reduced Cost X1 0.000000 0.7000000 X2 12.00000 0.000000 X3 0.000000 0.6166667 X4 30.00000 0.000000 X5 10.00000 0.000000 Row Slack or Surplus Dual Price 1 22.40000 -1.000000 2 0.000000 -0.5833333 3 4.100000 0.000000 4 0.000000 -4.166667 5 0.000000 0.8833333因此,每周每个动物的配料为饲料A2、A4
6、、A5分别为12、30 和10，合计为52，可使得饲养成本达到最小，最小成本为 22.4元；不选用饲料和的原因是因为这两种饲料的价格太 高了，没有竞争力。“Reduced Cost”分别等于0.7和0.617 ，说明当这两种饲料的价格分别降低0.7元和0.62元以上时 ，不仅选用这两种饲料而且使得饲养成本降低。从“Slack or Surplus”可以看出，蛋白质和维生素刚达到最低标准，矿 物质超过最低标准4.1；从“Dual Price”可以得到降低标准蛋 白质1单位可使饲养成本降低0.583元，降低标准维生素1单 位可使饲养成本降低4.167元，但降低矿物质的标准不会降低饲养成本，如果动物
7、的进食量减少，就必须选取精一些 的饲料但要增加成本，大约进食量降低1可使得饲养成本增 加0.88元。(1.5.1) 1.5.2 应用MATLAB求解线性规划 MATLAB（MATrix LABoratory）的基本含义是矩阵实验室 ，它是由美国MathWorks公司研制开发的一套高性能的集数 值计算、信息处理、图形显示等于一体的可视化数学工具软 件。它是建立在向量、数组和矩阵基础之上的，除了基本的 数值计算、数据处理、图形显示等功能之外，还包含功能强 大的多个“工具箱”，如优化工具箱（optimization toolbox）、 统计工具箱、样条函数工具箱和数据拟合工具箱等都是优化 计算的有力
8、工具。在这里仅介绍用MATLAB6.5优化工具箱 求解线性规划问题。 一般线性规划问题的数学模型为 其中C是目标函数的系数行向量（常数）, X 是n维列 向量（决策变量），A, A1是常数矩阵，b,b1是常数 向量，lb,ub是n维列向量分别表示决策变量X的下界 与上界。 在Matlab优化工具箱（Optimization Toolbox）中， 求解(1.5.1)的程序如下：x，fval，exitflag，output ，lambda = linprog (c，A，b，Aeq，beq，lb，ub ，x0，options) 说明：(1)A是不等式约束的系数矩阵，b是相应的常 数列向量，若没有不等
9、式约束，则均用 代替； (2) Aeq是等式约束的系数矩阵，beq是相应的常数 列向量，若没有等式约束，则均用代替； (3)如果某个变量无下界，则用-inf表示；如果某个变 量无上界，则用inf表示，若决策变量 无下界，则lb 用代替；若决策变量 无上界，则ub用代替； (4) x0是线性规划的初始解，这种设计仅对中规模算 法有效，通常可以缺省。 (5) 输出 是最优解，fval是最优值。 (6) 输出exitflag描述了程序的运行情况，若其值大 于零，表示程序收敛到最优解 ；若其值等于零， 表示计算达到了最大次数；若其值小于零，表示问 题无可行解，或程序运行失败。 (7)输出output表
10、示程序运行的某些信息，如迭代 次数(iterations)、所用算法(algorithm)、共轭梯度 (cgiterations)等。 (8)lambda表示解处的拉格朗日乘子，其中lower， upper，ineqlin，eqlin分别对应于下界、上界、不 等式约束与等式约束。例1.5.3 用MATLAB解线性规划问题（1.5.2 ） 解 Matlab程序如下: c=-2,-1,1; A=1,4,-1;2,-2,1; b=4;12; Aeq=1,1,2; beq=6; lb=0,0,-inf; ub=inf,inf,5; x,z=linprog(c,A,b,Aeq,beq,lb,ub) 运行
11、后得到输出 Optimization terminated successfully. x=4.66670.00000.6667 z=-8.6667例1.5.4 用MATLAB求解线性规划问题（1.5.3） 解 首先转化为求最小值问题Matlab程序如下 c=-2,-3,5; A=-2,5,-1; b=-10; Aeq=1,1,1; beq=7; lb=0,0,0; x,z=linprog(c,A,b,Ae q,beq,lb)运行后得到输出 x =6.42860.57140.0000 z =-14.5714 键入s=-z 运行后得到原问题的目标 函数最大值s=14.5714用MATLAB求解例
12、1.5.2的程序与输出结果为： c=0.2,0.7,0.4,0.3,0.5; A=-0.3,-2,-1,-0.6,-1.8;-0.1,-0.05,-0.02,-0.2,-0.05;-0.05,-0.1,-0.02, -0.2,-0.08;1,1,1,1,1; b=-60;-3;-8;52; lb=0,0,0,0,0; x,z=linprog(c,A,b,lb)Optimization terminated successfully. x =0.000012.00000.000030.000010.0000 z =22.4000习题1 1.建立下列线性规划问题的数学模型 （1）某工厂生产A、B、
13、C三种产品，三种产品对于材 料费用、劳动力和电力的单位消耗系数，资源限量和 单位产品价格如表1.1所示。问应如何确定生产计划可 使得总产值达到最大？建立线性规划问题的数学模型 。 表1.1 生产计划问题的数据产产品 资资源ABC资资源 限量材料费费用（元）22.54320劳动劳动 力（人天）618640 电电力（度）5510750 单单位价格（百元 ）6410（2）某疗养院营养师要为某类病人拟订一周的菜单。可供选择的蔬菜 及其费用和所含营养成分的数量以及这类病人每周所需各种营养成分 的最低数量如表1.2所示。另外，为了口味的需要，规定一周内所用卷 心菜不多于2份，其他蔬菜不多于4份。若病人每周
14、需要14份蔬菜，问 选用每种蔬菜各多少份，可使生活费用最小。建立线性规划问题的数 学模型。 表1.2 食谱问题的数据蔬菜每份蔬菜所含营营养成分费费用 （元/份）铁铁(mg)磷(mg)VA(单单位 )VC(mg )烟酸(mg)青豆0.451041580.31.5 胡萝萝卜0.4528906530.351.5 花菜1.05502550530.62.4 卷心菜0.42575270.150.6 甜菜0.5221550.251.8 土豆0.57523580.81.0 每周营营养 最低需求量6.0325175002455.0播种计划表亩产量表土地作物作物面积100600700500500800500850
15、 400400150300 土地面积200300500（4）某糖果厂用原料A、B、C加工成三种不同牌号的糖 果甲、乙、丙，已知各种牌号的糖果中A、B、C的含量 ，原料成本，各种原料的每月限制用量，三种牌号糖果 的单位加工费及售价如表1.4所示，问该厂每月应生产这 三种牌号的糖果各多少千克，使该厂获利最大？试建立 这个问题的线性规划数学模型。甲乙丙原料成本 (元/千克)每月限制用量 （千克）A8.002000B6.002500C4.001200加工费 （元/千克）2.01.61.2售价（元）13.611.49表1.4 糖果厂生产计划数据表0.00.51.01.51.92.53.03.54.04.5 1.00.90.71.52.02.43.22.02.73.5 5.05.56.06.67.07.68.59.010.0 1.04.03.62.75.74.66.06.87.3表1.5 随 变化的数据表求拟合以上数据的直线 ，目标为使y的各个观察值同按直线关系所预期的值的绝对偏差总和为
