Alphabeta Math
ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck 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.

An explicit weakly damped oscillator: the averaging constants and the energy-decay bound over five thousand seconds (about 796 undamped periods)

Physical scope

A worked numerical instance of the weak-damping averaging estimate: the model m x'' + epsilon gamma x' + k x = 0 with m = 1 kg, k = 1 kg/s^2, epsilon = 10^-4 and gamma = 1 kg/s, initial rotating-frame amplitude zeta_0 = 1, ball radii R = 1/2 and rho = 1/4 and horizon parameter A = 1/2, verifying the smallness and trapping hypotheses with the exact constants M = 3 sqrt2/2, L = sqrt2, K = 2 pi (3 sqrt2 + 6) exp(sqrt2/2) and concluding that the mechanical energy tracks the averaged exponential exp(-t/10^4) within 2.63 x 10^-2 on the interval [0,5000]. The numbers are exact-arithmetic evaluations of the model formulas, not measurements; no apparatus, no observed decay and no empirical statement is present.

Example

Adopt the damped oscillator model of The damped oscillator model: a harmonic restoring force with a phenomenological linear drag −cx˙, c≥0 with data

m=1 kg,k=1 kg/s2,ε=10−4,γ=1 kg/s,

so that c=εγ=10−4 kg/s and ω0=k/m=1 rad/s; the angles and frequencies are read under the convention of Dimensions and SI unit conventions of the non-relativistic classical mechanics framework. Fix amplitude unit a∗=1 m s−1/2 and time unit 1 s. Take ζ0=a∗, R=a∗/2, ρ=a∗/4 and A=(1 s)/2. Below these symbols, t and the constants in coefficient calculations mean their numerical values in these fixed units. In particular M is in a∗/s, L in s−1, B,K in a∗ and B1 is dimensionless. Every energy exponential below is multiplied by E0=1 J when read as a physical quantity, and the horizon [0,5000] means [0,5000 s]. For x>0, write ln⁡x:=log⁡x, where log⁡:(0,∞)→R is the inverse of the real exponential (The natural logarithm as the inverse of the exponential function). Then:

  1. the averaging data of Averaging for a weakly damped oscillator: the rotating-frame amplitude, the averaged exponential energy decay and the order-ε bound on the O(1/ε) time scale specialize to P=2π, F(t,ζ)=−ieit(cos⁡tIm⁡ζ−sin⁡tRe⁡ζ) and Fˉ(ζ)=−ζ/2, with averaged solution ζˉ(t)=e−εt/2=e−t/(2⋅104);
  2. the constants of the averaging estimate are M=322,L=2,B=62 π,B1=42 π,K=2π(32+6)e2/2, so that K≈1.30522×102,εK≈1.30522×10−2<0.25=ρ;
  3. the trapping hypothesis holds: ∣ζ0∣(1−e−γA/(2m))=1−e−1/4≈0.22120<R−ρ=0.25, and also 1/4=γA/(2m)≤ln⁡(4/3)≈0.28768;
  4. consequently the motion with datum x(0)=2 m, x′(0)=0 — which has ζ(0)=1 and mechanical energy E0=mω0∣ζ0∣2=1 J — satisfies ∣E(t)−e−t/104∣≤mω0εK(2∣ζ0∣+εK)≈2.62748×10−2 J(0≤t≤5000), the interval [0,5000] being [0,A/ε] and containing 5000/(2π)≈795.8 oscillation periods of the undamped frequency;
  5. at the end of the horizon the averaged energy is e−1/2≈0.60653 J, so the estimate bounds the deviation of the exact energy from the averaged exponential decay by about 4.3% of its value there.

Facts & Assumptions

Given: The model data m=1, k=1, ε=10−4, γ=1; the numbers ζ0=1, R=1/2, ρ=1/4, A=1/2; the constants M,L,B,B1,K of the Example.

[F1]

Averaging theorem (Averaging for a weakly damped oscillator: the rotating-frame amplitude, the averaged exponential energy decay and the order-ε bound on the O(1/ε) time scale clauses 1-4): for the model with data (m,εγ,k) the rotating-frame amplitude maps motions bijectively onto solutions of ζ′=εF(t,ζ); F is P=2π/ω0-periodic with average Fˉ(ζ)=−γζ/(2m); under the hypotheses εK<ρ and the trapping condition, the estimate ∥ζ(t)−ζˉ(t)∥2≤εK and the energy estimate ∣E(t)−E0e−εγt/m∣≤mω0εK(2∣ζ0∣+εK) hold on [0,A/ε], with M=(γ/m)2(∣ζ0∣+R), L=(γ/m)2, B=2MP, B1=2LP and K=eLA[B+(LB+B1M)A].

[F2]

Data specialization (Periodic averaging data: the average field, the oscillatory part and the periodic deviation, Real and imaginary parts, complex conjugation, and modulus): with ω0=1 and γ/m=1 the field of the theorem reads F(t,ζ)=−ieit(cos⁡tIm⁡ζ−sin⁡tRe⁡ζ) and P=2π; ∣ζ0∣=1.

[F3]

Elementary arithmetic (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, The derivatives of sine and cosine are cosine and minus sine): the displayed decimal values are evaluations of the exact expressions, obtained by direct arithmetic from 2≈1.41421356, π≈3.14159265 and e2/2≈2.02811498; in particular the datum with ζ0=1 corresponds to x(0)=2 m, x′(0)=0 and E0=mω0∣ζ0∣2=1 J.

[F4]

Exact numerical inequalities: 2<3/2 by squaring and 2<e<3 by The elementary numerical bound 2<e<3. For x∈[1/4,3/4], x(1−x)≥3/16 and 1+x2≤25/16, so x4(1−x)41+x2≥1625(316)4>0. Polynomial division gives x4(1−x)41+x2=x6−4x5+5x4−4x2+4−41+x2. By Principal arctangent: derivative, integral, power series, and the Gregory–Leibniz series, ∫01dt/(1+t2)=arctan⁡1. Also arctan⁡1=π/4: the principal inverse tangent of The principal inverse tangent arctan⁡:R→(−π/2,π/2) inverts tan⁡ on (−π/2,π/2); the cofunction identity of Cofunction, supplementary, quarter-turn, and reflection identities for the six trigonometric functions at x=π/4 gives sin⁡(π/4)=cos⁡(π/4), hence tan⁡(π/4)=1 by Tangent, cotangent, secant, and cosecant on their exact natural domains. Termwise integration of the polynomial using For a natural n≥1 the function x↦xn is differentiable everywhere with derivative ι(n) x n−1; for n=0 it is the constant 1, with derivative 0; for a natural n≥1 the function x↦x−n is differentiable at every x≠0 with derivative −ι(n) x−n−1; consequently every polynomial function is differentiable at every real, with the derivative computed term by term and Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G now gives ∫01x4(1−x)41+x2 dx=227−4arctan⁡1=227−π. The integrand is nonnegative on [0,1], and its integral on [1/4,3/4] is at least 12⋅1625(316)4>0 by If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a), so the full integral is positive and π<22/7. Also, for x>0, write ln⁡x:=log⁡x as in The natural logarithm as the inverse of the exponential function. Then ln⁡(4/3)=∫14/3dt/t≥34⋅13=14 by The natural logarithm has derivative 1/x and equals the integral from 1 to x of 1/t and integral monotonicity. For u>0, the fundamental theorem and exponential derivative give 1−e−u=∫0ue−s ds<u: on [0,u/2] the integrand is at most 1, and on [u/2,u] it is at most e−u/2<1 by strict increase of the exponential (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G, The exponential function is smooth and (exp⁡)′=exp⁡, The exponential function is strictly increasing, If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a)). Finally, e−1/2>1/3 follows from e1/2<e<3 and e−x=1/ex (The elementary numerical bound 2<e<3, The exponential function is strictly increasing, The exponential is positive and satisfies exp⁡(−x)=1/exp⁡(x)).

Verification

technique · direct; specialize the theorem's data, evaluate the constants, check the two explicit hypotheses, and state the conclusion
1.1F1F2algebra

Specialization: with m=k=γ=1 and ε=10−4 the quantities of [F1] become ω0=1, P=2π, Fˉ(ζ)=−ζ/2 and ζˉ(t)=e−εt/2=e−t/(2⋅104); the field is the one displayed in clause 1.

1.2F1F4givenalgebra

Constants: substituting ∣ζ0∣=1, R=1/2, γ/m=1 gives M=2⋅32=322 and L=2, hence B=2MP=62 π and B1=2LP=42 π. Then LB=2⋅62 π=12π and B1M=42 π⋅322=12π, so K=e2/2[62 π+12(12π+12π)]=2π(32+6)e2/2, and by [F4] K<198<200: 2<3/2 and e2/2<e<3 give K<2π⋅(3⋅32+6)⋅3=63π<63⋅22/7=198. This is clause 2.

2.1F1F4step 1.2givenalgebra

Hypotheses: εK<198/104<1/50<1/4=ρ by step 1.2, so the smallness condition holds; at the full horizon A/ε=5000 the trapping quantity is ∣ζ0∣(1−e−γA/(2m))=1−e−1/4<1/4=R−ρ by [F4]. The sufficient logarithmic condition also holds, since γA/(2m)=1/4≤ln⁡(4/3) by [F4]. Hence the averaged solution stays in the required ball. This is clause 3.

3.1F1F3F4step 1.1step 1.2step 2.1algebra

Conclusion: applying [F1] with the verified hypotheses and the horizon A/ε=5000 gives ∥ζ(t)−ζˉ(t)∥2≤εK<1/50 for 0≤t≤5000, hence the energy bound ∣E(t)−e−t/104∣≤εK(2+εK)<1012500<124 J; the datum corresponds to x(0)=2 m, x′(0)=0, E0=1 J, and the horizon contains more than 795 periods of the undamped frequency, because 5000>1590π and [F4] gives 1590π<1590⋅227=349807<5000; at t=5000 the averaged energy is e−1/2>1/3 by [F4], so the bound is below one eighth of it. This is clauses 4 and 5.

4.1step 1.1step 1.2step 2.1step 3.1∎

Clauses 1-5 are verified by exact computation from the formulas of the averaging theorem, with the rational and elementary estimates above in place of decimal arithmetic; no measured quantity, apparatus or empirical claim is present.

Remarks

  • Decimal orientation values. On the stated instance, K≈1.30522×102, εK≈1.30522×10−2 and the energy bound is approximately 2.63×10−2 J; the averaged energy at t=5000 is approximately 0.60653 J. These decimals are evaluations of the exact expressions of the proof and play no role in it.

  • Reading the constant. The dominating factor in K is eLA=e2/2, below 3, because the horizon parameter A=1/2 is small: the estimate is useful for horizons of order one in the slow time A, and K grows like AeLA with A. Taking A larger lengthens the interval [0,A/ε] but increases both K and the required smallness margin.

  • The finite check is arithmetic, not empirical. Every number above is computed from the model formulas; nothing in the example depends on an experiment, and the estimate is a conditional bound for the declared model.

Depends on

Dependency tree · 2 levels

Sources