Advances in Pure Mathematics
Vol.08 No.05(2018), Article ID:84750,22 pages
10.4236/apm.2018.85026

Existence and Stability Property of Almost Periodic Solutions in Discrete Almost Periodic Systems

Yoshihiro Hamaya

Department of Information Science, Okayama University of Science 1-1 Ridai-cho, Okayama, Japan

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: February 26, 2018; Accepted: May 21, 2018; Published: May 24, 2018

ABSTRACT

In this paper, we consider an almost periodic system which includes a system of the type x ( n + 1 ) = j 1 m a i j ( n ) x j k ( 1 i m ) , where k is a positive integer, a i j are almost periodic in n and satisfy a i j ( n ) 0 for i j , i = 1 m a i j ( n ) = 0 for 1 j m . In the special case where a i j ( n ) are constant functions, above system is a mathematical model of gas dynamics and was treated by T. Carleman and R. D. Jenks for differential systems. In the main theorem, we show that if the m × m matrix ( a i j ( n ) ) is irreducible, then there exists a positive almost periodic solution which is unique and has some stability. Moreover, we can see that this result gives R. D. Jenks’ result for differential model in the case where aij(n) are constant functions. In Section 3, we consider the linear system with variable cofficients x ( n + 1 ) = A ( n ) x ( n ) , x R m . Even in nonlinear problems, this linear system plays an important role, as their variational equations, and it is requested to determine the uniform asymptotically stability of the zero solution from the information about A ( n ) . In order to obtain the existence of almost periodic solutions of both linear and nonlinear almost periodic discrete systems: above linear system and x i ( n + 1 ) = j = 1 m a i j ( n ) g j ( x j ( n ) ) for 1 i m , respectively, we shall consider between certain stability properties, which are referred to as uniformly asymptotically stable, and the diagonal dominance matrix condition.

Keywords:

Almost Periodic Solutions, Linear and Nonlinear Almost Periodic Discrete Systems, Uniformly Asymptotically Stable, Diagonal Dominance Matrix Condition

1. Introduction

System of almost periodic difference equations has been studied to describe phenomena of oscillations in the natural and social sciences. The investigation of almost periodic systems has been developed quite widely during the twentieth century, since relationships with the stability theory have been found. A main interest of the subject is the existence theorem for almost periodic solutions. Obviously an almost periodic solution is a bounded solution, but the existence of bounded solutions does not necessarily imply the existence of almost periodic solutions. Therefore, in order to prove the existence of almost periodic solutions, we need some additional conditions to the existence of bounded solutions. A main subject of the investigation has been to find such additional conditions, and up to now, many conditions have been considered (for example, in the linear system, J. Favard’s separation condition [1] ).

In the Section 4, we consider the nonlinear almost periodic system of

x i ( n + 1 ) = j = 1 m a i j ( n ) x j ( n ) k , 1 i m , (1)

where k is a positive integer, a i j ( n ) are almost periodic in n and satisfy

i) i = 1 m a i j ( n ) = 0 , 1 j m ,

ii) a i j ( n ) 0 for i j .

In the special case where a i j ( n ) are constant functions, system (1) is a mathematical model of gas dynamics and was treated by T. Carleman [2] and R. D. Jenks [3] . In the main theorem, we show that if the m × m matrix ( a i j ( n ) ) is irreducible, then there exists a positive almost periodic solution which is unique and some stability. Moreover, we can see that this result gives R. D. Jenks’ result in the case where a i j ( n ) are constant functions. In the Section 5, we consider the linear almost periodic system with variable coefficients

x ( n + 1 ) = A ( n ) x ( n ) , n n 0 0, (2)

where x R m . Even in nonlinear problems, system (2) plays an important role, as their variational equations and moreover, it is requested to determine the uniformly asymptotic stability of the zero solution from the condition about A ( n ) . When A ( n ) is a constant matrix, it is well known that the stability is equivalent to the following condition (cf. [4] );

“Absolute values of all eigenvalues of A ( n ) A are less than one.”

However, it is not true in the case of variable coefficients, and hence we need additional conditions to (2). In the main theorem, we show that one of the such conditions is the diagonal dominance matrix condition on A ( n ) [5] , that is, A ( n ) satisfies

| a i i ( n ) | k = 1 , k i m | a k i ( n ) | , 1 i m .

This result improves a stability criterion based on results of F. Nakajima [6] for differential equations.

2. Preliminaries

We denote by Rm the real Euclidean m-space. Let R = ( , ) and R + = [ 0 , ) . Z is the set of integers, Z+ is the set of nonnegative integers. For x R m , let | x | be the Euclidean norm of x and x i be the i-th component. Let

D = { x R m | x i 0 for 1 i m } ,

Ω = { x D | i = 1 m x i = 1 }

and

Π = { x R m | i = 1 m x i = 0 } .

We introduce an almost periodic function f ( n , x ) : Z × U R m , where U is an open set in Rm.

Definition 1. f ( n , x ) is said to be almost periodic in n uniformly for x U , if for any ϵ > 0 and any compact set K in U there exists a positive integer L * ( ϵ , K ) such that any interval of length L * ( ϵ , K ) contains an integer τ for which

| f ( n + τ , x ) f ( n , x ) | ϵ

for all n Z and all x K . Such a number τ in above inequality is called an ò-translation number of f ( n , x ) .

In order to formulate a property of almost periodic functions, which is equivalent to the above definition, we discuss the concept of the normality of almost periodic functions. Namely, let f ( n , x ) be almost periodic in n uniformly for x U . Then, for any sequence { h k } Z , there exist a subsequence { h k } of { h k } and a function g ( n , x ) such that

f ( n + h k , x ) g ( n , x ) (3)

uniformly on Z × K as k , where K is a compact set in U. There are many properties of the discrete almost periodic functions [7] , which are corresponding properties of the continuous almost periodic functions f ( t , x ) C ( R × U , R m ) [cf. [8] [9] ]. We denote by T ( f ) the function space consisting of all translates of f, that is, f τ T ( f ) , where

f τ ( n , x ) = f ( n + τ , x ) , τ Z . (4)

Let H ( f ) denote the uniform closure of T ( f ) in the sense of (4). H ( f ) is called the hull of f. In particular, we denote by Ω ( f ) the set of all limit functions g H ( f ) such that for some sequence { n k } , n k as k and f ( n + n k , x ) g ( n , x ) uniformly on Z × S for any compact subset S in Rm. Specially, for a function f ( n ) on Z with values in Rm, H ( f ) denotes the set of all function g ( n ) such that for some sequence { n k } ,

f ( n + n k ) g ( n ) in Z as k ,

where the symbol “?” stands for the uniformly convergence on any compact set in Z (in short, “in Z”). Clearly, f H ( f ) .

By (3), if f : Z × U R m is almost periodic in n uniformly for x U , so is a function in Ω ( f ) .

We define the irreducible matrix to need after.

Definition 2. An m × m matrix A ( n ) = ( a i j ( n ) ) is said to be irreducible if for any two nonempty disjoint subsets I and J of the set of m integers { 1,2, , m } with I J = { 1,2, , m } , there exists an i in I and a j in J such that a i j ( n ) 0 . In the case where A ( n ) is scalar, A ( n ) is said to be irreducible if A ( n ) 0 . Otherwise, A ( n ) is said to be reducible, and we can assume that A ( n ) takes the form of

A ( n ) = ( 0 B ( n ) ) ,

where * is l × l ( 1 l m ) square matrix, *' is m l × l matrix, B ( n ) is m l × m l zero or a square irreducible matrix.

3. Linear Systems

We consider the system of linear difference equation

x ( n + 1 ) = A ( n ) x ( n ) , (5)

where x R m and the m × m matrix A ( n ) = ( a i j ( n ) ) is bounded on Z and almost periodic function in n. We state discretization of Jenks and Nakajima' results for differential equations [3] [10] .

Now we define stability properties with respect to the subset K in Rm. Here, we denote by x ( n , n 0 , x 0 ) the solution of system (5) with initial condition ( n 0 , x 0 ) .

Definition 3. The bounded solution u ( n ) of system (5) defined on Z is said to be;

i) uniformly stable (in short, U.S.) in K on Z+ if for any ϵ > 0 there exists a δ ( ϵ ) > 0 such that | u ( n ) x ( n , n 0 , x 0 ) | < ϵ for all n n 0 whenever x 0 K and | u ( n 0 ) x 0 | < δ ( ϵ ) at some n 0 in Z+.

ii) uniformly asymptotically stable (in short, U.A.S.) in K on Z+ if it is U.S. in K on Z+ and if there exists a δ 0 > 0 and, if for any ϵ > 0 there exists a T ( ϵ ) > 0 such that | u ( n ) x ( n , n 0 , x 0 ) | < ϵ for all n n 0 + T ( ϵ ) whenever x 0 K and | u ( n 0 ) x 0 | < δ 0 at some n 0 in Z+.

iii) uniformly asymptotically stable (in short, U.A.S.) in the whole K on Z+ if it is U.S. in K on Z+ and if for any ϵ > 0 and r > 0 there exists a T ( ϵ , r ) > 0 such that | u ( n ) x ( n , n 0 , x 0 ) | < ϵ for all n n 0 + T ( ϵ , r ) whenever x 0 K and | u ( n 0 ) x 0 | < r , at some n 0 in Z+.

When Z+ in the definitions (i), (ii) and (iii) is replaced by Z, we say that x ( n ) is U.S. in K on Z, U.A.S. in K on Z and U.A.S. in the whole K on Z, respectively. Clearly Definition 3 agrees with the definitions of the usual stability properties in the case where K = R m .

Throughout this paper, we suppose the following conditions;

i) i = 1 m a i j ( n ) = 0 for 1 j m ,

ii) a i j ( n ) 0 for i j

and

iii) each element in H ( A ) is irreducible.

First of all, we prove the following lemmas.

Lemma 1. Consider the m-equations x i ( n + 1 ) = f i ( n , x ( n ) ) , 1 i m , where f i ( n , x ) is continuous on second variable x in Rm, and assume that the initial value problem has a unique solution.

a) If i = 1 m f i ( n , x ) = 0 , then the set Π is invariant.

b) If f i ( n , x ) 0 for x i = 0 and x j 0 , then the set D is positively invariant, and in addition,

if i = 1 m f i ( n , x ) = 0 , then the set Ω is positively invariant.

In the case of differential equations, the proof of the similar lemma is obvious (for instance, see [11] ). We modify it to prove this lemma, but we omit it.

Lemma 2. If conditions (i) and (ii) are satisfied, then the trivial solution of system (5) is U.S. in P on Z and also it is U.S. on Z.

By modifying theorem in [5] , we can easily prove Lemma 2 at same technique.

Lemma 3. If each element in H ( A ) is irreducible, then the each element in H ( A ) , we say B H ( A ) and B ( n ) = ( b i j ( n ) ) , has the property that for any two nonempty disjoint subsets I and J of the set of m integers { 1,2, , m } with I J = { 1,2, , m } , there exists an i I and j J such that

l i m ¯ n | b i j ( n ) | 0.

Proof. Suppose not, Then there exists a B ( n ) = ( b i j ( n ) ) in H ( A ) and two nonempty disjoint subsets I and J of { 1,2, , m } with I J = { 1,2, , m } such that

lim ¯ n b i j ( n ) = 0 for all i I and j J .

Since B ( n ) is bounded on Z, there exists a subsequence n k , n k as k , such that

B ( n + n k ) C ( n ) in Z as k ,

where C ( n ) = ( c i j ( n ) ) H ( A ) . Clearly,

c i j ( n ) = l i m k b i j ( n + n k ) = 0 for n Z and i I , j J .

This show the reducibility of C ( n ) , which is a contradiction. This proves Lemma 3.

For system (5), we consider the system in H ( A ) of

x ( n + 1 ) = B ( n ) x ( n ) , (6)

where B H ( A ) .

Lemma 4. Assume that conditions (ii) and (iii) are satisfied for system (5), and let x ( n ) be a nontrivial solution of system (6) such that x ( n ) D on Z. Then there exists a constant c > 0 such that

x i ( n ) / | x ( n ) | c for n Z and 1 i m .

Proof. Let x ( n ) = ( x 1 ( n ) , x 2 ( n ) , , x m ( n ) ) be a solution of system (6) such that x ( n ) D on Z. First of all, we show that if x i ( n 0 ) = 0 at some n 0 Z , then

x i ( n ) = 0 for all n n 0 .

Since x i ( n ) satisfies the equation

x i ( n + 1 ) = b i i ( n ) x i ( n ) + j i m b i j ( n ) x j ( n ) ,

where ( b i j ( n ) ) = B ( n ) . Moreover, since j i m b i j ( n ) x j ( n ) 0 , we have

x i ( n + 1 ) b i i ( n ) x i ( n ) , (7)

which implies

x i ( n ) x i ( n 0 ) s = n 0 n 1 b i i ( s ) for n n 0 .

Thus, we obtain

x i ( n ) = 0 for n n 0 .

Because x i ( n 0 ) = 0 and x i ( n ) 0 on Z. Now suppose that Lemma 4 is not true. Then for some B in H ( A ) , the corresponding system (6) has a nontrivial solution x ( n ) , x ( n ) D on Z, such that for some sequence n k ,

x 1 ( n k ) / | x ( n k ) | 0 as k . (8)

Set ϕ k ( n ) = x ( n + n k ) / | x ( n k ) | . Then, ϕ k ( n ) satisfies

x ( n + 1 ) = B ( n + n k ) x (n)

and

ϕ k ( n ) D on Z , | ϕ k ( 0 ) | = 1.

Since the sequence { ϕ k ( 0 ) } is bounded, { ϕ k ( n ) } is uniformly bounded on any finite interval in Z, and hence there is a convergent subsequence of ϕ k , which is again denoted by ϕ k , such that

ϕ k ( n ) y ( n ) in Z for some function y ( n ) as k .

We can also assume that

B ( n + n k ) C ( n ) in Z as k ,

where C H ( A ) and C ( n ) = ( c i j ( n ) ) . Therefore, y ( n ) is the solution of the system

y ( n + 1 ) = C ( n ) y ( n ) , (9)

y ( n ) D on Z and | y ( 0 ) | = 1 . Moreover (8) implies that y 1 ( 0 ) = 0 . Thus, as was proved above, we have

y 1 ( n ) = 0 for n 0,

For this y ( n ) , we define two subsets I and J of { 1,2, , m } by I = { 1 i m | y i ( n ) 0 for n N i , where N i depends on y i ( n ) } and J = { 1 i m | y i ( n ) > 0 on Z } . Then I J = { 1,2, , m } , { 1 } I and J ϕ since y ( n ) 0 . By Lemma 3,

l i m ¯ n | c i 0 j 0 ( n ) | 0 for some i 0 I and some j 0 J . (10)

Now the i 0 -th equation of system (9) takes the form of

y i 0 ( n + 1 ) = k I c i 0 k ( n ) y k ( n ) + k J c i 0 k ( n ) y k ( n ) ,

and hence

k J c i 0 k ( n ) y k ( n ) = 0 for n + 1 m i n i I N i , (11)

because of the definition of the set I. Since each term in the left hand side of (11) is nonnegative, all of them are equal to zero. Therefore

c i 0 j 0 ( n ) y j 0 ( n ) = 0 for n + 1 min i I N i ,

which implies, by (10),

y j 0 ( n 0 ) = 0 at some n 0 .

This contradicts the definition of the set of J. The proof is completed.

The following proposition is an immediate result of Lemma 4.

Proposition 1. Under conditions (ii) and (iii), system (6) has no nontrivial solution x ( n ) such that

x ( n ) D on Z ,

where D = { x D | x i = 0 for some i ,1 i m } .

We next consider a non-homogeneous system corresponding to system (5)

x ( n + 1 ) = A ( n ) x ( n ) + f ( n ) (12)

and assume that A ( n ) satisfies conditions (i), (ii) and (iii).

Lemma 5. If f ( n ) is bounded on Z with values in Rm and i = 1 m f i ( n 1 ) is bounded on Z+, then all solutions of system (12) are bounded on Z+.

Proof. It is sufficient to show that (12) has at least one bounded solution on Z+, because the trivial solution of (5) is U.S. by Lemma 2. We consider the system with real parameter ò

x ( n + 1 ) = A ( n ) x ( n ) + ϵ f ( n ) (13)

and show that for a sufficiently small ò, system (13) has a bounded solution on Z+, which implies the existence of a bounded solution on Z+ for system (12) by replacing x in (13) with x / ϵ . For a 0 < δ < 1 / m and for the m-vector e each of whose components is 1, let D be a convex cone defined by

D = { x R m | e , x | e | | x | δ } ,

where , denotes the inner product and | x | 2 = x , x . Clearly, D D . Every solution x ( n ) of (13) satisfies

i = 1 m x i ( n + 1 ) = ϵ i = 1 m f i (n)

because of condition (i). By replacing n with n-1 of the above both sides,

i = 1 m x i ( n ) = ϵ i = 1 m f i ( n 1 ) M ˜ + ϵ M for x ( 0 ) Ω ,

where

M ˜ > 0 is sufficient small number and M = sup n > 0 i = 1 m f i ( n 1 ) .

When x ( n ) D , we have

i = 1 m x i ( n ) = e , x ( n ) | x ( n ) | | e | δ ,

and hence,

( M ˜ ϵ M ) / | e | | x ( n ) | ( M ˜ + ϵ M ) / ( | e | δ ) , (14)

Therefore, in order to show the boundedness of x ( n ) with x ( 0 ) in Ω, it is sufficient to prove that x ( n ) D on Z+ if ò is sufficiently small. Now suppose that for each solution x ϵ ( n ) of (13) with x ϵ ( 0 ) in Ω, there exists an l ϵ > 0 such that

x ϵ ( l ϵ ) D .

We can assume that

x ϵ ( n ϵ ) D at some n ϵ ,0 n ϵ < l ϵ

and

x ϵ ( n ) D D ¯ for n ϵ n l ϵ ,

where K and K ¯ denote the boundary and the closure of the set K, respectively. If we set y ϵ ( n ) = x ϵ ( n + n ϵ ) , y ϵ ( n ) is a solution of the system

y ( n + 1 ) = A ( n + n ϵ ) y ( n ) + ϵ f ( n + n ϵ )

such that y ϵ ( 0 ) D , y ϵ ( τ ϵ ) D at τ ϵ = l ϵ n ϵ > 0 and y ϵ ( n ) D D ¯ for 0 n τ ϵ . Thus, by (14),

( M ˜ ϵ M ) / | e | | y ϵ ( n ) | ( M ˜ + ϵ M ) / ( | e | δ ) for 0 n τ ϵ ,

The same argument in the proof of Lemma 4 enables us to assume that

y ϵ ( n ) w ( n ) in Z for some function w ( n ) as ϵ 0

and

A ( n + n ϵ ) B ( n ) in Z for some B H ( A ) as ϵ 0

Therefore, w ( n ) satisfies w ( n + 1 ) = B ( n ) w ( n ) and clearly, for τ = lim inf ϵ 0 τ ϵ ,

w ( n ) D D ¯ for 0 n < τ (15)

and

M ˜ / | e | | w ( n ) | M ˜ / ( | e | δ ) for 0 n τ . (16)

Moreover we have w ( 0 ) D , which implies by Lemma 1 that

w ( n ) D on Z + .

From this and (15) it follows that

w ( n ) D D D ¯ = D for 0 n < τ . (17)

Now we show that τ = . In fact, if τ < , we have

y ϵ ( τ ϵ ) w ( τ ) as ϵ 0.

Thus w ( τ ) D because y ϵ ( τ ϵ ) D , and hence

w ( τ ) D D = { 0 } ,

which contradicts (16). Therefore (16) and (17) hold for τ = . Moreover this enables us to assume that

M ˜ / | e | | w ( n ) | M ˜ / ( | e | δ ) for all n Z

and

w ( n ) D for all n Z .

Because H ( B ) is compact in the sense of the convergence. This contradicts the conclusion in Proposition 1. This proves that x ( n ) D on Z+ if ò is sufficiently small. The proof is completed.

Lemma 6. Under the assumptions (i) and (ii), if for each B in H ( A ) , the trivial solution of the system

x ( n + 1 ) = B ( n ) x (n)

is U.S. on Z and U.A.S. on Z+, then the trivial solution of system (5) is U.A.S. on Z.

Proof. Let x ( n , n 0 , x 0 ) be the solution of (5). Since the trivial solution of (5) is U.S. on Z by Lemma 2, as is seen from (ii) in Definition 3, it is sufficient to show that for any ϵ > 0 there exists a T ( ϵ ) > 0 such that

| x ( n , n 0 , x 0 ) | < δ ( ϵ ) for some n , n 0 n n 0 + T ( ϵ ) ,

whenever n 0 Z and | x 0 | < δ 0 = δ ( 1 ) , where δ ( ) is the number in (i) of Definition 3.

Now suppose that there exists an ϵ > 0 and sequences { n k } in Z and { x k } in Rm such that | x k | δ 0 and

| x ( n , n k , x k ) | δ ( ϵ ) for all n k , n k n n k + k .

Since | x k | < δ 0 = δ ( 1 ) ,

δ ( ϵ ) | x ( n , n k , x k ) | 1 for n k , n k n n k + k .

Set ϕ k ( n ) = x ( n + n k , n k , x k ) . Then, ϕ k ( n ) satisfies

x ( n + 1 ) = A ( n + n k ) x (n)

and

δ ( ϵ ) | ϕ k ( n ) | 1 for 0 n k .

We can assume that

ϕ k ( n ) y ( n ) in Z+ for some function y ( n ) as k

and

A ( n + n k ) B ( n ) in Z for some B H ( A ) as k .

Therefore y ( n ) is a solution of the system

y ( n + 1 ) = B ( n ) y ( n ) (18)

and

δ ( ϵ ) | y ( n ) | 1 on Z + .

On the other hand, we have

y ( n ) 0 as n ,

because the trivial solution of (18) is U.A.S. on Z+. Therefore there arises a contradiction. Thus the proof is completed.

We show the following theorem, before we will mention a definition of the exponential dichotomy of a linear system;

System (5) is said to possess an exponential dichotomy if there exists a projection matrix P and positive constants K 1 , K 2 , ρ 1 , ρ 2 , σ 1 and σ 2 such that

Φ ( n ) P Φ 1 ( l ) K 1 ρ 1 n l ( 0 < ρ 1 = e σ 1 ) , for n l ,

Φ ( n ) ( I P ) Φ 1 ( l ) K 2 ρ 2 n l ( 0 < ρ 2 = e σ 2 ) , for n l ,

where, I is a identical matrix and F is a fundamental matrix solution of system (5) (cf. [4] [5] [8] ).

Theorem 1. Assume that system (5) satisfies conditions (i), (ii) and (iii), Then the trivial solution of system (5) is U.A.S. in P on Z.

Proof. On the set Π which is invariant for system (5), the system is written as the ( m 1 ) -system

x ˜ ( n + 1 ) = A ˜ ( n ) x ˜ ( n ) (19)

where x ˜ = ( x 1 , x 2 , , x m 1 ) R m 1 and A ˜ ( n ) is an ( m 1 ) × ( m 1 ) matrix whose ( i , j ) element is given by a i j ( n ) a i m ( n ) for 1 i , j m 1 . First of all, we can show that for each B ˜ in H ( A ˜ ) , the system

x ˜ ( n + 1 ) = B ˜ ( n ) x ˜ ( n ) (20)

has an exponential dichotomy on Z+ since (20) has at least one bounded solution, and as is well known (cf. [12] ), it is equivalent to show the system

x ˜ ( n + 1 ) = B ˜ ( n ) x ˜ ( n ) + f ˜ ( n ) (21)

possesses at least one bounded solution on Z+ for any bounded function f ˜ ( n ) on Z+. For each B ˜ ( n ) in H ( A ˜ ) there corresponds some B ( n ) = ( b i j ( n ) ) in H ( A ) such that the ( i , j ) element of B ˜ ( n ) is equal to b i j ( n ) b i m ( n ) for 1 i , j m 1 . For f ˜ ( n ) = ( f 1 ( n ) , f 2 ( n ) , , f m 1 ( n ) ) , let g ( n ) be defined by

g i ( n ) = f i ( n ) for 1 i m 1,

g m ( n ) = i = 1 m 1 f i ( n ) .

Obviously g ( n ) and i = 1 m g i ( n 1 ) ( = 0 ) are bounded on Z+. Applying Lemma 5 to the m-system

x ( n + 1 ) = B ( n ) x ( n ) + g ( n ) ,

we obtain the bounded solution x ( n ) on Z+ with x ( 0 ) = 0 , and i = 1 m g i ( n 1 ) = 0 which yields

i = 1 m x i ( n ) = 0.

Hence we can verify that x ˜ ( n ) = ( x 1 ( n ) , x 2 ( n ) , , x m 1 ( n ) ) is a bounded solution on Z+ of system (21). The exponential dichotomy of (20) implies further that the trivial solution is U.A.S. on Z+, because the trivial solution is U.S. on Z by Lemma 2. Therefore it follows from Lemma 6 that the trivial solution of (19) is U.A.S. on Z, i.e., the trivial solution of (5) is U.A.S. in P on Z. The proof is completed.

4. Nonlinear Systems

We consider the nonlinear almost periodic system of

x i ( n + 1 ) = j = 1 m a i j ( n ) g j ( x j ( n ) ) for 1 i m , (22)

where A ( n ) = ( a i j ( n ) ) is almost periodic function of n with conditions

(iv) i = 1 m a i j ( n ) = 0 , 1 j m

and

(v) a i j ( n ) 0 for i j .

In addition, assume that g j ( u ) are continuously differentiable for u 0 , g j ( 0 ) = 0 and g ˙ j ( u ) > 0 for real number u > 0 , where g ˙ j ( u ) is the derivative of g j ( u ) at u.

We first consider the linear system

x ( n + 1 ) = L ( n ) x ( n ) (23)

and its perturbed system

x ( n + 1 ) = L ( n ) x ( n ) + f ( n , x ( n ) ) , (24)

where L ( n ) is an m × m matrix function, almost periodic function in n, f ( n , x ) is continuous with respect to its second argument and f ( n , x ) = o ( | x | ) uniformly for n Z . Assume that the set Π is invariant for both system (23) and (24).

First of all, we can prove the following lemmas.

Lemma 7. If the trivial solution of system (23) is U.A.S. in Π on Z, then the trivial solution of system (24) has also the same stability property.

Proof. Let x ˜ = ( x 1 , x 2 , , x m 1 ) for x = ( x 1 , x 2 , , x m ) R m . Then there are positive constants c 1 and c 2 such that

c 1 | x ˜ | | x | c 2 | x ˜ | for x Π , (25)

because x m = ( x 1 + x 2 + + x m 1 ) . On the set Π, systems (23) and (24) are written as

x ˜ ( n + 1 ) = L ˜ ( n ) x ˜ ( n ) (26)

and

x ˜ ( n + 1 ) = L ˜ ( n ) x ˜ ( n ) + g ( n , x ˜ ( n ) ) (27)

respectively, where the ( i , j ) element of L ˜ ( n ) ,1 i , j m 1 , is given by l i j ( n ) l i m ( n ) for L ( n ) = ( l i j ( n ) ) and g ( n , x ˜ ) = o ( | x ˜ | ) uniformly for n Z . Inequality (25) shows that the trivial solution of (23) is U.A.S. in P if and only if the trivial solution of (26) is U.A.S., and we have also the same equivalence between (24) and (27). As is well known, if the trivial solution of (26) is U.A.S., then the trivial solution of (27) has also the same stability property. Thus our assertion is proved.

The following lemma is obtained by the slight modification of the difference equation to Seifert’s result [13] . Then, we will omit the proof (cf. [9] ).

We consider the almost periodic nonlinear system

x ( n + 1 ) = f ( n , x ( n ) ) , n Z , x R m , (28)

where f ( n , x ) is almost periodic in n uniformly for x R m and for a constant L * > 0 , | f ( n , x ) f ( n , y ) | L * | x y | for n Z and x , y Ω .

Lemma 8. Assume that the set Ω is positively invariant for system (28) and all solutions in Ω on Z are U.A.S. in Ω on Z. Then the set of such solutions is finite and consists of only almost periodic solutions ϕ 1 , ϕ 2 , , ϕ m which satisfy

| ϕ i ( n ) ϕ j ( n ) | β on Z for i j and some constant β > 0 .

Now we can show the following theorem. Since the last statements of the following theorem are alternative, under each assumption of these statements we can prove the existence of almost periodic solutions in Ω and the module containment.

Theorem 2. Under the assumptions (iv) and (v), system (22) has a nontrivial almost periodic solution in Ω whose module is contained in the module of A ( n ) . In addition to the above assumptions, if A ( n ) is irreducible, then the almost periodic solution of (22) is unique in Ω, which remains in Ω 0 on Z, and it is U.A.S. in the whole Ω on Z, where Ω 0 = { x Ω | x i > 0 for all i , 1 i m } . Moreover, if A ( n ) is reducible, then at least one of the above almost periodic solutions p ( n ) satisfies that p ( n ) Ω on Z, where Ω = { x Ω | x i = 0 for some i , 1 i m } .

Proof. First of all, we consider the case where A ( n ) is irreducible. Since system (22) satisfies the conditions of Lemma 1, the set Ω is positively invariant, namely, y ( n ) Ω on Z+ for a solution y ( n ) of (22) with y ( 0 ) Ω , and furthermore we can assume that

y ( n ) Ω on Z

because of the almost periodicity of A ( n ) . We can show that this y ( n ) is U.A.S. in Ω on Z. If we set x = y ( n ) + w in system (22), then w Π for x in Ω and

w i ( n + 1 ) = j = 1 m a i j ( n ) g ˙ j ( y j ( n ) ) w j + o ( | w | ) for 1 i m , (29)

And Π is invariant for the above system. Considering the first approximation of system (29)

w ( n + 1 ) = L ( n ) w ( n ) , (30)

where L ( n ) = ( l i j ( n ) ) is defined by l i j ( n ) = a i j ( n ) g ˙ j ( y j ( n ) ) , condition (iv) implies that Π is also invariant for (30). Then, by Lemma 6, if the trivial solution of (30) is shown to be U.A.S. in Π on Z, then the trivial solution of (29) has the same stability, and consequently y ( n ) is U.A.S. in Ω on Z. Therefore it is sufficient to show that the trivial solution of (30) is U.A.S. in Π on Z. Clearly L ( n ) is bounded and we have

i = 1 m l i j ( n ) = i = 1 m a i j ( n ) g ˙ j ( y j ( n ) ) = 0 for 1 j m

and

l i j ( n ) = a i j ( n ) g ˙ j ( y j ( n ) ) 0 for i j (31)

because of conditions (iv) and (v), respectively. Thus L ( n ) satisfies conditions (i) and (ii). Condition (iii) will be verified in the following way. Applying the same argument as in the proof of Lemma 4 to system (22), we can see that there exists a constant c > 0 such that

1 y i ( n ) c for n Z and 1 i m , (32)

and hence there is a constant c ˜ > 0 such that

g ˙ j ( y j ( n ) ) c ˜ for n Z and 1 j m .

Therefore, (31) implies

l i j ( n ) c ˜ a i j ( n ) for i j ,

which guarantees that each element of H ( L ) is irreducible, because A ( n ) is irreducible and almost periodic. Thus it follows from Theorem 1 that the trivial solution of (30) is U.A.S. in Π on Z, i.e., all solutions of system (22) in Ω on Z are U.A.S. in Ω on Z. Therefore Lemma 8 concludes that system (22) possesses an almost periodic solution in Ω which remains in Ω 0 by (32), and the set of solutions in Ω on Z is finite and consists of only almost periodic solutions ϕ 1 , ϕ 2 , , ϕ m which satisfy

| ϕ i ( n ) ϕ j ( n ) | β on Z for i j and some constant β > 0 .

Next we can show that there exists a T > 0 such that each solution x ( n , n 0 , x 0 ) of (22) with x 0 Ω satisfies that for some ϕ j and the constant δ 0 of (ii) in Definition 3,

| x ( n , n 0 , x 0 ) ϕ j ( n ) | < δ 0 at some n , n 0 n n 0 + T ,

which implies

| x ( n , n 0 , x 0 ) ϕ j ( n ) | 0 as n , (33)

because ϕ j is U.A.S. in Ω. Suppose that this is not true. Then there exists a small constant α > 0 less than β and sequences { n k } in Z and { x k } in Ω such that

| x ( n , n k , x k ) ϕ j ( n ) | α for all n [ n k , n k + k ] and all j , 1 j m .

Since f ( n , x ) is almost periodic in n uniformly for x Ω , we can choose a sequence { τ k } , n k + [ k / 2 ] < τ k < n k + k , such that

f ( n + τ k , x ) f ( n , x ) in Z × Ω as k .

If we set ψ j ( k ) ( n ) = ϕ j ( n + τ k ) for 1 j m and ψ m + 1 ( k ) ( n ) = x ( n + τ k , n k , x k ) , these m + 1 functions ψ j ( k ) ( n ) satisfy

x ( n + 1 ) = f ( n + τ k , x (n)

and

ψ j ( k ) ( n ) Ω on Z for 1 j m ,

ψ m + 1 ( k ) ( n ) Ω for n [ k / 2 ] n k τ k ,

because ψ m + 1 ( k ) ( n k τ k ) = x k Ω . Moreover,

| ψ j ( k ) ( 0 ) ψ i ( k ) ( 0 ) | α for i j , 1 i , j m + 1.

We can assume that ψ j ( k ) ( n ) ψ j ( n ) in Z for some function ψ j ,1 j m + 1 , as k . Therefore ψ j ( j = 1 , 2 , , m ) are solutions of system (22), because f ( n + τ k , x ) f ( n , x ) in Z × Ω as k , and

ψ j ( n ) Ω on Z for 1 i m + 1

and

| ψ j ( 0 ) ψ i ( 0 ) | > α for i j , 1 i , j m + 1,

which shows that system (22) has m + 1 distinct solutions in Ω on Z. This is contradiction. Therefore, ψ j ( n ) is U.A.S. in the whole Ω on Z, if the uniqueness of ϕ j is shown.

Now we will prove the uniqueness of ϕ j . Suppose ϕ 1 ϕ j for j 2 and set

S 1 = { x 0 Ω | | x ( n ,0, x 0 ) ϕ 1 ( n ) | 0 as n }

and

S 2 = { x 0 Ω | | x ( n , 0 , x 0 ) ϕ j ( n ) | 0 for some ϕ j , j 2 , as n } .

Then S 1 and S 2 are open sets in Ω, and moreover these sets are nonempty and disjoint, because | ϕ 1 ( n ) ϕ j ( n ) | α on Z for j 2 . On the other hand , (33) shows that Ω = S 1 S 2 , which contradicts the connectedness of Ω. Thus the uniqueness of an almost periodic solution is proved, and moreover, as is seen from [14] , this uniqueness guarantees the module containment of the almost periodic solution.

Now consider the case where A ( n ) is reducible. We can assume A ( n ) takes the form of

A ( n ) = ( 0 B ( n ) ) ,

where B ( n ) is zero or a square irreducible matrix of order l ,2 l m 1 . If B ( n ) is zero, system (22) obviously has the constant solution p ( n ) in Ω such that p i ( n ) = 0 for 1 i m 1 and p m ( n ) = 1 . In the latter case, if we set in system (22)

x k = 0 for 1 k l m and y i = x m l + i for 1 i l ,

then system (22) is reduced to the lower dimensional system

y i ( n + 1 ) = j = 1 l b i j ( n ) g j ( y j ( n ) ) for 1 i l , (34)

where B ( n ) = ( b i j ( n ) ) . Since B ( n ) is irreducible, the above system (34) has an almost periodic solution y ( n ) such that

y i ( n ) > 0 for 1 i l and i = 1 l y i ( n ) = 1 (35)

and furthermore the module of y ( n ) is contained in the module of B ( n ) , i.e., of the module of A ( n ) . Thus, system (22) has an almost periodic solution p ( n ) in Ω on Z such that p i ( n ) = 0 for 1 i m l and p i ( n ) = y i m + l ( n ) for m l + 1 i m . The proof is completed.

Remark 1. As will be seen from the module containment, the above almost periodic solution is a critical point in the case where A ( n ) is a constant. Hence Theorem 2 is a discretization of Nakajimas’ result (Theorem 2 in [10] ).

5. A Stability Criteria of Linear Systems

We consider a stability criterion for solutions of a linear system with coefficient matrix of diagonal dominance type.

We again consider a linear system (5).

Let A ( n ) = ( a i j ( n ) ) , i , j = 1 , 2 , , m be an m × m matrix of functions for n Z . We assume the following conditions;

a j j ( n ) 0 for 1 j m and n Z , (36)

d e t A ( n ) 0, (37)

where, d e t A ( n ) denotes the determinant of matrix A ( n ) and

i = 1 , i j m | a i j ( n ) | | a j j ( n ) | for 1 j m and n Z . (38)

At first, we need the following lemmas for main theorem.

Lemma 9. If a square matrix A is irreducible and satisfies (38) and if for at least one j,

i = 1 , i j m | a i j | < | a j j | ,

then A is nonsingular.

For the proof, see [15] .

Lemma 10. If a nonsingular m × m matrix A = ( a i j ) satisfies (38), then all principal minors of A are nonsingular, namely,

det ( a j 1 j 1 a j 1 j 2 a j 1 j l a j 2 j 1 a j 2 j 2 a j 2 j l a j l j 1 a j l j 2 a j l j l ) 0 for 1 j 1 < j 2 < < j l m .

Proof. Let A 1 be an l × l ( l < m ) principal minor of A. Then, for a permutation matrix Q,

Q A Q T = ( A 1 A 2 ) ,

where A 2 has ( m l ) rows and l columns and Q T denotes the transposed matrix of Q. Moreover, from the definition of irreducibility, we can choose an l × l permutation matrix Q l such that

Q l A 1 Q l T = ( B 1 C 2 B 2 C p O B p ) (381)

where B i is an r i × r i irreducible matrix, i = 1 p r i = l , and C i has ( r 1 + r 2 + + r i 1 ) row and r i columns for 2 i p . In particular, in the case where A 1 is irreducible, B 1 must be A 1 itself, and the matrices B 2 , B 3 , , B p , C 2 , C 3 , , C p are not present. Setting P = ( Q l I ) Q for ( m l ) × ( m l ) unit matrix I, where Q l I is the direct sum of Q l and I, we have

B = P A P T = ( Q l A 1 Q l T A 2 Q l T ) = ( B 1 C 2 C p B 2 C 3 C p B p 1 C p O B p D 1 D 2 D p )

where A 2 Q l T = ( D 1 , D 2 , , D p ) and D i has m l rows and r i columns. Since the diagonal dominance condition (38) is invariant under the permutation of indexes, B also satisfies (38). Hence, letting

B i = ( b j k ) , C i = ( c j k ) and D i = ( d j k )

for a fixed i ,1 i p , we have

j = 1 , j k r i | b j k | + j = 1 , j k r i | c j k | + j = 1 , j k r i | d j k | | b k k | ( 1 k r i ) , (39)

where the summations on j are taken along columns and C 1 = O for convenience. If C i O or D i O , then

j = 1 , j k r i | c j k | + j = 1 , j k r i | d j k | > 0 for some k , 1 k r i ,

and hence for this k,

j = 1 , j k r i | b j k | < | b k k |

by (39). Therefore it follows from Lemma 9 that

d e t B i 0, (40)

since B i is irreducible. If C i = O and D i = O , then we have the form of

d e t B = d e t B i × (⋯)

which also implies (40), because d e t B 0 . In any case, we have d e t B i 0 . Since these are true for all i ,1 i p , it follows from (381) that

d e t A 1 0.

this proves Lemma 10.

Lemma 11. If system (5) satisfies conditions (36) and (38), then the norm of solution x ( n ) such that | x ( n ) | = i = 1 m | x i ( n ) | , is non-increasing, and consequently the zero solution is U.S..

For the proof, we can see (cf. [5] ).

In the following theorem, we can prove that the zero solution is U.A.S., if A ( n ) is bounded on Z and if condition (36), (37) and (38) are satisfied.

Theorem 3. In system (5), let A ( n ) be bounded on Z. Assume that conditions (36) and (38) are satisfied for all n Z and that there is a constant α > 0 such that

| d e t A ( n ) | α , on Z .

Then the zero solution is U.A.S..

Proof. As is stated in Lemma 11, the zero solution is U.S., and hence it is sufficient to show that for any ϵ > 0 there exists a T ( ϵ ) > 0 such that

| x ( n , n 0 , x 0 ) | < ϵ at some n , n 0 n n 0 + T ( ϵ ) ,

whenever | x 0 | 1 . Suppose that this is not true. Then there exists a constant ϵ > 0 , a sequence of solution { x ( k ) ( n ) } of (5) and a sequence n k such that

| x ( k ) ( n ) | ϵ on n k n n k + k 2

and

| x ( k ) ( n k ) | 1.

Since | x ( k ) ( n ) | is non-increasing, we have

ϵ | x ( k ) ( n ) | 1 , on n k n n k + k 2

and there exists a subinterval [ s k , s k + k ] of [ n k , n k + k 2 ] such that

| | x ( k ) ( n ) | | x ( k ) ( n ) | | < 1 k , for s k n , n s k + k .

Set φ ( k ) ( n ) = x ( k ) ( n + s k ) . Then, we obtain

φ ( k ) ( n + 1 ) = A ( n + s k ) φ ( k ) ( n ) , (41)

ϵ | φ ( k ) ( n ) | 1 , on 0 n k (42)

and

| | φ ( k ) ( n ) | | φ ( k ) ( n ) | | < 1 k , for 0 n , n k . (43)

Since A ( n ) is bounded, it follows from (41) and (42) that { φ ( k ) ( n ) } is uniformly bounded on any finite interval of Z, and thus, taking a subsequence, { φ ( k ) ( n ) } can be assumed to converge uniformly on any finite interval of Z. Defining y ( n ) by

y ( n ) = lim k φ ( k ) ( n ) ,

it follows from (42) and (43) that there is a constant β > 0 , ϵ β 1 , such that

| y ( n ) | = β , for n 0 .

Since y ( n ) = ( y 1 ( n ) , y 2 ( n ) , , y m ( n ) ) is defined on Z, we can choose an interval I = [ σ 1 , σ 2 ] (for some σ i Z ( i = 1 , 2 ) , 0 < σ 1 < σ 2 ) such that

y j 1 ( n ) , y j 2 ( n ) , , y j h ( n ) > 0 on I ,

y j h + 1 ( n ) , y j h + 2 ( n ) , , y j l ( n ) < 0 on I

and

y j l + 1 ( n ) y j l + 2 ( n ) y j m ( n ) 0 on I .

Here we note that { j 1 , j 2 , , j l } ϕ because y ( n ) 0 . Then

| y ( n ) | = p = 1 h y j p ( n ) p = h + 1 l y j p ( n ) = β on I .

Let

f ( k ) ( n ) = p = 1 h φ j p ( k ) ( n ) p = h + 1 l φ j p ( k ) ( n ) . (44)

Then, we have

lim k f ( k ) ( n ) = | y ( n ) | = β on I ,

and there is a sequence { θ k } I such that

lim k Δ f ( k ) ( θ k ) = 0. (45)

Moreover, since I is compact and A ( n ) is bounded on Z, we can assume that

lim k θ k = θ , for some θ I

and

lim k A ( s k + θ k ) = B , for some m × m matrix B .

Clearly B satisfies (36), (38) and | d e t B | α . Taking a difference of both sides of (44) at n = θ k and using relation (41), we find

Δ f ( k ) ( θ k ) = p = 1 h ( q = 1 m a j p j q ( s k + θ k ) φ j q ( k ) ( θ k , k ) ) p = h + 1 l ( q = 1 m a j p j q ( s k + θ k ) φ j q ( k ) ( θ k ) ) ,

and

lim k Δ f ( k ) ( θ k ) = p = 1 h ( q = 1 l b j p j q y j q ( θ ) ) p = h + 1 l ( q = 1 l b j p j q y j q ( θ ) ) ,

where B = ( b i j ) , since we have

lim k φ j ( k ) ( θ k ) = y j ( θ ) , for 1 j m

and

y j q ( θ ) = 0 for l + 1 q m .

By (45), we have

0 = q = 1 l ( p = 1 h b j p j q p = h + 1 l b j p j q ) y j q ( θ ) . (46)

Since B satisfies (36) and (38),

p = 1 h b j p j q p = h + 1 l b j p j q { 0 , for 1 q h , 0 , for h + 1 q l .

Therefore each term of right hand side of (46) is non-positive, because we have

y j q { > 0 , for 1 q h , < 0 , for h + 1 q l .

Then it follows from (46) that

( p = 1 h b j p j q p = h + 1 l b j p j q ) y j q ( θ ) = 0 , for 1 q l ,

which implies

p = 1 h b j p j q p = h + 1 l b j p j q = 0 , for 1 q l ,

since | y j q ( θ ) | 0 for 1 q l . Thus we have

det ( b j 1 j 1 b j 1 j 2 b j 1 j l b j 2 j 1 b j 2 j 2 b j 2 j l b j l j 1 b j l j 2 b j l j l ) = 0. (461)

On the other hand, B satisfies (38) and | det B | α > 0 , and thus it follows from Lemma 10 that all principal minors of B are nonsingular, which contradicts (461). This proves that the zero solution of system (5) is U.A.S..

Corollary 1. If system (5) is defined only for n 0 and all assumptions of Theorem 3 are satisfied for n 0 , then the zero solution is U.A.S. for n 0 .

Proof. We construct the system defined on Z by

x ( n + 1 ) = A 0 ( n ) x ( n ) , x R m , (47)

where

A 0 ( n ) = { A ( n ) , for n 0 , A ( 0 ) , for n < 0.

Since system (47) satisfies all assumptions of Theorem 3 on Z, the zero solution is U.A.S. on Z, and furthermore, since system (5) coincides with system (47) for n 0 , this prove our conclusion.

6. Application

Before Example 1, we state the following lemma is a special case of Theorem 3 in [13] .

In the nonlinear system

x ( n + 1 ) = F ( n , x ( n ) ) , (48)

let F ( n , x ) be almost periodic in n uniformly for x R m and for any r > 0 , let there exists a constant L * = L * ( r ) > 0 such that

| F ( n , x ) F ( n , y ) | L * | x y | for | x | , | y | r and n Z .

Lemma 12. If x ( n ) is a bounded solution of (48) on Z and if for any solution y ( n ) of (48), | x ( n ) y ( n ) | is monotone decreasing to zero as n , then x ( n ) is a unique almost periodic solution and its module is contained in the module of F ( n , x ) .

Example 1. Consider the variational linear difference equation

w ( n + 1 ) = A ( n ) w ( n ) , (49)

where

A ( n ) = ( p ( n ) α 2 p ( n ) 2 α q ( n ) α 2 ) .

We now assume that ( u 0 ( n ) , v 0 ( n ) ) is at least one bounded solution of (49) on Z and ( u ( n ) , v ( n ) ) is any solution of (49) on Z , and p ( n ) , q ( n ) are some bounded functions on Z such that

p ( n ) α / 2 and β 1 q ( n ) β 2 ,

for some positive constants α , β 1 and β 2 such that α 2 2 β 2 . We can verify that A ( n ) satisfies all assumptions in Corollary 1. First of all, A ( n ) is bounded in the future for n Z , because p ( n ) and q ( n ) are bounded function on Z . It is clear that the diagonal elements of A ( n ) are negative and

d e t A ( n ) = q ( n ) β 1 on Z .

The diagonal dominance condition (38) for A ( n ) requires that

p ( n ) + | p ( n ) 2 α q ( n ) | 0, and α 2 + | α 2 | 0,

which is equivalent to

q ( n ) α p (n)

and this is satisfied by

q ( n ) β 2 α 2 / 2 α p ( n ) .

Therefore, by Theorem 3, the zero solution of (49) is U.A.S. and

| u 0 ( n ) u ( n ) | + | v 0 ( n ) v ( n ) | 0 as n ,

where the convergence is monotone decreasing by Lemma 11. Thus, applying Lemma 12 to system (49), we find that there exists a unique almost periodic solution with the module contained in the module of A ( n ) .

7. Conclusion

In this paper, we obtain the existence and stability property of almost periodic solutions in discrete almost periodic systems. First, in Section 1, the research background is introduced. In Section 2, the fundamental concepts of the almost periodic solutions in discrete almost periodic systems is given. In Section 3, we are introduced to the several lemmas and have uniformly asymptotically stability theory of the linear system, and moreover, in Section 4, we consider the generalized gas almost periodic system, and if linear part is irreducible matrix, then we obtain the existence of almost periodic solutions of this system. Finally, in Sections 5 and 6, we consider and obtain an uniformly asymptotically stability criterion for solutions of a linear system with coefficient matrix of diagonal dominance conditions, and this result applies to meaningful example of a linear discrete system.

Cite this paper

Hamaya, Y. (2018) Existence and Stability Property of Almost Periodic Solutions in Discrete Almost Periodic Systems. Advances in Pure Mathematics, 8, 463-484. https://doi.org/10.4236/apm.2018.85026

References

  1. 1. Favard, J. (1933) Lecons les functions presque-periodiques. Gauthier-Villars, Paris.

  2. 2. Carleman, T. (1957) Problems Mathematiques dans la Theorie Cinetique des Gaz. Publ. Sci. Inst. Mittag-Leffler.

  3. 3. Jenks, R.D. (1968) Homogeneous Multidimensional Differential Systems for Mathematical Models. Journal of Differential Equations, 4, 549-565. https://doi.org/10.1016/0022-0396(68)90005-3

  4. 4. Elaydi, S. (2005) An Introduction to Difference Equations. 3rd Edition, Springer, Berlin.

  5. 5. Coppel, W.A. (1965) Stability and Asymptotic Behavior of Differential Equations. Heath Mathematical Monographs, Heath, Boston.

  6. 6. Nakajima, F. (1978) A Stability Criterion of Diagonal Dominance Type. SIAM Journal on Mathematical Analysis, 9, 815-824. https://doi.org/10.1137/0509063

  7. 7. Corduneanu, C. (1982) Almost Periodic Discrete Processes. Libertas Mathematica, 2, 159-169.

  8. 8. Fink, A.M. (1974) Almost Periodic Differential Equations. Lecture Notes in Mathematics 377, Springer-Verlag, Berlin.

  9. 9. Yoshizawa, T. (1975) Stability Theory and the Existence of Periodic Solutions and Almost Periodic Solutions. Applied Mathematical Sciences 14, Springer-Verlag, Berlin. https://doi.org/10.1007/978-1-4612-6376-0

  10. 10. Nakajima, F. (1976) Existence and Stability of Almost Periodic Solutions in Almost Periodic Systems. Publications of the Research Institute for Mathematical Sciences, Kyoto Univ., Kyoto, 12, 31-47. https://doi.org/10.2977/prims/1195190957

  11. 11. Krasnoselskii, M.A. (1964) Positive Solutions of Operator Equations. P. Noordhoff Ltd., Gronigen.

  12. 12. Massera, J.L. and Schaffer, J.J. (1958) Linear Differential Equations and Functional Analysis. Annals of Mathematics, 67, 517-573. https://doi.org/10.2307/1969871

  13. 13. Seifert, G. (1968) Almost Periodic Solutions and Asymptotic Stability. Journal of Mathematical Analysis and Applications, 21, 136-149. https://doi.org/10.1016/0022-247X(68)90248-5

  14. 14. Montandon, B. (1972) Almost Periodic Solutions and Integral Manifolds for Weakly Nonlinear Nonconservative Systems. Journal of Differential Equations, 12, 417-425. https://doi.org/10.1016/0022-0396(72)90015-0

  15. 15. Varga, R.S. (1962) Matrix Iterative Analysis. Prentice-Hall, Englewood Cliffs.