0-1整数规划.pptVIP

第一章 Matlab入门 第12章 0-1整数规划 1 预备知识：整数线性规划基本理论 2 整数线性规划MATLAB程序 3 0-1型整数线性规划 4 0-1型线性规划Matlab命令 5 建模与计算实验 12.1 整数线性规划基本理论 12.2 整数线性规划MATLAB程序 Matlab没有整数线性规划求解命令 算法：分支定界法 工具箱yalmip(可网上下载） 可求解全整数线性或混合整数线性规划 补充工具箱应用 yalmip用法 (1)定义变量： sqdvar（）实型；intvar（）整型；binvar（）0-1型； (2)设定目标函数: f=目标函数； (3)设定限定条件：F=set（限定条件）； 多个限定条件用加号相连：F=set（限定条件）+set（限定条件1）+set（限定条件2）……； (4) 求解： solvesdp（F,f）； 这里解得是F条件下目标函数f的最小值，要求最大值f前面加个负号。 求解之后查看数值: double（f） double（变量） Yalmip工具箱应用举例 目标函数：max z=4x1+6x2+2x3 ??? s.t.??????-x1+3x2=8 ???????????????-x2+3x3=10 ???????????????5x1-x3=8 ???????????????x1、x2、x3=0,x1、x2、x3为整数 ? 具体代码如下： ? x=intvar(1,3); ? f=[4 6 2]*x’; ? F=set(x0)+set([-1 3 0]*x’=8)+set([0 -1 3]*x’=10)+set([5 0 -1]*x’=8); ? solvesdp(F,-f) ? double(f) ? double(x) 非线性整数规划举例 Max z=x1^2+x2^2+3*x3^2+4*x4^2+2*x5^2-8*x1-2*x2-3*x3-x4-2*x5 s.t. x1+x2+x3+x4+x5=400 x1+2*x2+2*x3+x4+6*x5=800 2*x1+x2+6*x3=800 x3+x4+5*x5=200 0=xi=99 (i=1,2,...,5) xi为整数 在matlab中输入 x=intvar(1，5)； f=[1 1 3 4 2]*(x.^2)-[8 2 3 1 2]*x;F=set(0=x=99); F=F+set([1 1 1 1 1]*x=400)+set([1 2 2 1 6]*x =800)+set(2*x(1)+x(2)+6*x(3)=800); F=F+set(x(3)+x(4)+5*x(5)=200); solvesdp(F,-f) double(f) 80199 double(x) 53 99 99 99 0 12.3 0-1型整数线性规划 0-1型整数线性规划 12.4 0-1型线性规划Matlab命令 Matlab命令 x=bintprog(f, A, b)求解0-1型整数线性规划，用法类似于linprog。 x= bintprog (f, A, b, Aeq, beq)求解下面线性规划： min z = f ’?x, A?x ≤ b, Aeq?x = beq，x分量取值0或1。 x= bintprog(f, A, b, Aeq, beq, x0)指定迭代初值x0 如果没有不等式约束，可用[ ]替代A和b表示缺省，如果没有等式约束，可用[ ]替代Aeq和beq表示缺省；用[x, Fval]代替上述各命令行中左边的x，则可得到最优解处的函数值Fval。 例3 例3 求解下列0-1型整数规划 求解 首先化为标准形式， 转换求max为求min –f= 3x1-2x2+5x3 c=[3,-2,5]; a=[1,2,-1;1,4,1;1,1,0;0,4,1];b=[2;4;3;6]; [x,fval]=bintprog(c,a,b); xmax=x,zmax=-fval 12.5 建模与计算实验 例5 某厂拟用集装箱托运甲乙两种货物，每箱的体积、重量、可获利润以及托运所受限制如下表： 建模 设x1和x2分别为甲和乙两种货物的托运箱数 求解 化为min问题,用yalmip求解 例6 一架货运飞机，有效载重量为24吨，可运输物品的重量及运费收入如下表所示，其中各物品只有一件可供选择。 建模 求解 化为min问题,用bintprog求解 c=-[3,5,2,4,2,3];a=[8,13,6,9,5,7];b=24; [x,g]= bintprog(c,a,b) f=-g * 整数线性规划的提法： 在一般的线性规划中，