Theorem 6.7, p. 27 — min P_nmq ≤ U_nmq, with equality when P_nmq has an optimal solution with at most m non-zero entries
DisprovedSymPolyOpt.PowerSumUB.theorem_6_7Let with and , and let . Let be the optimal value of
and let be the value of the one-variable SDP (6.7): the infimum of over the monic real polynomials of degree with Newton sums for and (with ). Then:
- ;
- if has an optimal solution with at most non-zero entries, then
so that has the equivalent convex formulation (6.7).
The theorem complements the Hankel-matrix lower bound of Theorem 6.6 by an upper bound computed from an SDP in the single variable , whose size depends on and not on .
Formalization Note Both values are infima in the extended reals (EReal), so no optimal solution of (6.7) is assumed to exist. The hypothesis is the standing assumption of (6.4); with it forces . "An optimal solution" is a feasible point whose objective value equals ; "at most non-zero entries" counts the indices with .
import Mathlib import Definitions.Def_SymPolyOpt_PowerSumUB_Setting
namespace SymPolyOpt.PowerSumUB
theorem theorem_6_7 (n m q : ℕ) (γ : ℕ → ℝ) (hmn : m ≤ n) (hmq : m ≤ q)
(hq : q ≤ 2 * m - 2) :
SymPolyOpt.PowerSumLB.minP n m q γ ≤ valU m q γ ∧
((∃ x ∈ SymPolyOpt.PowerSumLB.feasP n m γ, ((SymPolyOpt.PowerSumLB.powerSum x q : ℝ) : EReal) = SymPolyOpt.PowerSumLB.minP n m q γ ∧
(Finset.univ.filter (fun i => x i ≠ 0)).card ≤ m) →
SymPolyOpt.PowerSumLB.minP n m q γ = valU m q γ) := by sorry
end SymPolyOpt.PowerSumUB
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.