实验二 线性规划的求解
学号:1 :何科
班级:2015级10班
一、 实验目的
1. 熟悉并掌握MATLAB的线性规划求解函数linprog()及其用法; 2. 熟悉并掌握LINGO软件求解线性规划的方法; 3. 能运用LINGO软件对线性规划问题进行灵敏度分析。
二、 实验任务
1. 对例1和例2,在MATLAB进行求解。
2. 对例3、4、5,在LINGO软件进行求解,并作灵敏度分析。 3. 对“3.3 投资的收益与风险”的模型I,在MATLAB中进行求解。 4. 对“习题5,6,7,8”进行建模与求解。
三、 实验过程与结果(对重要实验结果,截取全屏图,保存为JPG/PNG图片)
1. 例1: 代码:
f=[13 9 10 11 12 8]; A=[0.4 11 1 0 0 0; 0 0 0 0.5 1.2 1.3]; b=[800;900]; Aeq=[1 0 0 1 0 0; 0 1 0 0 1 0;
0 0 1 0 0 1]; beq=[400;600;500]; vlb=zeros(6,1); vub=[];
[x,fval]=linprog(f,A,b,Aeq,beq,vlb,vub)
结果:
x = 0.0000 600.0000 0.0000 400.0000 0.0000 500.0000 fval =1.3800e+04
例2:
代码: c=[40 36]; A=[-5 -3]; b=[-45]; Aeq=[]; beq=[]; vlb=zeros(2,1); vub=[9;15];
[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)
结果: x =
9.0000 0.0000
fval = 360
例3:
代码:
max=72*x1+64*x2; x1+x2<=50; 12*x1+8*x2<=480; 3*x1<=100;
结果:
Global optimal solution found.
Objective value: 3360.000 Infeasibilities: 0.000000 Total solver iterations: 2
Variable Value Reduced Cost X1 20.00000 0.000000 X2 30.00000 0.000000
Row Slack or Surplus Dual Price 1 3360.000 1.000000 2 0.000000 48.00000