Alphabeta Math
CorollaryStatement: AI-adaptedProof: AI-adaptedPipeline-generatedprecheck pass
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Variation of constants for z′=Az+b(t) via the matrix exponential

Statement

Let n≥1, equip Rn with the Euclidean inner product and norm of The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, let A∈Mn(R) (The vector space Mm×n(F):=F m×n of m by n matrices over a field, with entrywise operations), let t0∈R and z0∈Rn, and let b:R→Rn be continuous (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral). The Euclidean norm is a norm by Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation, and the standard coordinate vectors form an ordered basis by The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0. For A=(aij) and z∈Rn, define (Az)i:=∑j<naijzj; this is a linear map by Linear map between vector spaces over the same field, hence bounded by A linear map from a finite-dimensional normed space is bounded, and ∥A∥ denotes its induced Euclidean operator norm (The operator norm as the least bound and as the unit-sphere or unit-ball supremum). Matrix powers and products use Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes; their action agrees with composition of the corresponding linear maps by distributing the finite sums in the displayed coordinate formula. Define the matrix exponential

E:R⟶Mn(R),E(t):=∑k=0∞tkk!Ak.

Then:

  1. E is well defined and of class C1, with E(0)=In, E(u+v)=E(u)E(v), dduE(u)=AE(u)=E(u)A for all u,v∈R, and ∥E(u)∥≤e∥A∥∣u∣ for the operator norm (The operator norm as the least bound and as the unit-sphere or unit-ball supremum) subordinate to the Euclidean norm.
  2. The map z:R⟶Rn,z(t):=E(t−t0)z0+∫t0tE(t−s)b(s) ds, satisfies z(t0)=z0 and z′(t)=A z(t)+b(t) for every t∈R.
  3. Uniqueness: if J⊆R is an interval with nonempty interior containing t0 (Intervals of R: the nine order-convex forms, nondegeneracy, and length) and w:J→Rn is differentiable with w′=Aw+b on J and w(t0)=z0, then w=z∣J.

Facts & Assumptions

Given: n≥1, A∈Mn(R) identified with a linear endomorphism of Rn, z0∈Rn, a continuous b:R→Rn, and the matrix exponential E displayed above.

[F1]

The matrix action defined in the Statement is a bounded linear endomorphism A^ of normed Euclidean Rn: the Euclidean norm is a norm, and the standard coordinate vectors form a finite ordered basis, so the finite-dimensional boundedness theorem applies (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation, The standard list e:n→Fn with ei(i)=1F and ei(j)=0F for j≠i is an ordered basis of Fn; hence dim⁡FFn=n, and F0 is the zero space with basis ∅ and dimension 0, A linear map from a finite-dimensional normed space is bounded). Its induced metric agrees with d2 (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), which is complete (R and Rn for n≥1 with the Euclidean metric are complete, componentwise from the Cauchy criterion in R); hence this normed space is Banach. For matrices B=(bij) and C=(cjℓ), ((BC)z)i=∑ℓ<n(∑j<nbijcjℓ)zℓ=(B(Cz))i, so matrix multiplication represents composition (Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes). The exponential series lemma applies to A^ and gives absolute convergence in operator norm, E∈C1 with E′=AE=EA, E(0)=In, E(u+v)=E(u)E(v) and ∥E(u)∥≤e∥A∥∣u∣ (The exponential series of a bounded operator).

[F2]

The convolution differentiation lemma (Differentiation of a variable-limit convolution integral with a C-one matrix kernel) applies with interval I=R, open kernel domain U=R⊇I−I, g:=E:U→Mn(R) and f:=b:I→Rn: for every t, the map t↦∫t0tE(t−s)b(s) ds is differentiable with derivative E(0)b(t)+∫t0tE′(t−s)b(s) ds.

[F3]

Differentiation rules: the chain rule and the algebra of derivatives differentiate products and compositions; for a differentiable matrix curve M and constant vector v, t↦M(t)v is differentiable with derivative M′(t)v. [thm-chain-rule, thm-algebra-of-derivatives]

[F4]

Fundamental theorem for continuous derivatives: if u:J→Rn is differentiable with continuous derivative on an interval J containing t0, then ∫t0tu′=u(t)−u(t0) (If f:[a,b]→Rm is differentiable with integrable f′ then ∫abf′=f(b)−f(a); and a bounded derivative makes f Lipschitz); and ∥∫t0tϕ∥2≤∣∫t0t∥ϕ∥2∣ (For a≤b and f:[a,b]→Rm integrable when a<b, ∥∫abf∥2≤∫ab∥f∥2; for a<b, ∥f∥2 is integrable, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[F5]

Gronwall: if a continuous real u≥0 satisfies u(t)≤A+∫t0tB u for t≥t0 with constants A,B≥0, then u(t)≤AeB(t−t0); the time-reflected form with the integral from t to t0 applies for t≤t0 (Gronwall's integral inequality with variable and constant coefficients).

Proof

technique · direct; verify the explicit solution by the convolution lemma and the exponential identities, then prove uniqueness by Gronwall
1.1F1F3given

By [F1] the exponential E is C1 with E′=AE=EA, E(0)=In, E(u+v)=E(u)E(v) and ∥E(u)∥≤e∥A∥∣u∣ for all u,v; in particular t↦E(t−t0) is differentiable with derivative AE(t−t0), by the chain rule for the affine inner map t↦t−t0.

2.1step 1.1F2F3given

By step 1.1 and [F3], t↦E(t−t0)z0 has derivative AE(t−t0)z0. By [F2] the integral term has derivative E(0)b(t)+∫t0tE′(t−s)b(s) ds=b(t)+∫t0tAE(t−s)b(s) ds, where the second identity uses E′=AE of step 1.1 and linearity of the integral. Hence z′(t)=AE(t−t0)z0+b(t)+A∫t0tE(t−s)b(s) ds=Az(t)+b(t), by linearity of matrix-vector multiplication. Moreover z(t0)=E(0)z0+∫t0t0=z0, so clause 2 holds.

3.1step 2.1F4F5∎

For clause 3 let w solve w′=Aw+b on an interval J containing t0 with w(t0)=z0, and put u:=w−z∣J. By step 2.1, u′=Au and u(t0)=0, so by [F4] u(t)=∫t0tA u(s) ds(t∈J), and therefore ∥u(t)∥2≤∣∫t0t∥A∥ ∥u(s)∥2 ds∣≤∥A∥ ∣∫t0t∥u(s)∥2 ds∣ for t≥t0 and the reflected inequality for t≤t0. Applying [F5] with A0=0 and B=∥A∥ to the nonnegative continuous function ∥u∥2 gives ∥u(t)∥2≤0 on both sides of t0, hence u=0 and w=z∣J.

Depends on

Dependency tree · 2 levels

Sources