Hamiltonian and Lagrangian Convex Duality
📂Partial Differential EquationsHamiltonian and Lagrangian Convex Duality
Theorem
Legendre Transform
- L is a convex function.
- ∣v∣→∞lim∣v∣L(v)=+∞
Given the conditions, the Legendre transform L∗:Rn→R of L for the Lagrangian L:Rn→R is defined as follows:
L∗(p):=v∈Rnsup(p⋅v−L(v))∀ p∈Rn
Let’s assume the Lagrangian L satisfies the conditions for the Legendre transform to be defined. Let the Hamiltonian H be the Legendre transform of L.
H=L∗
Then, H also satisfies the two conditions necessary for the definition of the Legendre transform, and L=H∗ holds.
(a) H is a convex function.
λH(v1)+(1−λ)H(v2)≤H(λv1+(1−λ)v2)∀ v1,v2∈Rn,∀ 0≤λ≤1
(b) ∣p∣→∞lim∣p∣H(p)=+∞
(c) L=H∗
The converse also holds. Moreover, if H and L are differentiable at p and v∈Rn, then all the following statements are equivalent.
- p⋅v=L(v)+H(p)
- p=DL(v)
- v=DH(p)
Proof
(a)
Suppose p,q∈Rn and 0≤τ≤1. Then, the following holds.
H(τp+(1−τ)q)=v∈Rnsup((τp+(1−τ)q)⋅v−L(v))=v∈Rnsup(τ(p⋅v−L(v))+(1−τ)(q⋅v−L(v)))≤τvsup(p⋅v−L(v))+(1−τ)vsup(q⋅v−L(v))=τH(p)+(1−τ)H(q)
The second equality follows by adding τL(v)−τL(v). The third line holds because L is a convex function. The final equality is due to the definition of the Legendre transform and the assumption H=L∗.
■
(b)
Consider any fixed positive number λ>0 and suppose p=0∈Rn. Then, by assumption, the following holds.
H(p)=v∈Rnsup(p⋅v−L(v))
If we let v=λ∣p∣p, then by the definition of sup, if any v is substituted into the above equation, it results in a value less than or equal, thus the following is obtained.
H(p)=v∈Rnsup(p⋅v−L(v))≥λ∣p∣−L(λ∣p∣p)≥λ∣p∣−B(0,λ)maxL
Dividing both sides by ∣p∣ gives:
∣p∣H(p)≥λ−∣p∣maxL
This implies that the following equation holds.
∣p∣→∞liminf∣p∣H(p)≥λ
Since the above equation holds for any λ, we obtain:
∣p∣→∞liminf∣p∣H(p)=∞
Therefore, it is as follows:
∣p∣→∞lim∣p∣H(p)=∞
■
(c)
Since H=L∗ and L∗=sup(p⋅v−L(v)), for any v, H is always greater or equal.
H(p)≥p⋅v−L(v)∀ v∈Rn
Permuting H and L gives:
L(v)≥p⋅v−H(p)∀ v∈Rn
Taking sup for p∈Rn on both sides, the right side equals H∗ by the definition of the Legendre transform.
L(v)≥p∈Rnsup(p⋅v−H(p))=H∗(v)
Now, to prove that the opposite inequality L(v)≤H∗(v) holds, H∗ is defined as follows by definition.
H∗(v)=p∈Rnsup(p⋅v−H(p))
However, since we assumed H=L∗, we obtain:
H∗(v)=p∈Rnsup(p⋅v−r∈Rnsup(p⋅r−L(r)))=psuprinf(p⋅(v−r)+L(r) )
Lemma
If a function f:Rn→R is convex, for all x∈Rn, there exists r∈Rn satisfying the following:
f(y)≥f(x)+r⋅(y−x)∀ y∈Rn
Since L is convex, by the Lemma, for all v∈Rn, there exists s∈Rn satisfying the following condition.
L(r)≥L(v)+s⋅(r−v)∀ r∈Rn
Permuting s⋅(r−v) and substituting p=s gives (eq1) as follows.
H∗(v)≥rinf(s⋅(v−r))=L(v)
Therefore, H∗(v)≥L(v) and, H∗(v)=L(v) holds.
■