journal_list | How to participate | E-utilities
Bae: Nonlinear Behavior in Love Model with Discontinuous External Force

Abstract

This paper proposes nonlinear behavior in a love model for Romeo and Juliet with an external force of discontinuous time. We investigated the periodic motion and chaotic behavior in the love model by using time series and phase portraits with respect to some variable and fixed parameters. The computer simulation results confirmed that the proposed love model with an external force of discontinuous time shows periodic motion and chaotic behavior with respect to parameter variation.

1. Introduction

Over the last three decades, nonlinear dynamics including chaotic dynamics has been studied in many science fields, such as mathematics, physics, chemistry, engineering, and social science. In particular, studies of chaotic behavior in love models [15] have increased.
The dynamical love model has been defined in several forms, such as the Romeo and Juliet model [1, 37], the Laura and Petrarch model [8], and the Adam and Eve model [9]. Among these love models, the Romeo and Juliet model is mostly used in the study of nonlinear dynamics.
The love model based on Romeo and Juliet with a linear differential equation was first proposed by Sprott [1, 2], who described its linear and nonlinear behaviors. While the other love models have been the subject of some research [37], the Romeo and Juliet-based love model has been the most studied. Bae and his associates [37] proposed that the existence of periodic and chaotic behaviors based on the Romeo and Juliet love model can be represented through time series and phase portraits with identical and different time delays and an external force.
In this paper, we propose a love model for Romeo and Juliet with a discontinuous external force. The discontinuous external force closely reflects the real love between a man and woman. The discontinuous external force implies that the relation of love between Romeo and Juliet is not always the same. Sometimes they do not have thoughts about each other or sometimes they do not have concerns for each other. In addition, we investigated the periodic motion and chaotic behavior in the proposed Romeo and Juliet love model, with a discontinuous external force, by using time series and phase portraits.

2. Love Model

2.1 Romeo and Juliet Love Model

The basic love model for Romeo and Juliet was proposed by Sprott [1] and can be described as the following Eq. (1):
(1)
dRdt=aR+bJdJdt=cR+dJ,
where a and b specify Romeo’s romantic style, and c and d specify Juliet’s style.
According to parameters a, b, c, and d, we divided behavior into four types, as proposed by Strogatz [10]. When a > 0 and b > 0, Romeo is an eager beaver. That is, Romeo is encouraged by his and Juliet’s feelings. When a > 0 and b < 0, Romeo is a narcissistic nerd. This implies that Romeo wants more of what he feels but retreats from Juliet’s feelings. When a < 0 and b > 0, Romeo is a cautious or secure lover. That is, Romeo retreats from his own feelings but is encouraged by Juliet’s. Finally, when a < 0 and b < 0, Romeo is a hermit. That is, Romeo retreats from both his and Juliet’s feelings.
Eq. (1) can be rewritten as Eq. (2) by using an advanced form of Eq. (1), as follows:
(2)
dRdt=aR+bJ((1-J)dJdt=cR((1-R)+dJ.
From Eq. (2), we can describe a novel love model with a discontinuous external force, as in the following Eq. (3):
(3)
dRdt=aR+bJ((1-J)+f(t)dJdt=cR((1-R)+dJ,
where f (t)= 5 sinπt1 is a discontinuous external force. In this paper, we gives 2 seconds per every 10 seconds as t1 to have intermittence.

3. Nonlinear Behavior Analysis

For Eq. (3), we investigated the chaotic behavior by using time series and phase portraits according to parameter variation.

3.1 a = −7, b = −2, c = 1, and d = 1

With this parameter set, we obtained the time series and phase portrait shown in Figures 1 and 2, respectively.
Figures 1 and 2 confirm that the model has a 1-periodic motion.

3.2 a = −3.5, b = −2, c = 1, and d = 1

The time series and phase portrait obtained using this parameter set are shown in Figures 3 and 4, respectively, which again confirm 1-periodic motion.

3.3 a = −3, b = −2, c = 1, and d = 1

The time series and phase portrait obtained using this parameter set are shown in Figures 5 and 6, respectively, which display 2-periodic motion.

3.4 a = −2.5, b = −2, c = 1, and d = 1

The time series and phase portrait obtained using this parameter set are shown in Figures 7 and 8, respectively, which indicate that the model has quasi-periodic motion.

3.5 a = −2, b = −2, c = 1, and d = 1

The time series and phase portrait obtained using this parameter set are shown in Figures 9 and 10, respectively. The figures confirm quasi-periodic motion.

3.6 a = −1.5, b = −2, c = 1, and d = 1

The time series and phase portrait of the model obtained using this parameter set are shown in Figures 11 and 12, respectively. The figures show a quasi-periodic motion or quasi-attractor.

3.7 a = −1.4, b = −2, c = 1, and d = 1

The time series and phase portrait of the model obtained using this parameter set are shown in Figures 13 and 14, respectively, which show that the model has a chaotic attractor.

3.8 a = −1.3, b = −2, c = 1, and d = 1

The time series and phase portrait of the model obtained using this parameter set are shown in Figures 15 and 16, respectively. The figures confirm that the model has a chaotic motion.

3.9 a = −1.1, b = −2, c = 1, and d = 1

The time series and phase portrait of the model obtained using this parameter set are shown in Figures 17 and 18, respectively. These figures confirm a chaotic attractor with different type compare to other chaotic attractor such as Lorenz attractor, Chuas attractor and Chen attractor.

3.10 a = −1.0, b = −2, c = 1, and d = 1

The time series and phase portrait obtained using this parameter set are shown in Figures 19 and 20, respectively. These figures confirm a quasi-periodic motion.
Figures 120 confirm that the behavior of the love model with a discontinuous external force changes from periodic motion to quasi-attractor, chaotic attractor, and quasi-periodic motion as confirmed by the various time series and phase portraits. Thus, these results confirm that the system has a typical chaotic dynamics.

4. Conclusion

This paper proposed a love model for Romeo and Juliet that is organized by a discontinuous external force. From the proposed love model, we obtained periodic and chaotic motions by time series and phase portrait graphs by varying the parameters.

Biography

ijfis-16-064f21.tif
Youngchul Bae received his B.S., M.S., and Ph.D. degrees from the Department of Electrical Engineering, Kwangwoon University, Korea in 1984, 1985, and 1997, respectively. He worked at KEPCO (Korea Electric Power Company) during 1986–1991, and also worked at KISTI (Korea Institute of Science and Technology Information) during 1991–1997 as a senior researcher. He is currently a professor at the Division of Electrical, Electronic Communication and Computer Engineering, Chonnam National University, Yeosu, Korea. His research interests include nonlinear dynamics, chaos dynamics, robot control, intelligent systems, and motor control. He is a member of KIEE, KIECS, and KIIS. Tel: +82-61-659-7315, Fax: +82-61-659-7310

Notes

Conflict of Interest

No potential conflict of interest relevant to this article was reported.

References

1. Sprott JC. Dynamics of love and happiness. In : presented at the Chaos and Complex Systems Seminar; Madison, WI. 2001.
[Google Scholar]
2. Sprott JC. Dynamical models of happiness. Nonlinear Dynamics, Psychology, and Life Sciences. 9(1):23–36. 2005; PMID: 15629066.
[Google Scholar]
3. Kim SW, Shon YW, Bae Y. Mathematical modelling of love and its nonlinear analysis. Journal of the Korea Institute of Electronic Communication Sciences. 9(11):1297–1304. 2014; http://dx.doi.org/10.13067/JKIECS.2014.9.11.1297. DOI: 10.13067/JKIECS.2014.9.11.1297.
[CrossRef] [Google Scholar]
4. Bae Y. Behavior analysis of dynamic love model with time delay. Journal of the Korea Institute of Electronic Communication Sciences. 10(2):253–260. 2015; http://dx.doi.org/10.13067/JKIECS.2015.10.2.253. DOI: 10.13067/JKIECS.2015.10.2.253.
[CrossRef] [Google Scholar]
5. Bae Y. Modified mathematical modelling of love and its behaviour analysis. Journal of the Korea Institute of Electronic Communication Sciences. 9(12):1441–1446. 2014; http://dx.doi.org/10.13067/JKIECS.2014.9.12.1441. DOI: 10.13067/JKIECS.2014.9.12.1441.
[CrossRef] [Google Scholar]
6. Huang L, Bae Y. Behavior analysis in love model of Romeo and Juliet with time delay. Journal of Korean Institute of Intelligent Systems. 25(2):155–160. 2015; http://dx.doi.org/10.5391/JKIIS.2015.25.2.155. DOI: 10.5391/JKIIS.2015.25.2.155.
[CrossRef] [Google Scholar]
7. Huang L, Bae Y. Comparative behavior analysis in love model with same and different time delay. Journal of Korean Institute of Intelligent Systems. 25(3):210–216. 2015; http://dx.doi.org/10.5391/JKIIS.2015.25.3.210. DOI: 10.5391/JKIIS.2015.25.3.210.
[CrossRef] [Google Scholar]
8. Rinaldi S. Laura and Petrarch: An intriguing case of cyclical love dynamics. SIAM Journal on Applied Mathematics. 58(4):1205–1221. 1998; http://dx.doi.org/10.1137/S003613999630592X. DOI: 10.1137/S003613999630592X.
[CrossRef] [Google Scholar]
9. Cresswell C. Mathematics and Sex. Crows Nest, Australia: Allen & Unwin;2003.
10. Strogatz SH. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering. Reading, MA: Addison-Wesley;1994.

Figure 1
Time series of Romeo and Juliet model when a = −7, b = −2, c = 1, and d = 1.
ijfis-16-064f1.tif
Figure 2
Phase portrait of Romeo and Juliet model when a = −7, b = −2, c = 1, and d = 1.
ijfis-16-064f2.tif
Figure 3
Time series of Romeo and Juliet model when a = −3.5, b = −2, c = 1, and d = 1.
ijfis-16-064f3.tif
Figure 4
Phase portrait of Romeo and Juliet model when a = −3.5, b = −2, c = 1, and d = 1.
ijfis-16-064f4.tif
Figure 5
Time series of Romeo and Juliet model when a = −3, b = −2, c = 1, and d = 1.
ijfis-16-064f5.tif
Figure 6
Phase portrait of Romeo and Juliet model when a = −3, b = −2, c = 1, and d = 1.
ijfis-16-064f6.tif
Figure 7
Time series of Romeo and Juliet model when a = −2.5, b = −2, c = 1, and d = 1.
ijfis-16-064f7.tif
Figure 8
Phase portrait of Romeo and Juliet model when a = −2.5, b = −2, c = 1, and d = 1.
ijfis-16-064f8.tif
Figure 9
Time series of Romeo and Juliet model when a = −2, b = −2, c = 1, and d = 1.
ijfis-16-064f9.tif
Figure 10
Phase portrait of Romeo and Juliet model when a = −2, b = −2, c = 1, and d = 1.
ijfis-16-064f10.tif
Figure 11
Time series of Romeo and Juliet model when a = −1.5, b = −2, c = 1, and d = 1.
ijfis-16-064f11.tif
Figure 12
Phase portrait of Romeo and Juliet model when a = −1.5, b = −2, c = 1, and d = 1.
ijfis-16-064f12.tif
Figure 13
Time series of Romeo and Juliet model when a = −1.4, b = −2, c = 1, and d = 1.
ijfis-16-064f13.tif
Figure 14
Phase portrait of Romeo and Juliet model when a = −1.4, b = −2, c = 1, and d = 1.
ijfis-16-064f14.tif
Figure 15
Time series of Romeo and Juliet model when a = −1.3, b = −2, c = 1, and d = 1.
ijfis-16-064f15.tif
Figure 16
Phase portrait of Romeo and Juliet model when a = −1.3, b = −2, c = 1, and d = 1.
ijfis-16-064f16.tif
Figure 17
Time series of Romeo and Juliet model when a = −1.1, b = −2, c = 1, and d = 1.
ijfis-16-064f17.tif
Figure 18
Phase portrait of Romeo and Juliet model when a = −1.1, b = −2, c = 1, and d = 1.
ijfis-16-064f18.tif
Figure 19
Time series of Romeo and Juliet model when a = −1.0, b = −2, c = 1, and d = 1.
ijfis-16-064f19.tif
Figure 20
Phase portrait of Romeo and Juliet model when a = −1.0, b = −2, c = 1, and d = 1.
ijfis-16-064f20.tif
Formats:
Article | 
PDF LinksPDF(1.4M) | PubReaderPubReader | EpubePub | 
Download Citation
Share  |
         
METRICS
1,739
View
12
Save
4
Cited-By
In This Page: