应用Lingo与Matlab软件求解线性规划问题比较探究.docVIP

应用Lingo与Matlab软件求解线性规划问题比较探究 　　摘要：本文就一个给定的线性规划模型，通过介绍优化软件Lingo和科学计算软件Matlab中求解线性规划问题的命令和函数，指出Lingo软件在求解线性规划问题上占有一定优势。 关键词：线性规划 Lingo软件 Matlab软件 最优解 线性规划由前苏联经济学家康托洛维奇提出，它主要研究的是在线性等式（或不等式）约束条件下，使某一线性目标函数取得最大值（或最小值）的问题。随着计算机技术的发展，借助软件可以快速对线性规划问题进行求解和分析。目前，能够求解规划问题的数学软件比较多，常见的有优化软件Lingo和科学计算软件Matlab。 本文以如下线性规划为例，分别利用这二种软件来求解，并就它们在求解线性规划上的差异进行对比分析。 minz=10.8x11+10.95x12+11.1x13+11.25x14+11.1x22+ 11.25x23+11.4x24+11x33+11.15x34+11.3x44； s.t.x11+x12+x13+x14”“ 　　X（3，1） 0.000000 X（3，2） 0.000000 X（3，3） 25.00000 X（3，4） 5.000000 X（4，1） 0.000000 X（4，2） 0.000000 X（4，3） 0.000000 X（4，4） 10.00000 显然最优解同上，只是输出格式不同而已。 2 Matlab求解线性规划 2.1 Matlab软件简介 目前，Matlab提供了四十多个工具箱，这些工具箱专门针对某些具体应用领域。Matlab优化工具箱中提供了linprog函数来求解线性规划问题。 2.2 Matlab求解线性规划的命令介绍 Matlab中一般使用“[ ]”、 “，”或空格以及“；”来创建数组，“[ ]”中给出数组的所有元素，同行间的元素用“，”或者空格隔开，不同行之间用分号“；”隔开，并且用符号“■”置于矩阵右上角表示矩阵的转置运算。 Linprog函数的常见形式如下： 形式1：X=linprog（f，A，b） 用于求解目标函数为Minf′*x ，约束条件为A*x≤b的线性规划问题。其中X表示最优解，f 表示价值列向量，A表示约束不等式中的系数矩阵（二维数组），b（列向量）表示约束不等式中右端资源常数向量。 形式2：[X，fval]=Linprog（c，A，b，Aeq，beq） 相比较上面的问题，增加了等式约束，即Aeq*x=beq。其中X、c、A、b含义同上，fval表示最优解对应的目标函数值。若没有不等式存在，则令A=[]、b=[]. 形式3：[X，fval]=Linprog（c，A，b，Aeq，beq，vlb，vub） 增加了决策变量的上下界约束，即Vlb≤x≤vub，其中vlb、vub分别以列向量形式存储。如果没有不等式约束，令A=[]和b=[]；若没有等式约束，则令Aeq=[]、beq=[]. 2.3 Matlab求解上述线性规划的具体实现 Matlab程序如下： clear c=[10.8 10.95 11.1 11.25；0 11.1 11.25 11.4；0 0 11 11.15；0 0 0 11.3]； A=[1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1]； b=[25；35；30；10]； Aeq=[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1]； beq=[10；15；25；20]； f=c（：）； %把f变成列向量 vlb=zeros（16，1）； %确定决策变量的下界 vub=[inf 0 0 0 inf inf 0 0 inf inf inf 0 inf inf inf inf]； %通过取上下界值都为0，保证决策变量x21=x31=x41 =x32=x42=x43=0 [x，fval]=linprog（f，A，b，Aeq，beq，vlb，vub） 运行后，得结果： x =10.0000 0