A periodic function is defined by the property

f(t+T) =f(t)

If the function has an exponential-type power series expansion

f(t) =k=0 ck tkk!

then using the binomial theorem the periodic property translates to

f(t+T) =k=0 ck (t +T)k k! =k=0 ckk! l=0 k (kl) Tk-l tl f(t +T) =l=0 k=l ckl! (k-l)! Tk-l tl =l=0 m=0 cl+m m! Tm tll!

Comparing this to the original power series, one must have in general

ck =m=0 ck+m m! Tm m=1 ck+m m! Tm =0

where the second form removes the first term from both sides of the equivalence. Since the factorial and the period are both positive, the expansion coefficients must include negative values to satisfy these conditions.

For the fundamental periodic functions with expansions

sint= k=0 (1)k t2k+1 (2k+1)! cost= k=0 (1)k t2k (2k)!

one has immediately

k=0 (1)k (2π )2k+1 (2k+1)! =0 k=1 (1)k (2π )2k (2k)! =0

which can be confirmed with numerical approximations. To compare this with the general result above, write the expansion coefficients of the cosine function as

c2k =(1)k , c2k+1 =0

and insert in the general result:

m=1 c2k+m m! Tm =m=1 c2k+2m (2m)! T2m =m=1 (1 )k+m (2m)! (2π )2m =0

The cosine coefficients force the index here to be even as well, reproducing the immediate result just given apart from a minus sign. Oddly enough, repeating this procedure with the coefficients of the sine expansion

c2k =0 , c2k+1 =(1)k

produces the same result, since again the index is forced to be even. The immediate result above for the sine expansion has a different source.

As an aside, one can take either of the immediate results, cut off the expansion at some number of terms and treat it as a polynomial estimate for 2π that can always be solved numerically. This indicates a way to tell whether a given power series expansion represents a periodic function: successively increase the number of terms used in the process to determine that the resultant numerical value converges to a stable estimate.

Derivatives of a periodic function are also periodic, which can be seen by differentiating both sides of the defining equation:

f(n) (t+T) =f(n) (t)

Applying this to the exponential-type power series expansion gives

f(n) (t+T) =k=n ck (t+T )k-n (k-n)!

and comparing this to the original expansion at t=0 leads to a second general result:

cn =k=n ck Tk-n (k-n)! k=n+1 ck Tk-n (k-n)! =0

With n=1 and the cosine expansion coefficients above, one has

k=2 ck (k-1)! Tk-1 =k=1 c2k (2k-1)! T2k-1 =k=1 (1)k (2k-1)! (2π )2k-1 =0

which reproduces the immediate result from the sine expansion, apart from a shift in index and a subsequent minus sign. Repeating the procedure with the sine expansion coefficients apparently requires n=2 as the starting point.

Given that the second general result is a set of linear equations for the constants, one might think to stop the expansion at some larger number and solve backward for all of the coefficients of a general periodic function. First compare the equations to a summation formula for Bernoulli numbers, which have alternating signs:

m=1 N-1 (Nm) Bm =N! m=1 N-1 Bm (N-m)! m! =0

Equating the first factorial in the denominator to the one in the equations above gives the relation m=N+n-k  , so that it is not difficult to see that the recursive solution leading to n=0 is simply

ck=cN TN-k BN-k (N-k)!

Since the factorial increases faster than the Bernoulli number as N , the most general periodic function is a constant. Earth shattering!

The form of the coefficients can be confirmed with the following Mathematica code:

max = 10;
s[n_] := Sum[Subscript[c, n + m] T^m/m!, {m, 1, max - n}]
Solve[Table[s[i] == 0, {i, 0, max - 2}], 
 Table[Subscript[c, i], {i, 1, max - 1}]]

Given the expansion

xcothx =k=0 B2k (2k)! (2x )2k

one has immediately that the coefficients even in the period are identical to an expansion of T2 cothT2  . The linear term for cN-1 is oddly missing...

The nonpathological periodic functions used in physics can always be represented by Fourier series. If a function has the expansion

f(t) =n ancos( 2πntT ) +bnsin( 2πntT )

then consider it term by term:

f(t) =k=0 (1)k t2k (2k)! [ n an (2πnT )2k ] +(1)k t2k+1 (2k+1)! [ n bn (2πnT )2k+1 ]

This is potentially useful in correlating a known Fourier expansion with a power series expression of the same function. The coefficients in brackets are simple linear combinations of the Fourier coefficients.

Gravitational choreographies typically have constraints on the center of mass that leads to relationships of the form

f(t) +f(t+T) +f(t-T) =0

Again using an exponential-type power series expansion for the function, this relationship becomes

k=0 ckk! [tk +(t +T)k +(t -T)k] =0 k=0 ckk! [tk +l=0 k (kl) [1-(1 )k-l] Tk-l tl] =0 k=0 ckk! tk +l=0 k=l ckl! (k-l)! [1-(1 )k-l] Tk-l tl =0 k=0 ckk! tk+ l=0 m=0 cl+m m! [1-(1 )m] Tm tll! =0

Note that the double sum only contains even powers of the period. Repurposing the index, for gravitational choreographies one must have

ck +2 m=0 ck+2m (2m)! T2m =3ck +2 m=1 ck+2m (2m)! T2m =0

This statement holds separately for the coefficients of even and odd powers of the temporal variable.

If one again assumes that the expansion can be stopped at some large number and solved backward for the coefficients of a general gravitational choreography. In this case it is less easy to see the form of the final result, so Mathematica is useful. The code here is

max = 11;
s[k_] := 3 Subscript[c, k] + 
  2 Sum[Subscript[c, k + 2 m]/(2 m)! T^(2 m), {m, 1, (max - k)/2}]
Solve[Table[s[i] == 0, {i, 0, max - 2, 2}], 
 Table[Subscript[c, i], {i, 0, max - 2, 2}]]
Solve[Table[s[i] == 0, {i, 1, max - 2, 2}], 
 Table[Subscript[c, i], {i, 1, max - 2, 2}]]

Even and odd powers have the same numerical coefficients matching backwards from highest coefficients.

Picking any high numerator and searching on oeis.org reveals that the numerical coefficients have the same absolute value as A279120, which are the numerators in the expansion of 11 +2cosx  . The denominators are in the next sequence A279121 and have a common factor of three that can be removed by simple multiplication. This is an oddly specific function: what could be its significance for gravitational choreographies?

Uploaded 2021.04.30 — Updated 2022.11.14 analyticphysics.com