ํ•„ํ„ฐ ์ง€์šฐ๊ธฐ
ํ•„ํ„ฐ ์ง€์šฐ๊ธฐ

Dual problem and primal problem unbounded linear programming

์กฐํšŒ ์ˆ˜: 13 (์ตœ๊ทผ 30์ผ)
Klaus Hajdaraj
Klaus Hajdaraj 2021๋…„ 4์›” 14์ผ
Hello,
I need some help please to solve this problem I have . I have solved a max. problem , the primal problem , by using linprog.
After this , I have done the dual problem of the primal. After this proccess, I don't get the result that I should get , because the values of F0 should be the same in both cases(if I am not wrong).
I am sending the illustrated problem and code below. I hope someone can give me a help.
thank you !
Max 300๐‘ฅ1 + 400๐‘ฅ2
๐‘ฅ1 + ๐‘ฅ2 โ‰ค 320
4๐‘ฅ1 + 5๐‘ฅ2 โ‰ค 510
2๐‘ฅ1 + 3๐‘ฅ2 โ‰ค 430
3๐‘ฅ1 + ๐‘ฅ2 โ‰ค 300
๐‘ฅ1,๐‘ฅ2 โ‰ฅ 0
The Dual Problem
Min 320๐‘ฆ1 + 510๐‘ฆ2 + 430๐‘ฆ3 + 300๐‘ฆ4
๐‘ฆ1 + 4๐‘ฆ2 + 2๐‘ฆ3 + 3๐‘ฆ4 โ‰ฅ 300
4๐‘ฆ1 + 5๐‘ฆ2 + 3๐‘ฆ3 + ๐‘ฆ4 โ‰ฅ 400
๐‘ฆ1,๐‘ฆ2,๐‘ฆ3,๐‘ฆ4 โ‰ฅ 0
Matlab Code:
%Primary Problem
%Ax<=b
%coefficients of A
A = [1 4; 4 5; 2 3; 3 1];
%coefficients of B
b = [320 510 430 300];
%coefficients of the objective function
f = [300 400];
%The maximation of the linear function using the matlab linear
%programming function
[x0,F0] = linprog(-f,A,b,[],[],[0 0])
% disp('x0 = ')
% disp(x0);
%------------------------------------------------------------------
%Dual problem
%coefficients of AT
AT = [1 4 2 3;
4 5 3 1];
%coefficients of u
u = [300 400];
g = [320 510 430 300];
%setting a lower bounder
lb = zeros(1,4);
%setting an upper bounder
ub=[];
[u0,F0] = linprog(g,-AT,-u,[],[],lb,ub)
% disp('u0 = ');
% disp(u0);
%So,1st and 2nd constrains of the primary problem are active
%because u0 = [9.0909
% 72.7273
% 0
% 0]
%Solutions of Primary problem throug the dual problem
syms x1 x2
eqns = [x1 + x2 == 320, 4*x1 + 5*x2 == 510];
S = solve(eqns,[x1 x2]);
b = [S.x1 S.x2];
disp('Values of x ')
disp(b);
duality_gap = [x0(1)-u0(1); x0(2)-u0(2)]

๋‹ต๋ณ€ (0๊ฐœ)

์นดํ…Œ๊ณ ๋ฆฌ

Help Center ๋ฐ File Exchange์—์„œ Linear Programming and Mixed-Integer Linear Programming์— ๋Œ€ํ•ด ์ž์„ธํžˆ ์•Œ์•„๋ณด๊ธฐ

ํƒœ๊ทธ

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by