Next Article in Journal
Coincidence Point of Edelstein Type Mappings in Fuzzy Metric Spaces and Application to the Stability of Dynamic Markets
Previous Article in Journal
Enhanced Remora Optimization with Deep Learning Model for Intelligent PMSM Drives Temperature Prediction in Electric Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bogdanov–Takens Bifurcation Analysis of a Learning-Process Model

1
College of Mathematics and Data Science, Minjiang University, Fuzhou 350108, China
2
HNU-ASU Joint International Tourism College, Hainan University, Haikou 570228, China
3
School of Economics and Management, Southwest Jiaotong University, Chengdu 610031, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(9), 853; https://doi.org/10.3390/axioms12090853
Submission received: 8 June 2023 / Revised: 17 August 2023 / Accepted: 29 August 2023 / Published: 31 August 2023

Abstract

:
In this paper, as a complement to the works by Monterio and Notargiacomo, we analyze the dynamical behavior of a learning-process model in a case where the system admits a unique interior degenerate equilibrium. Meanwhile, we acquire the sufficient condition for the cusp of codimension 2 and verify that the system undergoes Bogdanov–Takens bifurcation around the cusp. Finally, we give a numerical simulation to support the theoretical results.

1. Introduction

With the development of human society, educational issues [1,2,3,4,5] are gradually attracting people’s attention and are not ignored. Doubts [6] are a universal phenomena during the process of learning knowledge. The influence and significance of doubt in experience of acquiring knowledge are focused on by many scholars [7,8,9]. As time goes by, there exists a dynamic process of mutual transformation between understanding and doubt. That is, understanding can be transformed into doubt, and doubt can also be transformed into understanding [10]. Therefore, Monterio and Notargiacomo [11] proposed that the learning process as the interplay between understanding and doubt can be studied by formulating and analyzing a dynamical system written in term of differential equations. In order to better study this, the entirety of knowledge is divided into two parts: one part is already understood, and the other is still doubted. The first “understanding–doubt” model was put forward as follows.
U ˙ = a { U ( U 1 ) ( α U ) f U D } { 1 ( U + D ) } , D ˙ = b { D ( β D ) + g U D } { 1 ( U + D ) } ,
U and D describe the level of understanding and doubt with time t during the learning process, respectively. 0 < α < 1 and 0 < β < 1 stand for the minimum background required to learn about a subject and the maximum level of doubt that learner can have about a subject that they did not learn about it, respectively. a and b represent the speed of the learning process. The parameters f and g describe the interaction between U ( t ) and D ( t ) and f g > 0 . At the standing point of physical meaning, the term 1 ( U + D ) restricts the dynamics to the right triangle domain given by 0 U 1 , 0 D 1 and 0 U + D 1 . Authors have studied the stability of the boundary equilibria and provided some numerical simulations to verify the theoretical results.
Liu, Ding and Chen performed deep studies of System (1) in [12]. They gave a complete analysis of the qualitative properties of the interior equilibria and a singular line segment. They also investigated some local bifurcation, including transcritical, pitchfork, and Hopf bifurcation [13,14] on isolated equilibrium, and transcritical bifurcation without parameters on non-isolated equilibrium.
However, dynamical behavior of the learning-process model in a case where the system admits a unique interior equilibrium can be studied further. Meanwhile, bifurcation of a higher codimension, which is the Bogdanov–Takens bifurcation of codimension 2, is worth analyzing since, compared to bifurcation of codimension 1, bifurcation of codimension 2 demonstrates more realistic dynamic behavior between understanding and doubt in the learning process.
Compared to the works by Liu, Ding and Chen in [12] we will give the analysis of Bogdanov–Takens bifurcation [15,16,17] of codimension 2 when System (1) admits a unique interior equilibrium (cusp).
This paper is arranged as follows. We obtain the sufficient conditions for the existence of the cusp [18] of codimension 2 in Section 2. In Section 3, we prove that System (1) undergoes Bogdanov–Takens bifurcation of codimension 2 and the theoretical results are verified by numerical simulation. A brief ends this paper in Section 4.

2. The Existence of the Cusp of Codimension 2

For simplicity, we make a time rescaling d τ = a d t , and System (1) becomes
U ˙ = { U ( U 1 ) ( α U ) f U D } { 1 ( U + D ) } , D ˙ = r { D ( β D ) + g U D } { 1 ( U + D ) } ,
where r = b a .
For physical meaning, we only consider the dynamic behavior of System (2) in the closure Ω : = { ( U , D ) : 0 U 1 , 0 D 1 , 0 U + D 1 } for all possibilities of ( α , β , f , g , r ) Σ : = { 0 < α < 1 , 0 < β < 1 , f g > 0 , r > 0 } .
The equilibria of System (2) are determined by the following equation:
{ U ( U 1 ) ( α U ) f U D } { 1 ( U + D ) } = 0 , r { D ( β D ) + g U D } { 1 ( U + D ) } = 0 .
There exists a singular line segment
L : = { ( U , D ) : U + D = 1 , 0 U 1 , 0 D 1 } ,
and three boundary equilibria E 0 ( 0 , 0 ) , E 1 ( 0 , β ) , E 2 ( α , 0 ) . The interior equilibria of System (2) are determined by the following equation:
( U 1 ) ( a U ) f D = 0 , β D + g U = 0 .
Substituting D = g U + β into Equation (4), we obtain
F ( U ) = U 2 + ( f g α 1 ) U + f β + α
and
F ( U ) = 2 U + ( f g α 1 ) U .
From F ( U ) = 0 , we have
α = U f g + U 2 + β f U U 1 .
The Jacobian matrix of System (2) at any equilibrium E ( U , D ) of System (2) takes the form
J E = a b c d ,
where
a = 4 U 3 + 3 D 3 α 6 U 2 + 2 f D + 2 α + 2 + 2 α 1 D U f D + α 1 D , b = U 3 + f 1 α U 2 + f D 1 + α + f D U , c = r r g D 2 + r g 1 U r g U + β D , d = 3 r D 2 2 r g U + β + 2 r 1 U D + r g U + β 1 U ,
and D e t ( J E ) = a d b c , T r ( J E ) = a + d .
Substituting (5) into D e t ( J E ) and F ( U ) , we can rewrite D e t ( J E ) and F ( U ) as
D e t ( J E ) = r U g + β U g + U + β 1 2 U F ( U )
and
F ( U ) = U 2 β f f g 2 U + 1 U 1 .
Therefore, from F ( 0 ) > 0 and F ( 1 ) > 0 , we consider quadratic function F ( U ) to have double zeros U = U * if F ( U * ) = 0 , which means that System (2) admits a degenerate equilibrium E * ( U * , β + g U * ) . Furthermore, if T r ( J E * ) = 0 and D e t ( J E * ) = 0 , then we have
α = f g + 2 U * 1 = : α 1 , β = ( U * 2 f g 2 U * + 1 ) f = : β 1 , r = g U * f 2 ( f g U * + U * 2 f g 2 U * + 1 ) = : r 1 .
Therefore, System (2) admits a degenerate equilibrium E * ( U * , β 1 + g U * ) .
In order to ensure 0 < α 1 < 1 , 0 < β 1 < 1 and r 1 > 0 , the parameters f, g and U * should satisfy max 1 2 1 2 f g , 1 f g + f < U * < min 1 1 2 f g , 1 f g , g > 0 , f > 0 and 0 < f g < 1 .
Additionally, equilibrium E * ( U * , β 1 + g U * ) should lie on the right triangle domain Ω , so we have 1 f g f < U * < 1 . On the other hand, 1 f g + f < 1 f g f , and we need to ensure 1 f g > 1 f g f , thus f > g ( g + 1 ) 2 holds.
Summing up the above, we have the following theorem:
Theorem 1.
Based on g > 0 , f > g ( g + 1 ) 2 , 0 < f g < 1 and
max 1 2 1 2 f g , 1 f g f < U * < min 1 1 2 f g , 1 f g ,
there exists a unique degenerate equilibrium E * ( U * , β 1 + g U * ) in system (2).
Combining above results with Theorem 2.1 in [12], we have the following theorem:
Theorem 2.
Depending on condition for existence of equilibrium E * , boundary equilibria of system (2) have the following qualitative properties:
( i ) The origin E 0 : ( 0 , 0 ) is a saddle;
( i i ) Equilibrium E 1 : ( 0 , β 1 ) is a stable node;
( i i i ) Equilibrium E 2 : ( α 1 , 0 ) is an unstable node.
In the following, we determine type of equilibrium E * .
Theorem 3.
From condition for existence of degenerate equilibrium E * and
f g ( f g + U * 2 1 ) + 2 U * ( U * 1 ) 2 f g + U * 1 0 ,
equilibrium E * is a cusp of codimension 2.
Proof. 
Making a transformation ( x , y ) = ( U U * , D β 1 g U * ) , then system (2) becomes
d x d t = a 10 x + a 01 y + a 20 x 2 + a 11 x y + a 02 y 2 + O ( | ( x , y ) | 3 ) , d y d t = b 10 x + b 01 y + b 20 x 2 + b 11 x y + b 02 y 2 + O ( | ( x , y ) | 3 ) ,
where
a 10 = g U * U * 1 f g + f + U * 1 , a 01 = U * U * 1 f g + f + U * 1 , a 20 = U * 3 + f 2 U * 2 + f 2 g 2 2 f 2 g + f g f + 1 U * + f 2 g 2 + f 2 g f g f , a 11 = f g + 2 U * f + U * 2 f 2 U * + 1 , a 02 = U * f , b 10 = U * 1 f g + f + U * 1 g 2 U * , b 01 = g U * U * 1 f g + f + U * 1 , b 20 = f g 2 U * , b 11 = 2 f g 2 + 2 g U * 2 g U * + 1 g U * f f g + U * 1 , b 02 = 2 f g + f + 2 U * 2 g U * f f g + U * 1 .
Since a 01 0 , we perform linear transformation x = X , y = a 10 a 01 X + 1 a 01 Y , then System (6) becomes
d X d t = Y + c 20 X 2 + c 11 X Y + c 02 Y 2 + O ( | ( X , Y ) | 3 ) , d Y d t = d 20 X 2 + d 11 X Y + d 02 Y 2 + O ( | ( X , Y ) | 3 ) ,
where
c 20 = U * U * 1 f g + f + U * 1 f , c 11 = 2 f g U * f g + 2 U * f + U * 2 f 2 U * + 1 U * U * 1 f g + f + U * 1 , c 02 = f U * U * 1 2 f g + f + U * 1 2 , d 20 = U * 2 U * 1 2 f g + f + U * 1 2 g f , d 11 = g f g + f + U * 1 f g U * 2 + 2 U * 1 f g + U * 1 , d 02 = f g U * 1 f g + U * 1 .
By Remark 1 of Section 2.13 in [19], we obtain an equivalent system of System (7) in the small neighborhood of ( 0 , 0 ) as follows
d X d t = Y + O ( | ( X , Y ) | 3 ) , d Y d t = d ¯ 20 X 2 + d ¯ 11 X Y + O ( | ( X , Y ) | 3 ) ,
where
d ¯ 20 = d 20 = U * 2 U * 1 2 f g + f + U * 1 2 g f < 0 , d ¯ 11 = 2 c 20 + d 11 = f g + f + U * 1 f g ( f g + U * 2 1 ) + 2 U * ( U * 1 ) 2 f g + U * 1 f .
If
f g ( f g + U * 2 1 ) + 2 U * ( U * 1 ) 2 f g + U * 1 0 ,
then d ¯ 20 d ¯ 11 0 , which means that degenerate equilibrium E * ( U * , g U * + β 1 ) is a cusp of codimension 2. □
Setting α = 3 25 , β = 13 60 , f = 3 5 , g = 1 5 , r = 18 95 , the phase portrait of cusp of codimension 2 of is described with the help of Matlab; see Figure 1.
Remark 1.
From [11], if doubt destroys comprehension, then f > 0 and g > 0 ; if doubt drives comprehension, then f < 0 and g < 0 . Therefore, according to sufficient condition for existence of cusp of codimension 2, doubt destroys comprehension.

3. Analysis of Bogdanov–Takens Bifurcation

From Section 2, we know that E * is a cusp of codimension 2. In this section, we investigate Bogdanov–Takens bifurcation around E * .
Theorem 4.
Under the condition for the existence of equilibrium E * , which is a cusp of codimension 2, we choose α and r as bifurcation parameters. As the parameters ( α , r ) disturb around ( α 1 , r 1 ) , system (2) undergoes Bogdanov–Takens bifurcation in a small neighborhood of E * .
Proof. 
We choose α and r as two bifurcation parameters and replace α and r with α 1 + ε 1 and r 1 + ε 2 , respectively, in System (2). Then, we acquire the unfolding system of System (2):
U ˙ = { U ( U 1 ) ( α 1 + ε 1 U ) f U D } { 1 ( U + D ) } , D ˙ = ( r 1 + ε 2 ) { D ( β 1 D ) + g U D } { 1 ( U + D ) } ,
where ε = ( ε 1 , ε 2 ) is a parameter vector in a small neighborhood of origin. Obviously, when ε 1 = ε 2 = 0 , System (9) has a unique positive equilibrium E * , which is a cusp of codimension 2. For simplicity, the universal unfolding will be obtained with the help of a series of transformation in the following.
Making a translation transformation u 1 = U U * , v 1 = D β 1 g U * , equilibrium E * will be translated into ( 0 , 0 ) , then System (9) becomes
d u 1 d t = a ^ 00 ( ε ) + a ^ 10 ( ε ) u 1 + a ^ 01 ( ε ) v 1 + a ^ 20 ( ε ) u 1 2 + a ^ 11 ( ε ) u 1 v 1 + a ^ 02 ( ε ) v 1 2 + O ( | ( u 1 , v 1 ) | 3 ) , d v 1 d t = b ^ 10 ( ε ) u 1 + b ^ 01 ( ε ) v 1 + b ^ 20 ( ε ) u 1 2 + b ^ 11 ( ε ) u 1 v 1 + b ^ 02 ( ε ) v 1 2 + O ( | ( u 1 , v 1 ) | 3 ) ,
where
a ^ 00 ( ε ) = U * ε 1 U * 1 2 f g + f + U * 1 f , a ^ 10 ( ε ) = U * 1 2 U * f g f g + 3 f U * + 2 U * 2 f 3 U * + 1 ε 1 f U * 1 f 2 g 2 U * + f 2 g U * + f g U * 2 U * f g f , a ^ 01 ( ε ) = U * U * 1 ε 1 + U * U * 1 f g + f + U * 1 , a ^ 20 ( ε ) = U * f g f g + 3 f U * + U * 2 2 f 2 U * + 1 ε 1 f U * 3 ( 2 f ) U * 2 ( f 2 g 2 + 2 f 2 g f g + f 1 ) + f 2 g 2 + f 2 g f g f , a ^ 11 ( ε ) = 1 2 U * ε 1 f g + 2 f U * + U * 2 f 2 U * + 1 , a ^ 02 ( ε ) = f U * , b ^ 10 ( ε ) = U * 1 f g + f + U * 1 g U * f g f g + U * 2 2 U * + 1 ε 2 f 2 U * 1 f g + f + U * 1 g 2 U * , b ^ 01 ( ε ) = U * 1 f g + f + U * 1 U * f g f g + U * 2 2 U * + 1 ε 2 f 2 + U * 1 f g + f + U * 1 g U * , b ^ 20 ( ε ) = g U * f g f g + U * 2 2 U * + 1 ε 2 f f g 2 U * , b ^ 11 ( ε ) = 2 f g 2 + 2 g U * 2 g U * + 1 U * f g f g + U * 2 2 U * + 1 ε 2 f g + U * 1 f 2 f g 2 + 2 g U * 2 g U * + 1 g U * f f g + U * 1 , b ^ 02 ( ε ) = 2 f g + f + 2 U * 2 U * f g f g + U * 2 2 U * + 1 ε 2 f g + U * 1 f + 2 f g + f + 2 U * 2 g U * f f g + U * 1 .
To save space, the current expression consists of a previous expression when System (10) undergoes each approximate identity transformation in the following.
Letting u 2 = u 1 , v 2 = d u 1 d t , System (10) can be written as
d u 2 d t = v 2 , d v 2 d t = c ^ 00 ( ε ) + c ^ 10 ( ε ) u 2 + c ^ 01 ( ε ) v 2 + c ^ 20 ( ε ) u 2 2 + c ^ 11 ( ε ) u 2 v 2 + c ^ 02 ( ε ) v 2 2 + O ( | ( u 2 , v 2 ) | 3 ) ,
where
c ^ 00 ( ε ) = 1 a 01 3 ( ε ) ( a ^ 00 3 ( ε ) a ^ 02 ( ε ) a ^ 11 ( ε ) a ^ 00 2 ( ε ) a ^ 01 2 ( ε ) b ^ 02 ( ε ) a ^ 00 2 ( ε ) a ^ 01 ( ε ) a ^ 02 ( ε ) b ^ 01 ( ε ) + 2 a ^ 00 2 ( ε ) a ^ 02 2 ( ε ) b ^ 00 ( ε ) + 2 a ^ 00 ( ε ) a ^ 01 2 ( ε ) a ^ 02 ( ε ) b ^ 00 ( ε ) a ^ 01 4 ( ε ) b ^ 00 ( ε ) + a ^ 00 ( ε ) a ^ 01 3 ( ε ) b ^ 01 ( ε ) ) ,
c ^ 10 ( ε ) = 1 a ^ 01 3 ( ε ) ( 2 a ^ 00 2 ( ε ) a ^ 02 ( ε ) a ^ 11 ( ε ) a ^ 00 ( ε ) a ^ 01 2 ( ε ) a ^ 11 ( ε ) 2 a ^ 00 ( ε ) a ^ 01 2 ( ε ) b ^ 02 ( ε ) 2 a ^ 00 ( ε ) a ^ 01 ( ε ) a ^ 02 ( ε ) b ^ 01 ( ε ) + 4 a ^ 00 ( ε ) a ^ 02 2 ( ε ) b ^ 00 ( ε ) + a ^ 01 3 ( ε ) a ^ 10 ( ε ) + a ^ 01 3 ( ε ) b ^ 01 ( ε ) + 2 a ^ 01 2 ( ε ) a ^ 02 ( ε ) b ^ 00 ( ε ) ) , c ^ 01 ( ε ) = 1 a ^ 01 3 ( ε ) ( a ^ 00 ( ε ) a ^ 02 ( ε ) a ^ 11 ( ε ) a ^ 01 2 ( ε ) a ^ 11 ( ε ) a ^ 01 2 ( ε ) b ^ 02 ( ε ) a ^ 01 ( ε ) a ^ 02 ( ε ) b ^ 01 ( ε ) + 2 a ^ 02 2 ( ε ) b ^ 00 ( ε ) ) , c ^ 20 ( ε ) = 1 a ^ 01 3 ( ε ) ( a ^ 00 2 ( ε ) a ^ 01 ( ε ) a ^ 11 2 ( ε ) 2 a ^ 00 2 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) a ^ 11 ( ε ) a ^ 00 ( ε ) a ^ 01 3 ( ε ) b ^ 11 ( ε ) 2 a ^ 00 ( ε ) a ^ 01 2 ( ε ) a ^ 02 ( ε ) b ^ 10 ( ε ) + 2 a ^ 00 ( ε ) a ^ 01 2 ( ε ) a ^ 10 ( ε ) b ^ 02 ( ε ) + 2 a ^ 00 ( ε ) a ^ 01 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) b ^ 01 ( ε ) + 2 a ^ 00 ( ε ) a ^ 01 ( ε ) a ^ 02 ( ε ) a ^ 11 ( ε ) b ^ 00 ( ε ) 4 a ^ 00 ( ε ) a ^ 02 2 ( ε ) a ^ 10 ( ε ) b ^ 00 ( ε ) + a ^ 01 4 ( ε ) b ^ 10 ( ε ) a ^ 01 3 ( ε ) a ^ 10 ( ε ) b ^ 01 ( ε ) + a ^ 01 3 ( ε ) a ^ 11 ( ε ) b ^ 00 ( ε ) 2 a ^ 01 2 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) b ^ 00 ( ε ) ) , c ^ 11 ( ε ) = 1 a ^ 01 3 ( ε ) ( a ^ 00 ( ε ) a ^ 01 ( ε ) a ^ 11 2 ( ε ) 2 a ^ 00 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) a ^ 11 ( ε ) 2 a ^ 01 3 ( ε ) a ^ 20 ( ε ) a ^ 01 3 ( ε ) b ^ 11 ( ε ) 2 a ^ 01 2 ( ε ) a ^ 02 ( ε ) b ^ 10 ( ε ) + a ^ 01 2 ( ε ) a ^ 10 ( ε ) a ^ 11 ( ε ) + 2 a ^ 01 2 ( ε ) a ^ 10 ( ε ) b ^ 02 ( ε ) 4 a ^ 02 2 ( ε ) a ^ 10 ( ε ) b ^ 00 ( ε ) + 2 a ^ 01 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) b ^ 01 ( ε ) + 2 a ^ 01 ( ε ) a ^ 02 ( ε ) a ^ 11 ( ε ) b ^ 00 ( ε ) ) , c ^ 02 ( ε ) = 1 a ^ 01 3 ( ε ) ( a ^ 00 ( ε ) a ^ 01 2 ( ε ) a ^ 11 ( ε ) a ^ 20 ( ε ) a ^ 00 ( ε ) a ^ 01 ( ε ) a ^ 10 ( ε ) a ^ 11 2 ( ε ) + a ^ 00 ( ε ) a ^ 02 ( ε ) a ^ 10 2 ( ε ) a ^ 11 ( ε ) a ^ 01 4 ( ε ) b ^ 20 ( ε ) + a ^ 01 3 ( ε ) a ^ 10 ( ε ) b ^ 11 ( ε ) a ^ 01 3 ( ε ) a ^ 11 ( ε ) b ^ 10 ( ε ) + a ^ 01 3 ( ε ) a ^ 20 ( ε ) b ^ 01 ( ε ) + 2 a ^ 01 2 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) b ^ 10 ( ε ) + 2 a ^ 01 2 ( ε ) a ^ 02 ( ε ) a ^ 20 ( ε ) b ^ 00 ( ε ) a ^ 01 2 ( ε ) a ^ 10 2 ( ε ) b ^ 02 ( ε ) a ^ 01 ( ε ) a ^ 02 ( ε ) a ^ 10 2 ( ε ) b ^ 01 ( ε ) 2 a ^ 01 ( ε ) a ^ 02 ( ε ) a ^ 10 ( ε ) a ^ 11 ( ε ) b ^ 00 ( ε ) + 2 a ^ 02 2 ( ε ) a ^ 10 2 ( ε ) b ^ 00 ( ε ) ) .
To remove v 2 2 from System (11), introduce a new time variable τ with d t = 1 c ^ 02 ( ε ) u 2 d τ and make u 3 = u 2 , v 3 = v 2 1 c ^ 02 ( ε ) u 2 , then System (11) becomes (rewriting τ as t):
d u 3 d t = v 3 , d v 3 d t = d ^ 00 ( ε ) + d ^ 10 ( ε ) u 3 + d ^ 01 ( ε ) v 3 + d ^ 20 ( ε ) u 3 2 + d ^ 11 ( ε ) u 3 v 3 + O ( | ( u 3 , v 3 ) | 3 ) ,
where
d ^ 00 ( ε ) = c ^ 00 ( ε ) , d ^ 10 ( ε ) = c ^ 10 ( ε ) 2 c ^ 00 ( ε ) c ^ 02 ( ε ) , d ^ 01 ( ε ) = c ^ 01 ( ε ) , d ^ 20 ( ε ) = c ^ 20 ( ε ) 2 c ^ 02 ( ε ) c ^ 10 ( ε ) + c ^ 00 ( ε ) c ^ 02 2 ( ε ) , d ^ 11 ( ε ) = c ^ 01 ( ε ) c ^ 02 ( ε ) + c ^ 11 ( ε ) .
From d ^ 20 ( 0 ) = U * 1 2 f g + f + U * 1 2 g U * 2 f , when ε i (i = 1, 2) are small, we have d ^ 20 ( ε ) < 0 . Therefore, d ^ 20 ( ε ) can be reduced into 1 with the following coordinates change:
u 4 = u 3 , v 4 = v 3 d ^ 20 ( ε ) , τ = d ^ 20 ( ε ) t ,
then System (12) can be transformed into (rewriting τ as t)
d u 4 d t = v 4 , d v 4 d t = e ^ 00 ( ε ) + e ^ 10 ( ε ) u 4 + e ^ 01 ( ε ) v 4 u 4 2 + e ^ 11 ( ε ) u 4 v 4 + O ( | ( u 4 , v 4 ) | 3 ) ,
where
e ^ 00 ( ε ) = d ^ 00 ( ε ) d ^ 20 ( ε ) , e ^ 10 ( ε ) = d ^ 10 ( ε ) d ^ 20 ( ε ) , e ^ 01 ( ε ) = d ^ 01 ( ε ) d ^ 20 ( ε ) , e ^ 11 ( ε ) = d ^ 11 ( ε ) d ^ 20 ( ε ) .
Making the transformation
u 5 = u 4 e ^ 10 ( ε ) 2 , v 5 = v 4 ,
u 4 in System (13) can be removed, then System (13) becomes
d u 5 d t = v 5 , d v 5 d t = f ^ 00 ( ε ) + f ^ 01 ( ε ) v 5 u 5 2 + f ^ 11 ( ε ) u 5 v 5 + O ( | ( u 5 , v 5 ) | 3 ) ,
where
f ^ 00 ( ε ) = e ^ 00 ( ε ) + 1 4 e ^ 10 2 ( ε ) , f ^ 01 ( ε ) = e ^ 01 ( ε ) + 1 2 e ^ 10 ( ε ) e ^ 11 ( ε ) , f ^ 11 ( ε ) = e ^ 11 ( ε ) .
From f ^ 11 ( 0 ) = f g + f + U * 1 f g ( f g + U * 2 1 ) + 2 U * ( U * 1 ) 2 f g + U * 1 f 0 , when ε i (i = 1, 2) are small, we have f ^ 11 ( ε ) 0 . Employing the change of variables again,
u 6 = f ^ 11 2 ( ε ) u 5 , v 6 = f ^ 11 3 ( ε ) v 5 , τ = 1 f ^ 11 ( ε ) t ,
we obtain the universal unfolding of System (9),
d u 6 d τ = v 6 , d v 6 d τ = g ^ 00 ( ε ) + g ^ 01 ( ε ) v 6 u 6 2 + u 6 v 6 + O ( | ( u 6 , v 6 ) | 3 ) ,
where
g ^ 00 ( ε ) = f ^ 00 ( ε ) f ^ 11 4 ( ε ) , g ^ 01 ( ε ) = f ^ 01 ( ε ) f ^ 11 ( ε ) .
The transversality condition
( g ^ 00 , g ^ 01 ) ( ε 1 , ε 2 ) ε 1 = ε 2 = 0 = f g ( f g + U * 2 1 ) + 2 U * ( U * 1 ) 2 5 f g + U * 1 4 U * 1 3 f 4 U * 6 g 3 0
holds. Therefore, from the Bogdanov–Takens bifurcation theorem [19,20], System (2) undergoes a Bogdanov–Takens bifurcation of codimension 2, which includes a sequence of bifurcations of codimension 1: saddle-node bifurcation, Hopf bifurcation and homoclinic bifurcation, when the parameters ε 1 and ε 2 vary in a small neighborhood of ( 0 , 0 ) . □
Turning to Matlab, we provide some data groups to simulate the dynamic process of Bogdanov–Takens bifurcation to support the theoretical results. Let U * = 1 2 , f = 3 5 and g = 1 5 , then α 1 = 3 25 , β 1 = 13 60 , r 1 = 18 95 and System (9) becomes
U ˙ = { U ( U 1 ) ( 3 25 + ε 1 U ) 3 5 U D } { 1 ( U + D ) } , D ˙ = ( 18 95 + ε 2 ) { D ( 13 60 D ) + 1 5 U D } { 1 ( U + D ) } .
System (16) will present different phase portraits according to following six examples for disturbance of ε 1 and ε 2 .
Example 1: System (16) admits a cusp of codimension 2 if ( ε 1 , ε 2 ) = ( 0 , 0 ) ; see Figure 1.
Example 2: there exists no equilibrium in System (16) if ( ε 1 , ε 2 ) = ( 0.001 , 0.001 ) ; see Figure 2a.
Example 3: set ( ε 1 , ε 2 ) = ( 0.001 , 0.074 ) , System (16) admits an unstable focus and a saddle; see Figure 2b.
Example 4: take ( ε 1 , ε 2 ) = ( 0.001 , 0.06 ) , System (16) has a stable focus surrounded by an unstable limit cycle; see Figure 2c.
Example 5: System (16) has an unstable homoclinic loop consisted of a saddle and an unstable homoclinic orbit if ( ε 1 , ε 2 ) = ( 0.001 , 0.05 ) ; see Figure 2d.
Example 6: if ( ε 1 , ε 2 ) = ( 0.001 , 0.04 ) , there is a stable focus and a saddle in System (16); see Figure 2e.

4. Conclusions

The dynamical behavior of System (1) around the unique positive equilibrium E * was investigated in this manuscript. The case that the Jacobian matrix J E * has two zero eigenvalues is focused. We proved that equilibrium E * is a cusp of codimension 2 through linear transformation and normal form theory [20].
From Theorems 1 and 3, in order to ensure that System (1) admits a cusp of codimension 2, parameters f and g must satisfy f > 0 and g > 0 . Therefore, doubt should be regarded as restraining force and destroys comprehension in the learning process if System (1) has a cusp of codimension 2 and undergoes Bogdanov–Takens bifurcation of codimension 2.
By Theorem 4, System (1) undergoes Bogdanov–Takens bifurcation of codimension 2, which leads to potentially dramatic changes in the system, hence bifurcation of codimension 2 should be not ignored. Due to a saddle-node bifurcation, the number of the interior equilibrium is zero, one, or two as parameters α and r are disturbed.Thus, there are some values of parameters such as that understanding and doubt co-exist in the form of interior equilibrium for different initial values. Through a Hopf bifurcation, there is a limit cycle appearing in System (1). As seen in Figure 2c, an unstable limit cycle is produced by the Hopf bifurcation. When all initial values lying inside the limit cycle, the orbit will tend to the interior equilibrium at last. That is, understanding and doubt will reach a steady state. In other words, knowledge can not be absorbed fully; once the learner decreases the time to understand, doubts regarding the knowledge will increase. From homoclinic bifurcation, a homoclinic loop appears in System (1). As seen in Figure 2d, an unstable homoclinic loop is consisted of a saddle and an unstable homoclinic orbit, which means that understanding and doubt co-exist in the form of an interior equilibrium for all initial values lying inside the homoclinic loop.
It is noteworthy that when all initial values lying outside the unstable limit cycle and unstable homoclinic loop, combining with Theorem 2, the orbit will tend to the equilibrium E 1 . That is, if enough time was taken, the learner will feel confident with the knowledge at last.
Above all, the occurrences of Bogdanov–Takens bifurcation yields to more complex dynamical behavior, including bifurcations of codimension one. It is very sensitive to the choice of coefficient of difficulty in studying material and the speed of the learning process, which indicates that the selection of learning and teaching methods should be put forward higher requirements for students and teachers in education.
However, the existence condition of Bogdanov–Takens of codimension 2 is sufficient and the Bogdanov–Takens bifurcation of higher codimensions in the distribution functions [21] is worth studying in the future.

Author Contributions

Formal analysis, Z.Z.; Writing—review and editing, Y.G., Z.Z. and Y.G. contributed equally to the writing of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Fujian Province grant number 2020J01499.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are thankful to Hebai Chen for their valuable suggestions to upgrade the manuscript significantly.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Monteiro, L.H. Population dynamics in educational institutions considering the student satisfaction. Commun. Nonlinear Sci. Numer. Simul. 2016, 30, 236–242. [Google Scholar] [CrossRef]
  2. Marton, F.; Säljö, R. On qualitative differences in learning: I-Outcome and process. Br. J. Educ. Psychol. 1976, 46, 4–11. [Google Scholar] [CrossRef]
  3. Bacdayan, W. A mathematical analysis of the learning production process and a model for determining what matters in education. Econ. Educ. Rev. 1997, 16, 25–37. [Google Scholar] [CrossRef]
  4. Pritchard, D.E.; Lee, Y.J.; Bao, L. Mathematical learning models that depend on prior knowledge and instructional strategies. Phys. Rev. Spec.-Top.-Phys. Educ. Res. 2008, 4, 010109. [Google Scholar] [CrossRef]
  5. Bordogna, M.; Albano, E.V. Phase transitions in a model for social learning via the Internet. Int. J. Mod. Phys. 2001, 12, 1241–1250. [Google Scholar] [CrossRef]
  6. Hecht, J.M. Doubt: A History; Harper Collins: New York, NY, USA, 2004. [Google Scholar]
  7. Schechter, C.; Ganon-Shilon, S. Reforming schools: The collective doubting perspective. Int. J. Educ. Manag. 2015, 29, 62–72. [Google Scholar] [CrossRef]
  8. Mills, C.M. Knowing when to doubt: Developing a critical stance when learning from others. Dev. Psychol. 2013, 49, 404–418. [Google Scholar] [CrossRef] [PubMed]
  9. Srikantia, P.; Pasmore, W. Conviction and doubt in organizational learning. J. Organ. Chang. Manag. 1996, 9, 42–53. [Google Scholar] [CrossRef]
  10. Burbules, N.C. Aporias, webs, and passages: Doubt as an opportunity to learn. Curric. Inq. 2000, 30, 171–187. [Google Scholar] [CrossRef]
  11. Monteiro, L.H.; Notargiacomo, P.C.S. Learning process as an interplay between understanding and doubt: A dynamical systems approach. Commun. Nonlinear Sci. Numer. Simul. 2017, 47, 416–420. [Google Scholar] [CrossRef]
  12. Liu, L.L.; Ding, K.W.; Chen, H.B. Dynamical analysis of a Lotka-Volterra Learning-Process Model. J. Appl. Anal. Comput. 2019, 9, 1855–1871. [Google Scholar] [CrossRef] [PubMed]
  13. Gou, W.; Jin, Z.; Wang, H. Hopf bifurcation for general network-organized reaction-diffusion systems and its application in a multi-patch predator-prey system. J. Differ. Equ. 2023, 346, 64–107. [Google Scholar] [CrossRef]
  14. Ma, T.; Meng, X.; Hayat, T.; Hobiny, A. Hopf bifurcation induced by time delay and influence of Allee effect in a diffusive predator-prey system with herd behavior and prey chemotaxis. Nonlinear Dyn. 2022, 108, 4581–4598. [Google Scholar] [CrossRef]
  15. Stróżyna, E.; Żoła̧dek, H. Analytic properties of the complete formal normal form for the Bogdanov-Takens singularity. Nonlinearity 2021, 34, 3046. [Google Scholar] [CrossRef]
  16. Lu, M.; Xiang, C.; Huang, J.C. Bogdanov-Takens bifurcation in a SIRS epidemic model with a generalized nonmonotone incidence rate. Discret. Contin. Dyn. Syst.-S 2020, 13, 3125–3138. [Google Scholar] [CrossRef]
  17. Zhu, Z.L.; Chen, Y.M.; Li, Z.; Chen, F.D. Stability and bifurcation in a Leslie-Gower predator-prey model with Allee effect. Int. J. Bifurc. Chaos. 2022, 32, 2250040. [Google Scholar] [CrossRef]
  18. Xiong, Y.Q.; Han, M.A. Limit cycles appearing from a generalized heteroclinic loop with a cusp and a nilpotent saddle. J. Differ. Equ. 2021, 303, 575–607. [Google Scholar] [CrossRef]
  19. Perko, L. Differential Equations and Dynamical Systems. In Texts in Applied Mathematics, 3rd ed.; Springer: New York, NY, USA, 2001; Volume 7. [Google Scholar]
  20. Chow, S.-N.; Li, C.Z.; Wang, D. Normal Forms and Bifurcation of Planar Vector Fields; Cambridge University Press: Cambridge, CA, USA, 1994. [Google Scholar]
  21. Wang, W.; Cherstvy, A.G.; Liu, X.; Metzler, R. Anomalous diffusion and nonergodicity for heterogeneous diffusion processes with fractional Gaussian noise. Phys. Rev. E 2020, 102, 012146. [Google Scholar] [CrossRef] [PubMed]
Figure 1. When α = 3 25 , β = 13 60 , f = 3 5 , g = 1 5 , r = 18 95 , equilibrium E * is a cusp of codimension 2 in System (2).
Figure 1. When α = 3 25 , β = 13 60 , f = 3 5 , g = 1 5 , r = 18 95 , equilibrium E * is a cusp of codimension 2 in System (2).
Axioms 12 00853 g001
Figure 2. Phase portrait of (16) under disturbance of ε 1 and ε 2 . (a) There exists no equilibrium in System (16) if ( ε 1 , ε 2 ) = ( 0.001 , 0.001 ) ; (b) System (16) admits an unstable focus and a saddle if ( ε 1 , ε 2 ) = ( 0.001 , 0.074 ) ; (c) System (16) has an unstable limit cycle if ( ε 1 , ε 2 ) = ( 0.001 , 0.06 ) ; (d) System (16) has an unstable homoclinic loop if ( ε 1 , ε 2 ) = ( 0.001 , 0.05 ) ; (e) there is a stable focus and a saddle in System (16) if ( ε 1 , ε 2 ) = ( 0.001 , 0.04 ) .
Figure 2. Phase portrait of (16) under disturbance of ε 1 and ε 2 . (a) There exists no equilibrium in System (16) if ( ε 1 , ε 2 ) = ( 0.001 , 0.001 ) ; (b) System (16) admits an unstable focus and a saddle if ( ε 1 , ε 2 ) = ( 0.001 , 0.074 ) ; (c) System (16) has an unstable limit cycle if ( ε 1 , ε 2 ) = ( 0.001 , 0.06 ) ; (d) System (16) has an unstable homoclinic loop if ( ε 1 , ε 2 ) = ( 0.001 , 0.05 ) ; (e) there is a stable focus and a saddle in System (16) if ( ε 1 , ε 2 ) = ( 0.001 , 0.04 ) .
Axioms 12 00853 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhu, Z.; Guan, Y. Bogdanov–Takens Bifurcation Analysis of a Learning-Process Model. Axioms 2023, 12, 853. https://doi.org/10.3390/axioms12090853

AMA Style

Zhu Z, Guan Y. Bogdanov–Takens Bifurcation Analysis of a Learning-Process Model. Axioms. 2023; 12(9):853. https://doi.org/10.3390/axioms12090853

Chicago/Turabian Style

Zhu, Zhenliang, and Yuxian Guan. 2023. "Bogdanov–Takens Bifurcation Analysis of a Learning-Process Model" Axioms 12, no. 9: 853. https://doi.org/10.3390/axioms12090853

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop