Discrete Dynamics in Nature and Society

Volume 2008, Article ID 595367, 18 pages

http://dx.doi.org/10.1155/2008/595367

## About the Stability and Positivity of a Class of Discrete Nonlinear Systems of Difference Equations

Department of Electricity and Electronics, Institute of Research and Development of Processes (IIDP), Faculty of Science and Technology, University of the Basque Country, P.O. Box 644, Leioa, 48080 Bilbao, Spain

Received 16 January 2008; Accepted 26 May 2008

Academic Editor: Antonia Vecchio

Copyright © 2008 M. De la Sen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper investigates stability conditions and positivity of the solutions of a coupled set of nonlinear difference equations under very generic conditions of the nonlinear real functions which are assumed to be bounded from below and nondecreasing. Furthermore, they are assumed to be linearly upper bounded for sufficiently large values of their arguments. These hypotheses have been stated in 2007 to study the conditions permanence.

#### 1. Introduction

There is a wide scientific literature devoted to investigate the properties of the solutions of nonlinear difference equations of several types [1–9]. Other equations of increasing interest are as follows:

(1)stochastic difference equations and systems (see, e.g., [10] and references therein);(2)nonstandard linear difference equations like, for instance, the case of time-varying coefficients possessing asymptotic limits and that when there are contributions of unmodeled terms to the difference equation (see, e.g., [11, 12]);(3)coupled differential and difference systems (e.g., the so-called hybrid systems of increasing interest in control theory and mathematical modeling of dynamic systems, [13–16] and the study of discretized models of differential systems which are computationally easier to deal with than differential systems; see, e.g., [17, 18]). In particular, the stability, positivity, and permanence of such equations are of increasing interest. In this paper, the following system of difference equations is considered [1]:

with , for all ; , , ; and , for all , under arbitrary initial conditions , for all . The identity allows the inclusion in a unified shortened notation via (1.1) of the dynamics:

as it follows by comparing (1.1) for with (1.2). The solution vector sequence of (1.1) will be denoted as , for all , under initial conditions , . The above difference system is very useful for modeling discrete neural networks which are very useful to describe certain engineering, computation, economics, robotics, and biological processes of populations evolution or genetics [1]. The study in [1] about the permanence of the above system is performed under very generic conditions on the functions , for all . It is only requested that the functions be bounded from below, nondecreasing, and linearly upper bounded for large values, exceeding a prescribed threshold, of their real arguments. In this paper, general conditions for the global stability and positivity of the solutions are investigated.

##### 1.1. Notation

, , . “” is the logic conjunction symbol. .
If ,
then is the transpose of *P*.

denote, respectively, *P* positive definite, semidefinite positive, negative definite, and
negative semidefinite. denote, respectively, *P* nonnegative (i.e., none of its entries is negative, also denoted as ), *P* positive (i.e., with at least one of its entries being positive),
and *P* strictly positive (i.e., all of its entries are positive). Thus, and ,
but the converses are not generically true. The same concepts and notation of
nonnegativity, positivity, and strict positivity will be used for real vectors.
Then, the solution vector sequence in of
(1.1) will be nonnegative in some interval *S*, denoted by (identical to ), for all ,
if all the components are nonnegative for .
If, in addition, at least one component is positive, then the solution vector is
said to be positive, denoted by (implying that ), for all .
If all of them are positive in *S*, then the solution vector is said to
be strictly on a discrete interval, denoted by (identical to and implying that and ), for all .

and are the and norms of vectors and induced norms of
matrices, respectively. is the *n*th identity matrix.

#### 2. Preliminaries

In order to characterize the properties of system (1.1), firstly define sets of nondecreasing and bounded-from-below functions in system (1.1) as follows irrespective of the initial conditions:

and sets of linearly upper bounded real functions:

for irrespective of the initial conditions as well. In a natural form, define also sets of nondecreasing, bounded-from-below, and linearly upper bounded real functions, again independent of the initial conditions, , that is,

for . The above definitions facilitate the potential restrictions on the functions , , required to derive the various results of the paper. The constraints on the functions , for all , used in the above definitions of sets, have been proposed by Stević for and then used to prove the conditions of permanence of (1.1) in [1] for some , , and , for all , subject to , , for all . The subsequent technical assumption will be then used in some of the forthcoming results.

*Assumption 2.1. * and .

The following two assertions are useful for the analysis of the difference system (1.1).

*Assertion 2.2. *For any given , , for all .

*Assertion 2.3. *(i) For any given , if , for all ,
for any real constants .
(ii) if , for all ,
for any real constants .(iii) if , for all ,
for any real constants .(iv) if Assumption 2.1 holds.

*Proof. *
Assertion 2.3
(i)–(iii) follow
directly from the definitions of and , for all .

Assertion 2.3
(iv) The proof is split into proving the two claims
below.*Claim 1. *.* Proof of Claim 1. * if if Assumption 2.1 holds.*Claim 2. *.*Proof of Claim 2. * if if Assumption 2.1 holds.

Then, Assertion 2.3
(iv)
has been proved from Claims 1-2.

The following result establishes that it is not possible to obtain equivalence classes from any collection of parts of the sets of functions in the definitions of , , and .

*Assertion 2.4. *For any ,
consider for some given 3-tuple in ,
and consider any discrete collection of distinct admissible triples () subject to the constraints and , for all ,
leading to the associated .
Define the binary relation in as . Then, is not an equivalence relation so that are not equivalence classes in with respect to .
Also, the sets and for any given respective collections , , , for all ,
are not equivalence classes, respectively, in and .

*Proof. * In view of Assertion 2.3(iv), can be all set equal to unity with no loss of
generality, which is done to simplify the notation in the proof. Note that for some .
Now, consider with such that . Then, . Since the equivalence classes with respect to any equivalence
relation are disjoint, in with respect to is not an equivalence class unless . Now, consider the linear function defined by so that . Thus, . Then, () are not equivalence relations, and there are
no equivalence classes in () with respect to (). The remaining part of the proof follows in
a similar way by using the definitions of the sets and ,
and it is omitted.

#### 3. Necessary Conditions for Stability and Positivity

Now, linear systems for system (1.1) with all the nonlinear functions in some specified class are investigated. Those auxiliary systems become relevant to derive necessary conditions for a given property to hold for all possible systems (1.1), whose functions are in some appropriate set , , or . This allows the characterization of the above properties under few sets of conditions on the nonlinear functions in the difference system (1.1). If , for all , then the auxiliary linear system to (1.1) is

If , for all , then the auxiliary linear system to (1.1) is

System (3.1)
may be equivalently rewritten as follows by defining the state vector sequence , for all ,
as the *k*th-order difference system:

with initial conditions for , where and

The one-step
delay may be removed by defining the following extended 2*k*th-order
system of state vector
satisfying

with and

Note that the extended system (3.8)-(3.9) is fully equivalent to system (3.3)–(3.7) since both have identical solutions for each given common set of initial conditions. Now, let be the -norm of real vectors of any order and associated induced norms of matrices (i.e., spectral norms of vectors and matrices). The following definitions are useful to investigate (1.1).

*Definition 3.1. *System (1.1) is said to be globally Lyapunov stable (or simply globally stable)
if any solution is bounded for any finite initial conditions.

*Definition 3.2. *System (1.1) is said to be permanent if any solution enters a compact set for for any bounded initial conditions with depending on the initial conditions.

*Definition 3.3. *System (1.1) is said to be positive if any solution is nonnegative for any
finite nonnegative initial conditions.

The system is
locally stable around an equilibrium point if any solution with initial
conditions in a neighborhood of such an equilibrium point remains bounded.
Local or global asymptotic stability to the equilibrium point occurs,
respectively, under local or global stability around a unique equilibrium point
if, furthermore, any solution tends asymptotically to such an equilibrium point
as .
Definition 3.2 is the definition of permanence in the sense used in [1],
which is compatible with global and local stability and with global or local
asymptotic stability according to Definition 3.1 and the above comments if .
However, it has to be pointed out that there are different definitions of
permanence, like, for instance, in [2], where vanishing solutions (related to
asymptotic stability to the equilibrium) or, even, negative solutions at
certain intervals are not allowed for permanence. On the other hand, note that
a continuous-time nonlinear differential system may be permanent without being
globally stable in the case that finite escape times *t* of the solution exist,
implying that because of unbounded discontinuities of the solution at finite
time *t*, that solution is unbounded in for some finite .
This phenomenon cannot occur for system (1.1) under the requirement , for all ,
which avoids the solution being infinity at
finite values of the discrete index *n* for any finite initial conditions.
The following result is concerned with necessary conditions of global Lyapunov stability of system (1.1) for all
the sets of functions , for all ,
since the linear system defined with , for all ,
in (1.1) has to be globally stable in order to keep global stability for any , for all .

Theorem 3.4. *System
(1.1) is globally stable and permanent for any given set of functions for any given and any given , for all ,
only if the subsequent properties hold. **
(i)
, for all . **
(ii)
,
equivalently, **where **,
and ,
with being the kth identity matrix. A
necessary condition is . *

*(iii) There exists*

*where , which is a solution to the matrix identity*

*for any given*

*Proof. * (i) Note that the identically zero
functions , for all , are all in for any , , , for all .
Proceed by contradiction by assuming that and for some ,
with the system being globally stable. Thus, if so that as ,
and then the system is unstable for some function .
Thus, the necessary condition for global stability has been proved, implying
also the permanence of all the solutions in some compact real interval .

(ii) Assume with everywhere in **R** so that , . Let the spectrum of *W* be ,
with each eigenvalue being repeated as many times as its multiplicity. Then, . It is first proved by complete
induction that if is an eigenvector of ,
then is an eigenvector of for any .
Assume that is an eigenvector of for some arbitrary and some eigenvalue .
Then, so that is also an eigenvector of for the same eigenvalue .
This property leads to

Proceed by
contradiction by assuming that system (1.1) is stable, for all ,
with . From (3.15), as ,
and then the system is unstable for a function for any real constant since it possesses an unbounded solution for some
finite initial conditions. Now, redefine the functions from the above , ,
as follows:

It is clear
by construction that if on an interval of infinite measure and if occurs on a real interval of finite measure,
then the above contradiction obtained for still applies for for any finite negative .
If occurs on an interval of finite measure and if occurs on an interval of infinite measure,
then the linear system resulting from (1.1) with the replacement is unstable so that any nontrivial solution is
unbounded. Furthermore, since as (function diverging to ) and being unbounded on **R** (implying that for being some monotonically increasing sequence of real numbers) are both
impossible situations for some since are all nondecreasing, it follows again that
the functions are bounded from below so that for some finite .
If the real subintervals within which equalizes or differs from are both of infinite measure, the result with some unbounded solution still applies
trivially for some finite .
Thus, system (1.1) is globally stable for any given set of functions for any and any , for all ,
only if the subsequent equivalent properties hold: , .
The necessary condition follows by inspecting the sum of entries of
the main diagonal of *W* which equalizes the sum of nonnegative real eigenvalues
of *W* (which are also the squares of the modules of the eigenvalues of ,
i.e., the squares of the singular values of ) which have to be all of modules not greater
than unity to guarantee global stability.

(iii) The
property derives directly from discrete Lyapunov global stability theorem and
its associate discrete Lyapunov matrix equation which has to possess a solution for any given . This
property is a necessary and sufficient condition for the global stability of the
extended linear system (3.8)-(3.9), and then
for that of system (3.3)–(3.7). The proof ends
by noting that system (3.8)-(3.9) has to be
stable in order to guarantee the global stability of system (1.1) for any set , for all ,
according to Property (ii).

Concerning positivity (Definition 3.3), it is well known that in the continuous-time and discrete-time linear and time-invariant cases, the positivity property may be established via a full characterization of the parameters (see, e.g., [2, 13, 17] as well as references therein). In particular, for a continuous-time linear time-invariant dynamic system to be positive, the matrix of dynamics has to be a Meztler matrix, while in a discrete-time one it has to be positive, where the control, output, and input-output interconnection matrices have to be positive in both (continuous-time and discrete-time) cases [2]. Under these conditions, each solution is always nonnegative all the time provided that all the components of the control and initial condition vectors are nonnegative [2, 13]. In general, in the nonlinear case, it is necessary to characterize the nonnegativity of the solutions over certain intervals and for certain values of initial conditions and parameters; that is, the positivity is not a general property associated with the differential system itself all the time but with some particular solutions on certain time intervals associated with certain constraints on the corresponding initial conditions. The positivity of (1.1) for linear functions is now invoked (in terms of necessary conditions) to guarantee the positivity of all the solutions of (1.1) for any set of nonnegative initial conditions and any potential set with for any given and any given , for all . This is formally addressed in the subsequent result.

Theorem 3.5. *System (1.1) is positive for any given set of nonnegative functions with for any given and any given , for all ,
only if , , , for all .*

*Outline of proof*

As argued in the proof of Theorem 3.4 for
stability, the linear system has to be positive in order to guarantee that it
is positive for any set with for any given and , for all .
The linear system (3.8)-(3.9) for is positive if and only if [3] since, in addition, this implies .
The proof follows since by direct inspection if and only if , , , for all .

Necessary joint conditions for stability, permanence, and positivity of (1.1) for any set with for any given and , for all , follow directly by combining Theorems 3.4 and 3.5.

#### 4. Main Stability Results

This section derives sufficiency-type conditions (easy to test) for global stability of the linear system (3.3)–(3.7) independently of the signs of the parameters , , and , (which are also allowed to take values out of the interval , but on their maximum sizes). It is allowed that be independent of the above parameters and negative, but fulfilling that their modules are less than unity. The mechanism of proof for the linear case is then extended directly to the general nonlinear system (1.1). The , , and , , are allowed to be negative but , , is required to formulate an auxiliary result for the main proof.

Theorem 4.1. *Assume
that , for all ,
and ** Then, the
linear system (3.3)–(3.7),
equivalently system (3.8)-(3.9), is
globally Lyapunov stable for any finite arbitrary initial conditions. It is
also permanent for any initial conditions: *

*Proof. *The successive use of the recursive
second identity in (3.3) with initial condition leads to

and taking -norms in (4.3) with ,
we get

where , ,
and since and

Note that
(4.4) is still valid if the term preceding the equality is any , for all ,
since they are all upper bounded by all the right-hand side upper bounds. Then,

for all , for all , for all , which implies directly that

If the
condition with holds, then the second term of the right-hand
side of (4.7) may be combined with the left-hand-side term to yield

for all , for all ,
depending on , which depends on , for any , for all .
Since is compact, it follows from (4.9) that any solution
sequence is bounded for any and any finite initial conditions. Thus, the linear
system (3.3)–(3.7) is globally
Lyapunov stable. Also, since is compact, it follows from (4.8) that any
solution sequence is permanent since it enters the prefixed compact set

for any and any finite initial conditions in .
Furthermore, is independent of each particular set of
initial conditions in . Thus,
the linear system (3.3)–(3.7) is
permanent.

The following technical result will be then useful as an auxiliary one to prove the stability of (1.1) under a set of sufficiency-type conditions based on extending the proof mechanism of Theorem 4.1 to the nonlinear case. Basically, it is proved that the functions , , grow at most linearly with their argument.

Lemma 4.2. *, for all .
In addition, is bounded for all .
The result also holds if , for all , for all .*

*Proof. * Now, it is
proved that (notation of “big Landau *O*” of *x*) for any , for all . First, note that for all for some , since for any , for all . The result
also holds if , for all ,
since . It is now proved by a contradiction argument
that if ,
then it is bounded, for all .
Assume with being arbitrarily large for some .
Thus, there exists being arbitrarily large so that for since so that it is monotonically nondecreasing.
This is a contradiction since is arbitrarily large. Thus, is bounded, for all .
Since it is bounded, then for some finite for as a result, so that on .
Again, the result still holds if , for all .

Theorem 4.3. *If , , , for all , for all , and , then system (1.1) is globally Lyapunov stable
for any finite arbitrary initial conditions. It is also permanent for any
initial conditions with the compact set defined in Theorem 4.1.*

*Proof. * If system (1.1) is taken, then
(4.4) is replaced with

where

The
description (4.6) is similar to (1.1) via an unforced linear system (3.3)–(3.7) with a
forcing sequence so that both solution sequences are identical
under identical initial conditions. One gets directly from (4.11) that

so that

Then by direct extension of (4.7) when using
(4.14),

with for some finite since , for all ,
from Lemma 4.2. Thus, may be regrouped in the left-hand side provided
that

Then, under
similar reasoning as that used to derive (4.8)-(4.9), one gets
from (4.15) that

for all , for all , depending on ,
which depends on ,
for any .
The solution sequences are all bounded under any finite initial conditions and
enter the compact set defined by

for all , for any set of initial conditions in the compact set .
Furthermore, is independent of each particular set of
initial conditions in .
Then, system (1.1) is globally Lyapunov stable and permanent.

Some simple properties concerning the instability of (1.1) based on simple constraints on the nonlinear functions, such as the stated boundedness from below of the strongest one of boundedness from above and below, are now established in the subsequent result.

Theorem 4.4. *The following properties hold. *

(i)*If and is bounded from above and below, then is bounded, for all . If and is bounded from above and below, then almost
all solution sequences for sufficiently large finite absolute values
of the initial conditions are unbounded. Thus, system (1.1) is unstable under
sufficiently large absolute values of the initial conditions for some .*(ii)*Assume
that and for some .
Then , for all ,
and as if ().
Thus, system (1.1) is unstable under such sufficiently large absolute values of
the initial conditions for some .*

*Proof. * (i) If , for all ,
for some , *j* = 1,2, and some ,
then

If ,
then the sequence is bounded so that the sequence may be unbounded only if .
If ,
then , and,
furthermore, if , then there
is a strictly monotonically increasing subsequence of ,
where is a countable subset of **N**, so that , for all ,
and as (i.e., it diverges).

If , then ,
for all and all ,
such that .

(ii) From (1.1), ,
and ,
it follows that if , for all ,
so that the absolute value of the solution sequence is monotonically increasing
so that it diverges. Less stringent condition for the initial conditions follows
by calculating the zeros of the convex function which are , which implies that if ,
and if .
This directly completes the proof.

#### 5. Positivity Results

Some positivity properties of the solution sequences of system (1.1) are now formulated in the subsequent formal result.

Theorem 5.1. * The following properties hold. *

(i)*Any solution vector sequence of (1.1) is nonnegative, for all ,
and any finite nonnegative , for all ,
if , for all ,
and*