Why Does the Optimal Solution Occur at Corner Points?

In the realm of linear programming, a fascinating and crucial concept is that the optimal solution, if it exists, will occur at one of the corner points (or vertices) of the feasible region. To comprehend why this happens, we need to delve into some fundamental principles of linear programming, geometry, and optimization.

What is Linear Programming?

Linear programming (LP) is a method to achieve the best outcome in a mathematical model whose requirements are represented by linear relationships. It’s widely used in various fields such as economics, business, engineering, and military applications. The primary objective of linear programming is to maximize or minimize a linear objective function, subject to a set of linear inequalities or equations, known as constraints.

Example Problem

Consider a simple example where a company wants to maximize its profit. The profit function could be represented as:

$P = 3x + 2y$

where $x$ and $y$ are quantities of two different products, and 3 and 2 are the profit contributions per unit of each product, respectively. The company has certain constraints, such as limited resources, which can be represented by inequalities like:

$x + y leq 4$

$x geq 0$

$y geq 0$

$x leq 3$

Feasible Region

The feasible region is the set of all possible points that satisfy all the constraints. In our example, the feasible region is a polygon on the $xy$-plane, bounded by the lines representing the constraints. Each corner point of this polygon is a potential candidate for the optimal solution.

Why Corner Points?

Geometric Interpretation

The objective function, $P = 3x + 2y$, represents a family of parallel lines. As we move these lines up or down (in the case of maximization), we look for the highest line that still touches the feasible region. The key insight here is that this highest line will touch the feasible region at one of its vertices or along an edge. However, in the case of a unique optimal solution, it will be at a vertex.

The Fundamental Theorem of Linear Programming

The Fundamental Theorem of Linear Programming states that if there is an optimal solution to a linear programming problem, it will occur at a vertex (or corner point) of the feasible region. This theorem is pivotal because it reduces the problem of finding the optimal solution to checking the vertices of the feasible region, which is a finite set, making the problem more manageable.

Algebraic Justification

From an algebraic perspective, each vertex of the feasible region is a solution to a system of linear equations formed by intersecting the constraint lines. In other words, at each vertex, at least two of the constraints are active (i.e., they hold as equalities). By solving these systems of equations, we can find the coordinates of the vertices and then evaluate the objective function at these points to determine the optimal solution.

Example Continued

Let’s go back to our example. The constraints form a polygon with vertices at (0,0), (3,0), (0,4), and (3,1). We evaluate the objective function at these vertices:

At (0,0): $P = 3(0) + 2(0) = 0$

At (3,0): $P = 3(3) + 2(0) = 9$

At (0,4): $P = 3(0) + 2(4) = 8$

At (3,1): $P = 3(3) + 2(1) = 11$

The maximum profit is 11, which occurs at the vertex (3,1).

Special Cases

Multiple Optimal Solutions

Sometimes, the objective function may be parallel to a constraint boundary, resulting in multiple optimal solutions along an edge of the feasible region. In such cases, any point along this edge, including the vertices, will be optimal.

Unbounded Feasible Region

In some cases, the feasible region may be unbounded. If the objective function can increase indefinitely, there might be no optimal solution. However, if an optimal solution exists, it will still occur at a vertex of the feasible region.

Practical Implications

Understanding that the optimal solution occurs at corner points simplifies the process of solving linear programming problems. Methods like the Simplex algorithm leverage this property by systematically moving from one vertex to another, improving the objective function value at each step, until the optimal solution is reached.

Conclusion

The principle that the optimal solution occurs at corner points of the feasible region is a cornerstone of linear programming. It stems from both geometric and algebraic properties of linear functions and constraints. By focusing on the vertices of the feasible region, we can efficiently find the optimal solution to complex problems in various fields. This understanding not only aids in solving linear programming problems but also provides deeper insights into the nature of optimization itself.

Citations

  1. 1. Khan Academy – Linear Programming
  2. 2. MIT OpenCourseWare – Linear Programming
  3. 3. Wolfram Alpha – Linear Programming

Related

(2) O3 + H → O2 + OH k2 = 1.78×10^-11 cm^3 s^-1 (3) O + OH → O2 + H k3 = 4.40×10^-11 cm^3 s^-1 (5) O + HO2 → O2 + OH k5 = 3.50×10^-11 cm^3 s^-1 (6) H + HO2 → O2 + H2 k6 = 5.40×10^-12 cm^3 s^-1 (9) OH + HO2 → O2 + H2O2 k9 = 4.00×10^-11 cm^3 s^-1 (10) HO2 + HO2 → O2 + H2O2 k10 = 2.50×10^-12 cm s^-1 (11) O + O2 + M → O3 + M k11 = 1.05×10^-34 cm^6 s^-1 (14) H + O2 + M → HO2 + M k14 = 8.08×10^-32 cm^6 s^-1 (15) H + H + M → H2O + M k15 = 3.31×10^-27 cm^6 s^-1 (16) O2 + hv → 2 O k16 = (1.26×10^-8 s^-1) φ (17) H2O + hv → H + OH k17 = (3.4×10^-6 s^-1) φ (18) O3 + hv → O2 + O k18 = (7.10×10^-5 s^-1) φ

Table 1 Reactions, rate constants and activation energies used in the model* No. Reaction kopt (M⁻¹ s⁻¹) 1 OH + H₂ → H + H₂O 3.74 x 10⁷ 2 OH + HO₂ → HO₂ + OH⁻ 5 x 10⁹ 3 OH + H₂O₂ → HO₂ + H₂O 3.8 x 10⁷ 4 OH + O₂ → O₂ + OH 9.96 x 10⁹ 5 OH + HO₂ → O₂ + H₂O 7.1 x 10⁹ 6 OH + OH → H₂O₂ 5.3 x 10⁹ 7 OH + e⁻aq → OH⁻ 3 x 10¹⁰ 8 H + O₂ → HO₂ 2.0 x 10¹⁰ 9 H + HO₂ → H₂O₂ 2.0 x 10¹⁰ 10 H + H₂O₂ → OH + H₂O 3.44 x 10⁷ 11 H + OH → H₂O 1.4 x 10¹⁰ 12 H + H → H₂ 1.94 x 10¹⁰ 13 e⁻aq + O₂ → O₂⁻ 1.9 x 10¹⁰ 14 e⁻aq + O₂ → HO₂⁻ + OH⁻ 1.3 x 10¹⁰ 15 e⁻aq + HO₂ 2.0 x 10¹⁰ 16 e⁻aq + H₂O₂ 1.1 x 10¹⁰ 17 e⁻aq + HO₂ → OH + OH⁻ 1.3 x 10¹⁰ 18 e⁻aq + H⁺ → H 2.3 x 10¹⁰ 19 e⁻aq + e⁻aq → H₂ + OH⁻ + OH⁻ 2.5 x 10⁹ 20 HO₂ + O₂ → O₂ + HO₂ 1.3 x 10⁹ 21 HO₂ + HO₂ → O₂ + H₂O₂ 8.3 x 10⁵ 22 HO₂ + HO₂ → O₂ + OH + H₂O 3.7 23 HO₂ + HO₂ → O₂ + O₂ + OH + H₂O 7 x 10⁵ s⁻¹ 24 H⁺ + O₂⁻ → HO₂ 4.5 x 10¹⁰ 25 H⁺ + O₂⁻ → O₂ 2.0 x 10¹⁰ 26 H⁺ + OH⁻ 1.4 x 10¹¹ 27 H⁺ + HO₂⁻ 2 x 10¹⁰ 28 H₂O₂ → HO₂ + H⁺ + OH⁻ 2.5 x 10⁻⁵ s⁻¹ 29 H₂O₂ → H⁺ + OH⁻ 1.4 x 10⁻⁷ s⁻¹ 30 O₂ + O₂ → O₂ + HO₂ + OH⁻ 0.3 31 O₂ + H₂O₂ → O₂ + OH + OH 16 32

(2) O3 + H → O2 + OH k2 = 1.78×10^-11 cm^3 s^-1 (3) O + OH → O2 + H k3 = 4.40×10^-11 cm^3 s^-1 (5) O + HO2 → O2 + OH k5 = 3.50×10^-11 cm^3 s^-1 (6) H2O + O → 2 OH k6 = 5.40×10^-12 cm^3 s^-1 (9) OH + HO2 → O2 + H2O k9 = 4.00×10^-11 cm^3 s^-1 (10) HO2 + HO2 → O2 + H2O2 k10 = 2.50×10^-12 cm s^-1 (11) O + O2 + M → O3 + M k11 = 1.05×10^-34 cm^6 s^-1 (14) H + O2 + M → HO2 + M k14 = 8.08×10^-32 cm^6 s^-1 (15) OH + H + M → H2O + M k15 = 3.31×10^-27 cm^6 s^-1 (16) O2 + hv → 2 O k16 = (1.26×10^-8 s^-1) φ (17) H2O + hv → H + OH k17 = (3.4×10^-6 s^-1) φ (18) O3 + hv → O2 + O k18 = (7.10×10^-8 s^-1) φ