Next Article in Journal
Fixed Points of g-Interpolative Ćirić–Reich–Rus-Type Contractions in b-Metric Spaces
Next Article in Special Issue
Singularly Perturbed Cauchy Problem for a Parabolic Equation with a Rational “Simple” Turning Point
Previous Article in Journal
Deep Reinforcement Learning Agent for S&P 500 Stock Selection
Previous Article in Special Issue
Axiomatic Approach in the Analytic Theory of Singular Perturbations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regularization Method for Singularly Perturbed Integro-Differential Equations with Rapidly Oscillating Coefficients and Rapidly Changing Kernels

by
Burkhan Kalimbetov
1,* and
Valeriy Safonov
2
1
Department of Mathematics, Akhmed Yassawi University, B. Sattarkhanov 29, Turkestan 161200, Kazakhstan
2
Department of Higher Mathematics, National Research University «MPEI», Krasnokazarmennaya 14, 111250 Moscow, Russia
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2020 / Revised: 9 November 2020 / Accepted: 10 November 2020 / Published: 13 November 2020

Abstract

:
In this paper, we consider a system with rapidly oscillating coefficients, which includes an integral operator with an exponentially varying kernel. The main goal of the work is to develop an algorithm for the regularization method for such systems and to identify the influence of the integral term on the asymptotic behavior of the solution of the original problem.

1. Introduction

In the study of various issues related to dynamic stability, with the properties of media with a periodic structure, in the study of other applied problems, one has to deal with differential equations with rapidly oscillating coefficients. Equations of this kind can describe some mechanical or electrical systems that are under the influence of high-frequency external forces, automatic control systems with a linear adjustable object, etc. As an example, we can cite the principle of operation of an oscillator with a small mass and a nonlinear restoring force, in which a high-frequency periodic force with a large amplitude acts. The presence of high-frequency terms creates serious problems for their direct numerical solutions. Therefore, asymptotic methods are usually applied to such equations first, the most famous of which are the Feshchenko–Shkil–Nikolenko splitting method [1,2,3,4,5] and the Lomov’s regularization method [6,7,8]. It should also be noted that singularly perturbed equations are the object of study by several Russian researchers, as well as other scientists (see, for example [9,10,11,12,13,14,15,16,17,18,19,20,21,22]).
In this paper, the Lomov’s regularization method is generalized to previously unexplored integro-differential equations with rapidly oscillating coefficients and with rapidly decreasing kernels of the form
ε d z d t a ( t ) z ε g ( t ) cos β ( t ) ε z t 0 t e 1 ε s t μ ( θ ) d θ K ( t , s ) z ( s , ε ) d s = h ( t ) , z ( t 0 , ε ) = z 0 , t [ t 0 , T ]
where z = z ( t , ε ) , h ( t ) , β ( t ) > 0 , a ( t ) > 0 , μ ( t ) < 0 , a ( t ) μ ( t ) t t 0 , T , g ( t ) are scalar functions, z 0 is a constant, ε > 0 is a small parameter. In the case β t = 2 γ t , and of the absence of an integral term, such a system was considered in [6,7,8].
The limit operator a ( t ) has a spectrum λ 1 t = a ( t ) , functions λ 2 t = i β t and λ 3 t = + i β t are associated with the presence in Equation (1) of a rapidly oscillating cos β ( t ) ε , and the function λ 4 t = μ ( t ) characterizes the rapid change in the kernel of the integral operator.
We introduce the following notations:
λ t = λ 1 t , , λ 4 t ,
m = m 1 , , m 4 is multi-index with non-negative components m j , j = 1 , 4 ¯ ,
m = j = 1 4 m j is multi-index height m ,
m , λ t = j = 1 4 m j λ j t .
Assume that the following conditions are met:
( 1 ) a ( t ) , β ( t ) , μ ( t ) C t 0 , T , R , g ( t ) , h ( t ) C t 0 , T , C , K ( t , s ) C t 0 s t T , C ;
(2) the relations m , λ t = 0 , m , λ t = λ j t , j 1 , , 4 for all multi-indices m with m 2 or are not fulfilled for any t t 0 , T , or are fulfilled identically on the whole segment t t 0 , T .
In other words, resonant multi-indices are exhausted by the following sets
Γ 0 = m : m , λ t 0 , m 2 , t t 0 , T , Γ j = m : m , λ t λ j t , m 2 , t t 0 , T , j = 1 , 4 ¯ .
Under these conditions, we will develop an algorithm for constructing a regularized [6] asymptotic solution of the problem (1).

2. Regularization of the Problem (1)

Denote by σ j = σ j ε independent of the t quantities σ 1 = e i ε β t 0 , σ 2 = e + i ε β t 0 , and rewrite the Equation (1) in the form
L z ( t , ε ) ε d z d t a ( t ) z ε g ( t ) 2 e i ε t 0 t β ( θ ) d θ σ 1 + e + i ε t 0 t β ( θ ) d θ σ 2 z t 0 t e 1 ε s t μ ( θ ) d θ K ( t , s ) z ( s , ε ) d s = h ( t ) , z ( t 0 , ε ) = z 0 , t [ t 0 , T ] .
We introduce regularizing variables
τ j = 1 ε t 0 t λ j ( θ ) d θ ψ j ( t ) ε , j = 1 , 4 ¯
and instead of problem (2) we consider the problem
L z ˜ ( t , τ , ε ) ε z ˜ t + j = 1 4 λ j ( t ) z ˜ τ j λ 1 ( t ) z ˜ ε g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z ˜ t 0 t e 1 ε s t λ 4 ( θ ) d θ K ( t , s ) z ˜ s , ψ ( s ) ε , ε d s = h ( t ) , z ˜ ( t , τ , ε ) | t = t 0 , τ = 0 = z 0 , t [ t 0 , T ]
for the function z ˜ = z ˜ t , τ , ε , where it is indicated (according to (3)): τ = τ 1 , , τ 4 , ψ = ψ 1 , , ψ 4 . It is clear that if z ˜ = z ˜ t , τ , ε is the solution of the problem (4), then the function z = z ˜ t , ψ ( t ) ε , ε is an exact solution of the problem (2), therefore, the problem (4) is an extension of the problem (2).
However, (4) cannot be considered completely regularized, since the integral term
J z ˜ = t 0 t e 1 ε s t λ 4 ( θ ) d θ K t , s z ˜ s , ψ ( s ) ε , ε d s
has not been regularized in it. To regularize J, we introduce a class M ε , asymptotically invariant with respect to the operator J z ˜ (see [6]; p. 62).
We first consider the space U of functions z t , τ , representable by sums
z t , τ , σ = z 0 t , σ + i = 1 4 z i t , σ e τ i + 2 m N z z m t , σ e m , τ , z 0 t , σ , z i t , σ , z m t , σ C t 0 , T , C , i = 1 , 4 ¯ , 2 m N z
where the asterisk ∗ above the sum sign indicates that in it the summation for m 2 occurs only over nonresonant multi-indices m = m 1 , , m 4 , i.e., over m i = 0 4 Γ i .
Note that in (5) the degree N z of the polynomial z t , τ , σ to exponentials e τ j depends on the element z . The elements of the space U depend on bounded in ε > 0 constants σ 1 = σ 1 ε and σ 2 = σ 2 ε , which do not affect the development of the algorithm described below, therefore in the notation of element (5) of this space U we omit the dependence on σ = σ 1 , σ 2 for brevity. We show that the class M ε = U | τ = ψ t / ε is asymptotically invariant with respect to the operator J .
The image of the operator J on the element (5) of the space U has the form:
J z t , τ = t 0 t e 1 ε s t λ 4 ( θ ) d θ K t , s z 0 s d s + i = 1 4 t 0 t e 1 ε s t λ 4 ( θ ) d θ K t , s z i s e 1 ε t 0 s λ i ( θ ) d θ d s +
+ 2 m N z t 0 t e 1 ε s t λ 4 ( θ ) d θ K t , s z m s e 1 ε t 0 s m , λ ( θ ) d θ d s =
= t 0 t e 1 ε s t λ 4 ( θ ) d θ K t , s z 0 s d s + e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z 4 s d s +
+ i = 1 , i 4 4 e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z i s e 1 ε t 0 s λ i θ λ 4 θ d θ d s +
+ 2 m N z e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z m s e 1 ε t 0 s m e 4 , λ θ d θ d s .
Integrating in parts, we have
J 0 t , ε = t 0 t K t , s z 0 s e 1 ε t 0 s λ 4 ( θ ) d θ d s = ε t 0 t K t , s z 0 s λ 4 s d e 1 ε t 0 s λ 4 ( θ ) d θ =
= ε K t , s z 0 s λ 4 s e 1 ε t 0 s λ 4 ( θ ) d θ s = t 0 s = t ε t 0 t s K t , s z 0 s λ 4 s e 1 ε t 0 s λ 4 ( θ ) d θ d s =
= ε K t , t z 0 t λ 4 t e 1 ε t 0 t λ 4 ( θ ) d θ K t , t 0 z 0 t 0 λ 4 t 0 ε t 0 t s K t , s z 0 s λ 4 s e 1 ε t 0 s λ 4 ( θ ) d θ d s .
Continuing this process further, we obtained the decomposition
J 0 t , ε = ν = 0 ( 1 ) ν ε ν + 1 I 0 ν K t , s z 0 s s = t e 1 ε t 0 t λ 4 ( θ ) d θ I 0 ν K t , s z 0 s s = t 0 , I 0 0 = 1 λ 4 s · , I 0 ν = 1 λ 4 s s I 0 ν 1 ( ν 1 ) .
Next, apply the same operation to the integrals:
J 4 , i t , ε = e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z i s e 1 ε t 0 s λ i θ λ 4 θ d θ d s = = ε e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z i s λ i s λ 4 s d e 1 ε t 0 s λ i θ λ 4 θ d θ =
= ε e 1 ε t 0 t λ 4 ( θ ) d θ K t , s z i s λ i s λ 4 s e 1 ε t 0 s λ i θ λ 4 θ d θ s = t 0 s = t
ε t 0 t s K t , s z i s λ i s λ 4 s e 1 ε t 0 s λ i θ λ 4 θ d θ d s =
= ε K t , t z i t λ i t λ 4 t e 1 ε t 0 t λ i θ d θ K t , t 0 z i t 0 λ i t 0 λ 4 t 0 e 1 ε t 0 t λ 4 θ d θ ε e 1 ε t 0 t λ 4 θ d θ t 0 t s K t , s z i s λ i s λ 4 s e 1 ε t 0 s λ i θ λ 4 θ d θ d s =
= ν = 0 ( 1 ) ν ε ν + 1 I i ν K t , s z i s s = t e 1 ε t 0 t λ i θ d θ I i ν K t , s z i s s = t 0 e 1 ε t 0 t λ 4 θ d θ , I i 0 = 1 λ i s λ 4 s · , I i ν = 1 λ i s λ 4 s s I i ν 1 , ν 1 , i = 1 , 3 ¯ .
Denote bay e 4 = ( 0 , 0 , 0 , 1 ) . Then
J m t , ε = e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z m s e 1 ε t 0 s m e 4 , λ θ d θ d s =
= ε e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z m s m e 4 , λ s d e 1 ε t 0 s m e 4 , λ θ d θ =
= ε e 1 ε t 0 t λ 4 ( θ ) d θ [ K t , s z m s m e 4 , λ s e 1 ε t 0 s m e 4 , λ θ d θ s = t 0 s = t t 0 t s K t , s z m s m e 4 , λ s e 1 ε t 0 s m e 4 , λ θ d θ d s ] =
= ν = 0 ( 1 ) ν ε ν + 1 [ I 4 , m ν K t , s z m s s = t e 1 ε t 0 t m , λ θ d θ I 4 , m ν K t , s z m s s = t 0 e 1 ε t 0 t λ 4 ( θ ) d θ ] , I 4 , m 0 = 1 m e 4 , λ s · , I 4 , m ν = 1 m e 4 , λ s s I 4 , m ν 1 , ν 1 , 2 m N z .
Here it is taken into account that m e 4 , λ s 0 , since by the definition of the space U multi-indices m Γ 4 . The image of the operator J on the space U element (5) is represented as a series
J z t , τ = e 1 ε t 0 t λ 4 ( θ ) d θ t 0 t K t , s z 4 s d s + ν = 0 ( 1 ) ν ε ν + 1 I 0 ν K t , s z 0 s s = t e 1 ε t 0 t λ 4 ( θ ) d θ
I 0 ν K t , s z 0 s s = t 0 + i = 1 , i 4 4 ν = 0 ( 1 ) ν ε ν + 1 I i ν K t , s z i s s = t e 1 ε t 0 t λ i ( θ ) d θ
I i ν K t , s z i s s = t 0 e 1 ε t 0 t λ 4 ( θ ) d θ +
+ 2 m N z ν = 0 ( 1 ) ν ε ν + 1 I 4 , m ν K t , s z m s s = t e 1 ε t 0 t m , λ θ d θ
I 4 , m ν K t , s z m s s = t 0 e 1 ε t 0 t λ 4 ( θ ) d θ , τ = ψ t / ε .
It is easy to show (see, for example, [23], pp. 291–294) that this series converges asymptotically for ε + 0 (uniformly in t [ t 0 , T ] ). This means that the class M ε is asymptotically invariant (for ε + 0 ) with respect to the operator J .
Let as introduce the operators R ν : U U , acting on each element z t , τ U of the form (5) according to the law:
R 0 z t , τ = e τ 4 t 0 t K t , s z 4 s d s ,
R 1 z t , τ = I 0 0 K t , s z 0 s s = t e τ 4 I 0 0 K t , s z 0 s s = t 0 +
+ i = 1 3 I i 0 K t , s z i s s = t e τ i I i 0 K t , s z i s s = t 0 e τ 4 +
+ 2 m N z I 4 , m 0 K t , s z m s s = t e ( m , τ ) I 4 , m 0 K t , s z m s s = t 0 e τ 4 ,
R ν + 1 z t , τ = I 0 ν K t , s z 0 s s = t e τ 4 I 0 ν K t , s z 0 s s = t 0 +
+ i = 1 3 ( 1 ) ν I i ν K t , s z i s s = t e τ i I i ν K t , s z i s s = t 0 e τ 4 +
+ 2 m N z I 4 , m ν K t , s z m s s = t e ( m , τ ) I 4 , m ν K t , s z m s s = t 0 e τ 4 , ν 1 .
Let now z ˜ t , τ , ε be an arbitrary continuous function in t , τ t 0 , T × τ : R e τ j 0 , j = 1 , 4 ¯ with the asymptotic expansion
z ˜ t , τ , ε = k = 0 ε k z k t , τ , z k t , τ U ,
converging as ε + 0 (uniformly in t , τ t 0 , T × τ : R e τ j 0 , j = 1 , 4 ¯ ). Then the image J z ˜ t , τ , ε of this function is expanded in the asymptotic series
J z ˜ t , τ , ε = k = 0 ε k J z k t , τ = r = 0 ε r s = 0 r R r s z s t , τ | τ = ψ t / ε .
This equality is the basis for introducing the extension of the operator J on the series type (7):
J ˜ z ˜ t , τ , ε J ˜ k = 0 ε k z k t , τ = d e f r = 0 ε r s = 0 r R r s z s t , τ .
Although the operator J ˜ is formally defined, its usefulness is obvious, since in practice they usually construct the N-th approximation of the asymptotic solution of problem (2), in which only the N-th partial sums of the series (7) will take part, which do not have a formal but true meaning. Now we can write down a problem that is completely regularized with respect to the original problem (2):
L z ˜ ( t , τ , ε ) ε z ˜ t + j = 1 4 λ j t z ˜ τ j λ 1 ( t ) z ˜ ε g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z ˜ J ˜ z ˜ = h ( t ) , z ˜ ( t , τ , ε ) | t = t 0 , τ = 0 = z 0 , t [ t 0 , T ] .

3. Iterative Problems and Their Solvability in the Space U

Substituting series (7) into (8) and equating the coefficients for the same powers ε , we obtain the following iterative problems:
L z 0 t , τ j = 1 4 λ j t z 0 τ j λ 1 ( t ) z 0 R 0 z 0 = h t , z 0 t 0 , 0 = z 0 ;
L z 1 t , τ = z 0 t + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z 0 + R 1 z 0 , z 1 t 0 , 0 = 0 ;
L z 2 t , τ = z 1 t + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z 1 + R 1 z 1 + R 2 z 0 , z 0 t 0 , 0 = 0 ;
L z k t , τ = z k 1 t + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z k 1 + R k z 0 + + R 1 z k 1 , z k t 0 , 0 = 0 , k 1 .
Each of the iterative problems can be written as
L z t , τ j = 1 4 λ j t z τ j λ 1 ( t ) z R 0 z = H t , τ , z t 0 , 0 = z ,
where H t , τ = H 0 t + i = 1 4 H i t e τ i + 2 m N H H m t e m , τ is the known function of the space U , z is the known number of complex the space C , and the operator R 0 has the form (see 6 0 )
R 0 z R 0 z 0 t + i = 1 4 z i t e τ i + 2 m N z z m t e m , τ = d e f e τ 4 t 0 t K t , s z 4 s d s .
We introduce the scalar product (for each t t 0 , T ) in the space U :
< z , w > < z 0 t + i = 1 4 z i t e τ i + 2 m N z z m t e m , τ ,
w 0 t + i = 1 4 w i t e τ i + 2 m N w w m t e m , τ > = d e f
= d e f z 0 t , w 0 t + i = 1 4 z i t , w i t + 2 m min N z , N w z m t , w m t ,
where (∗, ∗) we denote the ordinary scalar product in the complex space C : ( u , v ) = u · v ¯ . We prove the following statement.
Theorem 1.
Suppose that conditions (1) and (2) are satisfied and the right-hand side H t , τ = H 0 t + + i = 1 4 H i t e τ i + 2 m N H H m t e m , τ of the Equation (10) belongs to the space U . Then for the solvability of the Equation (10) in U it is necessary and sufficient that the identities
< H t , τ , e τ 1 > 0 , t t 0 , T
hold true.
Proof. 
We will determine the solution of the Equation (10) in the form of an element (5) of the space U:
z t , τ = z 0 t + i = 1 4 z i t e τ i + 2 m N H z m t e m , τ .
Substituting (12) into the Equation (10), we have
λ 1 t z 0 t + i = 1 4 λ i t λ 1 t z i t e τ i + 2 m N H m , λ t λ 1 t z m t e m , τ
e τ 4 t 0 t K t , s z 4 s d s = H 0 t + i = 1 4 H i t e τ i + 2 m N H H m t e m , τ .
Equating here separately the free terms and coefficients at the same exponents, we obtained the following equations:
λ 1 t z 0 t = H 0 t ,
λ i t λ 1 t z i t = H i t , i = 1 , 3 ¯ ;
λ 4 t λ 1 t z 4 t t 0 t K t , s z 4 s d s = H 4 t ;
m , λ t λ 1 t z m t = H m t , m Γ 1 , 2 m N H .
Since the function λ 1 t 0 t [ t 0 , T ] , the Equation (13 0 ) has a unique solution z 0 ( t ) = λ 1 1 t H 0 ( t ) . Since λ 4 t λ 1 t 0 t [ t 0 , T ] , then the Equation 13 4 can be written as
z 4 t = t 0 t λ 4 t λ 1 t 1 K t , s z 4 s d s λ 4 t λ 1 t 1 H 4 t .
Due to the smoothness of the kernel λ 4 t λ 1 t 1 K t , s and heterogeneity λ 4 t λ 1 t 1 H 4 t , this Volterra integral equation has a unique solution z 4 t C t 0 , T , C . The Equations 13 2 and 13 3 also have unique solutions
z i t = λ i t λ 1 t 1 H i t C t 0 , T , C , i = 2 , 3 ,
since λ i t λ 1 t , i = 2 , 3 . The Equation 13 1 is solvable in the space C t 0 , T , C if and only if identities H 1 t , e τ 1 0 t t 0 , T hold. It is easy to see that this identity coincides with identity (11).
Further, since m , λ t λ 1 t , 2 m N H ( m Γ 1 ) , then the Equation 13 m has a unique solution
z m t = m , λ t I A t 1 H m t C t 0 , T , C , 2 m N H .
Thus, condition (11) is necessary and sufficient for the solvability of the Equation (10) in the space U . The Theorem 1 is proved. □
Remark 1. 
If identity (11) holds, then under conditions (1) and (2), the Equation (10) has the following solution in the space U :
z t , τ = z 0 t + i = 1 4 z i t e τ i + 2 m N H z m t e m , τ z 0 t + α 1 t e τ 1 + + h 21 ( t ) e τ 2 + h 31 ( t ) e τ 3 + z 4 t e τ 4 + 2 m N H P m t e ( m , τ ) ,
where α 1 t C t 0 , T , C are arbitrary function, z 0 t = λ 1 1 ( t ) H 0 ( t ) , z 4 t is the solution of the integral Equation (14), and introduced notations
h 21 ( t ) H 2 t λ 2 t λ 1 t , h 31 ( t ) H 3 t λ 3 t λ 1 t , P m t m , λ t λ 1 t 1 H m t .

4. The Remainder Term Theorem

Along with problem (10), we consider the equation
L w t , τ = z t + g t 2 e τ 2 σ 1 + e τ 3 σ 2 z + R 1 z + Q t , τ ,
where z = z t , τ is the solution (15) of Equation (10), Q t , τ U is the known function of the space U (this form will have problems 9 k + 1 after calculating the solution of the problem 9 k in U). The right side of this equation:
G t , τ z t + g t 2 e τ 2 σ 1 + e τ 3 σ 2 z + R 1 z + Q t , τ =
= t z 0 t + i = 1 4 z i t e τ i + 2 m N H z m t e m , τ +
+ g t 2 e τ 2 σ 1 + e τ 3 σ 2 z 0 t + i = 1 4 z i t e τ i + 2 m N H z m t e m , τ + R 1 z + Q t , τ ,
may not belong to the space U , if z = z t , τ U . Indeed, taking into account the form (15) of function z = z t , τ U , we consider in G t , τ , for example, the terms
Z t , τ g t 2 e τ 2 σ 1 + e τ 3 σ 2 z 0 t + i = 1 4 z i t e τ i + 2 m N H z m t e m , τ =
= g t 2 z 0 t e τ 2 σ 1 + e τ 3 σ 2 + i = 1 4 g t 2 z i t e τ i + τ 2 σ 1 + e τ i + τ 3 σ 2 +
+ g t 2 e τ 2 σ 1 + e τ 3 σ 2 2 m N H P m t e m , τ .
Function Z t , τ U , since it contains resonant exponentials e τ 2 + τ 3 = e m , τ | m = 0 , 1 , 1 , 0 , e τ 2 + m , τ m 2 + 1 = m 3 , e τ 3 + m , τ m 3 + 1 = m 2 , and, therefore, the right-hand side G t , τ = Z t , τ of the Equation (16) also does not belong to the U . Then, according to the well-known theory (see [6], p. 234), we need to embed ∧: G t , τ G ^ t , τ the right-hand side G t , τ of the Equation (16) into the space U . This operation is defined as follows.
Let the function G t , τ = m = 0 N w m t e m , τ contain resonant exponentials, i.e., G t , τ , it has the form
G t , τ = w 0 t + i = 1 4 w i t e τ i + j = 0 4 m j = 2 : m j Γ j N w m j t e m j , τ + m = 2 , m m j , j = 0 , 4 ¯ N w m t e m , τ .
Then
G ^ t , τ = w 0 t + i = 1 4 w i t e τ i + j = 0 4 m j = 2 : m j Γ j N w m j t e τ j + m = 2 , m m j , j = 0 , 4 ¯ N w m t e m , τ .
Therefore, the embedding operation acts only on the resonant exponentials and replaces them with a unit or exponents e τ j of the first dimension according to the rule:
e m , τ | m Γ 0 = e 0 = 1 , e m , τ | m Γ j = e τ j , j = 1 , 4 ¯ .
Therefore, the right-hand sides of iterative problems 9 k (if they solve sequentially) may not belong to the space U . Then, according to [6] (p. 234), the right-hand sides of these problems must be embedded in U according to the above rule. As a result, we obtained the following problems:
L z 0 t , τ j = 1 4 λ j t z 0 τ j A ( t ) z 0 R 0 z 0 = h t , z 0 t 0 , 0 = z 0 ;
L z 1 t , τ = z 0 t + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z 0 + R 1 z 0 , z 1 t 0 , 0 = 0 ;
L z 2 t , τ = z 1 t + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z 1 + R 1 z 1 + R 2 z 0 , z 0 t 0 , 0 = 0 ;
L z k t , τ = z k 1 t + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z k 1 + R k z 0 + + R 1 z k 1 , z k t 0 , 0 = 0 , k 1
(images of linear operators t and R ν do not need to be embedded in the space U , since these operators act from U to U). Such a replacement will not affect the construction of an asymptotic solution to the original problem (1) (or its equivalent problem (2)), since the narrowing τ = ψ ( t ) ε of the series of problems 9 ¯ k will coincide with the series of problems 9 k (see [6], pp. 234–235).
It is easy to show that applying Theorem 1 to iterative problems 9 ¯ k , we can find their solutions uniquely in the space U. As a result, we can construct series (7) with coefficients z k ( t , τ ) U . As in [23] (pp. 303–308), we proved the following statement.
Theorem 2.
Suppose that conditions (1)–(2) are satisfied for the Equation (2). Then, when ε ( 0 , ε 0 ] ( ε 0 > 0 is sufficiently small) the Equation (2) has a unique solution z ( t , ε ) C 1 ( [ t 0 , T ] , C ) ; at the same time there is the estimate
| | z ( t , ε ) z ε N ( t ) | | C [ t 0 , T ] c N ε N + 1 , N = 0 , 1 , 2 , ,
where z ε N ( t ) is the narrowing (for τ = ψ ( t ) ε ) N -th partial sum of the series (7) (with coefficients z k t , τ U satisfying the iterative problems ( 9 ¯ k ) ), and the constant c N > 0 does not depend ε on ε ( 0 , ε 0 ] .

5. Construction of the Solution of the First Iteration Problem in the Space U

Using Theorem 1, we will tried to find a solution to the first iterative problem 9 ¯ 0 . Since the right-hand side h t of the equation 9 ¯ 0 satisfies condition (11), this equation has (according to (15)) a solution in the space U in the form
z 0 t , τ = z 0 0 t + α 1 0 t e τ 1 ,
where α 1 0 t C t 0 , T , C are arbitrary function, z 0 0 t = h ( t ) λ 1 ( t ) . Subordinating (17) to the initial condition z 0 t 0 , 0 = z 0 , we have
z 0 0 t 0 + α 1 0 t 0 = z 0 α 1 0 t 0 = z 0 + λ 1 1 t 0 h t 0 .
To fully calculate the function α 1 0 t , we pass to the next iterative problem 9 ¯ 1 . Substituting the solution (17) of the equation 9 ¯ 0 , into it, we arrived at the following equation:
L z 1 t , τ = d d t z 0 0 t d d t α 1 0 t e τ 1 + K t , t z 0 0 t λ 4 t e τ 4
K t , t 0 z 0 0 t 0 λ 4 t 0 + g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z 0 0 t + α 1 0 t e τ 1 +
+ K t , t α 1 0 t λ 1 t e τ 1 K t , t 0 α j 0 t 0 λ 1 t 0 ,
(here we used the expression 6 1 for R 1 z t , τ and took into account that when z t , τ = z 0 t , τ in the sum 6 1 only terms with e τ 1 and remain e τ 4 ). Let us calculate
M = g ( t ) 2 e τ 2 σ 1 + e τ 3 σ 2 z 0 0 t + α 1 0 t e τ 1 =
= 1 2 g t σ 1 α 1 0 t e τ 2 + τ 1 + σ 2 α 1 0 t e τ 3 + τ 1 + σ 1 z 0 0 t e τ 2 + σ 2 z 0 0 t e τ 3 .
Let us analyze the exponents of the second dimension included here for their resonance:
e τ 2 + τ 1 | τ = ψ t / ε = e 1 ε t 0 t i β ( θ ) + a ( θ ) d θ , i β + a = 0 , a , i β , + i β , μ , Ø ; e τ 3 + τ 1 | τ = ψ t / ε = e 1 ε t 0 t + i β ( θ ) + a ( θ ) d θ , + i β + a = 0 , a , i β , + i β , μ , Ø .
Thus, exponents e τ 2 + τ 1 ang e τ 3 + τ 1 are not resonant. Then, for solvability the Equation (18) it is necessary and sufficient that the condition
d d t α 1 0 t + K t , t α 1 0 t λ 1 t = 0
is satisfied. Attaching the initial condition α 1 0 t 0 = z 0 + λ 1 1 t 0 h t 0 , to this equation, we found uniquely the function
α 1 0 t = α 1 0 t 0 exp t 0 t K s , s λ 1 s d s ,
and therefore, we uniquely calculate the solution (17) of the problem 9 ¯ 0 in the space U . In this case, the leading term of the asymptotics of the solution to the problem (2) has the form
z ε 0 t = z 0 0 t + α 1 0 t 0 exp t 0 t K s , s λ 1 s d s e 1 ε t 0 t λ 1 θ d θ ,
where α 1 0 t 0 = z 0 + A 1 t 0 h t 0 , z 0 0 t = λ 1 1 t 0 h t .
Example. 
Consider a model problem
ε d z d t = z ε cos t 2 + t ε z t 0 t e 2 ( t s ) ε · t · s · z ( s , ε ) d s + h ( t ) , z ( t 0 , ε ) = z 0 , t [ t 0 , T ] ( t 0 0 ) ,
were a ( t ) 1 , μ ( t ) 2 , β ( t ) t 2 + t , K ( t , s ) t · s . The main term of the asymptotic solution of this problem has the form
z ε 0 t = h ( t ) + [ z 0 h ( t 0 ) ] exp [ t 0 3 t 3 3 ] e x p [ t t 0 ε ] .
For ε + 0 the function z ε 0 t tends to the solution of the degenerate equation z ¯ + h ( t ) = 0 uniformly on any interval [ t 0 + δ , T ] ( 0 < δ T t 0 ) and at the point t = t 0 takes on the value z ε 0 t 0 = z 0 . It is seen from (20) that the leading term of the asymptotics of the solution to problem (19) does not depend on cos t 2 + t ε and spectral value μ ( t ) 2 , but depends on the kernel K ( t , s ) t · s . Further calculations show that already the asymptotic solution z ε 1 t = z ε 0 t + ε z 1 t , ψ ( t ) ε of the first order will depend on both μ ( t ) 2 , and the frequency β ( t ) = 2 t + 1 of the rapidly oscillating cosine.

6. Conclusions

The function z ε 0 t shows that when passing from a differential equation of type (1) ( K ( t , s ) 0 ) to an integro-differential one ( K ( t , s ) 0 ), the main term of the asymptotic is influenced by the kernel K ( t , s ) of the integral operator. However, the main term of the asymptotics is not affected by the spectral values of the integral operator μ ( t ) and rapidly oscillating coefficients. Their effects are detected when constructing the next approximation z ε 1 t .

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grant No. AP05133858 of the Ministry of Education and Science of the Republic of Kazakhstan.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Feschenko, S.F.; Shkil, N.I.; Nikolenko, L.D. Asymptotic Methods in the Theory of Linear Differential Equations; Naukova Dumka: Kiev, Ukraine, 1966. [Google Scholar]
  2. Shkil, N.I. Asymptotic Methods in Differential Equations; Naukova Dumka: Kiev, Ukraine, 1971. [Google Scholar]
  3. Daletsky, Y.L.; Krein, S.G. On differential equations in Hilbert space. Ukr. Math. J. 1950, 2, 71–91. [Google Scholar]
  4. Daletsky, Y.L. The asymptotic method for some differential equations with oscillating coefficients. DAN USSR 1962, 143, 1026–1029. [Google Scholar]
  5. Daletsky, Y.L.; Krein, M.G. Stability of Solutions of Differential Equations in a Banach Space; Nauka: Moscow, Russia, 1970. [Google Scholar]
  6. Lomov, S.A. Introduction to General Theory of Singular Perturbations; Vol. 112 of Translations of Math; Monographs, American Math Society: Providence, RI, USA, 1992. [Google Scholar]
  7. Ryzhih, A.D. Asymptotic solution of a linear differential equation with a rapidly oscillating coefficient. Trudy MEI 1978, 357, 92–94. [Google Scholar]
  8. Ryzhih, A.D. Application of the regularization method for an equation with rapidly oscillating coefficients. Mater. All-Union. Conf. Asymptot. Methods 1979, 1, 64–66. [Google Scholar]
  9. Kalimbetov, B.T.; Temirbekov, M.A.; Khabibullaev, Z.O. Asymptotic solution of singular perturbed problems with an instable spectrum of the limiting operator. Abstr. Appl. Anal. 2012, 120192. [Google Scholar] [CrossRef]
  10. Imanbaev, N.S.; Kalimbetov, B.T.; Temirbekov, M.A. Asymptotics of solutions of singularly perturbed integro- differential equation with rapidly decreasing kernel. Bull.-Ksu-Math. 2013, 72, 55–63. [Google Scholar]
  11. Bobodzhanov, A.A.; Safonov, V.F. Singularly perturbed nonlinear integro-differential systems with rapidly varying kernels. Math. Notes. 2002, 72, 605–614. [Google Scholar] [CrossRef]
  12. Nefedov, N.N.; Nikitin, A.G. The Cauchy problem for a singularly perturbed integro-differential Fredholm equation. Comput. Math. Math. Phys. 2007, 47, 629–637. [Google Scholar] [CrossRef]
  13. Bobodzhanov, A.A.; Safonov, V.F. Singularly perturbed integro-differential equations with diagonal degeneration of the kernel in reverse time. J. Differ. Equ. 2004, 40, 120–127. [Google Scholar] [CrossRef]
  14. Bobodzhanov, A.A.; Safonov, V.F. Asymptotic analysis of integro-differential systems with an unstable spectral value of the integral operator’s kernel. Comput. Math. Math. Phys. 2007, 47, 65–79. [Google Scholar] [CrossRef]
  15. Bobodzhanov, A.A.; Safonov, V.F. The method of normal forms for singularly perturbed systems of Fredholm integro-differential equations with rapidly varying kernels. Sibir. Math. J. 2013, 204, 979–1002. [Google Scholar] [CrossRef]
  16. Vasil’eva, A.B.; Butuzov, V.F.; Nefedov, N.N. Singularly perturbed problems with boundary and internal layers. Proc. Steklov Inst. Math. 2010, 268, 258–273. [Google Scholar] [CrossRef]
  17. Imanaliev, M. Asymptotic Methods in the Theory of Singularly Perturbed Integro-Differential Equations; Ilim: Frunze, Kyrgyzstan, 1972. [Google Scholar]
  18. Bobodzhanov, A.; Safonov, V.; Kachalov, V. Asymptotic and pseudoholomorphic solutions of singularly perturbed differential and integral equations in the lomov’s regularization method. Axioms 2019, 8, 27. [Google Scholar] [CrossRef] [Green Version]
  19. Imanaliev, M. Methods for Solving Nonlinear Inverse Problems and Their Applications; Ilim: Frunze, Kyrgyzstan, 1977. [Google Scholar]
  20. Bobodzhanov, A.A.; Safonov, V.F. Asymptotic solutions of Fredholm integro-differential equations with rapidly changing kernels and irreversible limit operator. Russ. Math. 2015, 59, 1–15. [Google Scholar] [CrossRef]
  21. Bobodzhanov, A.A.; Safonov, V.F. A problem with inverse time for a singularly perturbed integro-differential equation with diagonal degeneration of the kernel of high order. Izv. Math. 2016, 80, 285–298. [Google Scholar] [CrossRef]
  22. Bobodzhanov, A.A.; Kalimbetov, B.T.; Safonov, V.F. Integro-differential problem about parametric amplification and its asymptotical integration. Int. J. Appl. Math. 2020, 33, 331–353. [Google Scholar] [CrossRef]
  23. Safonov, V.F.; Bobodzhanov, A.A. Course of Higher Mathematics. Singularly Perturbed Equations and the Regularization Method: Textbook; Publishing House of MPEI: Moscow, Russia, 2012. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kalimbetov, B.; Safonov, V. Regularization Method for Singularly Perturbed Integro-Differential Equations with Rapidly Oscillating Coefficients and Rapidly Changing Kernels. Axioms 2020, 9, 131. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms9040131

AMA Style

Kalimbetov B, Safonov V. Regularization Method for Singularly Perturbed Integro-Differential Equations with Rapidly Oscillating Coefficients and Rapidly Changing Kernels. Axioms. 2020; 9(4):131. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms9040131

Chicago/Turabian Style

Kalimbetov, Burkhan, and Valeriy Safonov. 2020. "Regularization Method for Singularly Perturbed Integro-Differential Equations with Rapidly Oscillating Coefficients and Rapidly Changing Kernels" Axioms 9, no. 4: 131. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms9040131

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