D1 5.4 - Two-stage Simplex and big-M methods

D1 5.4 - Two-stage Simplex and big-M methods

Lower-bound and equality constraints may need temporary artificial variables before a feasible basis is available. You will use two-stage Simplex and big-M for maximisation and minimisation, remove artificial variables correctly and interpret parameterised final tableaux (D1 5.4).

A feasible basis is a starting point, not the optimum

For x+y≥4x+y\ge4, x+2y≤6x+2y\le6, x,y≥0x,y\ge0, a first stage finds (x,y)=(4,0)(x,y)=(4,0) with every artificial variable zero. Has it already optimised 3x+2y3x+2y?

It has found a genuine feasible plan, but feasibility and optimisation answer different questions. This same plan will turn out to be neither the maximum nor the minimum of that expression.

You already know

Subtract surplus from a ≥\ge row and add an artificial variable if an initial basic column is needed; an equality may need an artificial variable without surplus. Make right-hand sides non-negative first. Original feasibility requires all artificials to vanish. Ordinary Simplex uses a canonical feasible basis and exact row operations.

For the opening constraints, write

x+y−s1+a1=4,x+2y+s2=6x+y-s_1+a_1=4,\qquad x+2y+s_2=6

Here s1s_1 is surplus, s2s_2 slack, and a1a_1 artificial. The computational starting basis is a1=4,s2=6a_1=4,s_2=6, with x=y=s1=0x=y=s_1=0. It is not an original feasible plan.

First-stage objective

In the two-stage convention used here, maximise W=−∑iaiW=-\sum_i a_i. Since each artificial variable is non-negative, W≤0W\le0. If its maximum is 00, every artificial variable is zero and an original feasible basis can be obtained. If its maximum is negative, the original constraints are infeasible.

For example, x≥3x\ge3, x≤1x\le1, x≥0x\ge0, gives x−s+a=3x-s+a=3. Since x≤1x\le1 and s≥0s\ge0, a=3−x+s≥2a=3-x+s\ge2. Thus the best W=−aW=-a is −2-2, attained at x=1,s=0x=1,s=0. No later optimisation can fix the incompatible original bounds.

An equality uses the same principle: x+y=4x+y=4 may start as x+y+a=4x+y+a=4. It has no surplus term, and reaching a=0a=0 restores the exact equality.

[DIAGRAM: asset_name: Two-stage Simplex separates feasibility from optimisation; asset_slug: 9fm0-d1-179-two-stage-flow; caption: A zero first-stage optimum permits the original objective to be restored; description: A left-to-right flow begins with adding artificial basic variables and then maximising W equal to minus the sum of artificial variables. A branch with maximum W below zero leads to original problem infeasible. A branch with maximum W zero leads to removing artificials and restoring the original objective, then to stage two optimising that objective. All arrows are full-opacity grey and the stage-two route is sage.]
A zero first-stage optimum permits the original objective to be restored

Key idea

Stage I asks whether the constraints can be satisfied without artificial help. Stage II optimises the original objective from a genuine feasible basis. A zero value of WW does not certify an original optimum.

Read the auxiliary objective as a feasibility certificate.

Find feasibility first, then restore the original objective

Use ordinary maximisation pivots on the auxiliary objective, then rebuild a canonical row for the original objective.

Two-stage Simplex

  1. Form non-negative slack, surplus and artificial variables with an initial feasible augmented basis.
  2. Set W=−∑aiW=-\sum a_i. Eliminate basic artificial variables from its row before pivoting.
  3. Maximise WW by Simplex. A negative optimum means the original problem is infeasible; an optimum zero permits Stage II.
  4. Remove artificials from the basis and delete their columns and the WW row. Restore the original objective, converting minimisation to maximisation if needed.
  5. Eliminate every current basic variable from that objective row, then continue Simplex and interpret the original result.

Not in the formula booklet: learn it.

An artificial variable can remain basic at value zero when Stage I ends. Pivot it out using a non-zero non-artificial non-basic coefficient in its zero-Value row. This keeps right-hand sides feasible. If that row has no such coefficient, it is redundant after artificials are set to zero and can be removed. For example, 2x+a=02x+a=0 permits a pivot on 22 to make x=0x=0 basic; a row consisting only of a=0a=0 is removed. Do not delete a non-redundant constraint merely because its artificial value is zero.

Worked example 1: a maximisation after feasibility

Question. Use two-stage Simplex to maximise P=3x+2yP=3x+2y, subject to x+y≥4x+y\ge4, x+2y≤6x+2y\le6, x,y≥0x,y\ge0.

Plan. First maximise W=−a1W=-a_1, then remove a1a_1 and restore a canonical row for PP.

The augmented constraints are those above. From W+a1=0W+a_1=0, subtract the first constraint to eliminate basic a1a_1:

W−x−y+s1=−4W-x-y+s_1=-4
b.v.xxyys1s_1s2s_2a1a_1Value
a1a_11111−1-1001144
s2s_2112200110066
WW−1-1−1-1110000−4-4

Choose xx from the tied entering columns. Ratios 4,64,6 select row 11; a1a_1 leaves. Use R2←R2−R1R_2\leftarrow R_2-R_1, RW←RW+R1R_W\leftarrow R_W+R_1:

b.v.xxyys1s_1s2s_2a1a_1Value
xx1111−1-1001144
s2s_200111111−1-122
WW000000001100

Stage I has optimum W=0W=0, with non-basic a1=0a_1=0. Delete its column and the WW row. Restore P−3x−2y=0P-3x-2y=0; add 3R13R_1 to eliminate basic xx:

b.v.xxyys1s_1s2s_2Value
xx1111−1-10044
s2s_20011111122
PP0011−3-3001212

Enter surplus s1s_1: only row 22 has a positive entry, giving ratio 22. Use R1←R1+R2R_1\leftarrow R_1+R_2, RP←RP+3R2R_P\leftarrow R_P+3R_2, and relabel row 22:

b.v.xxyys1s_1s2s_2Value
xx1122001166
s1s_10011111122
PP004400331818

Why this step? Restoring the raw objective is insufficient: the current basic xx must have zero coefficient in its row. A surplus variable can enter just like any other non-basic variable.

Answer. Set y=s2=0y=s_2=0: x=6x=6, maximum P=18P=18.

Check. The first constraint has surplus 22, the second is tight, and P=18−4y−3s2≤18P=18-4y-3s_2\le18. The Stage I point had objective 1212, so finding feasibility did not finish maximisation.

Worked example 2 changes only the required direction of optimisation of 3x+2y3x+2y, from maximum to minimum. Stage I is unchanged; Stage II must use a different objective.

Worked example 2: the same feasible start, a minimum instead

Question. Use two-stage Simplex to minimise C=3x+2yC=3x+2y under the same constraints.

Plan. Repeat the auxiliary feasibility stage, then maximise P=−CP=-C.

The first-stage initial tableau is

b.v.xxyys1s_1s2s_2a1a_1Value
a1a_11111−1-1001144
s2s_2112200110066
WW−1-1−1-1110000−4-4

Enter xx, using ratios 4,64,6, and apply R2←R2−R1R_2\leftarrow R_2-R_1, RW←RW+R1R_W\leftarrow R_W+R_1:

b.v.xxyys1s_1s2s_2a1a_1Value
xx1111−1-1001144
s2s_200111111−1-122
WW000000001100

Remove the non-basic artificial column and WW. The original converted row is P+3x+2y=0P+3x+2y=0. Apply RP←RP−3R1R_P\leftarrow R_P-3R_1:

b.v.xxyys1s_1s2s_2Value
xx1111−1-10044
s2s_20011111122
PP00−1-13300−12-12

Enter yy. Ratios 4,24,2 select row 22. Use R1←R1−R2R_1\leftarrow R_1-R_2, RP←RP+R2R_P\leftarrow R_P+R_2:

b.v.xxyys1s_1s2s_2Value
xx1100−2-2−1-122
yy0011111122
PP00004411−10-10

Why this step? The auxiliary objective does not depend on whether the original question asks for a maximum or minimum. The restored objective does, so its signs must be changed before Stage II.

Answer. x=y=2x=y=2, maximum P=−10P=-10, hence minimum C=10C=10.

Check. Both original constraints hold with equality. Since P=−10−4s1−s2≤−10P=-10-4s_1-s_2\le-10, the original cost satisfies C≥10C\ge10.

An optimum in a feasible Stage II tableau follows the ordinary Simplex tests, including the checks for unbounded directions or alternative optimal solutions. The extra work was obtaining a genuine basis and the correct restored row.

Watch out

Wrong: stop as soon as W=0W=0, or replace the auxiliary row by an unreduced original row containing basic-variable coefficients. Right: clean up the artificial basis, restore the original objective and eliminate every current basic variable before continuing.

Change only the original cost's xx coefficient in the next check.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account