数学建模算法大全线性规划

第一章 线性规划

§1 线性规划

在人们的生产实践中，经常会遇到如何利用现有资源来安排生产，以取得最大经济效益的问题。此类问题构成了运筹学的一个重要分支—数学规划，而线性规划(Linear Programming 简记LP)则是数学规划的一个重要分支。自从1947年G. B. Dantzig 提出求解线性规划的单纯形方法以来，线性规划在理论上趋向成熟，在实用中日益广泛与深入。特别是在计算机能处理成千上万个约束条件和决策变量的线性规划问题之后，线性规划的适用领域更为广泛了，已成为现代管理中经常采用的基本方法之一。

1.1 线性规划的实例与定义

例1某机床厂生产甲、乙两种机床，每台销售后的利润分别为4000元与3000元。生产甲机床需用B A 、机器加工，加工时间分别为每台2小时和1小时；生产乙机床需用C B A 、、三种机器加工，加工时间为每台各一小时。若每天可用于加工的机器时数分别为A 机器10小时、B 机器8小时和C 机器7小时，问该厂应生产甲、乙机床各几台，才能使总利润最大？

上述问题的数学模型：设该厂生产1x 台甲机床和2x 乙机床时总利润最大，则2

1,x x 应满足

（目标函数）2134m ax x x z += (1)

s.t.（约束条件）???????≥≤≤+≤+0

,781022122

121x x x x x x x （2）

这里变量21,x x 称之为决策变量，（1）式被称为问题的目标函数，（2）中的几个不等式

是问题的约束条件，记为s.t.(即subject to)。由于上面的目标函数及约束条件均为线性函数，故被称为线性规划问题。

总之，线性规划问题是在一组线性约束条件的限制下，求一线性目标函数最大或最小的问题。

在解决实际问题时，把问题归结成一个线性规划数学模型是很重要的一步，但往往也是困难的一步，模型建立得是否恰当，直接影响到求解。而选适当的决策变量，是我们建立有效模型的关键之一。

1.2 线性规划的Matlab 标准形式

线性规划的目标函数可以是求最大值，也可以是求最小值，约束条件的不等号可以是小于号也可以是大于号。为了避免这种形式多样性带来的不便，Matlab 中规定线性规划的标准形式为

b Ax x

c x

T ≤ that such min beq x Aeq =? ub x lb ≤≤

其中c 和x 为n 维列向量，A 、Aeq 为适当维数的矩阵，b 、beq 为适当维数的列向

量。

例如线性规划

b Ax x

c x

T ≥ that such max

的Matlab 标准型为

b Ax x

c x

T -≤-- that such min 1.3 线性规划问题的解的概念 一般线性规划问题的标准型为

∑==n

j j j x c z 1

max

(3)

s.t.

n

j x m

i b x a j n

j i

j ij

,,2,10

,,2,11

ΛΛ=≥==∑=(4)

可行解 满足约束条件（4）的解),,,(21n x x x x Λ=，称为线性规划问题的可行解，而使目标函数（3）达到最小值的可行解叫最优解。

可行域 所有可行解构成的集合称为问题的可行域，记为R 。

求解例1。对于每一固定的值z ，使目标函数值等于z 的点构成的直线称为目标函数等位线，当z 变动时，我们得到一族平行直线。对于例1，显然等位线越趋于右上方，其上的点具有越大的目标函数值。不难看出，本例的最优解为T

x )6,2(*=，最优目标值

26*=z 。

从上面的图解过程可以看出并不难证明以下断言：

（1）可行域R 可能会出现多种情况。R 可能是空集也可能是非空集合，当R 非空时，它必定是若干个半平面的交集（除非遇到空间维数的退化）。R 既可能是有界区域，也可能是无界区域。

（2）在R 非空时，线性规划既可以存在有限最优解，也可以不存在有限最优解（其目标函数值无界）。

（3）若线性规划存在有限最优解，则必可找到具有最优目标函数值的可行域R 的“顶点”。

上述论断可以推广到一般的线性规划问题，区别只在于空间的维数。在一般的n 维空间中，满足一线性等式

∑==n

i i

i b x

a 1

的点集被称为一个超平面，而满足一线性不等式

∑=≤n

i i

i b x

a 1

（或∑=≥n

i i i b x a 1

）的点集被称为一个半空间（其中),,(1n a a Λ为一n 维行

向量，b 为一实数）。若干个半空间的交集被称为多胞形，有界的多胞形又被称为多面体。易见，线性规划的可行域必为多胞形（为统一起见，空集Φ也被视为多胞形）。 在一般n 维空间中，要直接得出多胞形“顶点”概念还有一些困难。二维空间中的顶点可以看成为边界直线的交点，但这一几何概念的推广在一般n 维空间中的几何意义并不十分直观。为此，我们将采用另一途径来定义它。

定义 1 称n 维空间中的区域R 为一凸集，若R x x ∈?2

1

,及)1,0(∈?λ，有

R x x ∈-+21)1(λλ。

定义2 设R 为n 维空间中的一个凸集，R 中的点x 被称为R 的一个极点，若不

存在R x x ∈2

1、及)1,0(∈λ，使得21)1(x x x λλ-+=。

定义1 说明凸集中任意两点的连线必在此凸集中；而定义2 说明，若x 是凸集R 的一个极点，则x 不能位于R 中任意两点的连线上。不难证明，多胞形必为凸集。同样也不难证明，二维空间中可行域R 的顶点均为R 的极点（R 也没有其它的极点）。

1.5 求解线性规划的Matlab 解法

单纯形法是求解线性规划问题的最常用、最有效的算法之一。这里我们就不介绍单纯形法，有兴趣的读者可以参看其它线性规划书籍。下面我们介绍线性规划的Matlab 解法。

Matlab 中线性规划的标准型为

b Ax x

c T x

≤ such that min

beq x Aeq =?

ub x lb ≤≤

基本函数形式为linprog(c,A,b)，它的返回值是向量x 的值。还有其它的一些函数调用形式（在 Matlab 指令窗运行 help linprog 可以看到所有的函数调用形式），如：

[x,fval]=linprog(c,A,b,Aeq,beq,LB,UB,X 0,OPTIONS) 这里fval 返回目标函数的值，LB 和UB 分别是变量x 的下界和上界，0x 是x 的初始值，OPTIONS 是控制参数。

例2 求解下列线性规划问题 321532m ax x x x z -+=

???

??≥≥+-=++0,,105273

21321321x x x x x x x x x 解 （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 求解线性规划问题

32132 m in x x x z ++=

???

??≥≥+≥++0,,6

238243

2121321x x x x x x x x 解 编写Matlab 程序如下： c=[2;3;1];

a=[1,4,2;3,2,0]; b=[8;6];

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

1.6 可以转化为线性规划的问题

很多看起来不是线性规划的问题也可以通过变换变成线性规划的问题来解决。如： 例4 规划问题为

b

Ax x x x n ≤+++ t.s.|

|||||min 21Λ

其中T

n x x x ][1Λ=，A 和b 为相应维数的矩阵和向量。

要把上面的问题变换成线性规划问题，只要注意到事实：对任意的i x ，存在

0,>i i v u 满足

i i i v u x -=，i i i v u x +=||

事实上，我们只要取2||i i i x x u +=，2||i

i i x x v -=就可以满足上面的条件。

这样，记T n u u u ][1Λ=，T

n v v v ][1Λ=，从而我们可以把上面的问题

变成

∑=+n

i i i

v u

1

)(min

??

?≥≤-0

,)( t.s.v u b

v u A 例5 |}|max {min i y x i

i

ε

其中i i i y x -=ε。

对于这个问题，如果我们取||max 0i y i

x ε=，这样，上面的问题就变换成

0min

x

0011,, t.s.x y x x y x n n ≤-≤-Λ

此即我们通常的线性规划问题。 §2 运输问题(产销平衡)

例6 某商品有m 个产地、n 个销地，各产地的产量分别为m a a ,,1Λ，各销地的需求量分别为n b b ,,1Λ。若该商品由i 产地运到j 销地的单位运价为ij c ，问应该如何调运才能使总运费最省？

解：引入变量ij x ，其取值为由i 产地运往j 销地的该商品数量，数学模型为

∑∑==m i n

j ij

ij x

c 11

min

s.t. ????

?????≥====∑∑==0,,2,1,,,1,1

1ij m

i j ij n

j i ij x n j b x m i a x ΛΛ

显然是一个线性规划问题，当然可以用单纯形法求解。

对产销平衡的运输问题，由于有以下关系式存在：

∑∑∑∑∑∑=======??? ??=???? ??=m

i i n

j n j m i ij m

i n j ij j a x x b 1

11111

其约束条件的系数矩阵相当特殊，可用比较简单的计算方法，习惯上称为表上作业法（由

康托洛维奇和希奇柯克两人独立地提出，简称康—希表上作业法）。 §3 指派问题

3.1 指派问题的数学模型

例7拟分配n 人去干n 项工作，每人干且仅干一项工作，若分配第i 人去干第j 项工作，需花费ij c 单位时间，问应如何分配工作才能使工人花费的总时间最少？

容易看出，要给出一个指派问题的实例，只需给出矩阵)(ij c C =，C 被称为指派问题的系数矩阵。

引入变量ij x ，若分配i 干j 工作，则取1=ij x ，否则取0=ij x 。上述指派问题的数学模型为

∑∑==n i n

j ij

ij x

c 11

min

s.t.

1 0111

1????

?????===∑∑==或ij n

i ij n

j ij x x x (5)

(5)的可行解既可以用一个矩阵表示，其每行每列均有且只有一个元素为1，其余元素均为0，也可以用n ,,1Λ中的一个置换表示。

(5)的变量只能取0或1，从而是一个0-1规划问题。一般的0-1规划问题求解极为困难。但指派问题并不难解，其约束方程组的系数矩阵十分特殊（被称为全单位模矩阵，其各阶非零子式均为1±），其非负可行解的分量只能取0或1，故约束10或=ij x 可改写为0≥ij x 而不改变其解。此时，指派问题被转化为一个特殊的运输问题，其中n m =，

1==j i b a 。

3.2求解指派问题的匈牙利算法

由于指派问题的特殊性，又存在着由匈牙利数学家Konig 提出的更为简便的解法—匈牙利算法。算法主要依据以下事实：如果系数矩阵)(ij c C =一行（或一列）中每一元素都加上或减去同一个数，得到一个新矩阵)(ij b B = ，则以C 或B 为系数矩阵的指派问题具有相同的最优指派。

例8 求解指派问题，其系数矩阵为

?

?

??

?????

???=16221917171822241819211722191516C

解 将第一行元素减去此行中的最小元素15，同样，第二行元素减去17，第三行元素减去17，最后一行的元素减去16，得

?

?

??

?

????

???=0631

0157124074011B

再将第3列元素各减去1，得

?????

???????=**

**20531005711407301B 以2B 为系数矩阵的指派问题有最优指派 ???

? ??43124321 由等价性，它也是例7的最优指派。 有时问题会稍复杂一些。

例9 求解系数矩阵C 的指派问题

???????

?????????=61071041066141512141217766698979712C

解：先作等价变换如下

∨

∨

∨????????

?????

???→????????????????----- 2636040*08957510*00*0032202*05610710410661415121412177666989797124

6767 容易看出，从变换后的矩阵中只能选出四个位于不同行不同列的零元素，但5=n ，最优指派还无法看出。此时等价变换还可进行下去。步骤如下：

(1) 对未选出0元素的行打∨； (2) 对∨行中0元素所在列打∨；

(3) 对∨列中选中的0元素所在行打∨； 重复（2）、（3）直到无法再打∨为止。

可以证明，若用直线划没有打∨的行与打∨的列，就得到了能够覆盖住矩阵中所有零元素的最少条数的直线集合，找出未覆盖的元素中的最小者，令∨行元素减去此数，∨列元素加上此数，则原先选中的0元素不变，而未覆盖元素中至少有一个已转变为0，且新矩阵的指派问题与原问题也等价。上述过程可反复采用，直到能选取出足够的0元素为止。例如，对例5变换后的矩阵再变换，第三行、第五行元素减去2，第一列元素加上2，得

???

??

?

?

?

????????04140400811353800003420207

现在已可看出，最优指派为?

??

?

??5314254321。 §4 对偶理论与灵敏度分析

4.1原始问题和对偶问题

考虑下列一对线性规划模型：

x c T max

s.t. 0,≥≤x b Ax (P) 和

y b T min

s.t. 0,≥≥y c y A T (D)

称（P ）为原始问题，（D ）为它的对偶问题。

不太严谨地说，对偶问题可被看作是原始问题的“行列转置”：

（1） 原始问题中的第j 列系数与其对偶问题中的第j 行的系数相同； （2） 原始目标函数的各个系数行与其对偶问题右侧的各常数列相同； （3） 原始问题右侧的各常数列与其对偶目标函数的各个系数行相同； （4） 在这一对问题中，不等式方向和优化方向相反。 考虑线性规划：

0,s.t.m in

≥=x b Ax x

c T

把其中的等式约束变成不等式约束，可得

0, s.t. min ≥??

?

???-≥??????-x b b x A A x c T

它的对偶问题是

[

][

]

c y y A A y y b b T

T T

T

≤??

?

???-?

?

?

???-2121 s.t.max

其中1y 和2y 分别表示对应于约束b Ax ≥和b Ax -≥-的对偶变量组。令21y y y -=，

则上式又可写成

c y A y b T T ≤ s.t. m ax

原问题和对偶的对偶约束之间的关系：

min max

???

??≤≥无限制变量

00?????=≥≤000行约束 ??

?

??=≤≥000行约束

??

?

??≤≥无限制变量00 4.2 对偶问题的基本性质

1o 对称性：对偶问题的对偶是原问题。

2o 弱对偶性：若x 是原问题的可行解，y 是对偶问题的可行解。则存在

y b x c T T ≤。

3o 无界性：若原问题（对偶问题）为无界解，则其对偶问题（原问题）无可行解。

4o 可行解是最优解时的性质：设x

?是原问题的可行解，y ?是对偶问题的可行解，当y b x

c T

T ??=时，y x ?,?是最优解。 5o 对偶定理：若原问题有最优解，那么对偶问题也有最优解；且目标函数值相同。

6o 互补松弛性：若y x

?,?分别是原问题和对偶问题的最优解，则 0)?(?,0)?(?=-=-c y A x

b x A y T T T 例10 已知线性规划问题

5432132532m in

x x x x x ++++=ω

43254321≥++++x x x x x 33254321≥+++-x x x x x 5,,2,1,

0Λ=≥j x j

已知其对偶问题的最优解为5;5

3

,54*2*

1===z y y 。试用对偶理论找出原问题的最优解。

解 先写出它的对偶问题

2134m ax y y z +=

2221≤+y y ①

321≤-y y ② 53231≤+y y ③ 221≤+y y ④ 3321≤+y y ⑤ 0,21≥y y

将*

2*1,y y 的值代入约束条件，得②，③，④为严格不等式；由互补松弛性得

0*4*3*2===x x x 。因 0,*2*1>y y ；原问题的两个约束条件应取等式，故有

43*

5*1=+x x

32*5*1=+x x

求解后得到1,1*

5*1==x x ；故原问题的最优解为

5;]'10001[**==ωX 。

4.3 灵敏度分析

在以前讨论线性规划问题时，假定j i ij c b a ,,都是常数。但实际上这些系数往往是估计值和预测值。如市场条件一变，j c 值就会变化；ij a 往往是因工艺条件的改变而改变；

i b 是根据资源投入后的经济效果决定的一种决策选择。因此提出这样两个问题：当这

些系数有一个或几个发生变化时，已求得的线性规划问题的最优解会有什么变化；或者这些系数在什么范围内变化时，线性规划问题的最优解或最优基不变。这里我们就不讨论了。

4.4 参数线性规划

参数线性规划是研究j i ij c b a ,,这些参数中某一参数连续变化时，使最优解发生变化的各临界点的值。即把某一参数作为参变量，而目标函数在某区间内是这参变量的线性函数，含这参变量的约束条件是线性等式或不等式。因此仍可用单纯形法和对偶单纯形法进行分析参数线性规划问题。

习 题一

1．试将下述问题改写成线性规划问题：

????

????? ??∑∑∑===m

i i in m i i i m i i i x x a x a x a i

11211,,,min max Λ ???=≥=+++m i x x x x i

m ,,1,01st.21ΛΛ

2．试将下列问题改写成线性规划问题：

∑==n

j j j x c z 1

||max

?????==∑=取值无约束j

n

j i j ij x m i b x a 1),,2,1(st.Λ

3．线性回归是一种常用的数理统计方法，这个方法要求对图上的一系列点),(,),,(),,(2211n n y x y x y x Λ选配一条合适的直线拟合。方法通常是先定直线方程为

bx a y +=，然后按某种准则求定b a ,。通常这个准则为最小二乘法，但也可用其他准

则。试根据以下准则建立这个问题的线性规划模型：

∑=+-n

i i i

bx a y

1

|)(|min

4．某厂生产三种产品I ，II ，III 。每种产品要经过B A ,两道工序加工。设该厂有两种规格的设备能完成A 工序，它们以21,A A 表示；有三种规格的设备能完成B 工序，它们以321,,B B B 表示。产品I 可在B A ,任何一种规格设备上加工。产品II 可在任何规格的A 设备上加工，但完成B 工序时，只能在1B 设备上加工；产品III 只能在2A 与2B 设备上加工。已知在各种机床设备的单件工时，原材料费，产品销售价格，各种设备有效台时以及满负荷操作时机床设备的费用如下表，要求安排最优的生产计划，使该厂利

6．某战略轰炸机群奉命摧毁敌人军事目标。已知该目标有四个要害部位，只要摧毁其中之一即可达到目的。为完成此项任务的汽油消耗量限制为48000升、重型炸弹48枚、轻型炸弹32枚。飞机携带重型炸弹时每升汽油可飞行2千米，带轻型炸弹时每升汽油可飞行3千米。又知每架飞机每次只能装载一枚炸弹，每出发轰炸一次除来回路程汽油消耗（空载时每升汽油可飞行4千米）外，起飞和降落每次各消耗100升。有关

题的线性规划模型。

7．用Matlab 求解下列线性规划问题： 3213m ax x x x z --=

s.t. ???????≥=+-≥++-≤+-0

,,12324112321313

21321x x x x x x x x x x x

8．用Matlab 求解下列规划问题：

||4||3||2||m in

4421x x x x z +++=

s.t. ???

?

???

-

=+--=-+-=+--2132130432143214321x x x x x x x x x x x x

第二章 整数规划

§1 概论

1.1 定义 规划中的变量（部分或全部）限制为整数时，称为整数规划。若在线性规划模型中，变量限制为整数，则称为整数线性规划。目前所流行的求解整数规划的方法，往往只适用于整数线性规划。目前还没有一种方法能有效地求解一切整数规划。

1.2 整数规划的分类

如不加特殊说明，一般指整数线性规划。对于整数线性规划模型大致可分为两类： 1o 变量全限制为整数时，称纯（完全）整数规划。 2o 变量部分限制为整数的，称混合整数规划。 1.2 整数规划特点

（i ）原线性规划有最优解，当自变量限制为整数后，其整数规划解出现下述情况： ①原线性规划最优解全是整数，则整数规划最优解与线性规划最优解一致。 ②整数规划无可行解。 例1 原线性规划为 21m in x x z +=

0,0,5422121≥≥=+x x x x

其最优实数解为：4

5

min ,45,021===z x x 。

③有可行解（当然就存在最优解），但最优解值变差。

例2 原线性规划为 21m in x x z +=

0,0,6422121≥≥=+x x x x

其最优实数解为：2

3

min ,23,021===z x x 。

若限制整数得：2m in ,1,121===z x x 。

（ii ）整数规划最优解不能按照实数最优解简单取整而获得。 1.3 求解方法分类：

（i ）分枝定界法—可求纯或混合整数线性规划。 （ii ）割平面法—可求纯或混合整数线性规划。 （iii ）隐枚举法—求解“0-1”整数规划： ①过滤隐枚举法； ②分枝隐枚举法。

（iv ）匈牙利法—解决指派问题（“0-1”规划特殊情形）。 （v ）蒙特卡洛法—求解各种类型规划。

下面将简要介绍常用的几种求解整数规划的方法。 §2 分枝定界法

对有约束条件的最优化问题（其可行解为有限数）的所有可行解空间恰当地进行系统搜索，这就是分枝与定界内容。通常，把全部可行解空间反复地分割为越来越小的子集，称为分枝；并且对每个子集内的解集计算一个目标下界（对于最小值问题），这称为定界。在每次分枝后，凡是界限超出已知可行解集目标值的那些子集不再进一步分枝，这样，许多子集可不予考虑，这称剪枝。这就是分枝定界法的主要思路。

分枝定界法可用于解纯整数或混合的整数规划问题。在本世纪六十年代初由Land Doig 和Dakin 等人提出的。由于这方法灵活且便于用计算机求解，所以现在它已是解整数规划的重要方法。目前已成功地应用于求解生产进度问题、旅行推销员问题、工厂选址问题、背包问题及分配问题等。

设有最大化的整数规划问题A ，与它相应的线性规划为问题B ，从解问题B 开始，若其最优解不符合A 的整数条件，那么B 的最优目标函数必是A 的最优目标函数*

z 的上界，记作z ；而A 的任意可行解的目标函数值将是*

z 的一个下界z 。分枝定界法就是将B 的可行域分成子区域的方法。逐步减小z 和增大z ，最终求到*

z 。现用下例来说明：

例3 求解下述整数规划 219040Max x x z +=

???

??≥≤+≤+且为整数0,70

20756792

12121x x x x x x 解 （i ）先不考虑整数限制，即解相应的线性规划B ，得最优解为：

355.8779,8168.1,8092.421===z x x

可见它不符合整数条件。这时z 是问题A 的最优目标函数值*

z 的上界，记作z 。而

0,021==x x 显然是问题A 的一个整数可行解，这时0=z ，是*

z 的一个下界，记作z ，

即3560*

≤≤z 。

（ii ）因为21,x x 当前均为非整数，故不满足整数要求，任选一个进行分枝。设选1

x 进行分枝，把可行集分成2个子集：

44.8092][1=≤x ，514.8092][1=+≥x

因为4与5之间无整数，故这两个子集的整数解必与原可行集合整数解一致。这一步称为分枝。这两个子集的规划及求解如下：

问题1B ： 219040Max x x z +=

??

?

??≥≤≤≤+≤+0

,40702075679212121x x x x x x 最优解为：349,1.2,0.4121===z x x 。

问题2B ： 219040Max x x z += ???

??≥≥≤+≤+0

,570207567921

2121x x x x x x 最优解为：4.341,57.1,0.5121===z x x 。

再定界：3490*

≤≤z 。

（iii ）对问题1B 再进行分枝得问题11B 和12B ，它们的最优解为

340,2,4:112111===z x x B

327.14,00.3x 1.43,:122112===z x B

再定界：341340*

≤≤z ，并将12B 剪枝。

（iv ）对问题2B 再进行分枝得问题21B 和22B ，它们的最优解为

083,00.1x 5.44,:222121===z x B 22B 无可行解。

将2221,B B 剪枝。

于是可以断定原问题的最优解为：

340,2,4*21===z x x

从以上解题过程可得用分枝定界法求解整数规划（最大化）问题的步骤为：

开始，将要求解的整数规划问题称为问题A ，将与它相应的线性规划问题称为问题B 。

（i ）解问题B 可能得到以下情况之一：

（a ）B 没有可行解，这时A 也没有可行解，则停止．

（b ）B 有最优解，并符合问题A 的整数条件，B 的最优解即为A 的最优解，则停止。

（c ）B 有最优解，但不符合问题A 的整数条件，记它的目标函数值为z 。

（ii ）用观察法找问题A 的一个整数可行解，一般可取n j x j ,,1,0Λ==，试探，求得其目标函数值，并记作z 。以*

z 表示问题A 的最优目标函数值；这时有

z z z ≤≤*

进行迭代。

第一步：分枝，在B 的最优解中任选一个不符合整数条件的变量j x ，其值为j b ，以][j b 表示小于j b 的最大整数。构造两个约束条件

][j j b x ≤ 和 1][+≥j j b x

将这两个约束条件，分别加入问题B ，求两个后继规划问题1B 和2B 。不考虑整数条件求解这两个后继问题。

定界，以每个后继问题为一分枝标明求解的结果，与其它问题的解的结果中，找出最优目标函数值最大者作为新的上界z 。从已符合整数条件的各分支中，找出目标函数值为最大者作为新的下界z ，若无作用z 不变。

第二步：比较与剪枝，各分枝的最优目标函数中若有小于z 者，则剪掉这枝，即以后不再考虑了。若大于z ，且不符合整数条件，则重复第一步骤。一直到最后得到

z z =*为止。得最优整数解n j x j ,,1

,*Λ=。 §3 10-型整数规划

10-型整数规划是整数规划中的特殊情形，它的变量j x 仅取值0或1。这时j x 称为10-变量，或称二进制变量。j x 仅取值0或1这个条件可由下述约束条件：

10≤≤j x ，整数

所代替，是和一般整数规划的约束条件形式一致的。在实际问题中，如果引入 10-变量，就可以把有各种情况需要分别讨论的线性规划问题统一在一个问题中讨论了。我们先介绍引入10-变量的实际问题，再研究解法。

3.1 引入10-变量的实际问题

3.1.1投资场所的选定——相互排斥的计划

例4某公司拟在市东、西、南三区建立门市部。拟议中有7个位置（点）)7,,2,1(Λ=i A i 可供选择。规定

在东区。由321,,A A A 三个点中至多选两个； 在西区。由54,A A 两个点中至少选一个；

在南区，由76,A A 两个点中至少选一个。

如选用i A 点，设备投资估计为i b 元，每年可获利润估计为i c 元，但投资总额不能超过B 元。问应选择哪几个点可使年利润为最大？

解题时先引入10-变量)7,,2,1(Λ=i x i 令

?

??=.0,1点没被选中当点被选中当，，i A i A i x 7,,2,1Λ=i .

于是问题可列写成：

i i i x c z ∑==7

1

Max

?????

????=≥+≥+≤++≤∑=10,

11

2765

43217

1或i i i i x x x x x x x x B x b

3.1.2相互排斥的约束条件

有两个相互排斥的约束条件

244521≤+x x 或 453721≤+x x 。

为了统一在一个问题中，引入10-变量y ，则上述约束条件可改写为：

??

?

??=-+≤++≤+10)1(453724452121或y M y x x yM x x 其中M 是充分大的数。

约束条件

01=x 或 8005001≤≤x 可改写为

?

?

?=≤≤108005001或y y

x y 如果有m 个互相排斥的约束条件：

m i b x a x a i

n in i ,,2,111ΛΛ=≤++

为了保证这m 个约束条件只有一个起作用，我们引入m 个10-变量),,2,1(m i y i Λ=和一个充分大的常数M ，而下面这一组1+m 个约束条件

m i M y b x a x a i i n in i ,,2,111ΛΛ=+≤++（1）

11-=++m y y m Λ（2）

就合于上述的要求。这是因为，由于（2），m 个i y 中只有一个能取0值，设0*=i y ，代入（1），就只有*

i i =的约束条件起作用，而别的式子都是多余的。

3.1.3 关于固定费用的问题（Fixed Cost Problem ）

在讨论线性规划时，有些问题是要求使成本为最小。那时总设固定成本为常数，并在线性规划的模型中不必明显列出。但有些固定费用（固定成本）的问题不能用一般线性规划来描述，但可改变为混合整数规划来解决，见下例。

例5 某工厂为了生产某种产品，有几种不同的生产方式可供选择，如选定的生产方式投资高（选购自动化程度高的设备），由于产量大，因而分配到每件产品的变动成本就降低；反之，如选定的生产方式投资低，将来分配到每件产品的变动成本可能增加。所以必须全面考虑。今设有三种方式可供选择，令

j x 表示采用第j 种方式时的产量；

j c 表示采用第j 种方式时每件产品的变动成本；

j k 表示采用第j 种方式时的固定成本。

为了说明成本的特点，暂不考虑其它约束条件。采用各种生产方式的总成本分别为

?????=>+= 0

,00,j j j j j j x x x c k P 当当3,2,1=j .

在构成目标函数时，为了统一在一个问题中讨论，现引入10-变量j y ，令

??

??

?==>.00,0,1时种生产方式，即当不采用第时，

种生产方式，即当采用第j x j j x j j y （3） 于是目标函数

)()()(m in 333322221111x c y k x c y k x c y k z +++++=

（3）式这个规定可表为下述3个线性约束条件：

3,2,1,=≤≤j M y x y j j j ε（4）

其中ε是一个充分小的正常数，M 是个充分大的正常数。（4）式说明，当0>j x 时j y 必须为1；当0=j x 时只有j y 为0时才有意义，所以（4）式完全可以代替（3）式。 3.2 10-型整数规划解法之一（过滤隐枚举法）

解10-型整数规划最容易想到的方法，和一般整数规划的情形一样，就是穷举法，即检查变量取值为0或1的每一种组合，比较目标函数值以求得最优解，这就需要检查变量取值的n 2个组合。对于变量个数n 较大（例如10>n ），这几乎是不可能的。因此常设计一些方法，只检查变量取值的组合的一部分，就能求到问题的最优解。这样的方法称为隐枚举法（Implicit Enumeration ），分枝定界法也是一种隐枚举法。当然，对有些问题隐枚举法并不适用，所以有时穷举法还是必要的。

下面举例说明一种解10-型整数规划的隐枚举法。 例6321523Max x x x z +-=

?????

????=≤+≤+≤++≤-+1

0,,6

43

44223213

221321321或x x x x x x x x x x x x x 求解思路及改进措施：

（i ） 先试探性求一个可行解，易看出)0,0,1(),,(321=x x x 满足约束条件，故为一个可行解，且3=z 。

（ii ）因为是求极大值问题，故求最优解时，凡是目标值3

（iii ）改进过滤条件。

（iv ）由于对每个组合首先计算目标值以验证过滤条件，故应优先计算目标值z

大的组合，这样可提前抬高过滤门槛，以减少计算量。 §4 蒙特卡洛法（随机取样法）

前面介绍的常用的整数规划求解方法，主要是针对线性整数规划而言，而对于非线性整数规划目前尚未有一种成熟而准确的求解方法，因为非线性规划本身的通用有效解法尚未找到，更何况是非线性整数规划。

然而，尽管整数规划由于限制变量为整数而增加了难度；然而又由于整数解是有限个，于是为枚举法提供了方便。当然，当自变量维数很大和取值范围很宽情况下，企图用显枚举法（即穷举法）计算出最优值是不现实的，但是应用概率理论可以证明，在一定的计算量的情况下，完全可以得出一个满意解。

例7 已知非线性整数规划为：

5

43212

5

242322212328 243Max x x x x x x x x x x z -----++++=

?????

????≤++≤++≤++++≤++++=≤≤200

5200

62800622400)5,,1(9905433

215432154321x x x x x x x x x x x x x x x x i x i Λ 对该题，目前尚无有效方法求出准确解。如果用显枚举法试探，共需计算

10510)100(=个点，其计算量非常之大。然而应用蒙特卡洛去随机计算610个点，便可

找到满意解，那么这种方法的可信度究竟怎样呢？

下面就分析随机取样采集6

10个点计算时，应用概率理论来估计一下可信度。 不失一般性，假定一个整数规划的最优点不是孤立的奇点。

假设目标函数落在高值区的概率分别为0.01，0.00001，则当计算6

10个点后，有任一个点能落在高值区的概率分别为

多位）100(9999.099.011000000Λ≈-，

999954602.099999.011000000≈-。

解 （i ）首先编写M 文件mente.m 定义目标函数f 和约束向量函数g ，程序如下： function [f,g]=mengte(x);

f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)-8*x(1)-2*x(2)-3*x(3)... -x(4)-2*x(5); g(1)=sum(x)-400;

g(2)=x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800; g(3)=2*x(1)+x(2)+6*x(3)-200; g(4)=x(3)+x(4)+5*x(5)-200;

（ii ）编写如下程序求问题的解：

rand('state',sum(clock)); p0=0; tic

for i=1:10^5

x=99*rand(5,1);

x1=floor(x);x2=ceil(x); [f,g]=mengte(x1); if sum(g<=0)==4 if p0<=f

x0=x1;p0=f; end end

[f,g]=mengte(x2); if sum(g<=0)==4 if p0<=f

x0=x2;p0=f; end end end

x0,p0 toc

§5 整数规划的计算机解法

整数规划问题的求解可以使用Lingo 等专用软件。对于一般的整数规划问题，无法直接利用Matlab 的函数，必须利用Matlab 编程实现分枝定界解法和割平面解法。但对于指派问题等10-整数规划问题，可以直接利用Matlab 的函数bintprog 进行求解。

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

???????

????

??

???109

6

10

9532485724679278310

283 解：编写Matlab 程序如下：

c=[3 8 2 10 3;8 7 2 9 7;6 4 2 7 5 8 4 2 3 5;9 10 6 9 10]; c=c(:);

a=zeros(10,25); for i=1:5

a(i,(i-1)*5+1:5*i)=1; a(5+i,i:5:25)=1; end

b=ones(10,1);

[x,y]=bintprog(c,[],[],a,b); x=reshape(x,[5,5]),y

求得最优指派方案为151********=====x x x x x ，最优值为21。

习 题 二

1. 用分枝定界法解：

21Max x x z +=

???

?

?

?

???≥≤+-≤+整数21212121,,0,3121451149x x x x x x x x 2. 试将下述非线性的10-规划问题转换成线性的10-规划问题 3211m ax x x x x z -+=

??

?==≤++-)3,2,1(103

32321j x x x x j

，或 3. 某钻井队要从以下10个可供选择的井位中确定5个钻井探油，使总的钻探费

用为最小。若10个井位的代号为1021,,,s s s Λ，相应的钻探费用为1021,,,c c c Λ，并且井位选择上要满足下列限制条件：

（1） 或选择1s 和7s ，或选择钻探9s ；

（2） 选择了3s 或4s 就不能选5s ，或反过来也一样；

（3） 在8765,,,s s s s 中最多只能选两个；试建立这个问题的整数规划模型。