1. Introduction
The introduction of a relaxing parameter
in differential equations was found to provide stability properties for their corresponding solutions. This is a phenomenon well-known in numerical analysis where if the Ordinary Differential Equation (ODE)
is
stiff then one can try to use the
backward Euler method to obtain the sequence
by first considering the algebraic equations
for small time-steps
. If one can obtain
explicitly in terms of
then the iteration scheme often converges much faster, and for longer time intervals, than that provided by the
forward Euler method [
1], p. 349. That such a principle holds for ODEs as
was established through the study of Multiplicatively Advanced Differential Equations (MADEs) as
, and will be discussed further in this article. Part of our analysis of stability will require obtaining uniform apriori bounds. This will be achieved in a somewhat general setting, and the consequences will be presented in the form of examples of advanced differential equations.
1.1. Solutions of MADEs as Eigenfunctions
In [
2] solutions to equations of the form
were studied for
,
,
and
. In the case that
, with
, solutions
are referred to as
eigenfunctions since
as
. Specific asymptotic properties of solutions were obtained in Theorem 10 of [
3]. Here we only consider the case that
and
, however the derivatives may be of higher (integer) order than in Equation (
1). In addition, we extend solutions of these equations to all
so that the eigen equation, referred to as an eigen-MADE, has a solution
the Schwartz space of infinitely differentiable functions, with derivatives that decay faster than reciprocal polynomials (as defined in [
4] section V.3). An asymptotic theory near
can be developed indicating that an extension to
is quite natural. In this way the special functions that we study are eigenfunctions in
, although not in the traditional, local (
) sense. The significance of these functions will be demonstrated by examples, and convergence to familiar functions is obtained on compact subsets of
, as
.
1.2. Brief Overview
The study of multiply advanced differential equations falls within the area of functional differential equations, as is studied for instance in Fox, et al. [
2], Kato, et al. [
3] and Dung [
5]. There is also significant overlap with the area of
q-difference differential equations, where the multiplicative advancement
is referred to as a dilation and is denoted
. There is a rich and active study within the area of
q-difference differential equations with dilations involving
. These are highlighted by works of: L. Di Vizio [
6,
7,
8]; C. Hardouin [
7]; T. Dreyfus [
9,
10]; A. Lastra [
10,
11,
12,
13,
14,
15,
16,
17,
18,
19]; S. Malek [
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22]; J. Sanz [
17,
18,
19]; H. Tahara [
23]; and C. Zhang [
8,
24]; along with further references by these researchers and others. Often these studies in
q-difference differential equations overlap with the area of Gevrey asymptotics.
In the current work we continue by focusing on global solutions of a MADE on
. In particular, we discuss several techniques for starting with a given global solution to an original MADE and then generating solutions of new related MADEs. This theme will be developed as follows: In
Section 2, a known MADE solution first introduced in [
25], namely
, is used to produce a simple related solution
which is an eigensolution of a MADE in the sense of the Abstract. In turn,
is then used to obtain a new
q-advanced Airy function
satisfying a MADE analogue of the Airy differential equation. Then
itself is used along with convolution to generate families of functions
solving a
q-advanced PDE.
In
Section 3, a family of MADE solutions, under convolution and auto-correlation, are seen to produce related solutions of new MADEs. Furthermore, the least-element method in Poincare asymptotics is deployed to find natural extensions to related MADE solutions on the negative real line. A theory of asymptotic extensions to
is developed to clarify the notion that solutions to MADEs behave smoothly in a neigborhood of the origin. We also give conditions that ensures a natural extension to all of
, as is needed to even consider a Fourier transform. An investigation of the inhomogeneous MADEs that these solve is begun.
In
Section 4 we focus on considering solutions of MADEs as perturbations of classical solutions, and, mirroring a more direct convergence proof in
Section 2, we exhibit MADE solutions which converge to a classical solution of a damped-oscillation equation—the convergence being uniform on compact subsets of
.
In
Section 5, we return to the topics of convolution and auto-correlation to observe their impact when applied to MADE solutions. In this paper, we will discuss convolutions, correlations, and Fourier transforms for MADEs.
A table of Fourier transforms of global MADE solutions under study here is provided in
Section 6. These will be solutions of new MADEs, for which we obtain new elements in a table of Fourier transforms. This new table mimics what is often done for Laplace transforms, in the study of linear constant coefficient ODEs.
In various theories of differential equations, convolutions provide a useful tool since general solutions can be determined from fundamental solutions, as demonstrated here in Equation (
33). This is one motivation for obtaining solutions to homogenous equations, as appears in Proposition 2.
2. A Normalized Cosine Example and Extensions
From [
25], consider the following Schwartz functions, for
and all
,
where
There are several properties that we note. In particular, the function
is normalized, in that the uniform bound
holds, after some delicate work performed in [
25], for each
. It also solves the following eigen-MADE for all
and each
,
From (
6) we see that
satisfies an eigen-MADE in the sense of the Abstract, with
independently of the advancing parameter
. Note that
(as recorded in (
10) below) does not have an eigenvalue
independent of
q, thus we rely on
as the appropriate eigen-MADE solution.
Since
is not only
and bounded, but in fact Schwartz, we can obtain its Fourier transform, an operation defined for any
, as
In [
25] it was found that
where
was defined in Equation (
4) above, and the other normalizing constant is
To express the Fourier transform of linear, homogeneous MADEs, we found multiple uses of the Jacobi theta function
which allows the association that
, and which ensures that
for all
, due to the product formula. It will be of significance to note that the reciprocal
, for
, is Schwartz when extended to be identically 0 for
. Critical algebraic properties that we use are
A consequence is that the only zeros of
are for
for all
. This is obvious from the product definition of
in Equation (
8).
2.1. Uniform Convergence
Using Taylor series methods as an approach paralleling that in [
25] we show:
Proposition 1. On any compact subset of , approaches uniformly as .
Proof. A given compact set is contained in an interval for sufficiently large, so it suffices to prove the theorem on .
First, recall the following results shown in [
25]
From these, by induction on the even order derivatives of
, we obtain the higher order derivatives
and
We infer all derivatives of
via
and
Evaluating the derivatives of
at
yields
for all
.
Next computing
, the
degree Taylor polynomial for
expanded about
, gives
with remainder term
for appropriate
between 0 and
t. Using the sup norm
, along with the fact that
, to bound from above, we obtain
Let
and
denote the
degree Taylor polynomial and remainder terms for
respectively. Then, for each
and each
t with
, one has
Now, given any
choose
such that
. Then one has
. Next choose
with
. Then for all
one has
Next choose
such that
. Then for all
one has
For the given
, set
in (
19) and (20). Then for
and all
, applying the bounds (
21) and (
22) to (20) gives
verifying uniform convergence of
to
on
as
. □
Remark 1. Note that, alternatively, one can express Proposition 1 asA similar convergence proof is given in Section 4, with details related to the novelty of the result. 2.2. Application to PDE Example
We are now in a position to obtain
q-versions of various equations using
as a building block for relaxing equations. For example, define the Airy function (see page 570 in [
26])
Some properties of this function are that as , and . We now show:
Proposition 2. The q-advanced Airy function is defined here to befor . The functions and satisfy the homogeneous ODE and MADErespectively, for . Basic properties of for , are that is Schwartz with . Furthermore, for each , , and sufficiently large, so thatIn other words, uniformly for t in compact subsets of , as . Remark 2. Verifying convergence in Equation (27) may seem rather straight forward, due to the uniform convergence of to on compact sets. However, we need to use a careful argument, as demonstrated here. Proof. That
is Schwartz follows from the same for
, whereas the property
requires a manipulation of theta functions, and is shown in
Appendix A. We start with the second equation in (
26) since the first equation is known to hold [
26]. First define the function
Now compute, using
, and
, for
,
Next, to show convergence, consider any
and, without loss of generality, fix
. Let
t be in the interval
. Then, for any
, using integration by parts and boundedness of the sine function, we can write
Thus, for all
we can easily find
sufficiently large so that
The bound in Equation (
30) also holds if
is replaced with
since
and
for all
. Now, fix
sufficiently large so that the bounds in (
30), and also (
30) with cos replaced by
, are less than
. It is essential to note that this value of
R is independent of
.
Finally, for each
, define the function
The union of these
over
, is the interval
. From the uniform convergence in Equation (
24) we can choose
so that
for
,
, and
. This is now sufficient to verify the expression in Equation (
27). □
2.3. A q-Advanced PDE Example
The argument in the proof of Proposition 2 shows that knowledge of one MADE can help to generate and study other MADEs. In fact, this extends to Partial Differential Equations (PDEs). For example, consider the linear constant-coefficient Airy PDE [
27]
for
,
, and constant
. To obtain an advanced-type equation, consider the kernel function, defined for each
,
for appropriate
, to be determined. For any integrable
and any
, define,
(compare with Equation (2.2) of [
27]). Recall that the functional operation of convolution for integrable functions
gives a new function
defined by
where the last equality in Equation (
35) is the Convolution Theorem (see [
28] Theorem IX.4). To discover the PDE that
solves, first compute the
t-partial derivative of Equation (
34), to obtain
Now, taking three derivatives of Equation (
34) with respect to
x, gives
By replacing
in Equation (38), one can verify that the
q-advanced PDE, for
and
,
holds. To obtain consistency with the initial data
, first define the constant, for
,
which is finite since
is Schwartz (thus integrable) for
. Then we require,
The q-Airy Hypothesis: Given
, the expression in Equation (
40) does not vanish, i.e.,
.
In
Appendix B we show that the
q-Airy Hypothesis holds for all
. Then
where convergence in Equation (
41) is pointwise, and is shown in
Appendix C using a mollifier-type argument. If, in addition, we have
and
, then convergence in Equation (
41) becomes uniform.
4. Convergence of MADEs to Classical Solutions
In this section, we present another example where we can study convergence of a MADE solution to its classical analogue. This requires an apriori uniform bound in a fixed neighborhood of
for all
sufficiently small. Obtaining a uniform-in-
q bound for general
is rather deep, and complicated by the presence of the alternation
in Equation (
42). Here we study a series without this alternating factor, which defines a function that behaves like a damped oscillation. The details are more challenging than what appears in the proof of Proposition 1, so a full analysis is provided.
Consider the following linear third-order MADE
for
, on the interval
, satisfying the initial conditions
For small
, as
, Equations (
80) and (
81) can be considered to be a perturbation of the classical analogue, which is the ODE
with initial conditions
obtained by setting
in (
80) and (
81). One can check directly that (
82) and (
83) is solved uniquely by
Now, using techniques mirroring those of Theorem 3.2 of [
29], a particular solution to (
80) is
for
. Note that the expression in Equation (
85) does not have the alternation
, unlike the expression in Equation (
55) for
, and this will allow a sharp bound on
for all
, independent of
.
The first derivative of
is seen to be
where the fact that:
was used explicitly to obtain the last equality in (
86). Using this identity implicitly, we obtain:
and finally we verify:
A re-indexing
was used to move from (
88) to (89). Note that (89) gives that (
80) holds. From (
85)–(
87), one sees that
where the last equalities of (91) and (92) are obtained from (
8).
Normalizing
by
to obtain
one sees that
now satisfies the MADE (
80) along with the initial conditions (
81). The last initial condition follows from the fact that
where the last equality in (
94) follows from the next lemma.
Lemma 1. For the Jacobi theta function (8) satisfies Proof. For the first equality in (
95) one can write
where the second equality is obtained from Equation (
9) with
, and the last equality is the reciprocal identity in Equation (
9) with
. Dividing (
96) by
gives (
95). For the second equality in (
95), let
and
in Equation (
9). Then as above
. The lemma is shown. □
In addition to the last equality of (
94) being proven by the first equality in (
95) in Lemma 1, the second equality of (
95) proves that the second derivative of
at
equals
.
The following theta function bound will also be helpful.
Lemma 2. For the Jacobi theta function (8) satisfies Proof. Comparing each factor in the square brackets in (
98) with the corresponding factor in the square brackets of (
99) one sees that for all
from which one concludes that
giving the left inequality in (
97). The right inequality in (
97) holds via the assumption that
. □
Next we compute all derivatives of
at
and of
at
, in preparation for the computation of the Taylor series expansion at
for both
and
. From (
82) we immediately have that for
and
From (
102) and (
83) one concludes that for
The analogous results for are obtained in the following lemma.
Lemma 3. For and , let with given by (85). Then for and one has Furthermore, at one has Proof. We first establish (
104) for the case that
by induction on
k. So for
note that (
104) holds as a tautology for
, and for
it holds by (89). Assume that
for fixed
k. Then
where: the inductive hypothesis gives the rightmost equality in (
106), and that (89) along with the chain rule gives the first equality in (107). Thus, the
case holds for all
k. Now differentiate the expression
either
or
times to obtain (
104) in all remaining cases. Evaluating (
104) at
and relying on (
90)–(
94) gives (
105). □
Next, the
-degree Taylor polynomials
of
g and
f, respectively, expanded about
are given by
where (
108) follows from (
103), and (109) follows from (
105). For
, these have respective remainder terms
for some
and
. The goal of uniform convergence on compact subsets is now obtained in the following proposition.
Proposition 4. Let S be any compact set contained in . Then converges uniformly to on S as , where is given by both (93) and (85), while is given by (84). Proof. Without loss of generality, there is a
such that
, and it is sufficient to prove uniform convergence on
. For
, from the triangle inequality one has
Now for
and relying on (111), one starts with (114) to see
where: moving to (115) is obtained via (
93); (116) follows from (
85) and (91); the equality in (118) is obtained by (
8); and the inequality in (119) is given by (
97) in Lemma 2. Similarly, from (
110) and (
84), one has
Also, from (
108) and (109) if we let:
, then
Applying (118), (
121), and (
120) to (113) one has that for
Now, given
, choose
sufficiently large such that one has
. Then
Then for
one has
and
whence for
and
Applying (
125) to (
122) with
N taken to be
one has that for
So approaches uniformly on as , and the proposition is proven. □
6. Expanded Table of Fourier Transforms
In this final section we establish a short table of Fourier transforms for solutions of MADEs and their relations to Jacobi theta functions. Included are well-established results, along with new functions. The positive constants and are generic, but estimates are not presented here.
The introduction of new functions are as follows: For
see [
32] for decay constants
and
in
Table 1; The functions
and
are closely related to
and
, respectively, introducted in [
25], where constants
and
are obtained; The
q-Bessel functions, related to
, were introduced in [
33], along with decay constants
and
; Flat wavelets
have Fourier transforms that are averages of theta functions, first derived in [
29], along with constants
and
. The functions
and
, have Fourier transforms that involve theta functions, which can be used to obtain decay parameters
and
.
Note that similar tables for Laplace transforms are quite extensive, since applications only require control of function growth on . Here we are concerned with globally defined functions on for which a Fourier transform can be defined.