Next Article in Journal
Fixed Point Theorems through Modified ω-Distance and Application to Nontrivial Equations
Next Article in Special Issue
Complete Controllability Conditions for Linear Singularly Perturbed Time-Invariant Systems with Multiple Delays via Chang-Type Transformation
Previous Article in Journal
Generating Root-Finder Iterative Methods of Second Order: Convergence and Stability
Previous Article in Special Issue
Euclidean Space Controllability Conditions for Singularly Perturbed Linear Systems with Multiple State and Control Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Projector Approach to Constructing Asymptotic Solution of Initial Value Problems for Singularly Perturbed Systems in Critical Case

by
Galina Kurina
1,2,3,†
1
Voronezh State University, Voronezh 394018, Russia
2
Voronezh Institute of Law and Economics, Voronezh 394042, Russia
3
Federal Research Center “Computer Science and Control” of Russian Academy of Sciences, Moscow 119333, Russia
Voronezh State University, Universitetskaya pl.,1, Voronezh 394018, Russia.
Submission received: 13 February 2019 / Revised: 20 April 2019 / Accepted: 21 April 2019 / Published: 8 May 2019

Abstract

:
Under some conditions, an asymptotic solution containing boundary functions was constructed in a paper by Vasil’eva and Butuzov (Differ. Uravn. 1970, 6(4), 650–664 (in Russian); English transl.: Differential Equations 1971, 6, 499–510) for an initial value problem for weakly non-linear differential equations with a small parameter standing before the derivative, in the case of a singular matrix A ( t ) standing in front of the unknown function. In the present paper, the orthogonal projectors onto k e r A ( t ) and k e r A ( t ) (the prime denotes the transposition) are used for asymptotics construction. This approach essentially simplifies understanding of the algorithm of asymptotics construction.

1. Introduction

The bibliography of publications devoted to singularly perturbed problems is very extensive. Most of them deal with problems in which a degenerate equation, following from the original one where a small parameter is equal to zero, is resolvable with respect to a fast component of an unknown variable. If it is not so, then this more complicated case is known as critical [1], singular [2], nonstandard [3], or as a case where the unperturbed (degenerate) system is situated on the spectrum [4]. Numerous applications of singularly perturbed systems in the critical cases have been listed in [5].
Vasil’eva and Butuzov were the first to study initial value problems for singularly perturbed differential and difference systems in the critical case. Asymptotic solutions of boundary value problems for such systems have been obtained in [1,2,6]. Numerical methods for singularly perturbed systems in the critical case have been researched in [7] for initial value problems, and in [8] for boundary value problems.
An asymptotic solution containing boundary functions for the initial value problem of the weakly non-linear differential equation in a real m-dimensional space X:
ε d x d t = A ( t ) x + ε f ( x , t , ε ) , t [ 0 , T ] ,
x ( 0 , ε ) = x 0 ,
where x = x ( t , ε ) X and the matrix A ( t ) is singular, has been constructed in [4]. A discrete analogue of problem (1)-(2) was also considered. The results from this paper are also presented in [1,9]. In these publications, the purpose of studying equations of the last form is also explained. Here and further ε 0 means a small parameter, and the m × m matrix A ( t ) and the m-dimensional vector-function f ( x , t , ε ) are sufficiently smooth with respect to their arguments.
In contrast [4], the projector approach will be used in this paper for constructing an asymptotic solution of problem (1)-(2). It allows us to represent the algorithm of the boundary functions method for constructing an asymptotic solution of initial-value singularly perturbed problems in the critical case more clearly than in [4].
Note that the projector approach has been used in [10] for constructing the zero-order asymptotic solution for a singularly perturbed linear-quadratic control problem in the critical case.
We will assume the same assumptions as in [4] that the matrix A ( t ) has for each t [ 0 , T ] m eigenvalues λ 1 ( t ) , λ 2 ( t ) , . . . , λ m ( t ) , and that they satisfy the conditions:
Assumption 1.
λ j ( t ) = 0 for j = 1 , 2 , . . . , k , k < m .
Assumption 2.
All k eigenvectors v 1 ( t ) , v 2 ( t ) , . . . , v k ( t ) of the matrix A ( t ) , corresponding to λ j ( t ) = 0 , j = 1 , 2 , . . . , k , are linearly independent.
Following [4], we will here use eigenvectors having the same smoothness as the matrix A ( t ) . The existence of such eigenvectors has been proved in [11].
Furthermore, some assumptions will be yet added.
The transposition will be denoted by the prime. By I, as usual, we mean the identity operator. For the expansion of a function w ( ε ) into the series with respect to integer non-negative powers of ε w ( t , ε ) = j 0 ε j w j ( t ) , we introduce the notation [ w ( ε ) ] j = w j .
The paper is organized as follows. In Section 2, we present the standard decomposition of the original system (1) into systems with respect to functions from the asymptotic solution, depending on t, and with respect to so-called boundary functions, depending on the argument t / ε . In the next section, we introduce orthogonal projectors of the space X onto k e r A ( t ) and k e r A ( t ) . Based on these projectors, the algorithm of constructing the zero-order asymptotic approximation of a solution of problem (1)-(2) is given in Section 4, and the algorithm of constructing the n-th order asymptotic approximation, n 1 , is developed in Section 5. Table 1 and Table 2 in these two sections show the sequence of actions for finding asymptotics terms. In the sixth section, we present an example illustrating the projector approach for constructing the first-order asymptotic approximation. The last section presents our conclusions.

2. Problem Decomposition

In view of [4], we will seek the asymptotic solution of problem (1)-(2) in the form:
x ( t , ε ) = x ¯ ( t , ε ) + Π x ( τ , ε ) ,
where x ¯ ( t , ε ) = j 0 ε j x ¯ j ( t ) , Π x ( τ , ε ) = j 0 ε j Π j x ( τ ) , τ = t / ε . Functions Π j x ( τ ) will be found as in [4] with the help of the additional condition
Π j x ( τ ) 0 as τ + .
Following tradition (see, for instance, [1], p. 8), a series j 0 ε j x ¯ j ( t ) with terms depending on the original argument t is called regular series, in contrast with boundary series j 0 ε j Π j x ( τ ) consisting of so-called boundary functions depending on the argument τ 0 , which are essential only for arguments in some vicinities of points where additional conditions are prescribed (in a vicinity of zero in the considered case).
As usual in the theory of singular perturbations, the following representation will be used
f ( x ¯ ( t , ε ) + Π x ( τ , ε ) , t , ε ) f ¯ + Π f ,
where f ¯ = f ( x ¯ ( t , ε ) , t , ε ) = j 0 ε j f ¯ j ( t ) and Π f = f ( x ¯ ( ε τ , ε ) + Π x ( τ , ε ) , ε τ , ε ) - f ( x ¯ ( ε τ , ε ) , ε τ , ε ) = j 0 ε j Π j f ( τ ) .
Substituting expansion (3) into (1) and equating terms of the same order of ε separately depending on t and τ , we obtain the following equations for the terms of series (3):
d x ¯ j - 1 ( t ) d t = A ( t ) x ¯ j ( t ) + f ¯ j - 1 ( t ) ,
d Π j x ( τ ) d τ = A ( 0 ) Π j x ( τ ) + Π j - 1 f ( τ ) + [ ( A ( ε τ ) - A ( 0 ) ) Π x ( τ , ε ) ] j ,
where j = 0 , 1 , . . . ,
[ ( A ( ε τ ) - A ( 0 ) ) Π x ( τ , ε ) ] j = k = 0 j - 1 1 ( j - k ) ! d j - k A d t j - k ( 0 ) Π k x ( τ ) .
In order to write equations (5) and (6) in the same forms for the cases j = 0 and j > 0 , we suppose that terms of expansions with negative indices are equal to zero.
Substituting expansion (3) into (2) and equating terms of the same order of ε , we obtain the equalities:
x ¯ 0 ( 0 ) + Π 0 x ( 0 ) = x 0 ,
x ¯ j ( 0 ) + Π j x ( 0 ) = 0 , j > 0 .

3. Space Decomposition

Further, we will use the decompositions of the space X in the orthogonal sums (see, for instance, [12], p. 38)
X = k e r A ( t ) i m A ( t ) = k e r A ( t ) i m A ( t ) .
Orthogonal projectors P ( t ) and Q ( t ) of the space X onto the subspaces k e r A ( t ) and k e r A ( t ) , respectively, corresponding to the decompositions of the space X into two last orthogonal sums, will be applied. We can write the explicit form of these projectors. Namely, let V ( t ) = ( v 1 ( t ) , . . . , v k ( t ) ) and S ( t ) = ( s 1 ( t ) , . . . , s k ( t ) ) , where s 1 ( t ) , . . . , s k ( t ) are the eigenvectors of the matrix A ( t ) corresponding to eigenvalues λ j ( t ) = 0 , j = 1 , . . . , k . Following [9], we believe that the eigenvectors s i ( t ) have been chosen in such a way that V ( t ) S ( t ) is the k × k identity matrix. We explain that this is possible. The invertibility of the matrix V ( t ) S ( t ) is proved in [1]. If V ( t ) S ( t ) = B ( t ) I , then we take the columns of the matrix S ( t ) B ( t ) - 1 as s 1 ( t ) , . . . , s k ( t ) .
It easily follows from Assumption 2 that the k × k matrices V ( t ) V ( t ) and S ( t ) S ( t ) are invertible. It is not difficult to see that P ( t ) = V ( t ) ( V ( t ) V ( t ) ) - 1 V ( t ) and Q ( t ) = S ( t ) ( S ( t ) S ( t ) ) - 1 S ( t ) are orthogonal projectors of the space X onto the subspaces k e r A ( t ) and k e r A ( t ) , respectively, corresponding to the decompositions of the space X into the orthogonal sums.
The operator
A ( t ) = ( I - Q ( t ) ) A ( t ) ( I - P ( t ) ) : i m A ( t ) i m A ( t )
has the inverse operator. It will be denoted as A ( t ) + = ( I - P ( t ) ) A ( t ) + ( I - Q ( t ) ) .
The following condition is assumed.
Assumption 3.
For each t [ 0 , T ] the operator ( I - P ( t ) ) A ( t ) ( I - P ( t ) ) : i m A ( t ) i m A ( t ) is stable—that is, all eigenvalues of this operator have negative real parts.
It is not difficult to prove that the operator Q ( t ) P ( t ) : k e r A ( t ) k e r A ( t ) is invertible. Let us take a vector x from k e r A ( t ) . Then, x = V ( t ) c ( t ) , where c ( t ) = ( c 1 ( t ) , c 2 ( t ) , . . . , c k ( t ) ) and c i ( t ) , i = 1 , 2 , . . . , k , are some scalar functions. Consider the equation Q ( t ) P ( t ) x = 0 . It follows from this that V ( t ) S ( t ) ( S ( t ) S ( t ) ) - 1 S ( t ) V ( t ) c ( t ) = 0 . Since V ( t ) S ( t ) is a k × k identity matrix, then c ( t ) = 0 , which gives the provable invertibility.

4. Zero-Order Asymptotic Solution

From (5), we have the equation for x ¯ 0 ( t ) :
A ( t ) x ¯ 0 ( t ) = 0 .
Hence,
( I - P ( t ) ) x ¯ 0 ( t ) = 0 .
Using (9), we find from (7) the initial value
( I - P ( 0 ) ) Π 0 x ( 0 ) = ( I - P ( 0 ) ) x 0 .
From (6), we have the equation for Π 0 x ( τ )
d Π 0 x ( τ ) d τ = A ( 0 ) Π 0 x ( τ ) .
This equation is equivalent to two ones:
d ( I - P ( 0 ) ) Π 0 x ( τ ) d τ = ( I - P ( 0 ) ) A ( 0 ) ( I - P ( 0 ) ) Π 0 x ( τ ) ,
d ( P ( 0 ) Π 0 x ( τ ) ) d τ = P ( 0 ) A ( 0 ) ( I - P ( 0 ) ) Π 0 x ( τ ) .
In view of Assumption 3, we obtain a unique solution of initial problem (10)-(11) satisfying the inequality
( I - P ( 0 ) ) Π 0 x ( τ ) c exp ( - α τ ) , τ 0 ,
with some positive constants c and α independent of τ (see, for instance, [13], p. 106). In this estimate, any norm may be used, since all norms in a finite dimensional space are equivalent. Functions satisfying the last inequality are called exponential-type boundary functions.
From (12), we get the equality P ( 0 ) Π 0 x ( τ ) = P ( 0 ) Π 0 x ( 0 ) + 0 τ P ( 0 ) A ( 0 ) ( I - P ( 0 ) ) Π 0 x ( s ) d s . Since P ( 0 ) Π 0 x ( τ ) 0 as τ + , then P ( 0 ) Π 0 x ( 0 ) = - 0 + P ( 0 ) A ( 0 ) ( I - P ( 0 ) ) Π 0 x ( s ) d s . Using the exponential estimate for ( I - P ( 0 ) ) Π 0 x ( s ) , we uniquely define the exponential-type boundary function P ( 0 ) Π 0 x ( τ ) , namely,
P ( 0 ) Π 0 x ( τ ) = - τ + P ( 0 ) A ( 0 ) ( I - P ( 0 ) ) Π 0 x ( s ) d s .
Hence, the exponential-type boundary function Π 0 x ( τ ) has been found. Then, we can get the initial value from (7):
P ( 0 ) x ¯ 0 ( 0 ) = P ( 0 ) ( x 0 - Π 0 x ( 0 ) ) .
In view of (5), the equation for x ¯ 1 ( t ) has the form
A ( t ) x ¯ 1 ( t ) = - f ¯ 0 ( t ) + d x ¯ 0 ( t ) d t .
Taking into account (9), we can write the solvability condition for the last equation in the form
Q ( t ) d ( P ( t ) x ¯ 0 ( t ) ) d t = Q ( t ) f ( P ( t ) x ¯ 0 ( t ) , t , 0 ) .
Since
d ( P ( t ) x ( t ) ) d t = d ( P ( t ) 2 x ( t ) ) d t = d ( P ( t ) ) d t P ( t ) x ( t ) + P ( t ) d ( P ( t ) x ( t ) ) d t ,
we obtain the equation
d ( P ( t ) x ¯ 0 ( t ) ) d t = ( Q ( t ) P ( t ) ) - 1 Q ( t ) ( - d P ( t ) d t P ( t ) x ¯ 0 ( t ) + f ( P ( t ) x ¯ 0 ( t ) , t , 0 ) ) + ( I - P ( t ) ) d P ( t ) d t P ( t ) x ¯ 0 ( t ) .
If operator A ( t ) is constant, then projectors P ( t ) = P and Q ( t ) = Q are constant too, and the last equation has the form
d ( P x ¯ 0 ( t ) ) d t = ( Q P ) - 1 Q f ( P x ¯ 0 ( t ) , t , 0 ) .
We will yet assume the condition.
Assumption 4.
Problem (14)–(16) has a unique solution on the segment [ 0 , T ] .
A similar assumption regarding the solvability of some initial-value problem for a non-linear equation of the smaller dimension than the original one was presented in [1] (Assumption IV, p. 13).
Thus, the function x ¯ 0 ( t ) is defined. Hence, the zero-order asymptotics for a solution of problem (1)-(2) is found.
The following Table 1 shows the sequence of finding zero-order asymptotics terms.

5. Higher-Order Asymptotic Solutions

Suppose that the terms x ¯ j ( t ) and Π j x ( τ ) of expansion (3), j = 0 , 1 , . . . , n - 1 , n 1 , have been found.
From equation (5) with j = n , we obtain the relation
A ( t ) x ¯ n ( t ) = d x ¯ n - 1 ( t ) d t - f ¯ n - 1 ( t ) ,
where the right-hand side is known. Applying the operator ( I - Q ( t ) ) to this equation, we have
( I - Q ( t ) ) A ( t ) ( I - P ( t ) ) x ¯ n ( t ) = ( I - Q ( t ) ) ( d x ¯ n - 1 ( t ) d t - f ¯ n - 1 ( t ) ) .
From here, we find:
( I - P ( t ) ) x ¯ n ( t ) = A ( t ) + ( I - Q ( t ) ) ( d x ¯ n - 1 ( t ) d t - f ¯ n - 1 ( t ) ) .
Then, we can find from (8) with j = n the initial value
( I - P ( 0 ) ) Π n x ( 0 ) = - ( I - P ( 0 ) ) x ¯ n ( 0 ) .
The equation (6) with j = n has the form
d Π n x ( τ ) d τ = A ( 0 ) Π n x ( τ ) + Π n - 1 f ( τ ) + [ ( A ( ε τ ) - A ( 0 ) ) Π x ( τ , ε ) ] n .
This equation is equivalent to two ones.
d ( I - P ( 0 ) ) Π n x ( τ ) d τ = ( I - P ( 0 ) ) A ( 0 ) ( I - P ( 0 ) ) Π n x ( τ ) + ( I - P ( 0 ) ) Π n - 1 f ( τ ) + + ( I - P ( 0 ) ) [ ( A ( ε τ ) - A ( 0 ) ) Π x ( τ , ε ) ] n ,
d ( P ( 0 ) Π n x ( τ ) ) d τ = P ( 0 ) ( A ( 0 ) ( I - P ( 0 ) ) Π n x ( τ ) + Π n - 1 f ( τ ) + + [ ( A ( ε τ ) - A ( 0 ) ) Π x ( τ , ε ) ] n ) .
The sum of two last summands in the right-hand side in (20) is a known exponential-type boundary function. Therefore, in view of Assumption 3, we can find from (19) and (20) the exponential-type boundary function ( I - P ( 0 ) ) Π n x ( τ ) . Note that the proof of exponential estimates for boundary functions is given in detail in monograph [14].
As the function in the braces on the right-hand side in (21) is a known exponential type boundary function, we can get from (21) the exponential-type boundary function P ( 0 ) Π n x ( τ ) , namely
P ( 0 ) Π n x ( τ ) = - τ + P ( 0 ) ( A ( 0 ) ( I - P ( 0 ) ) Π n x ( s ) + Π n - 1 f ( s ) + [ ( A ( ε s ) - A ( 0 ) ) Π x ( s , ε ) ] n ) d s .
Hence, the exponential-type boundary function Π n x ( τ ) is defined. Then, we can find from (8) with j = n the initial value
P ( 0 ) x ¯ n ( 0 ) = - P ( 0 ) Π n x ( 0 ) .
Writing out equation (5) with j = n + 1 , we get
d x ¯ n ( t ) d t = A ( t ) x ¯ n + 1 ( t ) + f ¯ n ( t ) .
The solvability condition for this equation has the form
Q ( t ) d ( P ( t ) x ¯ n ( t ) ) d t = Q ( t ) ( f ¯ n ( t ) - d ( ( I - P ( t ) ) x ¯ n ( t ) ) d t ) .
In view of (15), we obtain from here the equation
d ( P ( t ) x ¯ n ( t ) ) d t = ( Q ( t ) P ( t ) ) - 1 Q ( t ) ( f ¯ n ( t ) + d P ( t ) d t P ( t ) x ¯ n ( t ) - - d ( ( I - P ( t ) ) x ¯ n ( t ) ) d t ) + ( I - P ( t ) ) d P ( t ) d t P ( t ) x ¯ n ( t ) .
If operator A ( t ) is constant, then this equation has the form:
d ( P x ¯ n ( t ) ) d t = ( Q P ) - 1 Q ( f ¯ n ( t ) - d ( ( I - P ) x ¯ n ( t ) ) d t ) .
It should be noted that equation (24) is linear with respect to P ( t ) x ¯ n ( t ) . As ( I - P ( t ) ) x ¯ n ( t ) has been found (see (18)), we can define the function P ( t ) x ¯ n ( t ) from (23) and (24).
Hence, we have found the terms of the n-th order in expansion (3).
The following Table 2 shows the sequence of finding the n-th order terms in expansion (3).
The previous arguments have, as a consequence, the following assertion.
Theorem 1.
Under Assumptions 1–4, the asymptotic solution of problem (1) and (2) in form (3) can be constructed with the help of orthogonal projectors onto k e r A ( t ) and k e r A ( t ) . The order of finding the asymptotics terms is the following: ( I - P ( t ) ) x ¯ j ( t ) , ( I - P ( 0 ) ) Π j x ( τ ) , P ( 0 ) Π j x ( τ ) , P ( t ) x ¯ j ( t ) , j 0 .

6. Illustrative Example

Consider the following initial value problem of form (1)-(2) on the segment [ 0 , T ] :
ε d y d t = - y + ε z 2 , ε d z d t = y ,
y ( 0 , ε ) = 1 , z ( 0 , ε ) = 1 .
Here, t [ 0 , 0 . 3 ] , x ( t , ε ) = ( y ( t , ε ) , z ( t , ε ) ) ; y = y ( t , ε ) , z = z ( t , ε ) I R ; x 0 = ( 1 , 1 ) ,
A ( t ) = A = - 1 0 1 0 , f ( x , t , ε ) = z 2 0 .
Hence,
λ 1 = 0 , λ 2 = - 1 , k e r A = { ( 0 , a ) } , V = ( 0 , 1 ) , k e r A = { ( a , a ) } ,
S = ( 1 , 1 ) , i m A = { ( - a , a ) } , i m A = { ( a , 0 ) } ,
P = 0 0 0 1 , I - P = 1 0 0 0 , Q = 1 / 2 1 / 2 1 / 2 1 / 2 ,
( I - P ) A ( I - P ) = - 1 0 0 0 : i m A i m A ,
( Q P ) - 1 = 0 0 0 2 : k e r A k e r A , A + = - 1 / 2 1 / 2 0 0 : i m A i m A .
We will construct the first-order approximation for the asymptotic solution of problem (26)-(27) using projectors P and Q.
Relation (9), in this case, has the form:
1 0 0 0 y ¯ 0 ( t ) z ¯ 0 ( t ) = y ¯ 0 ( t ) 0 = 0 0 .
Therefore, y ¯ 0 ( t ) = 0 .
From (10), we get Π 0 y ( 0 ) = 1 .
Equation (11) has the form:
d d τ Π 0 y ( τ ) 0 = - Π 0 y ( τ ) 0 .
Taking into account the initial value Π 0 y ( 0 ) = 1 found from (10), we obtain from the last equation Π 0 y ( τ ) = e - τ .
From (13), we find Π 0 z ( τ ) = - τ + Π 0 y ( s ) d s = - e - τ .
From (14) and (17) we derive, respectively,
0 0 0 1 y ¯ 0 ( 0 ) z ¯ 0 ( 0 ) = 0 z ¯ 0 ( 0 ) = 0 0 0 1 1 1 - 1 - 1 = 0 2 .
d d t 0 z ¯ 0 ( t ) = 0 0 0 2 1 / 2 1 / 2 1 / 2 1 / 2 ( z ¯ 0 ( t ) ) 2 0 = 0 ( z ¯ 0 ( t ) ) 2 ,
In view of the last two relations, we have z ¯ 0 ( t ) = 1 / ( 0.5 - t ) .
It is easy to verify that conditions 1–4 are satisfied for problem (26)-(27).
Thus, we have found the zero-order asymptotic solution of form (3) x ˜ 0 ( t , ε ) for the solution of problems (26)-(27). Namely, we have
y ˜ 0 ( t , ε ) = e - τ , z ˜ 0 ( t , ε ) = 1 / ( 0.5 - t ) - e - τ , τ = t / ε .
Now, we will seek for the first-order asymptotics.
Equation (18) for n = 1 has the form:
y ¯ 1 ( t ) 0 = - 1 / 2 1 / 2 0 0 1 / 2 - 1 / 2 - 1 / 2 1 / 2 d d t 0 z ¯ 0 ( t ) - ( z ¯ 0 ( t ) ) 2 0 = 1 / ( 0.5 - t ) 2 0 .
Therefore, y ¯ 1 ( t ) = 1 / ( 0.5 - t ) 2 .
From (19) with n = 1 , we get Π 1 y ( 0 ) = - 4 .
Equation (20) for n = 1 has the form:
d d τ Π 1 y ( τ ) 0 = - Π 1 y ( τ ) + 4 Π 0 z ( τ ) + ( Π 0 z ( τ ) ) 2 0 .
From the last two relations, we obtain Π 1 y ( τ ) = - ( 3 + 4 τ + e - τ ) e - τ .
From (22) with n = 1 , we find Π 1 z ( τ ) = - τ + Π 1 y ( s ) d s = ( 7 + 4 τ + e - τ / 2 ) e - τ .
From (23) and (25) with n = 1 , we derive, respectively,
d d t 0 z ¯ 1 ( t ) = 0 0 0 2 1 / 2 1 / 2 1 / 2 1 / 2 2 z ¯ 0 ( t ) z ¯ 1 ( t ) 0 - d d t y ¯ 1 ( t ) 0 = 0 2 z ¯ 0 ( t ) z ¯ 1 ( t ) - d y ¯ 1 ( t ) / d t ,
0 z ¯ 1 ( 0 ) = 0 - Π 1 z ( 0 ) .
In view of the last two relations, we have z ¯ 1 ( t ) = ( l n ( 0.5 - t ) 2 - 15 / 8 + l n 4 ) / ( 0.5 - t ) 2 .
Thus, we have found for problems (26)-(27) the first-order asymptotic solution of form (3) x ˜ 1 ( t , ε ) . Namely, we have
y ˜ 1 ( t , ε ) = y ˜ 0 ( t , ε ) + ε ( y ¯ 1 ( t ) + Π 1 y ( τ ) ) ,
z ˜ 1 ( t , ε ) = z ˜ 0 ( t , ε ) + ε ( z ¯ 1 ( t ) + Π 1 z ( τ ) ) .
Of course, these results can be obtained using the algorithm from [4], but we would like to demonstrate here the use of projectors for finding asymptotics terms. The results obtained by Maple 13 are given in Figure 1 and Figure 2. They have been presented for the completeness of the paper. The solid line represents the exact solution; the dash-dotted line—the solution of the degenerate problem, the line consisting of squares represents the zero-order approximation; and the dash line represents the the first-order approximation. These graphs show that an asymptotic solution is closer to the exact one if we use higher-order asymptotics. If we use the smaller value of ε , then it will result in an asymptotic solution more similar to the exact one. The graphs of the solution of the degenerate problem and the zero-order approximation illustrate the known property of boundary functions that are essential only for arguments in some vicinities of points where additional conditions are prescribed.

7. Conclusions

This paper dealt with a new approach to the algorithm of the method of boundary functions from [4] for asymptotic solving initial value problem of form (1)-(2) in the critical case. Namely, the algorithm was formulated with the help of orthogonal projectors of the space X onto k e r A ( t ) and k e r A ( t ) . Such an approach clearly shows the structure of the algorithm for finding asymptotics terms, given in Table 1 and Table 2 of the paper.

Funding

This research was funded by the RUSSIAN SCIENCE FOUNDATION, grant number 17-11-01220.

Acknowledgments

The author thank the anonymous referees for carefully reading the paper and for the numerous critical and constructive comments that have helped to improve the paper text. She thanks also N. T. Hoai for useful discussions and M. A. Kalashnikova for the help in printing the paper text.

Conflicts of Interest

The author declare no conflict 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. Vasil’eva, A.B.; Butuzov, V.F. Singularly Perturbed Equations in Critical Cases; Moscow University: Moscow, Russia, 1978; 107p, English transl.: University ofWisconsin-Madison: Wisconsin, WI, USA, 1980; 165p. [Google Scholar]
  2. O’Malley, R.E., Jr. A singular singularly-perturbed linear boundary value problem. SIAM J. Math. Anal. 1979, 10, 695–708. [Google Scholar] [CrossRef]
  3. Khalil, H.K. Feedback control of nonstandard singularly perturbed systems. IEEE Trans. Autom. Control 1989, 34, 1052–1060. [Google Scholar] [CrossRef]
  4. Butuzov, V.F.; Vasil’eva, A.B. Differential and difference systems of equations with a small parameter in the case when the unperturbed (degenerate) system is situated on the spectrum. Differ. Uravn. 1971, 6, 650–664. (In Russian); English transl.: Differ. Equ. 1971, 6, 499–510. [Google Scholar]
  5. Gu, Z.M.; Nefedov, N.N.; O’Malley, R.E., Jr. On singular singularly perturbed initial value problems. SIAM J. Appl. Math. 1989, 49, 1–25. [Google Scholar] [CrossRef]
  6. Schmeiser, C.; Wess, R. Asymptotic analysis of singular singularly perturbed boundary value problems. SIAM J. Math. Anal. 1986, 17, 560–579. [Google Scholar] [CrossRef]
  7. O’Malley, R.E., Jr.; Flaherty, J.E. Analytical and numerical methods for nonlinear singular singularly-perturbed initial value problems. SIAM J. Appl. Math. 1980, 38, 225–248. [Google Scholar] [CrossRef]
  8. Ascher, U. On some difference schemes for singular singularly-perturbed boundary value problems. Numer. Math. 1985, 46, 1–30. [Google Scholar] [CrossRef]
  9. Vasil’eva, A.B.; Butuzov, V.F.; Kalachev, L.V. The Boundary Function Method for Singular Perturbation Problems; SIAM: Philadelphia, PA, USA, 1995. [Google Scholar]
  10. Kurina, G.A.; Hoai, N.T. Projector approach for constructing the zero order asymptotic solution for the singularly perturbed linear-quadratic control problem in a critical case. AIP Conf. Proc. 2018, 1997, 020073-1–020073-7. [Google Scholar] [CrossRef]
  11. Sibuya, Y. Some global properties of matrices of functions of one variable. Math. Ann. 1965, 161, 67–77. [Google Scholar] [CrossRef]
  12. Kato, T. Perturbation Theory for Linear Operators; Mir: Moscow, Russia, 1972; 740p. [Google Scholar]
  13. Daletskii, Y.L.; Krein, M.G. Stability of Solutions of Differential Equations in Banach Space; Nauka: Moscow, Russia, 1970; 536p. [Google Scholar]
  14. Vasil’eva, A.B.; Butuzov, V.F. Asymptotic Expansions of Solutions of Singularly Perturbed Equations; Nauka: Moscow, Russia, 1973; 272p. [Google Scholar]
Figure 1. Trajectory y ( t , ε ) with ε = 0.01 and its approximations.
Figure 1. Trajectory y ( t , ε ) with ε = 0.01 and its approximations.
Axioms 08 00056 g001
Figure 2. Trajectory z ( t , ε ) with ε = 0.01 and its approximations.
Figure 2. Trajectory z ( t , ε ) with ε = 0.01 and its approximations.
Axioms 08 00056 g002
Table 1. The algorithm for finding the zero-order asymptotics terms.
Table 1. The algorithm for finding the zero-order asymptotics terms.
Asymptotics TermsFormulas
( I - P ( t ) ) x ¯ 0 ( t ) = 0 (9)
( I - P ( 0 ) ) Π 0 x ( τ ) (10), (11)
P ( 0 ) Π 0 x ( τ ) (13)
P ( t ) x ¯ 0 ( t ) (14), (16)
Table 2. The algorithm for finding the n-th order asymptotics terms, n 1 .
Table 2. The algorithm for finding the n-th order asymptotics terms, n 1 .
Asymptotics TermsFormulas
( I - P ( t ) ) x ¯ n ( t ) (18)
( I - P ( 0 ) ) Π n x ( τ ) (19), (20)
P ( 0 ) Π n x ( τ ) (22)
P ( t ) x ¯ n ( t ) (23),(24)

Share and Cite

MDPI and ACS Style

Kurina, G. Projector Approach to Constructing Asymptotic Solution of Initial Value Problems for Singularly Perturbed Systems in Critical Case. Axioms 2019, 8, 56. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms8020056

AMA Style

Kurina G. Projector Approach to Constructing Asymptotic Solution of Initial Value Problems for Singularly Perturbed Systems in Critical Case. Axioms. 2019; 8(2):56. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms8020056

Chicago/Turabian Style

Kurina, Galina. 2019. "Projector Approach to Constructing Asymptotic Solution of Initial Value Problems for Singularly Perturbed Systems in Critical Case" Axioms 8, no. 2: 56. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms8020056

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