整数规划Matlab实现

星辰之光.

已于 2022-05-19 20:05:17 修改

阅读量4.9k

点赞数 14

分类专栏： matlab实例文章标签：算法开发语言 matlab

于 2022-05-19 09:45:12 首次发布

本文链接：https://blog.csdn.net/m0_64548999/article/details/124855644

版权

matlab实例专栏收录该内容

11 篇文章 4 订阅

订阅专栏

$\max z=40x_1+90x_2$

$9x_1+7x_2\leq 56$

$7x_2+20x_2\leq 70$

$x_1,x_2\geq 0$ 且为整数

现实中很多问题都需要求整，下面我们用Matlab来解决这个问题，我们这里使用分支定界法。

先用线性规划解出大致范围。

%初步处理
clear all
clc
c=[40 90];%填目标函数系数
a=[9,7;7,20]; %就是要一致大于或小于，然后取条件函数系数
b=[56;70];%条件函数已知数部分
aeq=[];%若存在一个等式限制，则填写等式系数
beq=[];%等式已知数部分
lb=[0;0];%没有下限
ub=[inf;inf];%没有上限
[x,y]=linprog(-c,a,b,aeq,beq,lb,ub);
x
best=c*x%计算最优值

这是我们利用线性规划解出的。

分枝定界法

第一步

根据上一步结果，可以暂定z上限为 356。同样也可以看出当 $x_1$ , $x_2$ 分别为0时z最小，最小为0。故最大值z范围可定位 $0\leq z\leq 356$ 。

第二步

先对 $x_1$ 分枝

$x_1\leqslant [4.8092]= 4,x_1\geq [4.8092]=5$

4和5之间就没有整数，该方法叫做分枝

约束条件转变

$0\leq x_1<4$ $x_2>0$

%整数规划
clear all
clc
c=[40 90];%填目标函数系数
a=[9,7;7,20]; %就是要一致大于或小于，然后取条件函数系数
b=[56;70];%条件函数已知数部分
aeq=[];%若存在一个等式限制，则填写等式系数,没有等式约束，因此aeq,beq都为空
beq=[];%等式已知数部分
lb=[0;0];%下限依然都为0
ub=[4;inf];%x1上限为4，x2没有上限
[x,y]=linprog(-c,a,b,aeq,beq,lb,ub);
x    %获取对应x1,x2
best=c*x%计算最优值

运行结果：

我们得出 $x_1=4$ , $x_2=2.1$ , $z_1=349$

换分割的x范围，

$x_1\geq 5$ $x_2>0$

%整数规划
clear all
clc
c=[40 90];%填目标函数系数
a=[9,7;7,20]; %就是要一致大于或小于，然后取条件函数系数
b=[56;70];%条件函数已知数部分
aeq=[];%若存在一个等式限制，则填写等式系数,没有等式约束，因此aeq,beq都为空
beq=[];%等式已知数部分
lb=[5;0];
ub=[inf;inf];
[x,y]=linprog(-c,a,b,aeq,beq,lb,ub);
x    %获取对应x1,x2
best=c*x%计算最优值

运行结果：

我们得出 $x_1=5$ , $x_2=1.5714$ , $z=341.4286$

现在我们可以把最优解再定界，确定范围为： $0\leq z\leq 349$

第三步

对 $x_2$ 处理，这里不细写。

$0\leq x_1\leq 4$ $0\leq x_2\leq 2$

$0\leq x_1\leq 4$ $x_2 \geq 3$

因此我们可以再次确定z范围为 $340\leq z\leq 341$

第四步

对第二种情况分析

$x_1\geq 5$ $0\leq x_2\leq 2$

$x_1\geq 5$ $x_2 \geq 3$

可以看到，没有解

综合上面两个分枝，第一个仍有小数，不可取；第二个不可行解，也不可取。这样的操作称为剪枝。

结论

第二种情况，两分枝都不可取，第一种只有一个枝可取，即最优解为：

$x_1=4$ , $x_2=2$ , $z=340$

总结

开始，将要求解的整数规划问题称为问题A,将与它相应的线性规划问题称为问题B。求解情况如下:
1. B没有可行解，这时A也没有可行解，则停止
2. B 有最优解，并符合问题A的整数条件，B 的最优解即为A的最优解，则停止。
3. B有最优解，但不符合问题A的整数条件,记它的目标函数值为z,通过我上述所说的四个步骤使用分枝定界法进行分枝，剪枝,最后得出结果。

这里留一道题，你们可以练练手：