Next Article in Journal
Modified Advanced Encryption Standard Algorithm for Information Security
Next Article in Special Issue
Ulam’s Type Stability and Generalized Norms
Previous Article in Journal
Robust CFAR Detection for Multiple Targets in K-Distributed Sea Clutter Based on Machine Learning
Previous Article in Special Issue
Type and Cotype Constants and the Linear Stability of Wigner’s Symmetry Theorem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ulam Stability for a Class of Hill’s Equations

Department of Applied Mathematics, Okayama University of Science, Okayama 700-0005, Japan
*
Author to whom correspondence should be addressed.
Submission received: 28 October 2019 / Revised: 28 November 2019 / Accepted: 1 December 2019 / Published: 5 December 2019
(This article belongs to the Special Issue Ulam's Type Stability and Symmetry)

Abstract

:
This paper deals with Ulam’s type stability for a class of Hill’s equations. In the two assertions of the main theorem, we obtain Ulam stability constants that are symmetrical to each other. By combining the obtained results, a necessary and sufficient condition for Ulam stability of a Hill’s equation is established. The results are generalized to nonhomogeneous Hill’s equations, and then application examples are presented. In particular, it is shown that if the approximate solution is unbounded, then there is an unbounded exact solution.

1. Introduction

The concept of Ulam’s type stability was posed by Ulam, and its development is remarkable. Many researchers have studied this problem for functional equations. For an overview, see the book [1] written by Brzdęk, Popa, Raşa, and Xu. In 1998, Alsina and Ger [2] introduced this concept in the field of differential equations. After that, the study of Ulam’s type stability for differential equations continued to grow (see, [3,4,5,6,7,8,9,10]). Recently, Fukutaka and Onitsuka [11,12] dealt with Ulam’s type stability of the periodic linear differential equation
y p ( x ) y = 0
on R , where p : R R is a continuous periodic function. Throughout this paper, let I = ( a , b ) , a < b . We say that (1) has “Ulam stability (US)” on I if there exists a constant K > 0 with the following property: Let ε > 0 be a given arbitrary constant, and let η : I R be a continuously differentiable function. If | η ( x ) p ( x ) η ( x ) | ε holds for all x I , then there exists a solution y : R R of (1) such that | η ( x ) y ( x ) | K ε for all x I . We call such K a “US constant” for (1) on I. In 2020, Fukutaka and Onitsuka [12] established the following theorem of Ulam’s type stability and a necessary and sufficient condition.
Theorem 1.
Let P ( x ) and Q ( x ) be antiderivatives of p ( x ) and e P ( x ) on R , respectively. Suppose that p ( x ) is a periodic function with period ω > 0 on R . Then the following hold:
(i)
if b = and 0 ω p ( x ) d x > 0 , then (1) has US with minimum US constant
max x ( 0 , ω ] lim x Q ( x ) Q ( x ) e P ( x )
on I;
(ii)
if a = and 0 ω p ( x ) d x < 0 , then (1) has US with minimum US constant
max x ( 0 , ω ] Q ( x ) lim x Q ( x ) e P ( x )
on I.
Remark 1.
Let a = and b = . That is, I = R . Suppose that p ( x ) is a periodic function with period ω > 0 on R . Then it is known that the following facts:
(i)
if 0 ω p ( x ) d x > 0 , then lim x Q ( x ) < exists (see, (Lemma 3.1 [12]));
(ii)
if 0 ω p ( x ) d x < 0 , then lim x Q ( x ) > exists (see, (Lemma 3.1 [12]));
(iii)
lim x Q ( x ) Q ( x ) e P ( x ) and Q ( x ) lim x Q ( x ) e P ( x ) are periodic functions with period ω > 0 (see, (Lemma 2.3 [11])).
Theorem 2.
Suppose that p ( x ) is a periodic function with period ω > 0 on R . Then (1) has US on R if and only if 0 ω p ( x ) d x 0 .
In recent years, Ulam stability has been actively studied not only for the first-order linear differential equations but also for the second-order linear differential equations. In 2010, Li [13] dealt with Ulam’s type stability of the simple second-order linear differential equation
y = λ 2 y ,
where λ > 0 . This study extends to more general equations with constant coefficient. For example see the works of Li and Huang [14], Li and Shen [15], and Xue [16]. On the other hand, there are many studies on the second-order linear differential equations with variable coefficients (see, [17,18,19,20,21,22,23,24]). It is well known that the most commonly encountered variable coefficient second order differential equation is Hill’s equation
y + p ( x ) y = 0 ,
where p ( x ) is a periodic function. This equation briefly describes the behavior of a large number of physical systems. For example, we can find a pendulum with moving support, electrons in a periodic potential, and beam stabilization in alternating gradient proton synchrotron (see, (Chapter 7 [25])). However, there are no studies on the Ulam stability of the second-order linear differential equations with periodic coefficients. So, this paper focuses on the stability of differential equations with periodic coefficients. The main equation in this paper is
y λ 2 ( x ) λ ( x ) y = 0
on R , where λ : R R is a continuously differentiable, periodic function with period ω > 0 . Since
λ ( x + ω ) = lim h 0 λ ( x + ω + h ) λ ( x + ω ) h = lim h 0 λ ( x + h ) λ ( x ) h = λ ( x )
for all x R , λ 2 ( x ) λ ( x ) is also a periodic function with period ω > 0 on R , and thus, this equation is a member of Hill’s equations. If λ ( x ) λ then the above equation is reduced to the equation y = λ 2 y .
The main purpose of this study is to establish a necessary and sufficient condition for Ulam stability of Hill’s equations on R . In the second section, we will establish Ulam stability of nonhomogeneous equations. In Section 3, we will give the main theorem and its proof. Also, a instability theorem is given. By using the obtained results, a necessary and sufficient condition is established. In Section 4, we will disscus the minimal US constant for the case λ ( x ) λ . In addition, a comparison with previous results is also presented. In Section 5, the concept obtained in Section 3 is extended to nonhomogeneous Hill’s equations. As an application example, we will discuss Ulam stability of Hill’s equation whose coefficient is described by Fourier series. Finally, we will conclude that the unboundedness of the approximate solution implies the unboundedness of the exact solution.

2. Ulam Stability for Nonhomogeneous Equations

In this section, we consider the nonhomogeneous first-order linear differential equation
y p ( x ) y = f ( x )
on R , where p, f : R R are continuous. We say that (2) has “Ulam stability (US)” on I if there exists a constant K > 0 with the following property: Let ε > 0 be a given arbitrary constant, and let η : I R be a continuously differentiable function. If | η ( x ) p ( x ) η ( x ) f ( x ) | ε holds for all x I , then there exists a solution y : R R of (2) such that | η ( x ) y ( x ) | K ε for all x I . We say such K a “US constant” for (2) on I. Using the previous result Theorem 1, we find the following theorem.
Theorem 3.
Let P ( x ) and Q ( x ) be antiderivatives of p ( x ) and e P ( x ) on R , respectively. Suppose that p ( x ) is a periodic function with period ω > 0 on R . Then the following hold:
(i)
if b = and 0 ω p ( x ) d x > 0 , then (2) has US with US constant
max x ( 0 , ω ] lim x Q ( x ) Q ( x ) e P ( x )
on I;
(ii)
if a = and 0 ω p ( x ) d x < 0 , then (2) has US with US constant
max x ( 0 , ω ] Q ( x ) lim x Q ( x ) e P ( x )
on I.
Proof. 
First, we will prove (i). Suppose that b = and 0 ω p ( x ) d x > 0 . Let u ( x ) be a solution of (2) on R . Then we have
ε η ( x ) p ( x ) η ( x ) f ( x ) = η ( x ) p ( x ) η ( x ) u ( x ) p ( x ) u ( x ) = η ( x ) u ( x ) p ( x ) η ( x ) u ( x )
on I. Using Theorem 1 with the above inequality, the periodicity of p ( x ) and 0 ω p ( x ) d x > 0 , we see that there exists a solution v ( x ) of (1) such that
| η ( x ) u ( x ) v ( x ) | max x ( 0 , ω ] lim x Q ( x ) Q ( x ) e P ( x ) ε
on I. Now, we consider the function y ( x ) = u ( x ) + v ( x ) for all x R . Then
y ( x ) p ( x ) y ( x ) = u ( x ) p ( x ) u ( x ) + v ( x ) p ( x ) v ( x ) = f ( x ) .
This means that y ( x ) is a solution of (2) such that
| η ( x ) y ( x ) | = | η ( x ) u ( x ) v ( x ) | max x ( 0 , ω ] lim x Q ( x ) Q ( x ) e P ( x ) ε
on I. The proof of (ii) is omitted because it can be proved by the same way.  □

3. Main Results

We consider a Hill’s equation of the form
y λ 2 ( x ) λ ( x ) y = 0 ,
on R , where λ : R R is a continuously differentiable, periodic function with period ω > 0 . We say that (3) has “Ulam stability (US)” on I if there exists a constant K > 0 with the following property: Let ε > 0 be a given arbitrary constant, and let η : I R be a twice-continuously differentiable function. If η ( x ) λ 2 ( x ) λ ( x ) η ( x ) ε holds for all x I , then there exists a solution y : R R of (3) such that | η ( x ) y ( x ) | K ε for all x I . We call such K a “US constant” for (3) on I.
When I = R , we can establish the following theorem.
Theorem 4.
Let Λ ( x ) , Γ + ( x ) and Γ ( x ) be antiderivatives of λ ( x ) , e Λ ( x ) and e Λ ( x ) on R , respectively. Suppose that λ ( x ) is a periodic function with period ω > 0 on R . Then the following hold:
(i)
if 0 ω λ ( x ) d x > 0 , then (3) has US with US constant
max x ( 0 , ω ] Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x )
on R ;
(ii)
if 0 ω λ ( x ) d x < 0 , then (3) has US with US constant
max x ( 0 , ω ] lim x Γ + ( x ) Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] Γ ( x ) lim x Γ ( x ) e Λ ( x )
on R .
Proof. 
Let I = R . That is, a = , b = . Suppose that twice-continuously differentiable function η ( x ) satisfies η ( x ) λ 2 ( x ) λ ( x ) η ( x ) ε for all x R . Define ψ ( x ) = η ( x ) + λ ( x ) η ( x ) on R . Since λ ( x ) and η ( x ) are continuously differentiable functions on R , we see that ψ : R R is a continuously differentiable function. In addition, we have the inequality
| ψ ( x ) λ ( x ) ψ ( x ) | = ( η ( x ) + λ ( x ) η ( x ) ) λ ( x ) η ( x ) + λ ( x ) η ( x ) = η ( x ) + λ ( x ) η ( x ) + λ ( x ) η ( x ) λ ( x ) η ( x ) + λ ( x ) η ( x ) = η ( x ) λ 2 ( x ) λ ( x ) η ( x ) ε
for all x R .
First, we will prove (i). From Theorem 1 (i) with b = and 0 ω λ ( x ) d x > 0 , and | ψ ( x ) λ ( x ) ψ ( x ) | ε on R , there exists a solution w ( x ) of the differential equation
w λ ( x ) w = 0
such that
| η ( x ) λ ( x ) η ( x ) w ( x ) | = | ψ ( x ) w ( x ) | max x ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x ) ε
for all x R . From Theorem 3 (ii) with a = and 0 ω ( λ ( x ) ) d x < 0 , there exists a solution z ( x ) of the differential equation
z ( x ) λ ( x ) z ( x ) w ( x ) = 0
such that
| η ( x ) z ( x ) | max x ( 0 , ω ] Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x ) ε
for all x R . Since λ ( x ) and w ( x ) are differentiable function, z ( x ) is also differentiable. Then we have
z ( x ) λ 2 ( x ) λ ( x ) z ( x ) = z ( x ) + λ ( x ) z ( x ) + λ ( x ) z ( x ) λ ( x ) ( z ( x ) + λ ( x ) z ( x ) ) = ( z ( x ) + λ ( x ) z ( x ) ) λ ( x ) ( z ( x ) + λ ( x ) z ( x ) ) = w ( x ) λ ( x ) w ( x ) = 0
for all x R , and therefore, z ( x ) is a solution of (3).
Next, we will prove (ii). From Theorem 1 (ii) with a = and 0 ω λ ( x ) d x < 0 , and | ψ ( x ) λ ( x ) ψ ( x ) | ε on R , there exists a solution w ( x ) of (4) such that
| η ( x ) λ ( x ) η ( x ) w ( x ) | = | ψ ( x ) w ( x ) | max x ( 0 , ω ] Γ ( x ) lim x Γ ( x ) e Λ ( x ) ε
for all x R . From Theorem 3 (i) with b = and 0 ω ( λ ( x ) ) d x > 0 , there exists a solution z ( x ) of (5) such that
| η ( x ) z ( x ) | max x ( 0 , ω ] lim x Γ + ( x ) Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] Γ ( x ) lim x Γ ( x ) e Λ ( x ) ε
for all x R . By the same calculation as (6), we conclude that z ( x ) is a solution of (3) on R . This completes the proof of Theorem 4.  □
Remark 2.
If we compare US constants in (i) and (ii), respectively, it can be seen that they have symmetry. Now we will show this fact below. For i { 1 , 2 } , let Λ i ( x ) , Γ i + ( x ) and Γ i ( x ) be antiderivatives of λ i ( x ) , e Λ i ( x ) and e Λ i ( x ) on R , respectively. In (i), If λ ( x ) : = λ 1 ( x ) then we can find the assumption 0 ω λ 1 ( x ) d x > 0 and US constant
K 1 : = max x ( 0 , ω ] Γ 1 + ( x ) lim x Γ 1 + ( x ) e Λ 1 ( x ) max x ( 0 , ω ] lim x Γ 1 ( x ) Γ 1 ( x ) e Λ 1 ( x ) .
If λ 1 ( x ) : = λ 2 ( x ) then this assumption and US constant imply that 0 ω λ 2 ( x ) d x < 0 and
K 1 = max x ( 0 , ω ] Γ 2 ( x ) lim x Γ 2 ( x ) e Λ 2 ( x ) max x ( 0 , ω ] lim x Γ 2 + ( x ) Γ 2 + ( x ) e Λ 2 ( x )
since Λ 1 ( x ) = Λ 2 ( x ) + c 1 , Γ 1 + ( x ) = e c 1 Γ 2 ( x ) + c 2 and Γ 1 ( x ) = e c 1 Γ 2 + ( x ) + c 3 , where c 1 , c 2 and c 3 are real constants. That is, if λ 1 ( x ) : = λ 2 ( x ) then the assumption and US constant in (i) correspond to those in (ii). Therefore, it can be concluded that the assumptions and US constants in (i) and (ii) are symmetric with each other.
Lemma 1.
Let Λ ( x ) be an antiderivative of λ ( x ) on R . If λ ( x ) is a periodic function with period ω > 0 on R , then
Λ ( x + ω ) Λ ( x ) = 0 ω λ ( x ) d x
holds on R .
Proof. 
Since λ ( x ) is continuous, and Λ ( x ) is an antiderivative of λ ( x ) , there exists a constant c R such that
Λ ( x ) = 0 x λ ( s ) d s + c
for all x R . Using this fact, we can prove this lemma. For the reference, see (Lemma 2.1 [12]). □
Next, we will present an instability theorem.
Theorem 5.
Suppose that λ ( x ) is a periodic function with period ω > 0 on R . If 0 ω λ ( x ) d x = 0 , and a = or b = , then (3) does not have Ulam stability on I.
Proof. 
Set Λ ( x ) be an antiderivative of λ ( x ) on R . Using Lemma 1 and 0 ω λ ( x ) d x = 0 , we obtain
Λ ( x + ω ) Λ ( x ) = 0 .
This says that Λ ( x ) is also a periodic function with period ω , and thus, we have
min x ( 0 , ω ] Λ ( x ) Λ ( x ) max x ( 0 , ω ] Λ ( x ) .
For any ε > 0 , we define the function
η ( x ) = ε 0 x e 2 Λ ( s ) 0 s e Λ ( u ) d u d s + c 1 0 x e 2 Λ ( s ) d s + c 2 e Λ ( x ) ,
where c 1 , c 2 R . Since this function is a solution of the equation
η ( x ) λ 2 ( x ) λ ( x ) η ( x ) = ε ,
we find | η ( x ) λ 2 ( x ) λ ( x ) η ( x ) | = ε . Now, we consider the general solution of (3) which given by y ( x ) = c 3 0 x e 2 Λ ( s ) d s + c 4 e Λ ( x ) , where c 3 and c 4 are arbitrary constants. Namely, y ( x ) means all solutions of (3). From this, we obtain
| η ( x ) y ( x ) | = ε 0 x e 2 Λ ( s ) 0 s e Λ ( u ) d u d s + c 1 0 x e 2 Λ ( s ) d s + c 2 e Λ ( x ) c 3 0 x e 2 Λ ( s ) d s + c 4 e Λ ( x ) = ε 0 x e 2 Λ ( s ) 0 s e Λ ( u ) d u d s + ( c 1 c 3 ) 0 x e 2 Λ ( s ) d s + c 2 c 4 e Λ ( x ) = 0 x e 2 Λ ( s ) 0 s e Λ ( u ) d u + c 1 c 3 ε e 2 Λ ( s ) d s + c 2 c 4 ε e Λ ( x ) ε
for all x R .
Now, we consider the case b = . By (7), we have
0 x e Λ ( s ) d s e max x ( 0 , ω ] Λ ( x ) x 0 and e 2 max x ( 0 , ω ] Λ ( x ) e 2 Λ ( x ) e 2 min x ( 0 , ω ] Λ ( x ) > 0
for all x 0 . This implies,
e 2 Λ ( x ) 0 x e Λ ( s ) d s + c 1 c 3 ε e 2 Λ ( x ) e 2 min x ( 0 , ω ] Λ ( x ) max x ( 0 , ω ] Λ ( x ) x c 1 c 3 ε e 2 max x ( 0 , ω ] Λ ( x )
for all x 0 . From this, we see that
0 x e 2 Λ ( s ) 0 s e Λ ( u ) d u + c 1 c 3 ε e 2 Λ ( s ) d s + c 2 c 4 ε e 2 min x ( 0 , ω ] Λ ( x ) max x ( 0 , ω ] Λ ( x ) 2 x 2 c 1 c 3 ε e 2 max x ( 0 , ω ] Λ ( x ) x + c 2 c 4 ε
for all x 0 , and there exists a x 1 > max { 0 , a } such that
e min x ( 0 , ω ] Λ ( x ) 2 max x ( 0 , ω ] Λ ( x ) 2 x 2 c 1 c 3 ε e 2 max x ( 0 , ω ] Λ ( x ) x + c 2 c 4 ε > x
for all x x 1 . Therefore, using (7) and (8), we obtain
| η ( x ) y ( x ) | > x e Λ ( x ) ε x e max x ( 0 , ω ] Λ ( x ) ε
for all x x 1 , so that, lim x | η ( x ) y ( x ) | = . Hence, (3) does not have US on I = ( a , ) . The proof of the case a = is omitted as it can be proved in the same way. □
Remark 3.
Let J be an open interval satisfying J I . Then Ulam stability on I implies the same on J. On the other hand, instability on J implies the same on I.
Corollary 1.
Suppose that λ ( x ) is a periodic function with period ω > 0 on R . If 0 ω λ ( x ) d x = 0 , then (3) does not have Ulam stability on R .
Example 1.
Consider the case λ ( x ) = cos x for x I . If I = R , then (3) does not have US on R by using Corollary 1.
Theorem 4 and Corollary 1 give the following necessary and sufficient condition.
Theorem 6.
Suppose that λ ( x ) is a periodic function with period ω > 0 on R . Then (3) has US on R if and only if 0 ω λ ( x ) d x 0 .
Example 2.
Consider the case
λ ( x ) = α + β cos x
on R , where α , β R . Then (3) has US on R if and only if α 0 , by using Theorem 6.

4. Minimal US Constant for the Case λ ( x ) λ

Corollary 2.
If λ ( x ) λ R \ { 0 } , then (3) has US with US constant 1 / λ 2 on R .
Proof. 
Let Λ ( x ) = λ x , Γ + ( x ) = e λ x / λ , Γ ( x ) = e λ x / λ on R . Then we have
0 ω λ d x = λ ω > 0
if λ > 0 , and
0 ω λ d x = λ ω < 0
if λ < 0 . Thus, we have
max x ( 0 , ω ] Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) ] max x ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x ) = max x ( 0 , ω ] e λ x λ e λ x max x ( 0 , ω ] e λ x λ e λ x = 1 λ 2
if λ > 0 , and
max x ( 0 , ω ] lim x Γ + ( x ) Γ + ( x ) e Λ ( x ) ] max x ( 0 , ω ] Γ ( x ) lim x Γ ( x ) e Λ ( x ) = max x ( 0 , ω ] e λ x λ e λ x max x ( 0 , ω ] e λ x λ e λ x = 1 λ 2
if λ < 0 . From Theorem 4, we conclude that (3) has US with US constant 1 / λ 2 . □
Lemma 2.
Let λ ( x ) is a continuously differentiable function on R , and let b = . If
0 < sup x I λ 2 ( x ) λ ( x ) < ,
and (3) has US on I, then minimum US constant on I is at least 1 / sup x I λ 2 ( x ) λ ( x ) .
Proof. 
Let ε > 0 be a given arbitrary constant, and let η ( x ) be a solution of the differential equation
η ( x ) λ 2 ( x ) λ ( x ) η ( x ) = ε
on I. Since (3) has US, there exists a constant K 0 > 0 such that η ( x ) y ( x ) K 0 ε on I. Now, let y ( x ) be any solution of (3) such that the above inequality holds. We will show that K 0 1 / sup x I λ 2 ( x ) λ ( x ) . By way of contradiction, we assume that K 0 < 1 / sup x I λ 2 ( x ) λ ( x ) . From this assumption, we have
η ( x ) y ( x ) = λ 2 ( x ) λ ( x ) η ( x ) y ( x ) + ε λ 2 ( x ) λ ( x ) η ( x ) y ( x ) + ε λ 2 ( x ) λ ( x ) K 0 ε + ε sup x I λ 2 ( x ) λ ( x ) K 0 ε + ε = 1 sup x I λ 2 ( x ) λ ( x ) K 0 sup x I λ 2 ( x ) λ ( x ) ε
on I. Set
M = 1 sup x I λ 2 ( x ) λ ( x ) K 0 sup x I λ 2 ( x ) λ ( x ) ε > 0 ,
and set x 0 I . Then we obtain
η ( x ) y ( x ) x 0 x M d s + η ( x 0 ) y ( x 0 ) = M x M x 0 + η ( x 0 ) y ( x 0 )
for all x > x 0 . From this inequality, we have
η ( x ) y ( x ) x 0 x M s M x 0 + η ( x 0 ) y ( x 0 ) d s + η ( x 0 ) y ( x 0 ) = 1 2 M x 2 M x 0 η ( x 0 ) + y ( x 0 ) x 1 2 M x 0 2 + M x 0 η ( x 0 ) + y ( x 0 ) x 0 + η ( x 0 ) y ( x 0 )
for all x > x 0 . Using this inequality, we see that lim x η ( x ) y ( x ) = . This contradicts the inequality η ( x ) y ( x ) K 0 ε on I. □
Theorem 7.
If λ ( x ) = λ R \ { 0 } on R , then (3) has Ulam stability with minimum US constant 1 / λ 2 on R .
Proof. 
Let λ ( x ) = λ R \ { 0 } on R . Using Corollary 2, (3) has US and 1 / λ 2 is a US constant on R . From Lemma 2 with sup x R λ 2 ( x ) λ ( x ) = λ 2 , we see that the minimum US constant for (3) on R is at least 1 / λ 2 . This completes the proof of Theorem 7. □
Remark 4.
The decisive difference between the previous study described in the introduction and this study is whether the interval is finite or not. In the previous studies [13,14,15,16], they assume that < a < b < . If I = ( a , b ) = ( , ) = R , then Ulam stability cannot be judged by their results because the US constant in the proof of previous theorems depends on a or b, and, goes to ∞ when a or b . However, our study established Ulam stability on an infinite interval, and obtained the minimum US constant.

5. Nonhomogeneous Hill’s Equations and Applications

Consider the nonhomogeneous Hill’s equation
y λ 2 ( x ) λ ( x ) y = f ( x ) ,
where f : R R is continuous. Equation (9) has “Ulam stability (US)” on I if and only if there exists K > 0 such that the following holds: Let ε > 0 be given, and let η : I R be a twice-continuously differentiable function. If η ( x ) λ 2 ( x ) λ ( x ) η ( x ) f ( x ) ε on I, then there exists a solution y ( x ) of (9) such that | η ( x ) y ( x ) | K ε on I. K is called a “US constant” for (9) on I. The following result is obtained by using Theorem 4.
Theorem 8.
Let Λ ( x ) , Γ + ( x ) and Γ ( x ) be antiderivatives of λ ( x ) , e Λ ( x ) and e Λ ( x ) on R , respectively. Suppose that λ ( x ) is a periodic function with period ω > 0 on R . Then the following hold:
(i)
if 0 ω λ ( x ) d x > 0 , then (9) has US with US constant
max x ( 0 , ω ] Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x )
on R ;
(ii)
if 0 ω λ ( x ) d x < 0 , then (9) has US with US constant
max x ( 0 , ω ] lim x Γ + ( x ) Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] Γ ( x ) lim x Γ ( x ) e Λ ( x )
on R .
Proof. 
First, we will prove (i). Let u ( x ) be a solution of (9) on R . Suppose that 0 ω λ ( x ) d x > 0 and η ( x ) λ 2 ( x ) λ ( x ) η ( x ) f ( x ) ε on I. Then we have
ε η ( x ) λ 2 ( x ) λ ( x ) η ( x ) f ( x ) = η ( x ) λ 2 ( x ) λ ( x ) η ( x ) u ( x ) λ 2 ( x ) λ ( x ) u ( x ) = η ( x ) u ( x ) λ 2 ( x ) λ ( x ) η ( x ) u ( x )
on I. Using Theorem 4 with the above inequality, we see that there exists a solution v ( x ) of (3) such that
| η ( x ) u ( x ) v ( x ) | K 2 ε
on I, where
K 2 : = max x ( 0 , ω ] Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x ) .
Let y ( x ) = u ( x ) + v ( x ) for all x R . Then
y ( x ) λ 2 ( x ) λ ( x ) y ( x ) = u ( x ) λ 2 ( x ) λ ( x ) u ( x ) + v ( x ) λ 2 ( x ) λ ( x ) v ( x ) = f ( x ) .
Thus, y ( x ) is a solution of (9) such that
| η ( x ) y ( x ) | = | η ( x ) u ( x ) v ( x ) | K 2 ε
on I. The proof of (ii) is omitted because it can be proved by the same way. □
Theorem 8 and Corollary 1 imply the following result.
Theorem 9.
Suppose that λ ( x ) is a periodic function with period ω > 0 on R . Then (9) has US on R if and only if 0 ω λ ( x ) d x 0 .
US constants in Theorems 4 and 8 can be estimated as follows.
Lemma 3.
Let Λ ( x ) , Γ + ( x ) and Γ ( x ) be antiderivatives of λ ( x ) , e Λ ( x ) and e Λ ( x ) on R , respectively. Suppose that λ ( x ) is a periodic function with period ω > 0 on R . Then the following hold:
(i)
if λ ( x ) > 0 on R , then
1 max x ( 0 , ω ] λ ( x ) 2 max x ( 0 , ω ] Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) max t ( 0 , ω ] lim x Γ ( x ) Γ ( x ) e Λ ( x ) 1 min x ( 0 , ω ] λ ( x ) 2 ;
(ii)
if λ ( x ) < 0 on R , then
1 min x ( 0 , ω ] λ ( x ) 2 max x ( 0 , ω ] lim x Γ + ( x ) Γ + ( x ) e Λ ( x ) max x ( 0 , ω ] Γ ( x ) lim x Γ ( x ) e Λ ( x ) 1 max x ( 0 , ω ] λ ( x ) 2 .
Proof. 
We will prove (i). For simplicity, let
g ( x ) : = Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) ,
and
h ( x ) : = lim x Γ ( x ) Γ ( x ) e Λ ( x ) .
Then we have
g ( x ) = Γ + ( x ) lim x Γ + ( x ) e Λ ( x ) λ ( x ) g ( x ) = 1 λ ( x ) g ( x )
and
h ( x ) = lim x Γ ( x ) Γ ( x ) e Λ ( x ) + λ ( x ) h ( x ) = 1 + λ ( x ) h ( x ) .
From Remark 1, g ( x ) and h ( x ) are periodic functions with period ω > 0 , so that we can find x 1 ( 0 , ω ] such that g ( x 1 ) = max x ( 0 , ω ] g ( x ) . Since g ( x 1 ) is a maximum value of g ( x ) on ( 0 , ω ] , we see that
0 = g ( x 1 ) = 1 λ ( x 1 ) g ( x 1 ) = 1 λ ( x 1 ) max x ( 0 , ω ] g ( x ) .
That is, max x ( 0 , ω ] g ( x ) = 1 / λ ( x 1 ) . From this and λ ( x ) > 0 is a periodic function with period ω > 0 on R , we have
0 < 1 max x ( 0 , ω ] λ ( x ) max x ( 0 , ω ] g ( x ) = 1 λ ( x 1 ) 1 min x ( 0 , ω ] λ ( x )
on ( 0 , ω ] . Let x 2 R satisfy h ( x 2 ) = max x ( 0 , ω ] h ( x ) , then we have
0 = h ( x 2 ) = 1 + λ ( x 2 ) h ( x 2 ) = 1 + λ ( x 2 ) max x ( 0 , ω ] h ( x ) .
That is, max x ( 0 , ω ] h ( x ) = 1 / λ ( x 2 ) , and so that, we find
0 < 1 max x ( 0 , ω ] λ ( x ) max x ( 0 , ω ] h ( x ) = 1 λ ( x 2 ) 1 min x ( 0 , ω ] λ ( x )
on ( 0 , ω ] . Using (10) and (11), we obtain
1 max x ( 0 , ω ] λ ( x ) 2 max x ( 0 , ω ] g ( x ) max x ( 0 , ω ] h ( x ) 1 min x ( 0 , ω ] λ ( x ) 2 .
The proof of (ii) is omitted because it can be proved by the same way. □
Theorem 8 and Lemma 3 imply the following result.
Theorem 10.
Suppose that λ ( x ) is a periodic function with period ω > 0 on R . Then the following hold:
(i)
if λ ( x ) > 0 on R , then (9) has US with US constant 1 / min x ( 0 , ω ] λ ( x ) 2 ;
(ii)
if λ ( x ) < 0 on R , then (9) has US with US constant 1 / max x ( 0 , ω ] λ ( x ) 2 .
In many application fields, the coefficient p ( x ) of the Hill’s equation y + p ( x ) y = 0 is treated as a periodic function with period π . It is well known that p ( x ) can be rewritten by the Fourier series of p ( x ) . That is, Hill’s equation is given by
y + a 0 2 + n = 1 a n cos 2 n x + b n sin 2 n x y = 0
where a 0 , a n and b n are Fourier coefficients. Now, we consider the case
λ ( x ) = α + β cos 2 x + γ sin 2 x ,
in (3), where α , β and γ are real-valued constants. Then we have
λ 2 ( x ) λ ( x ) = ( α + β cos 2 x + γ sin 2 x ) 2 ( 2 β sin 2 x + 2 γ cos 2 x ) = α 2 + 2 ( α β γ ) cos 2 x + 2 ( α γ + β ) sin 2 x + ( β cos 2 x + γ sin 2 x ) 2 = α 2 + 2 ( α β γ ) cos 2 x + 2 ( α γ + β ) sin 2 x + β 2 1 + cos 4 x 2 + β γ sin 4 x + γ 2 γ 2 1 + cos 4 x 2 = 2 α 2 + β 2 + γ 2 2 + 2 ( α β γ ) cos 2 x + β 2 γ 2 2 cos 4 x + 2 ( α γ + β ) sin 2 x + β γ sin 4 x .
This implies that this case corresponds to the case:
a 0 = ( 2 α 2 + β 2 + γ 2 ) , a 1 = 2 ( α β γ ) , a 2 = β 2 γ 2 2 , a n = 0 ( i = 3 , 4 , 5 , ) ,
b 1 = 2 ( α γ + β ) , b 2 = β γ and b n = 0 ( i = 3 , 4 , 5 , ) .
Let η ( x ) be a solution of (12). Then it is an approximate solution of the nonhomogeneous Hill’s equation
y + a 0 2 + n = 1 a n cos 2 n x + b n sin 2 n x y = f ( x ) , | f ( x ) | ε
where a 0 , a n and b n are the same as above, and ε > 0 is given. In fact,
η ( x ) + a 0 2 + n = 1 a n cos 2 n x + b n sin 2 n x η ( x ) f ( x ) = | 0 f ( x ) | ε
holds on R . Since
0 π λ ( x ) d x = α π
holds, (13) has US on R if and only if α 0 by Theorem 9.
Next we consider the more simple case: α = 2 , β = γ = 1 . Then we have
λ ( x ) = 2 + cos 2 x + sin 2 x and λ 2 ( x ) λ ( x ) = 5 6 cos 2 t 2 sin 2 t + sin 4 t
and λ ( x ) < 0 for all x R . We will show that if the approximate solution is unbounded, then there is an unbounded exact solution. Using Theorem 10, we see that there exists a solution y ( x ) of (13) satisfying | η ( x ) y ( x ) | K 3 ε on R , where
K 3 : = 1 max x ( 0 , ω ] λ ( x ) 2 = 1 2 + 2 2 .
Since
η 1 ( t ) = e Λ ( x ) = exp 2 x 1 2 sin 2 x + 1 2 cos 2 x
is an unbounded solution of (12), we conclude that (13) has at least one unbounded solution.

6. Conclusions

In this paper, we have discussed the Ulam stability of periodic linear differential equations of the first order and second order, respectively. First, we have established a sufficient condition for Ulam stability of nonhomogeneous first order equations. Next, by using this result, Ulam stability of a class of Hill’s equations is established. In particular, the exact US constant is given in the theorem. If the coefficient is a constant, it turns out to be the minimum US constant. On the other hand, the instability theorem is also given. By combining the obtained results, we have established a necessary and sufficient condition for Ulam stability of a class of Hill’s equations. Using the idea of the results obtained first, this theory is extended to nonhomogeneous Hill’s equations. Finally, we have discussed the Ulam stability of Hill’s equations where the coefficient is described by a Fourier series. By using the obtained results, it can be concluded that the unboundedness of the approximate solution affects the unboundedness of the exact solution.

Author Contributions

These authors contributed equally to this work; conceptualization, R.F. and M.O.; methodology, R.F. and M.O.; investigation, R.F. and M.O.; resources, M.O.; writing–original draft preparation, R.F. and M.O.; writing–review and editing, R.F. and M.O.; supervision, R.F.; project administration, M.O.; funding acquisition, M.O.

Funding

The second author was supported by JSPS KAKENHI Grant Number JP17K14226.

Acknowledgments

These authors would like to thank the referee for carefully reading our manuscript and for giving useful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brzdęk, J.; Popa, D.; Raşa, I.; Xu, B. Ulam stability of Operators. Mathematical Analysis and Its Applications; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  2. Alsina, C.; Ger, R. On some inequalities and stability results related to the exponential function. J. Inequal. Appl. 1998, 2, 373–380. [Google Scholar] [CrossRef]
  3. Onitsuka, M. Hyers–Ulam stability of first order linear differential equations of Carathéodory type and its application. Appl. Math. Lett. 2019, 90, 61–68. [Google Scholar] [CrossRef]
  4. Onitsuka, M.; Shoji, T. Hyers–Ulam stability of first-order homogeneous linear differential equations with a real-valued coefficient. Appl. Math. Lett. 2017, 63, 102–108. [Google Scholar] [CrossRef] [Green Version]
  5. Popa, D.; Raşa, I. On the Hyers–Ulam stability of the linear differential equation. J. Math. Anal. Appl. 2011, 381, 530–537. [Google Scholar] [CrossRef] [Green Version]
  6. Popa, D.; Raşa, I. Hyers–Ulam stability of the linear differential operator with nonconstant coefficients. Appl. Math. Comput. 2012, 219, 1562–1568. [Google Scholar] [CrossRef]
  7. Rezaei, H.; Jung, S.-M.; Rassias, T.M. Laplace transform and Hyers–Ulam stability of linear differential equations. J. Math. Anal. Appl. 2013, 403, 244–251. [Google Scholar] [CrossRef]
  8. Wang, G.; Zhou, M.; Sun, L. Hyers–Ulam stability of linear differential equations of first order. Appl. Math. Lett. 2008, 21, 1024–1028. [Google Scholar] [CrossRef] [Green Version]
  9. Zada, B. Uniform exponential stability in the sense of Hyers and Ulam for periodic time varying linear systems. Differ. Equ. Appl. 2018, 10, 227–234. [Google Scholar] [CrossRef]
  10. Zada, A.; Wang, P.; Lassoued, D.; Li, T. Connections between Hyers–Ulam stability and uniform exponential stability of 2-periodic linear nonautonomous systems. Adv. Differ. Equ. 2017, 192, 7. [Google Scholar] [CrossRef] [Green Version]
  11. Fukutaka, R.; Onitsuka, M. Best constant in Hyers–Ulam stability of first-order homogeneous linear differential equations with a periodic coefficient. J. Math. Anal. Appl. 2019, 473, 1432–1446. [Google Scholar] [CrossRef]
  12. Fukutaka, R.; Onitsuka, M. A necessary and sufficient condition for Hyers–Ulam stability of first-order periodic linear differential equations. Appl. Math. Lett. 2020, 100, 106040. [Google Scholar] [CrossRef]
  13. Li, Y. Hyers–Ulam stability of linear differential equations y″ = λ2y. Thai J. Math. 2010, 8, 215–219. [Google Scholar]
  14. Li, Y.; Huang, J. Hyers–Ulam stability of linear second-order differential equations in complex Banach spaces. Electron. J. Differ. Equ. 2013, 184, 7. [Google Scholar]
  15. Li, Y.; Shen, Y. Hyers–Ulam stability of linear differential equations of second order. Appl. Math. Lett. 2010, 23, 306–309. [Google Scholar] [CrossRef] [Green Version]
  16. Xue, J. Hyers–Ulam stability of linear differential equations of second order with constant coefficient. Ital. J. Pure Appl. Math. 2014, 32, 419–424. [Google Scholar]
  17. Abdollahpour, M.R.; Park, C. Hyers–Ulam stability of a class of differential equations of second order. J. Comput. Anal. Appl. 2015, 18, 899–903. [Google Scholar]
  18. Alqifiary, Q.H.; Jung, S.-M. On the Hyers–Ulam stability of differential equations of second order. Abstr. Appl. Anal. 2014, 2014, 483707. [Google Scholar] [CrossRef]
  19. Anderson, D.R. Hyers–Ulam stability of higher-order Cauchy–Euler dynamic equations on time scales. Dyn. Syst. Appl. 2014, 23, 653–664. [Google Scholar]
  20. Ghaemi, M.B.; Gordji, M.E.; Madjid, E.; Alizadeh, B.; Park, C. Hyers–Ulam stability of exact second-order linear differential equations. Adv. Differ. Equ. 2012, 2012, 36. [Google Scholar] [CrossRef] [Green Version]
  21. Javadian, A.; Sorouri, E.; Kim, G.H.; Gordji, M.E. Generalized Hyers–Ulam stability of the second-order linear differential equations. J. Appl. Math. 2011, 2011, 813137. [Google Scholar] [CrossRef] [Green Version]
  22. Shen, Y. An integrating factor approach to the Hyers–Ulam stability of a class of exact differential equations of second order. J. Nonlinear Sci. Appl. 2016, 9, 2520–2526. [Google Scholar] [CrossRef] [Green Version]
  23. Shen, Y.; Chen, W.; Lan, Y. On the Ulam stability of a class of Banach space valued linear differential equations of second order. Adv. Differ. Equ. 2014, 2014, 294. [Google Scholar] [CrossRef] [Green Version]
  24. Zhao, X.; Zhang, X.; Ge, W. Hyers–Ulam stability of a class second differential equation y″(x) + p(x)y′(x) + q(x)y(x) = F(y(x)). Bull. Malays. Math. Sci. Soc. 2017, 40, 891–906. [Google Scholar] [CrossRef]
  25. Harris, C.J.; Miles, J.F. Stability of Linear Systems: Some Aspects of Kinematic Similarity. Mathematics in Science and Engineering; Academic Press, Inc.: London, UK, 1980. [Google Scholar]

Share and Cite

MDPI and ACS Style

Fukutaka, R.; Onitsuka, M. Ulam Stability for a Class of Hill’s Equations. Symmetry 2019, 11, 1483. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11121483

AMA Style

Fukutaka R, Onitsuka M. Ulam Stability for a Class of Hill’s Equations. Symmetry. 2019; 11(12):1483. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11121483

Chicago/Turabian Style

Fukutaka, Ryuma, and Masakazu Onitsuka. 2019. "Ulam Stability for a Class of Hill’s Equations" Symmetry 11, no. 12: 1483. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11121483

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop