• Views 393
• Citations 0
• ePub 17
• PDF 215
`The Scientific World JournalVolume 2013 (2013), Article ID 730736, 6 pageshttp://dx.doi.org/10.1155/2013/730736`
Research Article

## Analysis of a Fractional-Order Couple Model with Acceleration in Feelings

1Department of Mathematics, Faculty of Sciences, Gaziantep University, 27100 Gaziantep, Turkey
2Department of Mathematics, Faculty of Sciences, Ankara University, 06100 Ankara, Turkey

Received 1 October 2013; Accepted 26 November 2013

Academic Editors: A. Atangana and A. M. A. El-Sayed

Copyright © 2013 Ilknur Koca and Nuri Ozalp. 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

A fractional-order nonlinear dynamical model of couple has been introduced. Upper bounds are obtained for a fractional-order nonlinear dynamical model. Also different from other models, a model with the order 2α is discussed. We are expecting an acceleration in feelings; that is why we increase the order of the derivative between . Stability analysis of the fractional-order nonlinear dynamical model of involving two persons is studied using the fractional Routh-Hurwitz criteria. By using stability analysis on fractional-order system, we obtain sufficient condition on the parameters for the locally asymptotic stability of equilibrium points. Finally, numerical simulations are presented to verify the obtained results.

#### 1. Introduction

The first noninteger order differentiation and integration notion was considered in 1695 by Leibniz and L’Hôpital. In a letter to L’Hôpital in 1695, Leibniz raised the following question: “Can the meaning of derivatives with integer order be generalized to derivatives with noninteger orders?” L’Hôpital was somewhat curious about that question and replied by another question to Leibniz: “What if the order will be 1/2?” After the letter was answered by Leibniz, fractional order in the concept of derivative was formed [1].

There are lots of topics on fractional modeling, but in recent decades the study of interpersonal relationships has begun to be popular. Interpersonal relationships appear in many contexts, such as in family, kinship, acquaintance, work, and clubs [2]. Mathematical modeling in interpersonal relationships is very important for capturing the dynamics of people, but there are few models in this area and models have been limited to integer order differential equations. Another interesting dynamic is marriage. Marriage has been studied scientifically for the past sixty years [3]. Researchers are trying to understand why some couples divorce, but others do not, and why, among those who remain married, some are happy and some are miserable with one another [4]. Since experiments in these areas are difficult to generate, mathematical models may play a role in explanation of the dynamics of a couple and behavioral features.

Recently, a fractional-order system for the dynamics of love affair between a couple has been considered [5]. In this paper, different from [5], a model with the order is discussed. We are expecting an acceleration in feelings; that is why we increase the order of the derivative between . Also, upper bounds are discussed for the system.

We begin by giving the definitions and properties of fractional-order integrals and derivatives [6].

#### 2. Preliminaries and Definitions

The three most common definitions for fractional derivative can be given as the Grünwald-Letnikov definition, the Riemann-Liouville definition, and the Caputo definition.

Definition 1. The Riemann-Liouville type fractional integral of order of a function is defined by where is the gamma function.

Definition 2. The Grünwald-Letnikov definition is given as

Definition 3. The Riemann-Liouville type fractional derivative of order of a function is defined by where and is the integer part of .

Definition 4. The Caputo type fractional derivative of order of a function is defined by where and is the integer part of .

Some properties of the Caputo derivative and the Riemann-Liouville derivative formulas are given below: We see that, contrary to the Riemann-Liouville approach, in the case of the Caputo derivative, there are no restrictions on the values ().

#### 3. Equilibrium Points and Their Locally Asymptotic Stability

In this section, we consider a fractional-order nonlinear two-dimensional system as follows: where is the fractional derivative of order .  , , , and () are real constants. These parameters are oblivion, reaction, and attraction constants. In the equations above, we assume that feelings decay exponentially fast in the absence of partners. The parameters specify the romantic style of individuals 1 and 2. In the beginning of relationships, because they have no feelings towards each other, initial conditions are considered zero.

We note that, with zero initial conditions, the following equation is valid:

In that case, the system can be considered as follows:

Let us make the following changes of variables:

Then, transformed system is given below: with initial conditions where , , , , and () are real constants.

Let and consider the system with the initial values Here,

To evaluate the equilibrium points, let from which we can get the equilibrium points for and .

The Jacobian matrix for the system given in (14) is where To discuss the local stability of the equilibrium of the system given by (14), we consider the linearized system at . The characteristic equation of the linearized system is of the form If is taken as , we have the following reduced equation: where According to the fractional Routh-Hurwitz criteria, we have the following theorem.

Theorem 5. If and , then the equilibrium point is locally asymptotically stable for all .

Proof. equilibrium of the system given by (12) is asymptotically stable if all of the eigenvalues, , , of , satisfy the following condition (negative real part) [7, 8]: For , the Routh-Hurwitz criteria are just and . The characteristic polynomial satisfies eigenvalues as below: Now, suppose that and are positive. It is easy to see that if the roots are real, they are both negative, and if they are complex conjugates, they have a negative real part.
Next, to prove the converse, suppose that the roots are either negative or have a negative real part. Then, it follows that . If the roots are complex conjugates, , which implies that is also positive. If the roots are real, then since both of the roots are negative, it follows that .

Theorem 6. Let be as given in (20). If , then the positive equilibrium point of the system given in (12) is unstable.

Proof. If , from Descartes' rule of signs, it is clear that the characteristic equation has at least one positive real root. So, the equilibrium point of the system given in (12) is unstable.

#### 4. Analysis of a Model with Upper Bounds

In this section, we consider fractional-order system with the order between : A detailed analysis of this model is given in [5]. With the help of the following lemmas, upper bounds are discussed for the system.

Before giving our results, we give some useful lemmas [9, 10].

Lemma 7. Let , , , and be positive constants. Then, where , ve .

Lemma 8. Let , , and , , with be nondecreasing; let be nondecreasing in a variable for every fixed . If then, for , one has where , .

Let and consider the system with the initial conditions and . Here, and . Now, upper bounds for a fractional-order nonlinear system are discussed with the following theorem.

Theorem 9. Let and and satisfy the following inequality: where and ( and ,. Then, one has the following upper bounds for system of fractional order: for , where

Proof. Since and are assumed to be continuous functions, every solution of the initial value problem (IVP) given by (23) is also a solution of the following integral system for : Moreover, every solution of integral system is a solution of the IVP [11]. Now, we derive from (28) and (33) the following: where
Using Hölder's inequality for with in (32), we get the inequality below: where . Then, we have By Lemma 7, the following inequality is obtained for : By using the last relation in (35), we obtain where Since we have the following inequality: and the function is nondecreasing on . Now, with an application of Lemma 8 to (37) combining with (33), upper bounds for a fractional-order nonlinear system are obtained.

#### 5. Numerical Simulation

In nonlinear dynamic systems, predictability can be possible with stability. Also relationship development would be predictable given the right parameters. In this model, parameters provide the condition for the locally asymptotic stability of equilibrium points by using stability analysis on fractional-order transformed system.

In this paper, we focus on couple dynamics depending on the parameters as below. Many scenarios are possible. But in this model, secure or cautious lover (individual 1 retreats from his own feelings but is encouraged by that of individual 2 ( and )) and hermit (individual 1 retreats from his own feelings and that of individual 2 ( and )) are considered.

Let Now, we consider the system Let the initial conditions be

After the system is transformed, the following system is obtained with the order of : Let the initial conditions be Positive equilibrium point for the problem (44) and (45) is calculated as The approximate solutions and (resp., govern the feelings of to and the feelings of to are displayed in Figure 1 for with acceleration in feelings. Figure 2 shows the asymptotic approximation of to the equilibrium point for . For the numerical solution of the system, we use the predictor corrector method [12].

Figure 1: The graphs of (above) and (below) with the order .
Figure 2: Approximation to the equilibrium point with the order .

We have demonstrated via numerical simulations that the fractional-order nonlinear couple model (42) and (43) can exhibit asymptotic behavior in the presence of nonlinearity for an appropriate set of model parameters. We have observed that the model approaches the equilibrium points asymptotically.

#### 6. Conclusions

In this paper, stability analysis of the fractional-order nonlinear dynamical model of couple is studied by using the fractional Routh-Hurwitz criteria. By using stability analysis on fractional-order system, sufficient condition on the parameters for the locally asymptotic stability of equilibrium points is obtained. A fractional-order nonlinear dynamical model of couple with the order has been formulated and analyzed. In the discussed model, acceleration is observed in the solution. Also upper bounds for a system with the order have been obtained.

Finally, we have demonstrated via numerical simulations that a fractional-order nonlinear model of couple can exhibit asymptotic behavior in the presence of an appropriate set of model parameters.

#### References

1. J. T. Machado, V. Kiryakova, and F. Mainardi, “Recent history of fractional calculus,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1140–1153, 2011.
2. K. Barley and A. Cherif, “Stochastic nonlinear dynamics of interpersonal and romantic relationships,” Applied Mathematics and Computation, vol. 217, no. 13, pp. 6273–6281, 2011.
3. J. M. Gottman, J. D. Murray, C. C. Swanson, R. Tyson, and K. R. Swanson, The Mathematics of Marriage, The MIT Press, Cambridge, Mass, USA, 2002.
4. T. C. Martin and L. L. Bumpass, “Recent trends in marital disruption,” Demography, vol. 26, no. 1, pp. 37–51, 1989.
5. N. Ozalp and I. Koca, “A fractional order nonlinear dynamical model of interpersonal relationships,” Advances in Difference Equations, vol. 2012, article 189, 2012.
6. I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.
7. E. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models,” Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 542–553, 2007.
8. D. Matignon, “Stability results for fractional differential equations with applications to control processing,” Computational Engineering in Systems and Application Multiconference, vol. 2, pp. 963–968, 1996.
9. Q. -H. Ma and E. H. Yang, “Some new Gronwall-Bellman-Bihari type integral inequalities with delay,” Periodica Mathematica Hungarica, vol. 44, no. 2, pp. 225–238, 2002.
10. A. P. Prudnikov, Y. A. Brychkov, and O. L. Marichev, Integrals and Series of Elementary Functions, Nauka, Moscow, Russia, 1981.
11. V. Lakshmikantham and A. S. Vatsala, “Theory of fractional differential inequalities and applications,” Communications in Applied Analysis, vol. 11, no. 3-4, pp. 395–402, 2007.
12. K. Diethelm, N. J. Ford, and A. D. Freed, “A predictor-corrector approach for the numerical solution of fractional differential equations,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 3–22, 2002.