Distributionally robust chance constraint worst-case VaR bound, Eq. (1)
ProvedDRCVRP.RCI.chance_constraint_iff_worstCaseVaR_leLet be a nonempty set of probability distributions of the demand vector , let , , and let be a route, a list of distinct customers. Assume the values , , are bounded above. Then
where means .
A route satisfies its distributionally robust chance constraint exactly when the worst-case value-at-risk of its cumulative demand fits into one vehicle; this is what motivates the demand estimator (2).
Formalization Note The statement is the conjunction of the first-and-second and the first-and-third equivalences. The supremum is a real sSup; nonemptiness and boundedness above of the VaR values make it the paper's finite supremum (for an empty the paper's supremum is , which a real sSup cannot represent). The route is a List (Fin n) without duplicates, and the sum over the route is the sum over its underlying finite set.
import Mathlib import Definitions.Def_MultistageStochastic_RiskFunctional import Definitions.Def_DRCVRP_RCI_DemandEstimator import Definitions.Def_DRCVRP_RCI_Formulations open MeasureTheory
namespace DRCVRP.RCI
/-- Eq. (1), §3, p. 721: for a route `R_k` (a list of distinct customers) and an ambiguity set
`𝒫` (`Amb`) of probability distributions of the demand vector,
`ℙ[R_k ∈ ℛ(q̃)] ≥ 1 - ε ∀ ℙ ∈ 𝒫 ⟺ ℙ-VaR_{1-ε}[∑_{i ∈ R_k} q̃_i] ≤ Q ∀ ℙ ∈ 𝒫
⟺ sup_{ℙ ∈ 𝒫} ℙ-VaR_{1-ε}[∑_{i ∈ R_k} q̃_i] ≤ Q`.
The supremum is a real `sSup`; `Amb` nonempty and the VaR values bounded above make it the
paper's (finite) supremum. -/
theorem chance_constraint_iff_worstCaseVaR_le {n : ℕ} (Amb : Set (Measure (Fin n → ℝ)))
(hAmb : ∀ P ∈ Amb, IsProbabilityMeasure P) (hne : Amb.Nonempty)
(ε : ℝ) (hε0 : 0 < ε) (hε1 : ε < 1) (Q : ℝ) (r : List (Fin n)) (hr : r.Nodup)
(hbdd : BddAbove ((fun P => MultistageStochastic.valueAtRisk P
(fun q => ∑ i ∈ r.toFinset, q i) (1 - ε)) '' Amb)) :
(RouteChanceFeasible Amb ε Q r ↔
∀ P ∈ Amb, MultistageStochastic.valueAtRisk P (fun q => ∑ i ∈ r.toFinset, q i) (1 - ε) ≤ Q) ∧
(RouteChanceFeasible Amb ε Q r ↔ worstCaseVaR Amb ε r.toFinset ≤ Q) := by sorry
end DRCVRP.RCI
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.