Homepage Solution manuals Dimitris Bertsimas Introduction to Linear Optimization Exercise 2.3 (Basic feasible solutions in standard form polyhedra with upper bounds)

Exercise 2.3 (Basic feasible solutions in standard form polyhedra with upper bounds)

Consider a polyhedron defined by the constraints Ax = b and 0 x u, and assume that the matrix has linearly independent rows. Provide a procedure analogous to the one in Section 2.3 for constructing basic solutions, and prove an analog of Theorem 2.4.

Answers

We repeat Section 2.3 for a polyhedron with not only a lower bound 0 but also with an upper bound u:

P = {x nAx = b,0 x u}

for some matrix A m×n, vector u n and a vector b m.

Theorem 1. Consider the constraints Ax = b and 0 x u and assume that the m × n matrix A has linearly independent rows. A vector x n is a basic solution if and only if we have Ax = b, and there exist indices B(1),,B(m) such that

(a)
The columns AB(1),,AB(m) are linearly independent;
(b)
If iB(1),,B(m), then xi = 0 or xi = ui.

Proof. We demonstrate both directions.

  • Consider some x n : Ax = b and suppose that there are indices B(1),,B(m) that satisfy (a) and (b) in the statement of the theorem. Let I := {1,,n}{B(1),,B(m)} be the set of non-active indices, let U := I {i Ixi = ui} be the set of the indices where x hits the upper bound and let N := I {i Ixi = 0} be the set of the indices where x hits the lower bound. The active constraints and Ax = b imply that

    i=1mA B(i)xB(i) = i=1nA ixi0 iNAixi iUAixi = b iUAixi.

    Since the columns AB(i), i = 1,,m are linearly independent by the theorem hypothesis, by Theorem 1.2 (e,f), xB(1),,xB(m) are uniquely determined. In other words, the system of equations formed by the active constraints

    [ | | AB(1)AB(m) | | ] [xB(1) x B(m) ] = b iUAixi

    has a unique solution. By Theorem 2.2, this is equivalent to saying that there are n linearly independent active constraints. Coupled with the fact that all of the equality constraints are active by hypothesis, this implies that x is a basic solution.

  • For the converse, we assume that x is a basic solution of P = {x nAx = b,0 x u}; we will show that conditions (a) and (b) in the statement of the theorem are satisfied. Just as in the previous part, we partition the index into

    U := {1 i nxi = ui} N := {1 i nxi = 0}

    and denote the elements that are not indexed by either of these sets by xB(1),,xB(n) for some k 𝔑. We then have

    i=1kA B(i)xB(i) = i=1nA ixi0 iNAixi iUAixi = b iUAixi.

    Since x is a basic solution, there are n linearly independent constraints and by Theorem 2.2, the system of equations

    [ | | AB(1)AB(m) | | ] [xB(1) x B(m) ] = b iUAixi

    results in a unique solution. It follows that the columns AB(1),,AB(k) are linearly independent. [If they were not, we could find scalars λ1,,λk, no all of them zero, such that i=1kAB(i)λi = 0. This would imply that i=1kAB(i)(xB(i) + λi) = b iUAixi, contradicting the uniqueness of the solution.]
    We have shown that the columns AB(1),,AB(k) are linearly independent and this implies that k m. Since A has m linearly independent rows, it also has m linearly independent columns, which span m. It follows [cf. Theorem 1.3(b)] that we can find m k additional columns AB(k+1),,AB(m) so that the columns AB(i), i = 1,,m are linearly independent. In addition, if iB(1),,B(m), then iB(1),,B(k) (because k m), and we either have xi = 0 or xi = ui. Therefore, both conditions (a) and (b) in the statement of the theorem are satisfied.

In view of the above theorem, all basic solutions to a bounded form polynomial can be constructed according to the following procedure.

1.
Choose m linearly independent columns AB(1),,AB(m).
2.
Let xi = 0 or xi = ui for all iB(1),,B(m).
3.
Solve the system of m equations Ax = b for the unknowns xB(1),,xB(m).

If a basic solution constructed according to this procedure satisfies 0 x u, then it is feasible, and it is a basic feasible solution.

User profile picture
2021-11-08 18:51
Comments