Lemma 4.3 — the stochastic transportation problem is a convex program
OpenShorNonsmooth.Decomposition.stochastic_transport_convexconvex-analysisp2o-batch-b23ap2o-gran-per-chapterp2o-plan-bookp2o-v1stochastic-programmingtransportation-problem
In the stochastic transportation problem (4.163)–(4.165), let the demands have probability densities with finite mean (), and let the penalty coefficients satisfy . Then the problem is a convex programming problem: the feasible set
is convex, and the objective
is a convex function on it.
Convexity is what licenses solving the problem through its Lagrangian dual (4.166) by a subgradient-type method, as the book does next.
Formalization Note The book leaves implicit that (a penalty coefficient) and that the demands have finite mean (otherwise the expectation is infinite); both are hypotheses. The independence of the plays no role in this statement and is not assumed.
Preamble
import Mathlib import Definitions.Def_ShorNonsmooth_Decomposition_StochasticTransport
Formal statement
namespace ShorNonsmooth.Decomposition
open MeasureTheory
/-- Shor (1985), **Lemma 4.3** (p. 132): the stochastic transportation problem (4.163)–(4.165) is a
convex programming problem — its feasible set is convex and its objective
`Σ c_ij x_ij + Σ r_j E(ξ_j − Σ_i x_ij)⁺` is convex on it. The demands `ξ_j` have probability densities
`p_j` with finite mean (so that the expectations are finite), and the penalty coefficients `r_j` are
nonnegative. -/
theorem stochastic_transport_convex {m n : ℕ}
(c : Fin m → Fin n → ℝ) (a : Fin m → ℝ) (r : Fin n → ℝ) (hr : ∀ j, 0 ≤ r j)
(p : Fin n → ℝ → ℝ) (hp : ∀ j, IsDensity (p j))
(hmean : ∀ j, Integrable (fun z => z * p j z)) :
Convex ℝ (stochTransportFeasible (n := n) a) ∧
ConvexOn ℝ (stochTransportFeasible a) (stochTransportObjective c r p) := by sorry
end ShorNonsmooth.Decomposition
Source
Shor, Minimization Methods for Non-Differentiable Functions, Springer 1985, p. 132, Lemma 4.3 (problem (4.163)–(4.165), pp. 131–132)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.