1. Introduction
Periodic functions and their representation by Fourier series possess a great importance in different branches of fundamental and applied sciences, in particular, in optics. It is well known that the convergence of Fourier series of discontinuous functions show a peculiar property of oscillations in the neighborhoods of the discontinuities where they increase to higher becoming and remaining overshoots of oscillations with decreasing distances from the discontinuities about (or below) its level. This behavior is called the Gibbs phenomenon after J. W. Gibbs who described it in about 1898 [1]. The author of the short article Ul’yanov in the Mathematical Encyclopedia mentioned that the phenomenon was discovered already before Gibbs by H. Wilbraham in 1848 and gives the original citations1. The phenomenon is represented in detail, e.g., by Sommerfeld [2], Fichtenholz [3], Arfken [4] and surely by others. Körner [5] in his monograph also mentions Wilbraham as a precurser of the Gibbs phenomenon but without any references though with some pieces to history in the Appendices. More about Fourier series and the Gibbs phenomenon with many citations one finds in Weisstein’s MathWorld [6] under the corresponding notions.
In present paper we calculate the Fourier series of the 2π-periodic function
and of related functions which are not convergent in usual sense but are Generalized function or linear continuous functionals over a certain space of basis function (e.g., L. Schwartz [7], Gelfand and Shilov [8], Vilyenkin [9], Vladimirov [10]) as explained in the further text and illustrated by figures in Sections 2 - 4. The Tangent and related functions show new peculiarities of convergence in their Fourier series and it seems that they are worth to be described and illustrated. Apparently this phenomenon is present if the considered periodic function possesses singularities of, at least, first or higher order (but not logarithmic singularities). Whereas the Gibbs phenomenon is connected with discontinuities of the functions the here investigated phenomenon is connected with singularities of the kind of delta functions that arise from differentiation of functions with discontinuities but in our case from differentiation of logarithmic singularities. Insofar, the Gibbs phenomenon is connected with the phenomenon described but concerns functions with poles of first order for which in addition the Fourier series is not convergent in the ordinary sense. In Sections 5 and 6 we describe and illustrate by figures this phenomenon for two other related examples, the reciprocal Sine function
and the corresponding reciprocal Cosine function. In another form for identities with 2π-periodic delta functions and their Fourier series which are derived from the Poisson summation formula we investigate and illustrate the convergence phenomenon in Section 7.
In Section 8 we derive a modification of Fourier series with power-like expansions of the discussed 2π-periodic functions in (non-orthogonal) powers of Cosine and Sine functions. The considered examples of such expansions are convergent even in such cases when the Fourier series converge only in the sense of weak convergence of Generalized functions. A more general consideration of these power-like expansions with calculation of their coefficients is made in Sections 9 and 10 plus Appendices.
Another aspect is that the specially considered functions are connected with a kind of Higher Transcendental Functions called Polylogarithms (Apostol [11], Andrews, Askey, Roy [12], in last mainly Dilogarithms), where their notation
was introduced in more recent times (second half of last century) and which include many other Elementary and Special functions as special cases (e.g, Tangent and Logarithm function and also the Riemann zeta function
as the special case
of them although the last not as a function of z). We take into account the connection of the here-considered functions to the function class of Polylogarithms.
2. Fourier Series of
and
In this and in the following Sections we develop and discuss the Fourier series of the Tangent function and of related functions. Fourier series were established by Fourier in 1807 in the theory of heat and its propagation. We begin in this Section with the preliminary introduction of the Fourier series for the Tangent function and for its first integral.
We consider the 2π-periodic function
for which the following representations are important as starting point for further calculations
(2.1)
The function
is simply related to the function
by
(2.2)
This means that in transitions from formulae for
to formulae for
and related functions obtained by analogous differentiations or integrations one has to make the substitutions
, in particular
(2.3)
and it is not necessary to calculate them separately.
The next problem is to obtain the Fourier series of
. Generally, the Fourier series of a 2π-periodic function
(2.4)
can be represented in the form
(2.5)
with the well-known coefficients
(2.6)
Without paying attention to convergence we may formally obtain an expression for the Fourier series of
. In principle, it should be considered as known because in the monograph of Gel’fand and Shilov [8] and in the article of Vilenkin [9] and also in the monograph of Tolstov [13] we found the convergent Fourier series for
from which by differentiation follows the (in the usual sense) divergent Fourier series of
with the result2
(2.7)
Then using (2.3) we find the Fourier series of
with following result
(2.8)
Vice versa, from this follows using (2.3) the Fourier series (2.7) of
. In next section, we will give a more direct derivation of the Fourier series (2.8) of
.
For further investigations, we calculate the following definite integral over the function
in the sense of the Principal value concerning the singularities
(2.9)
with 2π-periodic continuation on the right-hand side3. Thus the integration constant of the indefinite integral is chosen here in such a way that the function on the right-hand side becomes zero for all its minima which are at once symmetry points of the function. The corresponding integral over the function
is
(2.10)
where the integrand is not locally integrable at the singularities
. The mentioned functions are illustrated in Figure 1.
Figure 1. Tangent and Cotangent functions and their definite first integral. The Tangent and Cotangent functions with poles of first order are the derivatives of the logarithmic functions with logarithmic singularities below them.
According to (2.9) from the series (2.8) with term-by-term integration and division by factor −2 follows the convergent Fourier series with exact agreement of function and its approximation for argument
(2.11)
and for (2.10), correspondingly with the exact agreement of function and its approximation at
(2.12)
The right-hand sides of the Fourier series (2.11) and (2.12) possess an uncommon form for Fourier series. Taking into account the following alternating series with conditional convergence with the well-known result in the natural order of summation of terms
(2.13)
and using this in (2.11) one finds
(2.14)
and, correspondingly, using (2.12) one finds the known Fourier series [8] [9]
(2.15)
In [3] is posed and dealt with the inverse problem to determine to a given Fourier series the corresponding function which it represents, in this case the series (2.15) with the function
as the result.
3. Derivation of the Fourier Series of
A derivation of the Fourier series can be made if one calculates the coefficients by the integrals in (2.6) in case of
(3.1)
The integrals
vanish due to the antisymmetry of the integrand. The integrals for
can be directly calculated or can be taken from an independent calculation of the Fourier series that provides then a result for these integrals. We will present here a derivation of the Fourier series by a method using complex representation of the involved trigonometric functions which is similar to a method used in [3] (pp. 469-471). In the way of this derivation, we meet soon an apparent problem with the complexity of the solution which we try to clarify.
As preparation to the derivation we first write down the following summation formula for the sum of the Geometric series.
(3.2)
which is uniformly convergent for
.
The representation of
by complex exponentials can be transformed then with
as follows
(3.3)
where
is only the case at the singularities
of the function
. Relation (3.3) is almost but not fully in agreement with (2.8) due to the presence of an additional part
. The additional part is imaginary-valued but the Tangent function on the left-hand side is real-valued. However, this apparent contradiction resolves since the additional imaginary part proves to be almost zero without the delta function due an identity (see (7.3)). The notation
on the left-hand side corresponds then to the Principal value of the function at the singularities4. First we write down the mentioned identity. Using well-known formulae of the theory of Generalized functions (see Section 7, (7.1) and (7.10)) from
(3.4)
which rests on the Poisson summation identity transforming a sum over a lattice into a sum over the reciprocal lattice as follows (see (7.3))
(3.5)
The singularities of the function
at
,
approaches the singularities of the functions
for which in the theory of Generalized functions is derived the well-known formula [7]-[10]
(3.6)
which means that one adds (or subtracts) in the argument x of the function an infinitesimal small imaginary
with positive
and lets them go
. A consequence of (3.6) is
(3.7)
with positive imaginary amplitude of the delta function. This can be compared now with the general complex extension of
(3.8)
One can see that this imaginary part for
indeed vanishes everywhere with the exception of the points of singularities of the function
. The Tangent function
and (3.3) can be written now more rigorously
(3.9)
In analogy to the above derivations for the Tangent function we obtain for the corresponding Cotangent function
(3.10)
The general complex extension of the function
is
(3.11)
In Section 5 we will come across the function
displaced on the real axis by
(3.12)
Its Fourier series follows from (3.9) by the substitutions
,
and is
(3.13)
In this and other formulae of this kind, in particular, in (3.9) and in (3.10) one has to check the sign of the imaginary part for general values of x and for infinitesimal positive y.
4. Convergence Behavior of the Fourier Series of
Illustrated
The Fourier series of the Tangent function (2.8) is not convergent in usual sense but it is convergent in the sense of a continuous linear functional over a space of appropriate basis functions
which has to be determined that means
(4.1)
In case of sufficiently well-behaved functions
such as the Sine functions on the right-hand side of (4.1) the linear continuous functional is simply the integral
(4.2)
In other cases such as for the delta “function” and its derivatives it has to be defined separately and it has to be proved then that it is a linear continuous functional.
The convergence behavior of Fourier series for the Tangent function (2.8) is illustrated for two different approximations of the number of first sum terms taken into account in Figure 2. It is seen that the approximations of the Fourier series for
oscillate to both sides of the curve above and below in a certain stripe. This stripe seems to be very stabile with two envelopes which do not depend on the number of sum terms in the approximation of the Fourier series but it is not very clear how to calculate these envelopes. The superposition of two different approximation in Figure 3 shows this behavior by its homogeneous color. Furthermore, this picture shows that at variable
the envelope functions
, apparently, are exactly
and for
are exactly
(in the limit of infinite sum terms). It seems that in a neighborhood of
the heuristically found envelope functions
![]()
Figure 2. Fourier series of Tangent function in different approximations. In the figure above we used 25 sum terms and in the figure below 50 sum terms of the Fourier series for the approximation. The Tangent function itself is shown embedded in a broad band of oscillations of the approximations by the Fourier series.
Figure 3. Superposition of two approximations with 25 and 50 sum terms in the Fourier series of
. Within the range
the conjectured approximate envelopes (4.3) are drawn. The homogeneous color shows that they seem to be very stabile and seem to be independent of the used number of sum terms in the approximation.
(4.3)
if periodically continued are “good” approximations for the two envelopes.
As already said the Tangent function
has to be considered as a Generalized function and its Fourier series converges in the sense of weak convergence to this function. As the space of basis functions can be chosen the space
of continuous and continuously differentiable functions which in infinity decrease faster than any power function
to which belong also all Gauss bell functions
with arbitrary width a and displacement x0 from the center. The more special class of finite functions of space
seems to be also appropriate (
) instead of
. The Tangent function
can be then considered as a Generalized function of the space
. When the Tangent function is combined in integrals with such functions then these integrals or linear functionals converge in ordinary sense.
The Fourier series of the first integral of the function
which is the function
is convergent in ordinary sense. This is illustrated in Figure 4 for two different (low) numbers of sum terms in the approximations. They do not show striking peculiarities. Fourier series within the frame of Generalized functions are discussed in all of the here cited works about Generalized functions [7]-[10] but similar to considered examples of such functions we found only in [8] [9] and additionally in [13] without support on Generalized functions.
Figure 4. Approximations of the Fourier series (2.11) to function
with 5 and 6 terms. The Fourier series converges in ordinary sense to mentioned function. The term
was taken as a whole and was not taken in the same approximation with only N sum terms of (2.13).
5. Fourier Series of
and
In analogy to Sections 2 and 3 we investigate now the Fourier series for the 2π-periodic function
(5.1)
which similar to the function
possess poles of first order at
and analogously the function
(5.2)
with poles of first order at
, the same kinds of poles as
and
. The two kinds of functions
and
on one side and
on the other side are connected by
(5.3)
The two functions
and
are only shifted to each other by
on the real axis and therefore we make derivations only for the first of these functions.
The first integrals are for (5.1) [14]
(5.4)
and for (5.2)
(5.5)
that can be checked by differentiation. The functions together with the Cosecant and Secant functions are illustrated in Figure 5. Since by (2.11) and (2.12) we know already the Fourier series of the two parts in (5.4) and (5.5) we have also their Fourier series and by term-wise differentiation the Fourier series of
and
. The result for the Fourier series of the function (5.5)
(5.6)
and for the function (5.5), correspondingly, with the variable substitution
in (2.10)
(5.7)
Figure 5. Cosecant and Secant functions and their definite first integral. The Cosecant and Secant functions are the derivatives of the functions below them which possess 2π-periodic logarithmic singularities.
They converge in ordinary way. By term-wise differentiation of (5.6) we find the Fourier series of
(5.8)
and by term-wise differentiation of (5.7) the Fourier series of
(5.9)
The last two Fourier series (5.8) and (5.9) can be only considered in the sense of Generalized functions and their convergence is weak convergence of Generalized functions.
Before we illustrate some previous results by figures in next Section we want to derive now the Fourier series (5.8) in another way in analogy to (3.3) by transition to the complex representation with application of the geometric series as follows
(5.10)
The function
is real-valued but on the right-hand side we find an
imaginary-valued additional term. Due to the identity (7.9) this additional term consists of 2π-periodic delta functions with positive and negative amplitudes and is zero in most of its range. In usual notation the left-hand side of (5.10) means the principal value of the function at the singularities which is real-valued. More correct are therefore the notations
(5.11)
and its principal value
(5.12)
In this connection we give the following complex extension of the function
(5.13)
In next Section we illustrate the convergence behavior of the Fourier series of
.
6. Convergence Behavior of the Fourier Series of
Illustrated
We illustrate now the convergence behavior of the Fourier series (5.12) for the function
(principal value) and the same for its first integral (5.6) which is the function
.
Figure 6 shows two approximations with different numbers of taken sum terms in the Fourier series of
. We see in this picture oscillations around the function which it should represent and the convergence is weak convergence in the sense of the theory of Generalized functions. This is very similar to Figure 2 for the function
. As in mentioned Figure 2 one believes to see two envelopes of the oscillations which by a suited definition not available up to now could be calculated. These envelopes apparently remain very stabile if we go to higher approximations and surely to the limiting case. One of the envelopes is obviously the x-axis
Figure 6. Fourier series of reciprocal Sine function in different approximations. In the figure above we used 15 sum terms and in the figure below 30 sum terms of their Fourier series for the approximation.
(6.1)
and for the other upper and lower envelopes we found heuristically the approximations
(6.2)
In Figure 7 is made a superposition of the two approximations of the Figure 6 with different numbers of taken sum terms and the uniform color of the oscillations shows that for these approximations is already reached a sufficient stability of the envelopes. These envelopes are drawn in Figure 7 in the range
and they show that the approximations are very well in the neighborhood of the minima or maxima, respectively but not very well in the neighborhood of the singularities.
Figure 7. Superposition of two approximations with 15 and 30 sum terms in Fourier series of
. This superposition shows that changing the number of taken sum terms of the Fourier series does not possess an essential influence to the envelope of the oscillations. Within the range,
the conjectured envelope functions (6.2) are drawn. In the neighborhood of the singularities, they are not satisfying approximations.
Figure 8 shows two approximations for the first integral of the function
in analogy to Figure 4. The Fourier series for this function converges
already in the ordinary sense of uniform convergence to the functions that they represent.
Figure 8. Approximations of the Fourier series to function
with 3 and 6 terms. The Fourier series converges fast in ordinary sense to the given function.
7. Derivation of Identities for 2π-Periodic Generalized Functions Resting on the Poisson Summation Formula
In the following, we derive identities for 2π-periodically distributed delta functions resting on the Poisson summation formula and investigate the convergence properties of their corresponding Fourier series. From the Poisson summation formula follows the well-known basic relation (e.g., [7]-[10])
(7.1)
On the left-hand side, we have a complex-valued function and on the right-hand side a real-valued function. Therefore the imaginary part of the left-hand side has to be identically zero
(7.2)
This is also immediately seen from the antisymmetry of the Sine function
.
Using
from (7.1) follows the identity
(7.3)
with 2π-periodic delta functions on the x-axis displaced by π on the right-hand side in comparison to (7.1).
We derive now further identities from (7.1). Due to (7.2) relation (7.1) can be written
(7.4)
and relation (7.3) correspondingly
(7.5)
where we used the symmetry of the Cosine function
. Approximations of the Fourier series (7.4) on the left-hand side to the Generalized functions on the right-hand side are illustrated in Figure 9 (Fourier series (7.5) is not illustrated).
From (7.4) by substitution
we find
(7.6)
where we used a known relation for delta functions with arguments multiplied by a real factor. Furthermore, using (7.4) and (7.6)
(7.7)
The identity (7.6) can be alternatively represented in the form
(7.8)
and the identity (7.7) correspondingly
(7.9)
which are more analogous to each other than the representations (7.6) and (7.7).
We mention here in addition the well-known formula for the delta function
with a function in its argument [7]
(7.10)
where the summation goes over all simple zeros of the function
. Therefore
Figure 9. Approximations of the Fourier series
to identity (7.4). The Fourier series of the 2π-periodic delta functions show oscillations in the neighborhood of the singularities which are comparable to the oscillations in the Gibbs phenomenon in the neighborhood of the discontinuities of the functions. The envelopes to these oscillations seem to be very stabile.
(7.11)
and the right-hand sides of the above formulae and of some formulae of Section 3 can be represented with it.
Approximations of the Fourier series for the considered 2π-periodic Generalized function (7.8) are represented in Figure 10 and for the Generalized function (7.9) in Figure 11. For comparison, all these Figures are represented in a frame with the same dimensions.
The delta function
becomes infinite at
but in the sense of Generalized functions, it is not correct to speak about a function value at the singularity. In the approximations of periodic delta functions by their Fourier series, they possess finite amplitudes at all singularities which depend very much on the number of taken sum terms that is to be seen in the figures. In the neighborhood of the singularities the oscillations in the truncated Fourier series increase but independently from the number of taken sum terms they form envelopes that never go exactly to zero as one would expect from the delta function off the singularities. They even go to negative values there if we consider delta functions with positive amplitudes and to positive values if we consider delta functions with negative amplitudes. This effect is somehow similar to the Gibbs phenomenon of Fourier series at points of discontinuities of the 2π-periodic functions which they represent.
8. Representation of 2π-Periodic Function by Power-Like Series of Cosine Functions Together with a Sine Function
Using the last of the representations of
given in (2.1) we find that this
function can be expanded in a Geometric series of powers of
(with a nonsingular multiplier function, here
) as follows
(8.1)
and analogously of
according to (2.2)
(8.2)
The powers
of
are nonorthogonal to each other. This is similar to Taylor series
where the power functions
of z are also nonorthogonal to each other. In the following we show that such power series as (8.1) and (8.2) possess a faster convergence behavior as the corresponding Fourier series but this must not be true for all Fourier series.
Figure 10. Approximations of Fourier series
to identity (7.8). Apparently, the oscillations of the approximations in the neighborhood of the delta functions seem to be very stabile to both sides of the abscissa and do not vanish in the limiing case.
Figure 11. Approximations of Fourier series
to identity (7.9). Again, the oscillations in the neighborhood of the delta functions are very stabile in their envelope.
If we bring the factor
from the right-hand to the left-hand side and
or
from the left-hand side to the right-hand side in (8.1) and (8.2) follows
(8.3)
and
(8.4)
Using furthermore (5.3) we find the expansions
(8.5)
and
(8.6)
Due to the similarity to Taylor series we call such expansions power-like series of Trigonometric functions.
The convergence behavior of (8.1) and (8.2) is in contrast to the convergence behavior of its Fourier series (see Figure 2 and Figure 3) and is uniform convergence on the whole real axis. This is illustrated for the function
in Figure 12 for two numbers of sum terms taken into account. Figure 13 illustrates the convergence behavior of the power-like series for the function
which is also in full contrast to the convergence behavior of its Fourier series (see Figure 6 and Figure 7).
The term-wise integration of the series for
using (2.9) provides
(8.7)
and from the series (2.10) for
using (2.9)
(8.8)
The convergence behavior of the power-like series of
which is the first integral of
is illustrated in Figure 14 for two numbers of sum terms taken into account. The Fourier series of this functions (see Figure 4) is convergent in ordinary sense and the convergence behavior of their power-like series is similar.
Figure 12. Two approximations of power-like series
. These series possess another convergence behavior than the corresponding Fourier series and converge uniformly on the whole range of definition of the function. The differences between the two approximations with different numbers of sum terms are rather small.
Figure 13. Two approximations of power-like series
. These series possess another convergence behavior than the corresponding Fourier series and converge uniformly on the whole range of definition of the function. The differences between the two approximations with different numbers of sum terms are rather small.
We come now to the other starlet of our investigations that is the function
and to its displaced function on the real axis
and expand them in power-like series of trigonometric functions. Since we have already its first integral in (5.4) which is closely related to the above derived expressions (8.7) and (8.8) we go the reversed way and derive the other results from its series. According to (5.4) and by applying (8.7) and (8.8) we find
Figure 14. Two approximations of power-like series
. These series possess another convergence behavior than the corresponding Fourier series and converge uniformly on the whole range of definition of the function.
(8.9)
The convergence behavior of this power-like series is illustrated in Figure 15 for two values of the number of sum terms taken into account in the approximation. It is widely analogous to that in Figure 14 for the power-like series of
and is convergence in ordinary sense.
By differentiation of the expression (8.9) one obtains together with (8.3) and (8.4) the third form of a power-like series for
(8.10)
The corresponding expansion power-like expansion for
is
(8.11)
It is immediately to see that from (8.10) and (8.11) follow the power-like expansions
Figure 15. Two approximations of power-like series
. These series possess another convergence behavior than the corresponding Fourier series and converge uniformly on the whole range of definition of the function.
(8.12)
The power-like series (8.10) - (8.12) in trigonometric functions converge in ordinary sense.
In next section, we make a more general consideration of the power-like expansions in Sine and Cosine function and derive a general formula for the coefficients in these expansions.
9. Remark about Power-Like Series of Cosine and Sine Functions
In last section, we expanded some functions into power series of Sine and Cosine functions about which we make some more general remarks. Every function
can be separated in an even part
and an odd part
according to
(9.1)
In the Fourier series (2.5) with the coefficients (2.6) these two parts are separated by the even functions
and the odd functions
. This is not similar for power-like series in trigonometric functions. Whereas the powers
are even functions for arbitrary n the powers
are even functions for even
but are odd functions for odd
. Furthermore, the different powers of
and the different powers of
taken alone are, in general, no more mutually orthogonal. The last is similar to the powers in Taylor series. A complete set of symmetric and antisymmetric functions for the power-like expansion of 2π-periodic functions is formed by the two sets
. One way to get power-like expansions in Cosine and Sine functions for 2π-periodic functions is to use within the formulae for Fourier series (2.5) the following transition formulae from Cosine and Sine with multiple arguments to expansions in powers of these functions as follows
(9.2)
where
and
denotes the Chebyshev polynomials of first and second kind, correspondingly. These relations are known (e.g., Gradshteyn and Ryzhik [14]) and can be proved by complete induction. The inversion of these formulae is
(9.3)
A general procedure of the transition from Fourier series to power-like series can be made using the transition relations (9.2) and the completeness relation for expansions of 2π-periodic functions into Fourier series as explained in the next Section. The final result is an expansion of the form
(9.4)
in the two sets of nonorthogonal (symmetric) functions
and nonorthogonal (antisymmetric) functions
.
The power-like expansions in powers of Cosine and Sine functions from (8.1) - (8.10) obtained in the previous Section are of the kind (9.4) and (8.11) is a modification of them. As the cases of expansions of
and of
show where their Fourier series converge only in the sense of weak convergence their power-like series converge in the usual sense of uniform convergence and one needs only a relatively “low” number of sum term to have sufficient approximations (Figure 14 and Figure 15). Nevertheless, one has to keep in mind as already said that in the derivation of the coefficients
and
is made a substitution of the summation indices which causes a change of the order of summation in the double sums and which is only allowed if these series are absolutely convergent.
10. Completeness Relations to Fourier Series
The whole information contained in the three formulae for Fourier series (2.4), (2.5) and (2.6) can be compactly written in the following basic relation which from its kind is a completeness relation for the involved basic trigonometric functions in the basic interval
or in any equivalent interval
(10.1)
and the proof of the formulae (2.5) and (2.6) under the assumption (2.4) is equivalent to the proof of this relation. For example, the Fourier coefficients in a Fourier series of the function
follow by multiplication of (10.1) with
and integration over the basic interval
(10.2)
Also the orthogonality relations
(10.3)
follow from (10.1) inserting
and
for
in (10.2) where the last in (10.3) follows from antisymmetry of the integrand.
Inserting (9.2) into the completeness relation (10.1) leads to
(10.4)
We transform and calculate the right-hand side of this expression in detail in Appendix A. The result can be written as expansion in powers
which are symmetric functions and in functions
which are antisymmetric functions
(10.5)
with the abbreviations
(10.6)
In Section 8 we obtained expansions of the kind (9.4) in special cases more directly. The coefficients of these expansions (10.6) are related to the Generating functions of the Chebyshev polynomials with even and odd indices as the following representations show
(10.7)
In Appendix B we derive basic Generating functions for Chebyshev polynomials with even and odd indices which show the kind of functions with which we, principally, have to do it.
11. Representations of Power-Like Arithmetic Expansions by Polylogarithms
In Section 8 we found in explicit form expansions in functions
and
. These expansions can be expressed by Special Higher Transcendental functions which are called Polylogarithms [11] [12].
Polylogarithms with the notation
are defined by the series
(11.1)
The special case
is related to the Geometric series with absent first sum term
(11.2)
The next special case
is the logarithmic function
(11.3)
The special case
is called Dilogarithm and it is the following definite integral of the logarithmic function (11.3)
(11.4)
It is already not representable in a finite way using only the basic operations of addition, multiplication and division of elementary functions.
The Polylogarithms can be differentiated in a simple way according to
(11.5)
but their integration is not very simple.
In the following we provide the obtained representations of power-like expansions by Polylogarithm functions. On the right-hand side we add in addition the Fourier series of the functions. From (8.1) and (8.2) follow
(11.6)
For the power-like expansions of the function
we found in Section 8
the three different power-like representations (8.1), (8.2) and (8.10) depending on the factors in front of the relations on the right-hand side
(11.7)
They converge in an ordinary sense that is illustrated in comparison to each other in Figure 16. The last partial picture of Figure 16 represents the convergence behavior of the power-like series of the function
. Together with
which is displaced by π relatively to
one finds
Figure 16. Three different power-like series for
and for
the last in lower row with approximations by N sum terms. In the last power-like series we see a noticeable difference for approximations with even and odd sum terms. The even approximations remain on the concave side and the odd approximations on the convex side of the curves with stabile height of the hump and only their widths decrease with higher approximations.
(11.8)
An interesting phenomenon of the power-like representations which converge in ordinary sense is here that the approximations with even and odd sum terms show a different behavior. For even numbers of sum terms taken they remain at x = 0 on the concave side of the curves and for odd numbers of sum terms taken they are fully on the convex side of the curves. Only the widths of the humps change when the even or odd number of sum terms taken is changed and go to zero in the limit. The last given representations are the Fourier series which converge in the sense of the theory of Generalized functions similar to what is
Figure 17. Fourier series of
for two approximations by N sum terms. It converges in the sense of weak convergence of Generalized functions and shows that the envelope of the oscillations is here not stabile changing the number of taken sum terms. The last is likely caused because we have here in addition a linear term in “k” in its Fourier series.
Figure 18. Superposition of two approximation to Fourier series of
with
and
sum terms. It shows that the envelopes for this function depend on the sum terms in the approximation.
shown in Figure 2 and Figure 3 and Figure 6 and Figure 7. They are obtained by term-wise differentiation of the Fourier series for
and
, correspondingly. This is represented for
for two numbers of taken sum terms of the Fourier series in Figure 17. In Figure 18 is represented the superposition of two approximations with different numbers of taken sum terms in the Fourier series of
where the inhomogeneous color shows the dependence of the envelopes from its number.
One of the analogous representations to (11.7) for
is
(11.9)
With the representation of the Fourier series and of the power-like series by Polylogarithms we connected the hope that we can also determine the higher than first-order integrals of
and of
. Indeed, the program “Mathematica 10.4” provides also a result for the second-order integrals but it is so complicated and lengthy and with complex arguments of the Polylogarithms that it is hardly to use and we do not give it.
12. Conclusion
The convergence behavior of Fourier series of 2π-periodic functions which converge in the sense of weak convergence of Generalized functions was investigated and represented in figures that seems to be new. In particular, this concerns the basic functions
from which some further related functions such as their first integrals may be derived. It was found that the oscillations of the approximations are distributed in broad bands to both sides of the functions. The Fourier series of functions with weak convergence can be applied in connection with the multiplication by “sufficiently well-behaved” functions which rapidly decrease in infinity such as the Gaussian bell functions and integration over such products. In the course of calculations we came across to expansions which were called power-like expansions and which converge in ordinary sense also in cases where the Fourier series are only weakly convergent. These are expansions in powers
for even parts of functions and
for odd parts of functions. The derivation of general formulae for the coefficients of the terms in such expansion was not fully successful up to now. The obtained examples were found by application of the Geometric series. Furthermore, we investigated distributed 2π-periodic delta function with origin from the Poisson identity. In the neighborhood of the singularities (the delta functions) on the x-axis they show oscillations which also in the limiting case of their Fourier series do not vanish. This is very similar to the Gibbs phenomenon of Fourier series in the neighborhood of discontinuities of the functions. It is represented in some of the figures.
Appendix A
Calculation of the completeness relation for the power-like expansion in trigonometric functions
Starting from the completeness relation for Fourier series of 2π-periodic functions (10.1) we insert the transition relations (9.2) that leads to
(A.1)
In this relation, we make a reordering of the sum terms of the double sum by the substitution
using the formula
(A.2)
with
in our case and with the result
(A.3)
If we express in this formula
and
by the Chebyshev polynomials of first and second kind (see (9.2)) we find
(A.4)
This representation by the Chebyshev polynomials shows that we may introduce two functions which depend only on
as follows
(A.5)
with the abbreviations
(A.6)
Expression (A.5) is the new completeness relation for power-like expansions of 2π-periodic function in the sets of functions
. It is asymmetric with respect to functions of variables x and y (some duality of these functions). Formula (10.6) does not uniquely determine
which has to be set
. If we multiply (A.5) by an arbitrary function
and integrate over a basis range of length 2π one finds
(A.7)
with the coefficients
(A.8)
Altogether, the general formulae for the determination of the coefficients in power-like expansions are very difficult. All power-like expansions of 2π-periodic functions explicitly calculated in present paper were simpler determined with the help of the Geometric series than by the here-discussed general procedure.
Appendix B
Generating functions for Chebyshev polynomials with even or odd indices
The further treatment of relations (10.6) is connected with Generating functions for the Chebyshev polynomials with summation over even or odd indices, respectively to which we make here first steps. From the basic well-known Generating functions for Chebyshev polynomials
(B.1)
we derive now Generating functions for Chebyshev polynomials with even or odd indices and Cosine functions in the argument.
For Chebyshev polynomials of first kind
with even or odd indices only one finds
(B.2)
Analogously, for Chebyshev polynomials of second kind
with even or odd indices only one finds
(B.3)
One may simplify these results by the substitutions
and finds finally from (B.2)
(B.4)
and from (B.3)
(B.5)
It is easy to write these Generating functions by variables
but it is not so easy to transform them to indices with substitution
or
with arbitrary integers n where a certain number of initial sum terms has to be excluded and which would fully finish the explicit calculation of (A.6).
NOTES
1Wilbraham H., “Cambridge and Dublin Ma. J.”, 1848, vol. 3, pp. 198-201, and Gibbs J. W., “Nature”, 1898, v. 59, p. 200.
2With respect to this series Gel’fand and Shilov cite the monograph of Fikhtengol’ts, vol. III [3] (pp. 469-472), see also Vilyenkin [9], p. 463, Tolstov [13] pp. 117-121.
3The modulus of the trigonometric function in (2.9) and similar formulae is important since the Logarithm as real-valued function of real variable x is only defined for
. This means that for Trigonometric functions without modulus in the argument of the Logarithm would arise gaps of real-valued definition everywhere where these Trigonometric functions possess negative values.
4The Principal value is often denoted by the symbols “
” or “V. p.” prepended to the function.