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 ∈ ℜn∣Ax = 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 i≠B(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 ∈ I∣xi = ui} be the set of the indices where x hits the upper bound and let N := I ∖{i ∈ I∣xi = 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 ixi−0∑ i∈NAixi−∑ i∈UAixi = b−∑ i∈UAixi.

    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−∑ i∈UAixi

    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 ∈ ℜn∣Ax = 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 ≤ n∣xi = ui} N := {1 ≤ i ≤ n∣xi = 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 ixi−0∑ i∈NAixi−∑ i∈UAixi = b−∑ i∈UAixi.

    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−∑ i∈UAixi

    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 −∑ i∈UAixi, 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 i≠B(1),…,B(m), then i≠B(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 i≠B(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