TheoremBase

The class and the derivatives are computed from the coordinate expression of the quadratic form; translation is handled by noting that the difference quotients of the translated function are those of the original; and the convexity and semiconvexity claims follow from an algebraic identity for the quadratic form along a segment.

Proof

Throughout we use the notation of the statement. We first record an identity used twice below.

(P) Let x,y∈Rnx,y\in\mathbb{R}^{n}, let t∈Rt\in\mathbb{R}, put s=1−ts=1-t and w=t x+s yw=t\,x+s\,y. Then

t (x⋅(Mx))+s (y⋅(My))−w⋅(Mw)=t s ((x−y)⋅(M(x−y))).t\,\bigl(x\cdot(Mx)\bigr)+s\,\bigl(y\cdot(My)\bigr)-w\cdot(Mw)=t\,s\,\bigl((x-y)\cdot\bigl(M(x-y)\bigr)\bigr).

Indeed, since M⊤=MM^{\top}=M, claim 5 of Elementary Properties of the Transpose of a Real Matrix together with claim 1 of The Dot Product on Rn\mathbb{R}^n Is Symmetric, Bilinear and Positive Definite gives x⋅(My)=y⋅(Mx)x\cdot(My)=y\cdot(Mx). By claim 3 of Linearity of the Matrix-Vector Product and the Quadratic Form as a Double Sum we have Mw=t(Mx)+s(My)Mw=t(Mx)+s(My), so expanding with claims 2, 3, 4 and 5 of The Dot Product on Rn\mathbb{R}^n Is Symmetric, Bilinear and Positive Definite,

w⋅(Mw)=t2 x⋅(Mx)+2ts x⋅(My)+s2 y⋅(My).w\cdot(Mw)=t^{2}\,x\cdot(Mx)+2ts\,x\cdot(My)+s^{2}\,y\cdot(My).

Hence the left-hand side of the display equals (t−t2) x⋅(Mx)+(s−s2) y⋅(My)−2ts x⋅(My)(t-t^{2})\,x\cdot(Mx)+(s-s^{2})\,y\cdot(My)-2ts\,x\cdot(My). Since t+s=1t+s=1 we have t−t2=t(1−t)=tst-t^{2}=t(1-t)=ts and s−s2=s(1−s)=sts-s^{2}=s(1-s)=st, so this equals

ts(x⋅(Mx)−2 x⋅(My)+y⋅(My))=ts ((x−y)⋅(M(x−y))),ts\Bigl(x\cdot(Mx)-2\,x\cdot(My)+y\cdot(My)\Bigr)=ts\,\bigl((x-y)\cdot\bigl(M(x-y)\bigr)\bigr),

the last equality by expanding (x−y)⋅(M(x−y))(x-y)\cdot(M(x-y)) in the same way and using x⋅(My)=y⋅(Mx)x\cdot(My)=y\cdot(Mx).

Claim 1. By claim 4 of Linearity of the Matrix-Vector Product and the Quadratic Form as a Double Sum and the coordinate description of the dot product,

Q(z)=12∑i=1n∑j=1nMijzizj+∑i=1nqizi+c.Q(z)=\tfrac{1}{2}\sum_{i=1}^{n}\sum_{j=1}^{n}M_{ij}z_{i}z_{j}+\sum_{i=1}^{n}q_{i}z_{i}+c .

For i∈[n]i\in[n] let πi:Rn→R\pi_{i}:\mathbb{R}^{n}\to\mathbb{R} be the iith coordinate function, πi(z)=zi\pi_{i}(z)=z_{i}. The set Rn\mathbb{R}^{n} is open, every point being the centre of a ball contained in it. By claim 2 of Constants, Affine Functions, Sums, Products and Quotients of Cᵏ Functions the functions πi\pi_{i} and the constant functions are smooth on Rn\mathbb{R}^{n}, hence of class C2C^{2}, and by claim 3 of that lemma sums, scalar multiples and products of functions of class C2C^{2} are again of class C2C^{2}. Since a finite sum is defined recursively, induction on the upper limit shows that a finite sum of functions of class C2C^{2} is of class C2C^{2}. Applying this to the display shows that QQ is of class C2C^{2} on Rn\mathbb{R}^{n}.

We next compute the partial derivatives. Directly from Partial Derivative of a Real-Valued Function on an Open Subset of Rn\mathbb{R}^n, for k,i∈[n]k,i\in[n] and a∈Rna\in\mathbb{R}^{n} the difference quotient of πi\pi_{i} at aa with respect to the kkth variable is constantly 11 if k=ik=i and constantly 00 otherwise, so ∂kπi(a)=1\partial_{k}\pi_{i}(a)=1 if k=ik=i and ∂kπi(a)=0\partial_{k}\pi_{i}(a)=0 otherwise. By claim 1 of Constants, Affine Functions, Sums, Products and Quotients of Cᵏ Functions, extended to finite sums by the same induction as above, and by the product rule there,

∂kQ(z)=12∑i=1n∑j=1nMij(∂kπi(z) zj+zi ∂kπj(z))+∑i=1nqi ∂kπi(z).\partial_{k}Q(z)=\tfrac{1}{2}\sum_{i=1}^{n}\sum_{j=1}^{n}M_{ij}\Bigl(\partial_{k}\pi_{i}(z)\,z_{j}+z_{i}\,\partial_{k}\pi_{j}(z)\Bigr)+\sum_{i=1}^{n}q_{i}\,\partial_{k}\pi_{i}(z).

By claim 7 of Properties of Finite Sums the last sum equals qkq_{k}. Splitting the double sum by claim 2 of Properties of Finite Sums, applied to the inner and then the outer sum, gives two double sums. In the first, for fixed ii claim 3 of Properties of Finite Sums gives ∑jMij∂kπi(z)zj=∂kπi(z)∑jMijzj=∂kπi(z) (Mz)i\sum_{j}M_{ij}\partial_{k}\pi_{i}(z)z_{j}=\partial_{k}\pi_{i}(z)\sum_{j}M_{ij}z_{j}=\partial_{k}\pi_{i}(z)\,(Mz)_{i} by Matrix-Vector Product, and then claim 7 of Properties of Finite Sums over ii leaves only the term i=ki=k, giving (Mz)k(Mz)_{k}. In the second, for fixed ii only the term j=kj=k survives, again by claim 7, giving MikziM_{ik}z_{i}, and summing over ii and using Mik=MkiM_{ik}=M_{ki} gives ∑iMkizi=(Mz)k\sum_{i}M_{ki}z_{i}=(Mz)_{k}. Hence

∂kQ(z)=12((Mz)k+(Mz)k)+qk=(Mz)k+qk=(Mz+q)k,\partial_{k}Q(z)=\tfrac{1}{2}\bigl((Mz)_{k}+(Mz)_{k}\bigr)+q_{k}=(Mz)_{k}+q_{k}=(Mz+q)_{k},

so DQ(z)=Mz+qDQ(z)=Mz+q by Gradient of a Real-Valued Function on an Open Subset of Rn\mathbb{R}^n. Consequently ∂lQ\partial_{l}Q is the function z↦∑j=1nMljzj+qlz\mapsto\sum_{j=1}^{n}M_{lj}z_{j}+q_{l}, whose partial derivative with respect to the kkth variable is, by the same computation, ∂k∂lQ(z)=∑j=1nMlj∂kπj(z)=Mlk\partial_{k}\partial_{l}Q(z)=\sum_{j=1}^{n}M_{lj}\partial_{k}\pi_{j}(z)=M_{lk}. By The Hessian Matrix of a Function of Class C² the entry of D2Q(z)D^{2}Q(z) in row kk and column ll is ∂k∂lQ(z)=Mlk=Mkl\partial_{k}\partial_{l}Q(z)=M_{lk}=M_{kl}, so D2Q(z)=MD^{2}Q(z)=M.

Finally, for open V⊆RnV\subseteq\mathbb{R}^{n} the restriction of QQ to VV is of class C2C^{2} on VV by claim 3 of Restriction of a Cᵏ Function to an Open Subset, and its partial derivatives at points of VV agree with those of QQ by claim 1 of that lemma, so its gradient and Hessian at z∈Vz\in V are Mz+qMz+q and MM.

Claim 2. Let z0∈V−bz_{0}\in V-b, so z0+b∈Vz_{0}+b\in V. Since VV is open there is a positive ρ\rho with {z′:dE(z′,z0+b)<ρ}⊆V\{z':d_{E}(z',z_{0}+b)<\rho\}\subseteq V. For z∈Rnz\in\mathbb{R}^{n} we have (z+b)−(z0+b)=z−z0(z+b)-(z_{0}+b)=z-z_{0}, so dE(z+b,z0+b)=∥z−z0∥=dE(z,z0)d_{E}(z+b,z_{0}+b)=\lVert z-z_{0}\rVert=d_{E}(z,z_{0}) by claim 2 of Properties of the Euclidean Norm on Rn\mathbb{R}^n: Square, Distance, Vanishing, Coordinate Bounds, Homogeneity, Triangle Inequalities and Expansion of a Square; hence dE(z,z0)<ρd_{E}(z,z_{0})<\rho implies z+b∈Vz+b\in V, that is z∈V−bz\in V-b. So V−bV-b is open.

Let ψ\psi be of class C2C^{2} on VV and let i∈[n]i\in[n] and z∈V−bz\in V-b. For h∈Rh\in\mathbb{R} the point obtained from zz by adding hh to its iith coordinate, translated by bb, is the point obtained from z+bz+b by adding hh to its iith coordinate; therefore the difference quotients of ψb\psi_{b} at zz and of ψ\psi at z+bz+b with respect to the iith variable coincide, and Partial Derivative of a Real-Valued Function on an Open Subset of Rn\mathbb{R}^n gives that ∂iψb(z)\partial_{i}\psi_{b}(z) exists and equals ∂iψ(z+b)\partial_{i}\psi(z+b). In other words ∂iψb=(∂iψ)b\partial_{i}\psi_{b}=(\partial_{i}\psi)_{b}, the translate of ∂iψ\partial_{i}\psi in the same sense.

The translation map T:V−b→VT:V-b\to V, T(z)=z+bT(z)=z+b, satisfies dE(T(z),T(z′))=dE(z,z′)d_{E}(T(z),T(z'))=d_{E}(z,z') as computed above, hence is Lipschitz and so continuous by A Lipschitz Map is Uniformly Continuous and Continuous. Therefore, by claim 3 of Semicontinuity and Continuity Under Composition with a Continuous Map, the translate of any function continuous on VV is continuous on V−bV-b. Since ψ\psi is of class C2C^{2} on VV, it is of class C1C^{1} and each ∂iψ\partial_{i}\psi is of class C1C^{1} on VV, by clauses 1 and 2 of Functions of Class Cᵏ on an Open Subset of ℝⁿ. Applying the two previous paragraphs, ψb\psi_{b} is continuous on V−bV-b, its partial derivatives exist and equal the translates of those of ψ\psi and are continuous, so ψb\psi_{b} is of class C1C^{1}; and the same argument applied to each ∂iψ\partial_{i}\psi shows that ∂iψb=(∂iψ)b\partial_{i}\psi_{b}=(\partial_{i}\psi)_{b} is of class C1C^{1} on V−bV-b. By clause 2 of Functions of Class Cᵏ on an Open Subset of ℝⁿ, ψb\psi_{b} is of class C2C^{2} on V−bV-b. Finally Dψb(z)=(∂1ψ(z+b),…,∂nψ(z+b))=Dψ(z+b)D\psi_{b}(z)=\bigl(\partial_{1}\psi(z+b),\dots,\partial_{n}\psi(z+b)\bigr)=D\psi(z+b), and the entry of D2ψb(z)D^{2}\psi_{b}(z) in row ii and column jj is ∂i∂jψb(z)=∂i((∂jψ)b)(z)=∂i∂jψ(z+b)\partial_{i}\partial_{j}\psi_{b}(z)=\partial_{i}\bigl((\partial_{j}\psi)_{b}\bigr)(z)=\partial_{i}\partial_{j}\psi(z+b), so D2ψb(z)=D2ψ(z+b)D^{2}\psi_{b}(z)=D^{2}\psi(z+b).

Claim 3. Let x,y∈Cx,y\in C and let t∈Rt\in\mathbb{R} with 0≤t≤10\le t\le1; put s=1−ts=1-t, so 0≤s0\le s, and w=t x+s yw=t\,x+s\,y, which lies in CC because CC is convex. Since t+s=1t+s=1, claims 2, 4 and 5 of The Dot Product on Rn\mathbb{R}^n Is Symmetric, Bilinear and Positive Definite give t (q⋅x+c)+s (q⋅y+c)=q⋅w+ct\,(q\cdot x+c)+s\,(q\cdot y+c)=q\cdot w+c. Hence, by (P),

t Q(x)+s Q(y)−Q(w)=12(t x⋅(Mx)+s y⋅(My)−w⋅(Mw))=12 ts ((x−y)⋅(M(x−y))).t\,Q(x)+s\,Q(y)-Q(w)=\tfrac{1}{2}\Bigl(t\,x\cdot(Mx)+s\,y\cdot(My)-w\cdot(Mw)\Bigr)=\tfrac{1}{2}\,ts\,\bigl((x-y)\cdot\bigl(M(x-y)\bigr)\bigr).

By hypothesis 0n⪯M0_{n}\preceq M, which by The Positive Semidefinite Ordering on Symmetric Matrices means ζ⋅(0nζ)≤ζ⋅(Mζ)\zeta\cdot(0_{n}\zeta)\le\zeta\cdot(M\zeta) for every ζ\zeta; since ζ⋅(0nζ)=∑i=1n∑j=1n0⋅ζiζj=0\zeta\cdot(0_{n}\zeta)=\sum_{i=1}^{n}\sum_{j=1}^{n}0\cdot\zeta_{i}\zeta_{j}=0 by claim 4 of Linearity of the Matrix-Vector Product and the Quadratic Form as a Double Sum together with two applications of claim 7 of Properties of Finite Sums, every summand being 00, the quadratic form (x−y)⋅(M(x−y))(x-y)\cdot(M(x-y)) is nonnegative. As tsts is nonnegative, the right-hand side is nonnegative, so Q(w)≤t Q(x)+s Q(y)Q(w)\le t\,Q(x)+s\,Q(y). By Convex Real-Valued Function on a Convex Subset of Rn\mathbb{R}^n the restriction of QQ to CC is convex on CC.

Claim 4. Write L=∥M∥L=\lVert M\rVert, which is nonnegative by claim 1 of Properties of the Norm of a Symmetric Real Matrix; hence 0≤λ+L0\le\lambda+L. Let x,y∈Cx,y\in C and t∈Rt\in\mathbb{R} with 0≤t≤10\le t\le1, put s=1−ts=1-t and w=t x+s y∈Cw=t\,x+s\,y\in C. Since ff is semiconvex on CC with constant λ\lambda, Quadratic Increment Characterisation of Semiconvexity gives

f(w)≤t f(x)+s f(y)+λ2 ts ∥x−y∥2.f(w)\le t\,f(x)+s\,f(y)+\tfrac{\lambda}{2}\,ts\,\lVert x-y\rVert^{2}.

Using t (q⋅x+c)+s (q⋅y+c)=q⋅w+ct\,(q\cdot x+c)+s\,(q\cdot y+c)=q\cdot w+c as in claim 3, we obtain

g(w)−(t g(x)+s g(y))=(f(w)−tf(x)−sf(y))+12(t x⋅(Mx)+s y⋅(My)−w⋅(Mw)),g(w)-\bigl(t\,g(x)+s\,g(y)\bigr)=\Bigl(f(w)-t f(x)-s f(y)\Bigr)+\tfrac{1}{2}\Bigl(t\,x\cdot(Mx)+s\,y\cdot(My)-w\cdot(Mw)\Bigr),

so by the previous display and (P),

g(w)−(t g(x)+s g(y))≤λ2 ts ∥x−y∥2+12 ts ((x−y)⋅(M(x−y))).g(w)-\bigl(t\,g(x)+s\,g(y)\bigr)\le\tfrac{\lambda}{2}\,ts\,\lVert x-y\rVert^{2}+\tfrac{1}{2}\,ts\,\bigl((x-y)\cdot\bigl(M(x-y)\bigr)\bigr).

By claim 2 of Properties of the Norm of a Symmetric Real Matrix and claim 3 of Properties of the Absolute Value in an Ordered Field, (x−y)⋅(M(x−y))≤L ∥x−y∥2(x-y)\cdot(M(x-y))\le L\,\lVert x-y\rVert^{2}, and tsts is nonnegative, so

g(w)≤t g(x)+s g(y)+λ+L2 ts ∥x−y∥2.g(w)\le t\,g(x)+s\,g(y)+\frac{\lambda+L}{2}\,ts\,\lVert x-y\rVert^{2}.

As x,y∈Cx,y\in C and tt were arbitrary, Quadratic Increment Characterisation of Semiconvexity shows that gg is semiconvex on CC with constant λ+L=λ+∥M∥\lambda+L=\lambda+\lVert M\rVert.

Citations

Loading…

Dependencies

Uses0

Loading…

Comments

Log in to comment.

Loading…