对偶单纯形法 - 360文档中心

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

1.对偶单纯形法

2.F(x)=3x1+4x2+5x3

X1+2x2+3x3>=5

2x1+2x2+x3>=6

xi>=0

f=[3;4;5];

A=[-1 -2 -3

-2 -2 -1];

b=[-5;-6];

lb=zeros(3,1);

[x,fval,exitflag,output,lambda] = linprog(f,A,b,[],[],lb) x =

1.0000

2.0000

0.0000

fval =

11.0000

exitflag =

1

output =

iterations: 8

algorithm: 'large-scale: interior point' cgiterations: 0

message: 'Optimization terminated.' lambda =

ineqlin: [2x1 double]

eqlin: [0x1 double]

upper: [3x1 double]

lower: [3x1 double]

x, lambda.ineqlin, lambda.lower

x =

1.0000

2.0000

0.0000

ans =

1.0000

ans =

0.0000

1.0000

2.单纯形法

f(x)=-9x1+-16x2

x1+4x2+x3=80

2x1+3x2+x4=90

xi>=0

f=[-9;-16;0;0];

Aeq=[1 4 1 0

2 3 0 1];

beq=[80;90];

lb=zeros(4,1);

[x,fval,exitflag,output,lambda] = linprog(f,[],[],Aeq,beq,lb) Optimization terminated.

x =

24.0000

14.0000

0.0000

fval =

-440.0000

exitflag =

1

output =

iterations: 5

algorithm: 'large-scale: interior point'

cgiterations: 0

message: 'Optimization terminated.' lambda =

ineqlin: [0x1 double]

eqlin: [2x1 double]

upper: [4x1 double]

lower: [4x1 double]

x, lambda.ineqlin, lambda.lower

x =

24.0000

14.0000

0.0000

ans =

Empty matrix: 0-by-1

ans =

0.0000

1.0000 4.0000