利用Matlab解决数学问题.pdf

利用Matlab解决数学问题

一、线性规划

求解线性规划的Matlab解法

单纯形法是求解线性规划问题的最常用、最有效的算法之一。单纯形法是首先由

GeorgeDantzig于1947年提出的，近60年来，虽有许多变形体已被开发，但却保持着同样

的基本观念。由于有如下结论：若线性规划问题有有限最优解，则一定有某个最优解是可行

区域的一个极点。基于此，单纯形法的基本思路是：先找出可行域的一个极点，据一定规则

判断其是否最优；若否，则转换到与之相邻的另一极点，并使目标函数值更优；如此下去，

直到找到某一最优解为止。这里我们不再详细介绍单纯形法，有兴趣的读者可以参看其它线

性规划书籍。下面我们介绍线性规划的Matlab解法。

Matlab5.3中线性规划的标准型为

T

mincxsuchthatAxb

x

基本函数形式为linprog(c,A,b)，它的返回值是向量的值。还有其它的一些函数调用形式（在x

Matlab指令窗运行helplinprog可以看到所有的函数调用形式），如：

[x,fval]=linprog(c,A,b,Aeq,beq,LB,UB,X,OPTIONS)

0

Aeq*xbeq

这里fval返回目标函数的值，Aeq和beq对应等式约束，LB和UB分别是变

x

x0x

量的下界和上界，是的初始值，OPTION是控制参数。

例2求解下列线性规划问题

maxz2x3x5x

123

xxx7



123

2x5xx10

123

x,x,x0



123

解（i）编写M文件

c=[2;3;-5];

a=[-2,5,-1];b=-10;

aeq=[1,1,1];

beq=7;

x=linprog(-c,a,b,aeq,beq,zeros(3,1))

value=c*x

（ii）将M文件存盘，并命名为example1.m。

（iii）在Matlab指令窗运行example1即可得所求结果。

例3求解线性规划问题

minz2x3xx

123

x4x2x8

123

3x2x6

12

x,x,x0

123

解编写Matlab程序如下：

c=[2;3;1];

a=[1,4,2;3,2,0];

b=[8;6];

[x,y]=linprog(c,-a,-b,[],[],zeros(3,1))

二、整数规划

整数规划问题的求解可以使用Lingo等专用软件。对于一般的整数规划规划问题，无法

直接利用Matlab的函数，必须利用Matlab编程实现分枝定界解法和割平面解法。但对于指

01A

派问题等特殊的整数规划问题或约束矩阵是幺模矩阵时，有时可以直接利用Matlab

的函数linprog。

例8求解下列指派问题，已知指派矩阵为

382103



87297

64275



84235

9106910

解：编写Matlab程序如下：

c=[382103;87297;64275

84235;9106910];

c=c(:)