Linear Programming Lecture Notes - Penn State Personal Web Server
Linear Programming Lecture Notes - Penn State Personal Web Server
Linear Programming Lecture Notes - Penn State Personal Web Server
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
This leads to the following linear programming problem:<br />
⎧<br />
min x1 + 1.5x2<br />
⎪⎨ s.t. 2x1 + 3x2 ≥ 20<br />
(3.60)<br />
x1 + 2x2 ≥ 12<br />
⎪⎩<br />
x1, x2 ≥ 0<br />
Before turning to Matlab, let’s investigate this problem in standard form. To transform<br />
the problem to standard form, we introduce surplus variables s1 and s2 and our problem<br />
becomes:<br />
⎧<br />
min x1 + 1.5x2<br />
⎪⎨ s.t. 2x1 + 3x2 − s1 = 20<br />
(3.61)<br />
x1 + 2x2 − s2 = 12<br />
⎪⎩<br />
x1, x2, s1, s2 ≥ 0<br />
This leads to a set of two linear equations with four variables:<br />
2x1 + 3x2 − s1 = 20x1 + 2x2 − s2 = 12<br />
We can look at the various results from using Expression 3.51 and Definition 3.47. Let:<br />
<br />
2 3 −1 0 20<br />
(3.62) A =<br />
b =<br />
1 2 0 −1 12<br />
Then our vector of decision variables is:<br />
⎡ ⎤<br />
⎢<br />
(3.63) x = ⎢<br />
⎣<br />
x1<br />
x2<br />
s1<br />
s2<br />
⎥<br />
⎦<br />
We can use Gauss-Jordan elimination on the augmented matrix:<br />
<br />
2 3 −1 0 20<br />
1 2 0 −1 12<br />
Suppose we wish to have xB = [s1 s2] T . Then we would transform the matrix as:<br />
−2 −3 1 0 −20<br />
−1 −2 0 1 −12<br />
<br />
This would result in x1 = 0, x2 = 0 (because xN = [x1 x2] T and s1 = −20 and s2 = −12.<br />
Unfortunately, this is not a feasible solution to our linear programming problem because<br />
we require s1, s2 ≥ 0. Alternatively we could look at the case when xB = [x1 x2] T and<br />
xN = [s1 s2] T . Then we would perform Gauss-Jordan elimination on the augmented matrix<br />
to obtain: 1 0 −2 3 4<br />
0 1 1 −2 4<br />
<br />
That is, x1 = 4, x2 = 4 and of course s1 = s2 = 0. Notice something interesting:<br />
−1 −2 3<br />
1 −2<br />
= −<br />
2 3<br />
1 2<br />
48