Exercise 4.19 (Null variables)

Let P = {x nAx = b,x 0} be a nonempty polyhedron, and let m be the dimension of the vector b. We call xj a null variable if xj = 0 whenever x P.

(a)
Suppose that there exists some p m for which pA 0,pb = 0, and such that the j th component of pA is positive. Prove that xj is a null variable.
(b)
Prove the converse of (a): if xj is a null variable, then there exists some p m with the properties stated in part (a).
(c)
If xj is not a null variable, then by definition, there exists some y P for which yj > 0. Use the results in parts (a) and (b) to prove that there exist x P and p m such that:

pA 0,pb = 0,x + Ap > 0.

Answers

(a)

Proof. Consider the following pair of problems that are duals of each other:

minimize xj subject toAx b x 0, maximize pb subject topA ej,
(1)

where ej denotes the unit vector in jth direction. Let c denote pA’s jth direction; by assumption, c > 0. Define q := 1 cp. Then it is obvious that q is feasible for the dual problem as well, and the corresponding cost is qb = 1 cpb = 0. By weak duality (Theorem 4.3), the optimal cost of the primal problem cannot go lower than 0. But it cannot be positive either; thus, we conclude that xj will always be a null variable. □

(b)

Proof. Since xj is a null variable, the optimal cost of the primal problem in (1) will be 0. By strong duality (Theorem 4.4), the optimal cost of the dual problem is 0 as well. Let q denote the corresponding optimal solution. Then it is easy obvious that p := q satisfies pA ej0, (pA)j 1 > 0 and pb = 0. □

(c)

Proof. Let I denote the indices of non-null variables, and let N denote the indices of null variables.
By definition, for each j I, we can find a vector x(j) such that xj(j) > 0. By setting x := jIx(j), we ensure that x is positive at all nonnull variables.
Similarly, by part (a), for each j N, we can find a vector p(j) such that its jth component (p(j)A)j is positive. Setting p := jNp(j), we ensure that pA is strictly positive at every null variable.
Finally, since both x 0 and pA 0, adding x + pA results in a strictly positive vector. □

User profile picture
2022-03-21 00:20
Comments

c. I think it should be set x := 1 |I| jIx(j) to satisfy Ax = b.

User profile picture
2022-06-21 13:13
Comments