Documentation

Books.ConvexAnalysis_Rockafellar_1970.Chapters.Chap03.section12_part9

theorem quadraticHalfLinear_translate_sub_of_dotProduct_symmetric {n : } (Q : (Fin n) →ₗ[] Fin n) (hQsymm : ∀ (x y : Fin n), Q x ⬝ᵥ y = x ⬝ᵥ Q y) (a x : Fin n) :
quadraticHalfLinear n Q (x - a) = quadraticHalfLinear n Q x + ↑(x ⬝ᵥ -Q a) + ↑(1 / 2 * a ⬝ᵥ Q a)

Translation identity for the shifted argument: rewrite quadraticHalfLinear n Q (x - a) in terms of quadraticHalfLinear n Q x, plus a linear term and a constant, assuming dot-product symmetry of Q.

theorem sq_segment_le (a b t : ) (ht0 : 0 t) (ht1 : t 1) :
((1 - t) * a + t * b) ^ 2 (1 - t) * a ^ 2 + t * b ^ 2

Convexity of the square function: t ↦ ((1 - t)a + t b)^2 lies below the chord between a^2 and b^2 for t ∈ [0,1].

theorem properConvexFunctionOn_diagonalQuadratic {n : } (d : Fin n) (hd : ∀ (i : Fin n), 0 d i) :
ProperConvexFunctionOn Set.univ fun (x : Fin n) => ↑(1 / 2 * i : Fin n, d i * x i ^ 2)

A weighted diagonal quadratic function with nonnegative weights is a proper convex function on ℝ^n.

Elementary partial quadratic convex functions are proper convex on Set.univ.

theorem properConvexFunctionOn_precomp_linearEquiv {n : } (A : (Fin n) ≃ₗ[] Fin n) {f : (Fin n)EReal} (hf : ProperConvexFunctionOn Set.univ f) :
ProperConvexFunctionOn Set.univ fun (x : Fin n) => f (A x)

Precomposition by a linear equivalence preserves proper convexity on Set.univ.

theorem properConvexFunctionOn_translate {n : } (a : Fin n) {f : (Fin n)EReal} (hf : ProperConvexFunctionOn Set.univ f) :
ProperConvexFunctionOn Set.univ fun (x : Fin n) => f (x - a)

Translating the input preserves proper convexity on Set.univ.

theorem properConvexFunctionOn_affine_dotProduct {n : } (v : Fin n) (α : ) :
ProperConvexFunctionOn Set.univ fun (x : Fin n) => ↑(x ⬝ᵥ v + α)

The affine function x ↦ ⟪x, v⟫ + α is a proper convex function on Set.univ.

theorem isPartialQuadraticConvexFunction_iff_exists_elementary_affine (n : ) (f : (Fin n)EReal) :
IsPartialQuadraticConvexFunction n f ∃ (h : (Fin n)EReal) (A : (Fin n) ≃ₗ[] Fin n) (a : Fin n) (aStar : Fin n) (α : ), IsElementaryPartialQuadraticConvexFunction n h f = fun (x : Fin n) => h (A (x - a)) + ↑(x ⬝ᵥ aStar) + α

Text 12.3.3: Characterization of partial quadratic convex functions by affine change of coordinates to an elementary partial quadratic convex function plus a linear term.

theorem ereal_iSup_neg_eq_neg_iInf {ι : Sort u_1} (g : ιEReal) :
⨆ (i : ι), -g i = -⨅ (i : ι), g i

Negation turns iInf into iSup in EReal.

theorem fenchelConjugate_zero_eq_neg_iInf (n : ) (f : (Fin n)EReal) :
fenchelConjugate n f 0 = -⨅ (x : Fin n), f x

The value of the Fenchel conjugate at 0 is the negative of the infimum of f.

theorem iInf_fenchelConjugate_eq_neg_eval_zero (n : ) (f : (Fin n)EReal) (hf_closed : LowerSemicontinuous f) (hf_convex : ConvexFunction f) (hf_proper : ProperConvexFunctionOn Set.univ f) :
⨅ (xStar : Fin n), fenchelConjugate n f xStar = -f 0

For closed proper convex f, the infimum of the Fenchel conjugate is -f 0.

theorem inf_eq_at_zero_iff_conjugate_inf_eq_at_zero (n : ) (f : (Fin n)EReal) (hf_closed : LowerSemicontinuous f) (hf_convex : ConvexFunction f) (hf_proper : ProperConvexFunctionOn Set.univ f) :
⨅ (x : Fin n), f x = f 0 f 0 = 0 ⨅ (xStar : Fin n), fenchelConjugate n f xStar = fenchelConjugate n f 0 fenchelConjugate n f 0 = 0

Text 12.3.4: Let f be a closed proper convex function on ℝ^n. Then inf_x f x = f 0 = 0 if and only if inf_{x*} f* x* = f* 0 = 0, where f* is the Fenchel conjugate (here f* = fenchelConjugate n f).

@[reducible, inline]
abbrev IsSymmetricFunction {α : Type u_1} (n : ) (f : (Fin n)α) :

Defn 12.3: A function f on ℝ^n is symmetric (or even) if it satisfies f (-x) = f x for all x ∈ ℝ^n.

Equations
Instances For
    theorem iSup_comp_neg {α : Type u_1} {β : Type u_2} [AddGroup α] [CompleteLattice β] (g : αβ) :
    ⨆ (x : α), g (-x) = iSup g

    Reindex an iSup by the involution x ↦ -x.

    If f is even, then its Fenchel conjugate f* is even.

    theorem even_clConv_of_even_fenchelConjugate (n : ) {f : (Fin n)EReal} (hf_convex : ConvexFunction f) (hEvenConj : Function.Even (fenchelConjugate n f)) :

    If the Fenchel conjugate of f is even, then the closed convex envelope clConv n f is even.

    theorem even_of_even_fenchelConjugate_of_closedProperConvex (n : ) {f : (Fin n)EReal} (hf_closed : LowerSemicontinuous f) (hf_convex : ConvexFunction f) (hf_proper : ProperConvexFunctionOn Set.univ f) (hEvenConj : Function.Even (fenchelConjugate n f)) :

    For closed proper convex f, evenness of its Fenchel conjugate implies evenness of f.

    Text 12.3.5: Let f be a closed proper convex function on ℝ^n. Then f is symmetric (i.e. f (-x) = f x for all x) if and only if its conjugate f* is symmetric (here f* = fenchelConjugate n f).

    def IsRotationallySymmetric {α : Type u_1} (n : ) (f : (Fin n)α) :

    Text 12.3.5: A function f on ℝ^n is symmetric with respect to all orthogonal transformations (rotationally symmetric) if and only if it has the form f x = g (‖x‖), for some function g on [0, +∞), where ‖·‖ is the Euclidean norm.

    Furthermore, for such an f, it is a closed proper convex function on ℝ^n if and only if g satisfies:

    1. g is convex on [0, +∞);
    2. g is non-decreasing;
    3. g is lower semicontinuous;
    4. g 0 is finite.
    Equations
    • One or more equations did not get rendered due to their size.
    Instances For
      noncomputable def euclidNorm (n : ) (x : Fin n) :

      The Euclidean norm on Fin n → ℝ, defined via the equivalence with EuclideanSpace.

      Equations
      Instances For
        def nonnegRay :
        Set (Fin 1)

        The nonnegative ray in Fin 1 → ℝ used to model [0, +∞).

        Equations
        Instances For

          The nonnegative ray nonnegRay is convex.

          In Euclidean space, vectors of equal norm are related by some linear isometry.

          theorem isRotationallySymmetric_iff_exists_norm_comp {α : Type u_1} (n : ) (f : (Fin n)α) :
          IsRotationallySymmetric n f ∃ (g : (Fin 1)α), ∀ (x : Fin n), f x = g fun (x_1 : Fin 1) => euclidNorm n x

          Text 12.3.5: Rotational symmetry is equivalent to depending only on the Euclidean norm.