Theorem 3.1, (15) — the RLS problem is the SOCP minimize λ s.t. ‖Ax − b‖ ≤ λ − τ, ‖[x; 1]‖ ≤ τ
ProvedRobustLS.Unstructured.rls_socp_formulationp2o-batch-p200ap2o-gran-per-chapterp2o-plan-paperp2o-v1robust-least-squaressecond-order-cone-programming
Let , and , and consider the second-order cone program in the variables , :
Then:
- for every , the worst-case residual (with ) is the smallest for which some makes feasible for (15);
- minimizes over if and only if there are such that is an optimal solution of (15).
So the robust least-squares problem and the SOCP (15) have the same optimal value and the same optimal .
Formalization Note Optimality in (15) is stated as: is feasible and for every feasible . The hypothesis is the one of Theorem 3.1.
Preamble
import Mathlib import Definitions.Def_RobustLS_Unstructured_Core open Matrix
Formal statement
namespace RobustLS.Unstructured
/-- El Ghaoui & Lebret (1997), Theorem 3.1, p. 1040 (PDF p. 6), SOCP (15): "This problem can be
formulated as the SOCP minimize λ subject to ‖Ax − b‖ ≤ λ − τ, ‖[x; 1]‖ ≤ τ." For `n ≥ 1` and
`ρ = 1`: (i) for every `x`, `r(A, b, x)` is the least `λ` for which some `τ` makes
`(x, λ, τ)` feasible for (15); (ii) `x` minimizes `r(A, b, ·)` over `ℝ^m` if and only if
`x` is the `x`-part of an optimal solution of (15). -/
theorem rls_socp_formulation {n m : ℕ} (hn : 0 < n) (A : Matrix (Fin n) (Fin m) ℝ)
(b : Fin n → ℝ) :
(∀ x : Fin m → ℝ,
IsLeast {lam : ℝ | ∃ τ : ℝ, SocpFeasible A b x lam τ} (worstCaseResidual A b 1 x)) ∧
(∀ x : Fin m → ℝ,
(∀ y : Fin m → ℝ, worstCaseResidual A b 1 x ≤ worstCaseResidual A b 1 y) ↔
∃ lam τ : ℝ, SocpFeasible A b x lam τ ∧
∀ (y : Fin m → ℝ) (lam' τ' : ℝ), SocpFeasible A b y lam' τ' → lam ≤ lam') := by sorry
end RobustLS.Unstructured
Source
El Ghaoui & Lebret, Robust Solutions to Least-Squares Problems with Uncertain Data, SIAM J. Matrix Anal. Appl. 18(4):1035–1064 (1997), p. 1040, Theorem 3.1, Eq. (15)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.