整数规划Matlab.pptVIP

* 例：用分枝定界法求解整数规划问题（用图解法计算） 记为（IP） 解：首先去掉整数约束，变成一般线性规划问题(原整数规划问题的松驰问题) 记为（LP） f=[-1 -5]; a=[-1 1;5 6;1 0]; b=[2 30 4]; lb=[0 0]; [x, fv]=linprog(f,a,b,[],[],lb) % x1 =1.6364, x2=3.6364, fv = -19.8182 　　　　　　　　　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　 function [a1,a2,b1,b2]=intp(k,n,a,b,p) 　 % k为分支变量的编号,n为矩阵a的列数,p是要分支变量值　 d1=zeros(1,n); d1(k)=1; d2=d1; d2(k)=-1; a1=[a;d1] a2=[a;d2] c=fix(x) b1=[b;c] b2=[b;-1-c] 对于x1＝18/11≈1.64， 取值x1 ≤1， x1 ≥2 将（LP）划分为（LP1） 和（LP2）. 用图解法求（LP）的最优解，如图所示。 x1 x2 ⑴ ⑵ 3 (18/11,40/11) ⑶ x1＝18/11, x2 =40/11 Z(0) =－218/11≈(－19.8) 2 1 1 2 3 有下式： 现在只要求出（P1）和（P2）的最优解即可。 p=x(1) %取x1分枝,x1≤1，x1≥2 [a1,a2,b1,b2]=intp(1,2,a,b,p) %形成分支后的约束矩阵. [x1,fv1,e1]=linprog(f,a1,b1,[],[],lb) %x1 =(1, 3), fv1 = -16, e1 = 1,整数解不再分支. [x2,fv2,e2]=linprog(f,a2,b2,[],[],lb) % x2 = (2,10/3), fv2 = -18.6667-16, e2 =1,继续分支 x1 x2 ⑴ ⑵ 3 3 (18/11,40/11) ⑶ 先求（P1）,如图所示。此时在B点取得最优解。 x1＝1, x2 =3, Z(1)＝－16(上界) 找到整数解，问题已探明，此枝停止计算。 1 1 同理求（P2） ,如图所示。 在C 点取得最优解。 即x1＝2, x2 =10/3, Z(2) ＝－56/3≈－18.7 ∵Z2 Z1＝－16 ∴原问题有比（－16）更小的最优解，但 x2 不是整数，故继续分支。 B A C －１６． 加入条件： x2≤3, x2≥4 有下式： 只要求出（P3）和（P4）的最优解即可。 p=x2(2) % [a3,a4,b3,b4]=intp(2,2,a2,b2,p) %形成分支后的约束矩阵 [x3,fv3,e3]=linprog(f,a3,b3,[],[],lb) %x3 =(2.4, 3), fv3 = -17.4-16, e3=1,继续分支. [x4,fv4,e4]=linprog(f,a4,b4,[],[],lb) % e40, 无可行解 x1 x2 ⑴ ⑵ 3 3 (18/11,40/11) ⑶ 1 1 B A C 先求（LP3）,如图所示。此时D 在点取得最优解。 即 x1＝12/5≈2.4, x2 =3, Z(21)＝-87/5≈-17.4 Z(1)=-16 但x1＝12/5不是整数，可继续分枝。即 2≤x1≤3。 D 求（LP4），如图所示。 无可行解，不再分枝。 在（P3）的基础上继续分枝。加入条件2≤x1≤3有下式： 只要求出（P5）和（P6）的最优解即可。 且为整数 且为整数