There is perhaps no better way starting the discussion than quoting directly from Iserles and Nørsett  : “By any yardstick, Fourier series are one of the greatest and most influential concepts of contemporary mathematics. … It is thus with a measure of trepidation and humility that we wish to pursue a variation upon the Fourier theme in this paper.” Since trigonometric series was first used by d’Alembert in 1747, the full formation of Fourier theories surprisingly took more than a century of endeavors highlighted by the famous d’Alembert-Euler-Bernoulli controversy and many important and/or pioneering contributions from Euler, Dirichlet, Lagrange, Lebesgue and other leading mathematicians of the time. Nevertheless, it was not long before mathematicians and scientists came to appreciate the power and far-reaching implications of Fourier’s claim that any function could be expanded into a trigonometric series. Fourier’s discovery is easily ranked in the “top ten” mathematical advances of all time.
Despite what has been said, the Fourier series will lose much of its luster when used to expand a sufficiently smooth nonperiodic function defined on a compact interval. It is well known that a continuous function can always be expanded into a Fourier series inside the interval (the word “inside” is highlighted to emphasize the fact that the two end points shall not be automatically included). This is actually the primary reason for the inefficiency of the Fourier series in approximating a nonperiodic function, and, understandably, in solving various boundary value problems. This work is aimed at overcoming the said difficulties associated with the conventional Fourier series.
It is known that a continuous function f(x) defined on the interval [−π, π] can always be expanded into a Fourier series
where the expansion coefficients are calculated from
The Fourier series, (1.1), reduces to
if f(x) is an odd function;
if f(x) is an even function.
The convergence of the Fourier series, (1.1), is well understood through the following theorems.
THEOREM 1. If is an absolutely integrable piecewise smooth function of period of 2π, then the Fourier series of converges to at points of continui-
ty and to at points of discontinuity. If is continuous everywhere, then the series converges absolutely and uniformly.
Proof. Pages 75-78 of Ref.  .
THEOREM 2. For any absolutely integrable function, its Fourier coefficients satisfy
Proof. Pages 70-71 of Ref.  .
THEOREM 3. Let be a continuous function of period 2π, which has n derivatives, where n − 1 derivatives are continuous and the n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points). Then, the Fourier series of all n derivatives can be obtained by term-by-term differentiation of the Fourier series of, where all the series, except possibly the last, converge uniformly to the corresponding derivatives. Moreover, the Fourier coefficients of the function satisfy the relations
Proof. Pages 84, 130, and 131 of Ref.  .
As a matter of fact, (1.7) can be replaced by more explicit expressions  
where is the partial sum of the Fourier series defined as
The aforementioned convergence theorems are established based on the condition that f(x) is a periodic function of period 2π. It is known that the Fourier series of an analytic 2π-periodic function can actually converge at an exponential rate  . However, once the periodicity condition is removed, the convergence of a series expansion can be seriously deteriorated or even there is no convergence in the maximum norm. When f(x) is defined only on a compact interval [−π, π], it can be viewed as the part of the 2π-periodic function which is the periodic extension of f(x) onto the whole x-axis. Thus, even f(x) is sufficiently smooth on [−π, π], the extended periodic function may only be piece-wise smooth due to the potential discontinuities at (). As a consequence, the series expansion of f(x) converges to f(x) for every x Î (−π, π), and to at x = ±π. Understandably, such a Fourier expansion converges very slowly.
Assume, for example, that is continuous on [−π, π] with an absolutely integrable derivative (which may not exist at certain points). Then we have
where am and bm are the Fourier coefficients of and. Since
the Fourier coefficients of an absolutely integrable function tend to zero as (Theorem 2), it is obvious that
If, then we have
which recovers the convergence rate for a continuous 2π-periodic function. Unfortunately, the condition, , is generally not true for an arbitrary function.
In recognizing this slow convergence problem, the subtraction of polynomials has been developed to remove the Gibbs phenomenon with (or its related derivatives) and to thus accelerate the convergence of resulting Fourier expansions  -  . In polynomial subtraction schemes, a new (or corrected) function will be created with a desired smoothness through removing the potential jumps, such as, at the end points
where is a polynomial of degree 2K + 1, satisfying
The polynomials can be easily constructed, for example, using the Lanczos’s system of polynomials:
Lanczos polynomials of even (odd) degrees are obviously even (odd) functions. It should be noted that Lanczos polynomials are closely related to Bernoulli polynomials which are also widely used in the methods of polynomial subtraction.
The first few Lanczos polynomials can be explicitly expressed as
For complete Fourier expansion of,
For the sine expansion of,
For the cosine expansion of:
Assume that is Cn−1 continuous on [−π, π] and its n-th derivative is absolutely integrable. Then the corrected function can be periodically extended into a 2π-periodic function of: a) CK continuity for K ≤ n − 1 in (1.12); b) C2K+1 continuity for 2K + 1 < n in (1.13) and c) C2K+2 continuity for 2K + 2 < n in (1.14).
By recognizing the slower convergence of sine series than its cosine counterpart, a modified Fourier series was proposed as 
If is differentiable and its derivative has bounded variation, the expansion coefficients, am and bm, in (1.15) will both decay like   , which is still considered relatively slow in many cases.
2. An Alternative Form of Fourier Cosine Series
For a sufficiently smooth function defined on a compact interval [0, p], it can always be expanded into
It is known that expansion coefficients am decay like.
To accelerate the convergence and maintain a close similarity to classical Fourier series, an alternative trigonometric expansion of is here sought in the form of
and coefficients bp are to be determined as described below.
THEOREM 4. Let have Cn−1continuity on the interval [0, π] and its n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points). If n ³ 2, then the Fourier coefficient am, as defined in (2.4), decays at a polynomial rate as
Proof. By integrating by part, we have
for sufficiently large m.
Substracting (2.9) from (2.8) leads to
The first two terms in (2.10) vanish if
In order to have a unique and smallest set of coefficients, bp, we set Q = P in (2.11) and (2.12), or equivalently, in (2.6) and (2.7). The convergence estimate, (2.5), becomes evident from (2.10) according to Theorem 2. W
Remark. If, the relationship, (2.5), in Theorem 4 can be modified to
or, more explicitly,
Alternatively, (2.4) can be expressed as
Equations (2.6) and (2.7) can be rewritten in matrix form as
Determination of the coefficients, B1 and B2, involves the inversion of a Vandermonde-like matrix
which is always invertable if xk ¹ xj for j ¹ k.
Consider a polynomial of degree 2P − 1
Then it is obvious that
where δij is Kronecker’s symbol.
According to (2.28), matrix C = [cik] is actually the inverse of matrix X.
To find an explicit expression for matrix C, let
Thus, we have
Comparing (2.32) with (2.27) leads to
In light of (2.34), the coefficients bp () can be obtained as
By making use of (2.21), the first few coefficients, for example, are readily found as:
for P = 2;
for P = 3.
EXAMPLE 1. Consider function. Its conventional Fourier expansions are easily obtained as
Under the current framework, this function can be expanded as:
for P = 1;
for P = 2;
for P = 3.
A graphic display of the results, (2.48), (2.49), (2.51), (2.53) and (2.55), is given in Figure 1 for A = 4, B = 2, and C = 1. The corresponding truncation errors are plotted in Figure 2 and Figure 3.
Figure 1. Decays of expansion coefficients: bm in (2.48), am in (2.49), am in (2.51), am in (2.53), and am in (2.55).
3. An Alternative Form of Fourier Sine Series
Similarly, can also be expanded into sine series:
where bp are the expansion coefficients to be determined, and
THEOREM 5. Let have Cn−1 continuity on the interval [0, π] and the n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points). Then for the Fourier coefficients of defined in (3.1) decay at a polynomial rate as
Proof. By integrating by part, we have
The first two terms in (3.6) will both vanish if
Figure 2. Truncation errors, , for the series expansions: (2.48), (2.49), (2.51), (2.53), and (2.55). M = 20.
Figure 3. Errors, , for series expansion (2.51): M = 10, M = 20 and M = 40.
In order to have a unique and smallest set of coefficients, bp, we set Q = P in (3.7) and (3.8), or equivalently, in (3.4) and (3.5). The convergence estimate, (3.3), then becomes evident according to Theorem 2. W
Remark. If, the relationship (3.3) can be modified to
which can be further written in a shaper form as
The expansion coefficients, am, can be alternatively expressed as
(see (2.16)). (3.13)
We can rewrite (3.4) and (3.5) in matrix form as
Following the same procedures as described earlier, coefficients bp can be obtained from
Using this formula, the first several coefficients are easily determined as:
for P = 1;
for P = 2;
for P = 3.
EXAMPLE 2. Consider function. The classical sine series expansion of this function is easily found as
In the context of the current framework, this function can be expressed as
for P = 1;
for P = 2.
4. An Alternative Form of Fourier Series Expansion
Let be defined on the interval [−π, π]. It can also be expanded into a complete trigonometric series as
where am and bm are the expansion coefficients to be calculated from
and and are to be determined as follows.
THEOREM 6. Let have Cn−1 continuity on the interval [−π, π] and the n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points). Then the Fourier coefficients of, as defined in (4.2)and (4.3), decay at a polynomial rate as
Proof. Function can be considered as the superposition of an even function and an odd function. Theorems 4 is then directly applicable to g(x) on [0, π]. Thus, (4.4) holds if
(4.10) and (4.11) can be rewritten as (4.6) and (4.7), respectively.
Similarly, relationship (4.5) is readily obtained from applying Theorem 5 to the odd function h(x) on interval [0, π] by recognizing that
The expansion coefficients of g(x) are determined from
Similarly, the expansion coefficients of h(x) are determined from
The even (odd) extension of () onto [−π, 0) will lead to an even (odd) function
()on [−π, π]. Expression (4.1) will then become evident. W
Alternatively, (4.2) and (4.3) can be expressed as
where is given in (2.16).
The coefficients and are readily calculated from (2.35) and (3.19), respectively, by letting
EXAMPLE 3. Consider function. Its classical Fourier expansion is easily found as
By setting P = 0 and Q = 1 in (4.1), we have
It is seen from (4.25) that the sine series now converges at a rate of m−3 which is faster than m2 for its cosine counterpart. If desired, the convergence of the series expansion in the form of (4.1) can be further accelerated by setting P = Q = 1. Accordingly, in addition to (4.23), we have
The series expansion given in (4.27) will converge at a rate of m−3 in comparison with m−2 for that in (4.24).
COROLLARY 1. Let have Cn−1continuity on the interval [0, π] and the n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points). Assume n ³ 2. Then can be expanded as
where and are calculated from(2.4) and (2.35), respectively.
Proof. For, denote
or, alternatively, (5.3) and (5.4).
Then expansion (5.1) follows immediately from (2.3) in view that
Since for, (5.2) is evident from (2.5). W
COROLLARY 2. Let have Cn−1 continuity on the interval [0, π] and the n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points). Then for, can be expanded as
provided that Am and satisfy (5.3) and (5.4), and and are calculated from (3.2) and (3.19), respectively.
Proof. By (5.5) and (5.6), we have
Thus, (5.7) and (5.8) become obvious from (3.1) and (3.3), respectively. W
COROLLARY 3. Let have Cn−1 continuity on the interval [−π, π] and the n-th derivative is absolutely integrable (the n-th derivative may not exist at certain points. Then for and, can be expanded as
where, , and are defined in the same way as those in (4.1).
Since Corollary 3 is obvious from Theorem 6 and Corollaries 1 and 2, its proof will not be given here.
Notice that in (5.1)
and (superposed * indicates taking complex conjugate).
Remark. In Corollary 1, (5.1) can be alternatively written as
Similarly, (5.7) in Corollary 2 can be written as
And (5.9) in Corollary 3 as
Alternative Fourier series expansions have been presented in an effort of better representing a sufficiently smooth function in a compact interval. The series expansions can take various forms, resulting in different rates of convergence. When one of the series expansions, for example, is used to solve a boundary value problem, its convergence rate needs to be compatible with the smoothness of the solution “physically” dictated by the problem. Thus, there may exist the best form for any given problem. Among other important applications, the new Fourier series will potentially lead to a new path for solving differential equations and boundary value problems.
 Iserles, A. and Norsett, S.P. (2008) From High Oscillation to Rapid Approximation I: Modified Fourier Expansions. IMA Journal of Numerical Analysis, 28, 862-887.
 Gottlieb, D. and Orszag, S.A. (1977) Numerical Analysis of Spectral Methods: Theory and Applications. SIAM, Philadelphia.
 Tadmor, E. (1986) The Exponential Accuracy of Fourier and Chebyshev Differencing Methods. SIAM Journal on Numerical Analysis, 23, 1-10.
 Jones, W.B. and Hardy, G. (1970) Accelerating Convergence of Trigonometric Approximations. Mathematics of Computation, 24, 547-560.
 Baszenski, G., Delvos, F.J. and Tasche, M. (1995) A United Approach to Accelerating Trigonometric Expansions. Computers & Mathematics with Applications, 30, 33-49.
 Lyness, J.N. (1974) Computational Techniques Based on the Lanczos Representation. Mathematics of Computation, 28, 81-123.
 Eckhoff, K.S. (1993) Accurate and Efficient Reconstruction of Discontinuous Functions from Truncated Series Expansions. Mathematics of Computation, 61, 745-763.
 Eckhoff, K.S. (1995) Accurate Reconstructions of Functions of Finite Regularity from Truncated Fourier Series Expansions. Mathematics of Computation, 64, 671-690.
 Eckhoff, K.S. (1998) On a High Order Numerical Method for Functions with Singularities. Mathematics of Computation, 67, 1063-1087.
 Olver, S. (2009) On the Convergence Rate of a Modified Fourier Series. Mathematics of Computation, 78, 1629-1645.
 Adcock, B. (2011) Convergence Acceleration of Modified Fourier Series in One or More Dimensions. Mathematics of Computation, 80, 225-261.