Dynamics of Nonlinear Interactions among Forces with Lagged Effect ()
1. Introduction
Differential equations models that incorporate the history of the phenomenum into the model have been extensively considered as model for many problems because of a commum feature to them, which is, the appearance of oscillations, a fact which is very important in the model problems coming from ecology and mathematical economics [1] [2] [3] [4] [5] . Since delay-differential equations share many properties with ordinary differential equations, we have use methods and techniques of geometric dynamical systems theory that have been implemented in functional differential equations to describe the dynamics of flow associated with the system of equations
(1)
where
are sufficiently smooth functions,
and
are non negative constants. Assume that f, h and g vanish at
,
,
,
,
,
and
, near an equilibrium point (see [2] [4] [6] [7] ). The model represented by system (1) has been considered because it appears as generator of self-sustaining cycles with time lags incorporated in the nonlinearities. Recently, dynamical systems with delays have been found in biology and economics, and delay dependence on the state of the system cannot be eliminated by any change of variables (see [2] [8] [9] [10] references cited there).
When
, the authors in [11] studied the stability and Hopf bifurcation of the ratio-dependent predator-prey system similar to (1). Yuan in [12] exploited the behaviuor of plankton-population models, one with instantaneous predation, another with delayed predation and has obtained conditions to guarantee the coexistence of two species, and addresses the stability and bifurcation under different density of fish, with or without maturation time delay. We have decided to keep to the model represented by system 1 with
to improve understanding of the combined effects of functional response and time delays on the dynamics of predator-prey systems. Point
that is equilibrium point for the system 1, therefore also equilibrium of (2). Taking a delay as a paramenter (
), our purpose is to relate the dynamics of the two systems in the neighbourhood of
.
The linear problem around
, associated to System (1), is given by
(2)
and its characteristic equation is given by
(3)
By assuming that
(4)
the Equation (3) becomes
(5)
where
,
,
and
and
and
. If
, we can choose
for the Equation (5) as the characteristic Equation (3.5) examined [2] .
For a moment assume that either
and
or
and
, so it is not true that all roots of the Equation (5) have negative real part neither, since on the real axis
and
as
, then system 1 will be unbounded solutions.
Because of the two delays, the analysis on how the roots of (5) locate with respect to the imaginary axis is a classical problem that, besides being important in itself, plays an important role in the study of asymptotic behavior in the theory of delayed differential equations (see [6] ). In view of the difficulty to find all the values of the parameters for which all the roots
in (5) have negative real parts, only the special situations of (5) are considered [1] [2] [6] [11] [13] . To analyze the behavior dynamics of System (1) close to equilibrium solution it is necessary to have detailed information about the behavior of the eigenvalue of the linear equation associated to it, and so this problem can be reduced to the fact that all roots of Equation (5) have negative real part. There have been many attempts to describe the region K in the parameters space in which the solutions of (5) approach equilibrium solution (see [7] section 11.2). In fact, each point that belongs to the boundary of K,
, is of primary interest because it represents the point where the equilibrium point of System (1) can undergo a bifurcation from stability to instability. In the analysis, we have described the region K for System (1) in Theorem 1 and found parameters belonging to
, which correspond to the point where the equilibrium of System (1) switches from being stable to unstable with periodic oscillations as the parameter
crosses the value
.
2. Stability
Let
be a solution of Equation (5). Separating real and imaginary part in (5) one obtains the following equations system for x and y:
(6)
In order to simplify notation we set
(7)
Theorem 1. Assume that (4) is satisfied and
(7a)
(7b)
(7c)
Then, there is
so that if
, all roots of Equation (5) have negative real part.
Proof: Observe that
and
are positive,
is negative and
(see (4)). Suppose that
in (6). The system has been reduced to
. If
then
. If
, the first inequality in (7b) results in
. From this it follows that
and so has no positive solution with null imaginary part. The roots of the system (6) with null real part must be roots of the system
(8)
We will show that system (6) has no purely imaginary solution. Suppose
is a solution of (8), then one must have
(9)
It is clear that
, thus we will consider only
. A direct calculation shows that
and that
(10)
is a minimum global point of
on the interval
.
Since
and
, then
. From the first inequality in (7a) and the inequality of (7c) we notice that
and (5) have no purely imaginary root.
Let consider the sets
for
given by
(11)
For each
, the first equation in (6) becomes equivalent to
(12)
where
. We notice that
(13)
If
, then it follows from (13) that
. If
, we can verify that
and so, Equation (12) has no solution belonging to
. In fact, function
is increasing in
and yeilds the positive solution for the second-degree equation at
given by
.
The second equation in (6) is equivalent to
(14)
Since
, there is no solution of the Equation (14) that belongs to
. It is clear that
for
. Moreover, for all
, we have
. Thus, there is no solution of (12) that belongs to
. Let consider
, then
is positive because
(see (7a)). Hence, there is no solution of (12) that belongs to
. For each
, we define
. We can check directly that
and
. Then,
is positive if
(see (7)), hence (14) has no solution at
. Let
. Because,
and
, then
is positive if
(see (7)) and so, has no solution at
. Assume
, since
,
,
is positive if
(see (7a)). Hence, there is no solution for (12) that belongs to
.
By general arguments on the compactness of the interval
and the continuity of
, we arrive at the existence of a
so that System (6) has no solution belonging to
. If
, conditions (7a), (7b), (7c) and
imply that
, where
(see (11)). So, it can be proven that System (6) has no solution that lies at
and so the proof of the Theorem is complete.
Corollary 1. Let
be given by (5), satisfing
, (7b) and (7c). Then, there is
so that if
the equilibrium solution of System (1) is asymptotically stable for all
and
(see (7)).
Proof: Let consider
and
given in Theorem 1. It is easy to see that if
then
for
and in this case
is indepentent from
and
(see (11)).
3. Hopf Bifurcation
A search for purely imaginary solutions of Equation (6) plays a key role in the analysis of stability and bifurcation of periodic solution of System (1). We will find solutions for Equation (6) with real part null, which actully solutions of the System (8).
Theorem 2. Assume (4),
,
. Given
, consider
so that
. For
and
there exists
, solution of the second equation of System (8). If
, then
is the solution of System 8. Moreover,
.
Proof: Consider
in the second equation of the System (8). We set
and
. It is easy to see that
and
is concave with a maximum at
,
is near
. Clearly,
is an increasing function for
and
. If
with
, there exists
, so that
. Thus, given
so that
, consider
and for this
consider
given by
. It is
possible to find
so that
is the solution of
, which is the first equation of System 8 or else consider
. Therefore
is the solution of System 8 for
and
given above.
Now, by using implicit derivative in (5) we can show that
(15)
is positive because
and
.
A simple observation yields the following important corollary.
Corollary 2. Let
satisfy the hypotheses of Theorem 2. Then for each
there corresponds
given by,
and
so that
solves System 8. Moreover,
and the equilibrium solution of System (1) undergoes a Hopf bifurcation at the points
.
4. Two-Delayed Kaldor-Kalecki Model
Consider a Kaldor-Kalecki model of business cycle with two delays where memory-dependence may be not only due to regulation, but also to the inertia of institutional, technical systems including the deployment of research and development results described by
(16)
where
,
,
,
is the gross product,
the capital stock,
the adjustment coefficient in the goods market,
the depreciation rate of capital stock,
and
investment and saving are nonlinear and we assume are smooth enough fuctions, and
are time lags representing delays for the investment in the capital stock due to past investment decision, which we assume are in both delays the gross product and the human or physical capital stockrespectively. Let
,
and
be real numbers such as
,
,
. If
, we define
(17)
Then,
is a solution of the non-linear system
(18)
which is also a steady state for the system 16. The characteristic equation of the linear problem around
, associated to System 16, is given by the Equation (5) where
If we assume
,
,
, and
, then a, b, and c are positive real numbers and
. If
and
satisfy all conditions of theorem 1 the steady state
of the Kaldor- Kalecki system 16 is asymptocaly stable. If the conditions of corollary 1 are satisfied the stability is delay-independent. If the conditions of theorem 2 are satisfied then the steady state
becomes unstable and undergoes a Hopf bifurcation. This Hopf bifurcation is the first one of a cascade taking place as the delays go to infinity. The structure of this cascade is described by corollary 2.
5. Concluding Remarks
This paper studies dynamic characteristics of feed-back effects that incorporate the memory of the phenomenum into the model because of a commum feature to them, which is, the appearance of oscillations, a fact which is very important in the model problems coming from mathematical economics or ecology (see 1). We have chosen an approach that makes it clear how the delay phenomena and anti-damping force dominate, lead either asymptotic stability or fluctuations to system 1, and it can also make trajectories converge to the other limit cycle which may be stable or not. We apply our results to describe Kaldor-Kaleck model dynamics 16, which can be used to model decision-making dynamics in business cycle. Ours results describe special circumstances which show that equilibrium point of system 1 can either have a stable fixed point surrounded by an unstable cycle that appears because Hopf bifurcation.