Dynamics of Nonlinear Interactions among Forces with Lagged Effect

Abstract

This paper is concerned with the dynamics of the steady state of a two-delayed nonlinear system of functional differential equations. The stability of the steady state together with its dependence on the magnitude of time delays has been examined by means of characteristic equation corresponding to the nonlinear equation. General criteria for stability involving the two-delay equations have been obtained.

Share and Cite:

Santos, J. and Azevedo, K. (2017) Dynamics of Nonlinear Interactions among Forces with Lagged Effect. Applied Mathematics, 8, 1437-1444. doi: 10.4236/am.2017.810104.

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

{ u ˙ ( t ) = f ( u ( t ) , v ( t ) ) v ˙ ( t ) = h ( u ( t ) , v ( t ) ) + g ( u ( t r ) , v ( t σ ) ) (1)

where f , h , g : R 2 R are sufficiently smooth functions, r and σ are non negative constants. Assume that f, h and g vanish at P 0 = ( u 0 , v 0 ) , f u ( P 0 ) = κ 11 , f v ( P 0 ) = κ 12 , h u ( P 0 ) = κ 21 , h v ( P 0 ) = κ 22 , g u ( P 0 ) = δ and g v ( P 0 ) = κ , 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 σ = r , 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 σ r to improve understanding of the combined effects of functional response and time delays on the dynamics of predator-prey systems. Point P 0 that is equilibrium point for the system 1, therefore also equilibrium of (2). Taking a delay as a paramenter ( σ > 0 ), our purpose is to relate the dynamics of the two systems in the neighbourhood of P 0 .

The linear problem around P 0 , associated to System (1), is given by

{ u ˙ ( t ) = κ 11 u ( t ) + κ 12 v ( t ) v ˙ ( t ) = κ 21 u ( t ) + κ 22 v ( t ) + δ u ( t r ) + κ v ( t σ ) (2)

and its characteristic equation is given by

p ( λ ) = ( κ 11 λ ) [ κ 22 λ + κ e λ σ ] κ 12 ( κ 21 + δ e λ r ) = 0. (3)

By assuming that

κ 11 = 0 , κ 12 > 0 , κ 22 < 0 , κ 21 < δ < 0 and κ > 0 , (4)

the Equation (3) becomes

p ( λ ) = λ 2 d λ b λ e λ σ a e λ r + c = 0 , (5)

where a = δ κ 12 , b = κ , c = κ 12 κ 21 and d = κ 22 and a , b , c and d > 0 . If σ = r , we can choose a , b , c , d for the Equation (5) as the characteristic Equation (3.5) examined [2] .

For a moment assume that either d = 0 and c < a or a > 0 and c < 0 , so it is not true that all roots of the Equation (5) have negative real part neither, since on the real axis p ( 0 ) = c a < 0 and p ( λ ) 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, ( 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 ( K ) , 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 .

2. Stability

Let λ = x + i y be a solution of Equation (5). Separating real and imaginary part in (5) one obtains the following equations system for x and y:

{ x 2 ( y 2 c ) d x b e σ x [ x c o s y σ + y s i n y σ ] a e r x c o s y r = 0 2 x y d y + b e σ x [ x s i n y σ y c o s y σ ] + a e r x s i n y r = 0 (6)

In order to simplify notation we set

x ¯ = b + b 2 + 2 a 2 , y ¯ = x ¯ 2 + b 2 d 2 + ( x ¯ 2 + b 2 d 2 ) 2 + a + c . (7)

Theorem 1. Assume that (4) is satisfied and

d < min { 2 b / 2 , 2 c } , 2 r < 2 σ < π , (7a)

4 a + ( b + d ) 2 < 4 c , y ¯ < 2 , (7b)

b < 2 2 ϒ ( a , c ) ( d 2 2 c ) , where ( 8 c 2 + a 2 ) ϒ ( a , c ) = c 2 a 2 , (7c)

Then, there is ϵ 0 > 0 so that if b < 2 ϵ 0 , all roots of Equation (5) have negative real part.

Proof: Observe that a , b and c are positive, d is negative and c > a (see (4)). Suppose that y = 0 in (6). The system has been reduced to x 2 ( b e σ x + d ) x a e r x + c = 0 . If ϑ ( x ) = x 2 ( b e σ x + d ) x a e r x + c then ϑ ( 0 ) > 0 . If x > 0 , the first inequality in (7b) results in ϑ ( x ) > x 2 ( b + d ) x a + c > 0 . From this it follows that p ( λ ) = 0 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

{ b y sin y σ + a cos y r = ( y 2 c ) b y cos y σ + a sin y r = d y (8)

We will show that system (6) has no purely imaginary solution. Suppose y 0 is a solution of (8), then one must have

sin ( ( σ r ) y ) = ( c y 2 ) 2 + ( d 2 b 2 ) y 2 a 2 2 a b y : = ρ ( y ) . (9)

It is clear that ρ ( y ) = ρ ( y ) , thus we will consider only y > 0 . A direct calculation shows that 2 b a y 2 ρ ˙ ( y ) = 3 y 4 + ( d 2 b 2 2 c ) y 2 ( c 2 a 2 ) and that

y 0 6 = [ ( d 2 ( 2 c + b 2 ) ) 2 + 12 ( c 2 a 2 ) ( d 2 ( 2 c + b 2 ) ) ] 1 2 , (10)

is a minimum global point of ρ on the interval ( 0, ) .

Since a b y 2 ρ ˙ ( y ) = a b y ρ ( y ) + y 4 ( c 2 a 2 ) and ρ ˙ ( y 0 ) = 0 , then 3 a b y 0 ρ ( y 0 ) = 2 ( c 2 a 2 ) + ( d 2 2 c b 2 ) y 0 2 . From the first inequality in (7a) and the inequality of (7c) we notice that ρ ( y 0 ) > 1 and (5) have no purely imaginary root.

Let consider the sets S i for i { 0,1,2, ,6 } given by

S 0 = { ( x , y ) 2 : 0 < x x ¯ and 0 < y ¯ y } ; S 1 = { ( x , y ) 2 : 0 < 4 σ y π and 2 x > b } ; S 2 = { ( x , y ) 2 : π 4 σ y 2 π and 0 < y x } ; S 3 = { ( x , y ) 2 : π 2 σ y 2 π and x > 0 } ; S 4 = { ( x , y ) 2 : 0 < y x and y x ¯ } ; S 5 = { ( x , y ) 2 : 0 < x y and x x ¯ } ; S 6 = { ( x , y ) 2 : π 4 σ y 2 π and 0 < x y } . (11)

For each ( x , y ) R + 2 , the first equation in (6) becomes equivalent to

Γ 1 ( x , y ) = : x y d x c x + y b e x σ ζ 1 ( x , y ) e x r x + y a cos y r = 0 , (12)

where ζ 1 ( x , y ) = ( x + y ) 1 [ x cos y σ + y sin y σ ] . We notice that

| ζ 1 ( x , y ) | 1 and a e r x x + y a e r x c o s r y x + y a e r x x + y a e r x y . (13)

If ϵ * ( x , y ) = y 2 + ( b + x d ) y + a + c , then it follows from (13) that y Γ 1 ( x , y ) < ϵ * ( x , y ) . If ( x , y ) S 0 , we can verify that max { ϵ * ( x , y ) : ( x , y ) S 0 } < 0 and so, Equation (12) has no solution belonging to S 0 . In fact, function 2 y ( x ) = ( b + x d ) + ( b + x d ) 2 + 4 ( a + c ) is increasing in x and yeilds the positive solution for the second-degree equation at y given by ϵ * ( x , y ) = 0 .

The second equation in (6) is equivalent to

Γ 2 ( x , y ) = 2 x y d y + b e σ x [ x sin σ y y cos σ y ] + a e r x sin r y = 0. (14)

Since Γ 2 ( x , y ) > y ( 2 x b ) > 0 for ( x , y ) S 1 , there is no solution of the Equation (14) that belongs to S 1 . It is clear that x s i n σ y y c o s σ y 0 for ( x , y ) S 2 . Moreover, for all ( x , y ) S 2 , we have Γ 2 ( x , y ) 2 x y d y + a e r x sin r y > 0 . Thus, there is no solution of (12) that belongs to S 2 . Let consider ( x , y ) S 3 , then Γ ( x , y ) 2 x y d y is positive because d < 0 (see (7a)). Hence, there is no solution of (12) that belongs to S 3 . For each ( x , y ) S 4 , we define ζ 2 ( x , y ) = ( x y ) 1 [ x sin y σ y cos y σ ] . We can check directly that y ζ 2 ( x , y ) 2 and y s i n y r x . Then, Γ 2 ( x , y ) > ( x 1 y 3 ) [ 2 y 2 b y a ] is positive if y x ¯ (see (7)), hence (14) has no solution at S 4 . Let ( x , y ) S 5 . Because, x ζ 2 ( x , y ) 2 and ( x y ) 1 sin ( y r ) > x 2 , then x 2 Γ 2 ( x , y ) > x 1 y 1 ( 2 x 2 2 b x a ) is positive if x > x ¯ (see (7)) and so, has no solution at S 5 . Assume ( x , y ) S 6 , since x y x 2 , b y 2 / 2 < b y e σ x cos σ y , Γ 2 ( x , y ) > ( d y b y 2 / 2 ) is positive if d < b 2 / 2 (see (7a)). Hence, there is no solution for (12) that belongs to S 6 .

By general arguments on the compactness of the interval [ 0, π 4 ] and the continuity of ρ , we arrive at the existence of a ϵ 0 > 0 so that System (6) has no solution belonging to S ϵ 0 = { ( x , y ) R 2 : 0 x ϵ 0 , 0 σ y π 4 } . If R + 2 = { ( x , y ) R 2 , x > 0 , y > 0 } , conditions (7a), (7b), (7c) and b < 2 ϵ 0 imply that R + 2 S S ϵ 0 , where S = i = 0 6 S i (see (11)). So, it can be proven that System (6) has no solution that lies at C l ( R + 2 ) and so the proof of the Theorem is complete.

Corollary 1. Let a , b , c , d be given by (5), satisfing d < min { 2 b / 2 , 2 c } , (7b) and (7c). Then, there is ϵ 0 > 0 so that if y ¯ < ϵ 0 the equilibrium solution of System (1) is asymptotically stable for all σ and r (see (7)).

Proof: Let consider ϵ 0 and S ϵ 0 given in Theorem 1. It is easy to see that if y ¯ < ϵ 0 then S i S ϵ 0 for i { 1,2,3,4 } and in this case C l ( R + 2 ) S ϵ 0 S 0 S 5 S 6 is indepentent from r 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), d = 2 b / 2 , σ = π / 2 . Given a , b > 0 , consider r so that r < b a ( 1 2 2 ) < π 2 . For a , b and r there exists 0 < y < 1 / 2 , solution of the second equation of System (8). If c = b y sin ( π y / 2 ) + a cos ( r y ) + ( y ) 2 , then ( π / 2 , r , y ) is the solution of System 8. Moreover, ( λ ˙ ( π 2 ) | x = 0 ) > 0 .

Proof: Consider d = 2 b / 2 in the second equation of the System (8). We set H ( y ) = a sin ( r y ) and G ( y ) = y b ( 2 2 + cos ( π y / 2 ) ) . It is easy to see that G ( 0 ) = G ( 1 / 2 ) = 0 and G is concave with a maximum at y ¯ , y ¯ is near y = 0.285281 . Clearly, H ( y ) is an increasing function for 0 < y < 1 / 2 and H ( 0 ) = 0 . If H ( 0 ) = a r < b ( 1 2 2 ) = G ( 0 ) with r < π / 2 , there exists y , so that H ( y ) = G ( y ) . Thus, given a , b so that b a ( 1 2 2 ) < π 2 , consider r < b a ( 1 2 2 ) and for this r consider y given by H ( y ) = G ( y ) . It is

possible to find c so that y is the solution of b y s i n ( π y / 2 ) + a c o s ( r y ) + y 2 = c , which is the first equation of System 8 or else consider b y sin ( π y / 2 ) + a cos ( r y ) + ( y ) 2 = c . Therefore ( π / 2 , r , y ) is the solution of System 8 for a , b and c given above.

Now, by using implicit derivative in (5) we can show that

( λ ˙ ( σ ) | x = 0 ) = b y 2 [ 2 y sin σ y d cos σ y ] + b a r cos ( σ r ) y . (15)

( λ ˙ ( π / 2 ) | x = 0 ) = b ( y ) 2 [ 2 y sin π / 2 y d cos π / 2 y ] + b a r cos ( π / 2 r ) y > 0 is positive because 0 < y < 1 / 2 and a , b , d > 0 .

A simple observation yields the following important corollary.

Corollary 2. Let a , b , c , d satisfy the hypotheses of Theorem 2. Then for each ( k , m ) Z × Z there corresponds ( σ k , r m ) given by, y σ k = y π / 2 + 2 k π and y r m = y r + 2 m π so that ( σ k , r m , y ) solves System 8. Moreover, ( λ ˙ ( σ k ) | x = 0 ) > 0 and the equilibrium solution of System (1) undergoes a Hopf bifurcation at the points σ k , k = 0 , 1 , 2 , .

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

{ Y ˙ ( t ) = α [ I ( Y ( t ) , K ( t ) ) S ( Y ( t ) , K ( t ) ) ] K ˙ ( t ) = β I ( m ˜ Y ( t ) + ( 1 m ˜ ) Y ( t η ) , m K ( t τ ) + ( 1 m ) K ( t ) ) q K b 1 d 0 ( t ) (16)

where 0 < d 0 < b 1 , m ˜ 0 , 0 m 1 , Y is the gross product, K the capital stock, α > 0 the adjustment coefficient in the goods market, q ( 0,1 ) the depreciation rate of capital stock, I ( Y , K ) and S ( Y , K ) investment and saving are nonlinear and we assume are smooth enough fuctions, and η , τ > 0 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 a 0 , b 0 and c 0 be real numbers such as 0 < a 0 < c 0 , b 0 b 1 , a 0 + b 0 c 0 + d 0 . If γ * = a 0 d 0 b 0 c 0 ( b 1 d 0 ) ( a 0 c 0 ) 1 , we define

I ( Y , K ) = Y a 0 ( t ) K b 0 ( t ) q 1 K ( b 0 d 0 ) and S ( Y , K ) = Y c 0 ( t ) K d 0 ( t ) q 1 Y ( c 0 a 0 ) . (17)

Then, P 0 = ( Y 0 , K 0 ) = ( ( 1 + q β 1 ) 1 b 0 b 1 ; ( 1 + q β 1 ) d 0 b 1 b 0 ) is a solution of the non-linear system

{ I ( Y ( t ) , K ( t ) ) S ( Y ( t ) , K ( t ) ) = 0 , I ( Y ( t ) , K ( t ) ) q β 1 K ( b 1 d 0 ) ( t ) = 0 I Y ( Y ( t ) , K ( t ) ) S Y ( Y ( t ) , K ( t ) ) = 0. (18)

which is also a steady state for the system 16. The characteristic equation of the linear problem around P 0 , associated to System 16, is given by the Equation (5) where

a = ( I K ( P 0 ) S K ( P 0 ) ) ( 1 m ˜ ) I Y ( P 0 ) , b = m I K ( P 0 ) , c = ( I K ( P 0 ) S K ( P 0 ) ) m ˜ I Y ( P 0 ) and d = ( 1 m ) I K ( P 0 ) q .

If we assume 0 < ( 1 m ) I K ( P 0 ) < q , I K ( P 0 ) S K ( P 0 ) < 0 , I Y ( P 0 ) > 0 , and I K ( P 0 ) > 0 , then a, b, and c are positive real numbers and d < 0 . If a , b , c , d , r and σ satisfy all conditions of theorem 1 the steady state P 0 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 P 0 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.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Huang, W. (2013) Information Lag and Dynamic Stability. Journal of Mathematical Economics, 44, 513-529.
[2] Pal, P.J. and Mandal, P.K. (2014) Bifurcation Analysis of a Madified Leslie-Gower Predator-Prey Model with Beddington-DeAngelis Functional Response and Strog Allee Effect. Mathematics and Computers in Simulation, 97, 123-146. https://doi.org/10.1016/j.matcom.2013.08.007
[3] Rustichini, A. (1989) Hopf Bifurcation for Functional Differential Equations of Mixed. Journal of Dynamics and Differential Equations, 1, 145-177. https://doi.org/10.1007/BF01047829
[4] Santos, J.S. and Dos Reis, J.G. (1996) Unbounded Continum of Periodic Solutions for a Autonmous Delay Equation. Journal of Mathematical Analysis and Applications, 203, 840-849. https://doi.org/10.1006/jmaa.1996.0415
[5] Wu, X.P. (2012) Zero-Hopf Bifurcation Analysis of a Kaldor-Kalecki Model of Business Cycle with Delay. Nonlinear Analysis, Real World Appiled, 13, 736-754.
https://doi.org/10.1016/j.nonrwa.2011.08.013
[6] Aki, S.M.T. and Godoy, S.M.S. (1981) Permanence of Stability for a Class of System of Differential Equations with Two Delays. Nonliner Analysis, 1, 172-184.
[7] Hale, J.K. and Lunel, S.M.V. (1993) Introduction to Functional Differential Equations. Springer Verlag-New York. https://doi.org/10.1007/978-1-4612-4342-7
[8] Matsunaga, H. (2008) Delay-Dependent and Delay-Independent Stability Criteria for a Delay Differential System. Proceedings of the American Mathematical Society, 136, 4305-4312. https://doi.org/10.1090/S0002-9939-08-09396-9
[9] Ruan, S. and Wei, J. (1999) Periodic Solutions of Planar Systems with Two Delays. Physica D, 129-A, 1017-1032.
[10] van den Driessche, P. and Wu, J. (1981) Stabilization Role of Inhibitory Self-Connections in a Delayed Neuralnetwork. Physica D, 150, 84-90.
[11] Liu, Z. and Yuan, R. (2004) Stability and Bifurcation in a Delayed Predator-Prey System with Beddington-Deangelis Functional Response. Journal of Mathematical Analysis and Applications, 296, 521-527. https://doi.org/10.1016/j.jmaa.2004.04.051
[12] Yuan, Y. (2012) A Coupled Plankton System with Instantaneous and Delayed Predation. Journal of Biological Dynamics, 6, 148-165. https://doi.org/10.1080/17513758.2010.544409
[13] Faria, T. (2000) On a Planar System Modelling a Neuron Network with Memory. Journal of Differential Equations, 168, 129-49. https://doi.org/10.1006/jdeq.2000.3881

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.