Processing math: 100%
Research article Special Issues

New classes of graphs with edge δ graceful labeling

  • Graph labeling is a source of valuable mathematical models for an extensive range of applications in technologies (communication networks, cryptography, astronomy, data security, various coding theory problems). An edge δ graceful labeling of a graph G with p vertices and q edges, for any positive integer δ, is a bijective f from the set of edge E(G) to the set of positive integers {δ,2δ,3δ,,qδ} such that all the vertex labels f[V(G)], given by: f(u)=(uvE(G)f(uv))mod(δk), where k=max(p,q), are pairwise distinct. In this paper, we show the existence of an edge δ graceful labeling, for any positive integer δ, for the following graphs: the splitting graphs of the cycle, fan, and crown, the shadow graphs of the path, cycle, and fan graph, the middle graphs and the total graphs of the path, cycle, and crown. Finally, we display the existence of an edge δ graceful labeling, for the twig and snail graphs.

    Citation: Mohamed R. Zeen El Deen, Ghada Elmahdy. New classes of graphs with edge δ graceful labeling[J]. AIMS Mathematics, 2022, 7(3): 3554-3589. doi: 10.3934/math.2022197

    Related Papers:

    [1] E Laxmi Lydia, Chukka Santhaiah, Mohammed Altaf Ahmed, K. Vijaya Kumar, Gyanendra Prasad Joshi, Woong Cho . An equilibrium optimizer with deep recurrent neural networks enabled intrusion detection in secure cyber-physical systems. AIMS Mathematics, 2024, 9(5): 11718-11734. doi: 10.3934/math.2024574
    [2] Majdy M. Eltahir, Ghadah Aldehim, Nabil Sharaf Almalki, Mrim M. Alnfiai, Azza Elneil Osman . Reinforced concrete bridge damage detection using arithmetic optimization algorithm with deep feature fusion. AIMS Mathematics, 2023, 8(12): 29290-29306. doi: 10.3934/math.20231499
    [3] Wahida Mansouri, Amal Alshardan, Nazir Ahmad, Nuha Alruwais . Deepfake image detection and classification model using Bayesian deep learning with coronavirus herd immunity optimizer. AIMS Mathematics, 2024, 9(10): 29107-29134. doi: 10.3934/math.20241412
    [4] Hend Khalid Alkahtani, Nuha Alruwais, Asma Alshuhail, Nadhem NEMRI, Achraf Ben Miled, Ahmed Mahmud . Election-based optimization algorithm with deep learning-enabled false data injection attack detection in cyber-physical systems. AIMS Mathematics, 2024, 9(6): 15076-15096. doi: 10.3934/math.2024731
    [5] Thavavel Vaiyapuri, M. Sivakumar, Shridevi S, Velmurugan Subbiah Parvathy, Janjhyam Venkata Naga Ramesh, Khasim Syed, Sachi Nandan Mohanty . An intelligent water drop algorithm with deep learning driven vehicle detection and classification. AIMS Mathematics, 2024, 9(5): 11352-11371. doi: 10.3934/math.2024557
    [6] Amina Turki, Omar Kahouli, Saleh Albadran, Mohamed Ksantini, Ali Aloui, Mouldi Ben Amara . A sophisticated Drowsiness Detection System via Deep Transfer Learning for real time scenarios. AIMS Mathematics, 2024, 9(2): 3211-3234. doi: 10.3934/math.2024156
    [7] Maha M. Althobaiti, José Escorcia-Gutierrez . Weighted salp swarm algorithm with deep learning-powered cyber-threat detection for robust network security. AIMS Mathematics, 2024, 9(7): 17676-17695. doi: 10.3934/math.2024859
    [8] Sultanah M. Alshammari, Nofe A. Alganmi, Mohammed H. Ba-Aoum, Sami Saeed Binyamin, Abdullah AL-Malaise AL-Ghamdi, Mahmoud Ragab . Hybrid arithmetic optimization algorithm with deep learning model for secure Unmanned Aerial Vehicle networks. AIMS Mathematics, 2024, 9(3): 7131-7151. doi: 10.3934/math.2024348
    [9] Nuha Alruwais, Hayam Alamro, Majdy M. Eltahir, Ahmed S. Salama, Mohammed Assiri, Noura Abdelaziz Ahmed . Modified arithmetic optimization algorithm with Deep Learning based data analytics for depression detection. AIMS Mathematics, 2023, 8(12): 30335-30352. doi: 10.3934/math.20231549
    [10] Muhammad Suleman, Muhammad Ilyas, M. Ikram Ullah Lali, Hafiz Tayyab Rauf, Seifedine Kadry . A review of different deep learning techniques for sperm fertility prediction. AIMS Mathematics, 2023, 8(7): 16360-16416. doi: 10.3934/math.2023838
  • Graph labeling is a source of valuable mathematical models for an extensive range of applications in technologies (communication networks, cryptography, astronomy, data security, various coding theory problems). An edge δ graceful labeling of a graph G with p vertices and q edges, for any positive integer δ, is a bijective f from the set of edge E(G) to the set of positive integers {δ,2δ,3δ,,qδ} such that all the vertex labels f[V(G)], given by: f(u)=(uvE(G)f(uv))mod(δk), where k=max(p,q), are pairwise distinct. In this paper, we show the existence of an edge δ graceful labeling, for any positive integer δ, for the following graphs: the splitting graphs of the cycle, fan, and crown, the shadow graphs of the path, cycle, and fan graph, the middle graphs and the total graphs of the path, cycle, and crown. Finally, we display the existence of an edge δ graceful labeling, for the twig and snail graphs.



    The financial crisis will have a great impact on the economic order of China and even the whole world. During the past several decades, the international financial crises have continually burst out and spread to many countries or regions quickly. For example, "tequila crisis" of Latin American countries in 1994, "Russian virus" in 1998, and the financial crises of Southeast Asian in 1997, etc. [1]. The financial crises is very harmful and will lead to great disorder of economic development. Thus it is an important task for us to deal with the various financial models to reveal their inherent change law in order to control vicious economic development and serve human beings. At present, there are many valuable works on all kinds of finance models. For example, Yu et al. [2] reported the bifurcation and its control issue for a hyperchaotic finance model, Cao [3] investigated the chaos control of a hyperchaotic finance model, Liao et al. [4] revealed the impact of policy lag on the Hopf bifurcation and chaos for a macroeconomic model. For more related studies, one can see [5-11]

    Hopf bifurcation is an important dynamical phenomenon in delayed systems. In particular, Hopf bifurcation and its control in economic systems plays a vital role in maintaining economic stability and virtuous circle of development. Thus it is important for us to explore this topic in economic or financial models.

    In 2011, Chen and Ying [1] investigated the following financial crises contagions model:

    {du1(t)dt=αu1(t)u22(t),du2(t)dt=u2(t)(β+u1(t)u2(t)), (1.1)

    where u1,u2 denote the stock return rates of country I and country II, respectively, α>0 represents the increasing rate of the average stock returns of country I under the normal situation, and β>0 represents the decreasing rate of the stock returns of country II. In details, one can see [1]. By virtue of the stability theory of ordinary differential equation, Chen and Ying [1] systematically analyzed the stability of different equilibrium points of model (1.1).

    Considering that the stock return rate of country I is affected by the stock return rate of country II during the past time and the stock return rate of country II is affected by the stock return rate of country I during the past time, we think that it is more suitable for us to introduce the time delay into model (1.1), then we can establish the following delayed financial crises contagions model:

    {du1(t)dt=αu1(t)u22(tσ),du2(t)dt=u2(t)(β+u1(tσ)u2(t)), (1.2)

    where u1,u2 denote the stock return rates of country I and country II, respectively, α>0 represents the increasing rate of the average stock returns of country I under the normal situation, and β>0 represents the decreasing rate of the stock returns of country II, σ is a delay.

    From a mathematical point of view, fractional-order dynamical model is more efficient instrument to describe the real financial phenomenon in economics than integer-order ones since fractional-order dynamical model possesses the memory trait and hereditary peculiarity for all kinds of economic variables and inherent development process [12,17-23], Inspired by this idea, we modify the delayed financial crises contagions model (1.2) as the following fractional-order form:

    {duμ1(t)dtμ=αu1(t)u22(tσ),duμ2(t)dtμ=u2(t)(β+u1(tσ)u2(t)), (1.3)

    where 0<μ<1 is a constant, u1,u2 denote the stock return rates of country I and country II, respectively, α>0 represents the increasing rate of the average stock returns of country I under the normal situation, and β>0 represents the decreasing rate of the stock returns of country II, σ is a delay. The fractional-order financial crises contagions model (1.3) owns greater advantages in describing economic laws than the integer-order ones and Hopf bifurcation property can effectively depict the stock return rates of country I and country II. Motivated by this idea, we think that it is necessary to deal with the Hopf bifurcation and its control issue for model (1.3). In particular, the key object is to discuss the stability and Hopf bifurcation of system (1.3) and analyze the Hopf bifurcation control issue of system (1.3). In addition, we still reveal the effect of delay on Hopf bifurcation of system (1.3).

    The key contributions of this work are as follow: (1) A novel fractional-order delayed financial crises contagions model is built. (2) A delay-independent sufficient condition guaranteeing the stability and the creation of Hopf bifurcation for the involved fractional-order delayed financial crises contagions model is obtained. (3) A suitable delayed feedback controller is successfully designed to control the Hopf bifurcation of the involved fractional-order delayed financial crises contagions model.

    The work is arranged as follows. The requisite theory about fractional-order differential system is prepared in Section 2. The delay-independent stability and bifurcation criteria remaining the stability and the onset of Hopf bifurcation for fractional-order delayed financial crises contagions model are built in Section 3. The delay-independent stability and bifurcation criteria maintaining the stability and the onset of Hopf bifurcation for fractional-order delayed controlled financial crises contagions model are built in Section 4. The computer simulations substantiating the studied key results are performed in Section 5. The conclusion is drawn in Section 6.

    In this section, some necessary important definitions and lemmas about fractional-order dynamical system are given.

    Definition 2.1. [12] The fractional integral of order μ of the function g(η) is given by

    Iμg(η)=1Γ(μ)ηη0(ηs)μ1g(s)ds,

    where ηη0,μ>0, and Γ(s)=0ηs1eηdη denotes Gamma function.

    Definition 2.2. [12] Let g(η)C([η0,),R). Define the Caputo fractional-order derivative of order μ of g(η) as follows:

    Dμg(η)=1Γ(lμ)ηη0g(m)(s)(ηs)μm+1ds,

    where ηη0 and m denotes a positive integer which satisfies m1μ<m. Furthermore, when 0<μ<1, then

    Dμg(η)=1Γ(1μ)ηη0g(s)(ηs)μds.

    Definition 2.3. [13] For the given system:

    Dμxl(t)=hl(xl(t)),l=1,2,,k, (2.1)

    where μ(0,1],xl(t)=(x1(t),x2(t),,xk(t)),hl(t)=(h1(t),h2(t),,hk(t)). If hl(xl)=0, then (x1,x2,,xk) is said to be the equilibrium point of system (2.1).

    Lemma 2.1. [14] For the given fractional order system Dμy=Ly,y(0)=y0 where 0<μ<1,yRk,LRk×k. Assume that λh(h=1,2,,k) is the root of the characteristic equation of Dμy=Ly. Then system Dμy=Ly is said to be asymptotically stable |arg(λh)|>μπ2(h=1,2,,k). Besides, this system is said to be stable |arg(λh)|>μπ2(h=1,2,,k) and every critical eigenvalue that satisfies |arg(λh)|=μπ2(h=1,2,,k) has geometric multiplicity one.

    Lemma 2.2. [15] For the given fractional order system Dμw(t)=T1w(t)+T2w(tσ), where w(t)=ϕ(t),t[σ,0],μ(0,1],wRn,T1,T2Rn×n,μR+(n×n). The characteristic equation of the system can be expressed as det|sμIT1T2esσ|=0. Then the zero solution of the system is asymptotically stable if each root of the equation det|sμIT1T2esσ|=0 owns negative real part.

    In this section, we are to analyze the influence of time delay σ on Hopf bifurcation for the fractional-order delayed financial crises contagions model (1.3).

    Let (u1,u2) be the equilibrium point of model (1.3), then

    {αu1u22=0,u2(β+u1u2)=0 (3.1)

    It follows from (3.1) that system (1.3) has the unique positive equilibrium point U(u1,u2) where u1=β2α,u2=αβ.

    let

    {˜u1(t)=u1(t)u1,˜u2(t)=u2(t)u2, (3.2)

    then system (1.3) is expressed as the following form:

    {d˜uμ1(t)dtμ=α(˜u1(t)+u1)(˜u2(tσ)+u2)2,d˜uμ2(t)dtμ=(˜u2(t)+u2)[β+(˜u1(tσ)+u1)(˜u2(t)+u2)]. (3.3)

    The linear system of (3.3) near (0,0) owns the expression:

    {d˜uμ1(t)dtμ=u22˜u1(t)2u1u2˜u2(tσ),d˜uμ2(t)dtμ=u22˜u1(tσ)+(2u1u2β)˜u2(t). (3.4)

    Let ui denote ˜ui(i=1,2), then system (3.4) becomes

    {duμ1(t)dtμ=u22u1(t)2u1u2u2(tσ),duμ2(t)dtμ=u22u1(tσ)+(2u1u2β)u2(t). (3.5)

    The characteristic equation of Eq (3.5) is given by

    det[sμ+u222u1u2esσu22esσsμ(2u1u2β)]=0, (3.6)

    which leads to

    s2μ+a1sμ+a2+b1e2sσ=0, (3.7)

    where

    {a1=u222u1u2+β,a2=u2(β2u1u2),b1=2u1u32. (3.8)

    Assume that

    (K1)a1>0,a2+b1>0

    holds.

    Lemma 3.1. For system (1.3), the positive equilibrium point U(u1,u2) is locally asymptotically stable provided that (K1) holds true.

    Proof. If σ=0, then (3.7) becomes

    λ2+a1λ+a2+b1=0. (3.9)

    It follows from (K1) that every root λh of (3.7) satisfies |arg(λh)|>μπ2(h=1,2). By Lemma 3.1, one knows that the positive equilibrium point U(u1,u2) is locally asymptotically stable. The proof completes.

    Assume that s=iγ=γ(cosπ2+isinπ2) is a root of Eq. (3.7), then

    {b1cos2γσ=γ2μcosγπa1γμcosγπ2a2,b1sin2γσ=γ2μsinγπa1γμsinγπ2. (3.10)

    According to (3.10), we have

    {cos2γσ=1b1[γ2μcosγπa1γμcosγπ2a2],sin2γσ=1b1[γ2μsinγπa1γμsinγπ2]. (3.11)

    and

    b21=[γ2μcosγπ+a1γμcosγπ2+a2]2+[γ2μsinγπ+a1γμsinγπ2]2, (3.12)

    which leads to

    γ4μ+ϵ1γ3μ+ϵ2γ2μ+ϵ3γμ+ϵ4=0 (3.13)

    where

    {ϵ1=2a1(cosγπcosγπ2+sinγπsinγπ2),ϵ2=2a2cosγπ,ϵ3=2a1a2cosγπ2,ϵ4=a22b21. (3.14)

    Denote

    Ψ(γ)=γ4μ+ϵ1γ3μ+ϵ2γ2μ+ϵ3γμ+ϵ4. (3.15)

    Suppose that

    (K2)a22<b21.

    By (K2), one derives ϵ4<0. Notice that dΨ(γ)dγ>0, for each γ>0, then Eq (3.13) has at least one positive real root. So, Eq (3.7) has at least a pair of purely roots.

    Suppose that Eq (3.15) has four real roots (say γh>0(h=1,2,3,4). By (3.11), we have

    σlh=12γh[arccos(γ2μcosγπa1γμcosγπ2a2b1)+2lπ], (3.16)

    where l=0,1,2,,h=1,2,3,4. Let

    γ0=minh=1,2,3,4{γ0h},γ0=γ|σ=σ0. (3.17)

    Assume that

    (K3)Q1S1+Q2S2>0,

    where

    {Q1=2μγ2μ10cos(2μ1)π2+μa1γμ10cos(μ1)π2,Q2=2μγ2μ10sin(2μ1)π2+μa1γμ10sin(μ1)π2,S1=2bγ0sin2γ0σ0,S2=2bγ0cos2γ0σ0. (3.18)

    Lemma 3.2. Suppose that s(σ)=ρ1(σ)+iρ2(σ) is the root of (3.7) near σ=σ0 such that ρ1(σ0)=0,ρ2(σ0)=γ0, then Re[dsdσ]σ=σ0,γ=γ0>0.

    Proof. It follows from (3.7) that

    (2μs2μ1+μa1sμ1)dsdσ2b1e2sσ(dsdσσ+s)=0. (3.19)

    It follows from (3.19) that

    (dsdσ)1=2μs2μ1+μa1sμ12sb1e2sσσs, (3.20)

    which leads to

    Re[(dsdσ)1]=Re[2μs2μ1+μa1sμ12sb1e2sσ]. (3.21)

    Thus

    Re[(dsdσ)1]σ=σ0,γ=γ0=Q1S1+Q2S2S21+S22. (3.22)

    Applying (K3), we have

    Re[(dsdσ)1]σ=σ0,γ=γ0>0, (3.23)

    which ends the proof.

    By means of the analysis above, the following assertion holds.

    Theorem 3.1. If (K1)(K3) are satisfied, then the positive equilibrium point (u1,u2) of system (1.3) is locally asymptotically stable if 0σ<σ0 and system (1.3) generates Hopf bifurcation around the positive equilibrium point (u1,u2) when σ passes through the delay value σ0.

    In this section, we are to analyze the influence of time delay σ on Hopf bifurcation for the fractional-order delayed controlled financial crises contagions model. we design a time delay feedback controller [16] which takes the form:

    ξ(t)=θ[u1(tσ)u1(t)], (4.1)

    where θ is feedback gain coefficient.

    {duμ1(t)dtμ=αu1(t)u22(tσ)+θ[u1(tσ)u1(t)],duμ2(t)dtμ=u2(t)(β+u1(tσ)u2(t)), (4.2)

    Clearly, system (4.2) has the unique positive equilibrium point U(u1,u2) where u1=β2α,u2=αβ.

    let

    {ˉu1(t)=u1(t)u1,ˉu2(t)=u2(t)u2, (4.3)

    then system (4.2) is expressed as the following form:

    {dˉuμ1(t)dtμ=α(ˉu1(t)+u1)(ˉu2(tσ)+u2)2+θ[ˉu1(tσ)ˉu1(t)],dˉuμ2(t)dtμ=(ˉu2(t)+u2)[β+(ˉu1(tσ)+u1)(ˉu2(t)+u2)]. (4.4)

    The linear system of (4.4) near (0,0) owns the expression:

    {dˉuμ1(t)dtμ=(u22+θ)ˉu1(t)+θˉu1(tσ)2u1u2ˉu2(tσ),dˉuμ2(t)dtμ=u22ˉu1(tσ)+(2u1u2β)ˉu2(t). (4.5)

    Let ui denote ˉui(i=1,2), then system (4.5) becomes

    {duμ1(t)dtμ=(u22+θ)u1(t)+θu1(tσ)2u1u2u2(tσ),duμ2(t)dtμ=u22u1(tσ)+(2u1u2β)u2(t). (4.6)

    The characteristic equation of Eq (4.6) is given by

    det[sμ+(u22+θ)θesσ2u1u2esσu22esσsμ(2u1u2β)]=0, (4.7)

    which leads to

    s2μ+c1sμ+c2(sμ+d1)esσ+d2e2sσ=0, (4.8)

    where

    {c1=u222u1u2+β+θ,c2=(u22θ)(β2u1u2),d1=β2u1u2,d2=2u1u32. (4.9)

    Assume that

    (K4)c1>1,c2d1+d2>0

    holds.

    Lemma 4.1. For system (4.2), the positive equilibrium point U(u1,u2) is locally asymptotically stable provided that (K4) holds true.

    Proof. If σ=0, then (4.8) becomes

    λ2+(c11)λ+c2d1+d2=0. (4.10)

    It follows from (K4) that every root λj of (4.8) satisfies |arg(λj)|>μπ2(j=1,2). By Lemma 3.1, one knows that the positive equilibrium point U(u1,u2) is locally asymptotically stable. The proof completes.

    By (4.8), we have

    (s2μ+c1sμ+c2)esσ(sμ+d1)+d2esσ=0. (4.11)

    Let s=iϱ=ϱ(cosπ2+isinπ2) be the root of Eq. (4.11), then

    {m1cosϱσm2sinϱσ=m3,n1cosϱσ+n2sinϱσ=n3, (4.12)

    where

    {m1=ϱ2μcosμπ+c1ϱμcosμπ2+c2+d2,m2=ϱ2μsinμπ+c1ϱμsinμπ2+d2,m3=ϱμcosμπ2+d1,n1=ϱ2μsinμπ+c1ϱμsinμπ2,n2=ϱ2μcosμπ+c1ϱμcosμπ2+c2d2,n3=ϱμsinμπ2. (4.13)

    It follows from (4.12) that

    {cosϱσ=m3n2+n3m2m1n2+n1m2,sinϱσ=m1n3n1m3m1n2+n1m2, (4.14)

    which leads to

    (m1n2+n1m2)2=(m3n2+n3m2)2+(m1n3n1m3)2. (4.15)

    For the convenience of calculation, we deal with this formula (4.13) properly. Let

    {ν1=cosμπ,ν2=c1cosμπ2,ν3=c2+d2,ν4=sinμπ,ν5=c1sinμπ2,ν6=d2,ν7=cosμπ2,ν8=d1,ν9=sinμπ,ν10=c1sinμπ2,ν11=cosμπ,ν12=c1cosμπ2,ν13=c2d2,ν14=sinμπ2. (4.16)

    then (4.13) becomes

    {m1=ν1ϱ2μ+ν2ϱμ+ν3,m2=ν4ϱ2μ+ν5ϱμ+ν6,m3=ν7ϱμ+ν8,n1=ν9ϱ2μ+ν10ϱμ,n2=ν11ϱ2μ+ν12ϱμ+ν13,n3=ν14ϱμ. (4.17)

    Notice that

    {(m1n2+n1m2)2=ς1ϱ8μ+ς2ϱ7μ+ς3ϱ6μ+ς4ϱ5μ+ς5ϱ4μ+ς6ϱ3μ+ς7ϱ2μ+ς8ϱμ+ς9,(m3n2+n3m2)2=ς10ϱ6μ+ς11ϱ5μ+ς12ϱ4μ+ς13ϱ3μ+ς14ϱ2μ+ς15ϱμ+ς16,m1n3n1m3)2=ς17ϱ6μ+ς18ϱ5μ+ς19ϱ4μ+ς20ϱ3μ+ς21ϱ2μ, (4.18)

    where

    {ς1=(ν1ν11+ν4ν9)2,ς2=2(ν1ν11+ν4ν9)(ν1ν12+ν2ν11+ν5ν9+ν4ν10),ς3=(ν1ν12+ν2ν11+ν5ν9+ν4ν10)2+2(ν1ν11+ν4ν9)(ν1ν13+ν2ν12+ν3ν11+ν6ν9+ν5ν10),ς4=2(ν1ν11+ν4ν9)(ν2ν13+ν3ν12+ν6ν10)+2(ν1ν12+ν2ν11+ν5ν9+ν4ν10)×(ν1ν13+ν2ν12+ν3ν11+ν6ν9+ν5ν10),ς5=(ν1ν13+ν2ν12+ν3ν11+ν6ν9+ν5ν10)2+2(ν1ν12+ν2ν11+ν5ν9+ν4ν10)×(ν2ν13+ν3ν12+ν6ν10)+2ν3ν13(ν1ν11+ν4ν9),ς6=2(ν1ν13+ν2ν12+ν3ν11+ν6ν9+ν5ν10)×(ν2ν13+ν3ν12+ν6ν10)+2ν3ν13(ν1ν12+ν2ν11+ν5ν9+ν4ν10),ς7=2ν3ν13(ν1ν13+ν2ν12+ν3ν11+ν6ν9+ν5ν10)+(ν2ν13+ν3ν12+ν6ν10)2,ς8=2ν3ν13(ν2ν13+ν3ν12+ν6ν10),ς9=(ν3ν13)2,ς10=(ν9ν11+ν4ν14)2,ς11=2(ν9ν11+ν4ν14)(ν9ν12+ν8ν11+ν5ν14),ς12=(ν9ν11+ν4ν14)2+2(ν9ν11+ν4ν14)×(ν9ν13+ν8ν12+ν6ν14),ς13=ν8ν13(ν9ν11+ν4ν14)+2(ν9ν13+ν8ν12+ν6ν14)×(ν9ν12+ν8ν11+ν5ν14),ς14=(ν9ν13+ν8ν12+ν6ν14)2+ν8ν13(ν9ν12+ν8ν11+ν5ν14),ς15=ν8ν13(ν9ν12+ν8ν11+ν5ν14),ς16=(ν8ν13)2,ς17=(ν1ν14ν7ν9)2,ς18=2(ν1ν14ν7ν9)(ν2ν14ν8ν9ν7ν10),ς19=2(ν1ν14ν7ν9)(ν3ν14ν8ν10)+(ν2ν14ν8ν9ν7ν10)2,ς20=2(ν3ν14ν8ν10)(ν2ν14ν8ν9ν7ν10),ς21=(ν3ν14ν8ν10)2. (4.19)

    By (4.15) and (4.18), we get

    υ1ϱ8μ+υ2ϱ7μ+υ3ϱ6μ+υ4ϱ5μ+υ5ϱ4μ+υ6ϱ3μ+υ7ϱ2μ+υ8ϱμ+υ9=0, (4.20)

    where

    {υ1=ς1,υ2=ς2,υ3=ς3ς10ς17,υ4=ς4ς11ς18,υ5=ς5ς12ς19,υ6=ς6ς13ς20,υ7=ς7ς14ς21,υ8=ς8ς15,υ9=ς9ς16. (4.21)

    Denote

    Φ(ϱ)=υ1ϱ8μ+υ2ϱ7μ+υ3ϱ6μ+υ4ϱ5μ+υ5ϱ4μ+υ6ϱ3μ+υ7ϱ2μ+υ8ϱμ+υ9. (4.22)

    Suppose that

    (K5)ς9<ς16.

    By (K5), one derives υ9<0. Notice that dΦ(ϱ)dϱ>0, for each ϱ>0, then Eq (4.20) has at least one positive real root. So, Eq (4.11) has at least a pair of purely roots.

    Suppose that Eq (4.11) owns eight real roots (say ϱj>0(j=1,2,,8). By (4.14), we have

    σij=1ϱj[arccos(m3n2+n3m2m1n2+n1m2)+2iπ], (4.23)

    where i=0,1,2,,j=1,2,,8. Let

    σ0=minj=1,2,,8{σ0j},ϱ0=ϱ|σ=σ0. (4.24)

    Assume that

    (K6)R1V1+R2V2>0,

    where

    {R1=[2μϱ2μ10cos(2μ1)π2+μc1ϱμ10cos(μ1)π2]cosϱ0σ0[2μϱ2μ10sin(2μ1)π2+μc1ϱμ10sin(μ1)π2]sinϱ0σ0μϱμ10cos(μ1)π2,R2=[2μϱ2μ10cos(2μ1)π2+μc1ϱμ10cos(μ1)π2]sinϱ0σ0+[2μϱ2μ10sin(2μ1)π2+μc1ϱμ10sin(μ1)π2]cosϱ0σ0μϱμ10sin(μ1)π2,V1=(ϱ2μ0cosμπ+c1ϱμ0cosμπ2+c2)ϱ0sinϱ0σ0+d2ϱ0sinϱ0σ0+(ϱ2μ0sinμπ+c1ϱμ0sinμπ2+c2)ϱ0cosϱ0σ0,V2=(ϱ2μ0cosμπ+c1ϱμ0cosμπ2+c2)ϱ0cosϱ0σ0+d2ϱ0cosϱ0σ0+(ϱ2μ0sinμπ+c1ϱμ0sinμπ2+c2)ϱ0sinϱ0σ0. (4.25)

    Lemma 4.2. Suppose that s(σ)=η1(σ)+iη2(σ) is the root of (4.11) near σ=σ0 such that η1(σ0)=0,η2(σ0)=ϱ0, then Re[dsdσ]σ=σ0,ϱ=ϱ0>0.

    Proof. It follows from (4.11) that

    (2μs2μ1+μc1sμ1)dsdσesσ+esσ(dsdσσ+s)(s2μ+c1sμ+c2)μsμ1dsdσd2esσ(dsdσσ+s)=0. (4.26)

    It follows from (4.26) that

    (dsdσ)1=RVσs, (4.27)

    where

    {R=(2μs2μ1+μc1sμ1)esσμsμ1,V=esσs(s2μ+c1sμ+c2)+d2sesσ. (4.28)

    Then

    Re[(dsdσ)1]=Re[(RV)1]. (4.29)

    Thus

    Re[(dsdσ)1]σ=σ0,ϱ=ϱ0=R1V1+R2V2V21+V22. (4.30)

    Applying (K6), we have

    Re[(dsdσ)1]σ=σ0,ϱ=ϱ0>0, (4.31)

    which ends the proof.

    By means of the analysis above, the following assertion holds.

    Theorem 4.1. If (K4)-(K6) are satisfied, then the positive equilibrium point (u1,u2) of system (4.2) is locally asymptotically stable if 0σ<σ0 and system (4.2) generates Hopf bifurcation around the positive equilibrium point (u1,u2) when σ passes through the delay value σ0.

    Example 5.1 Consider the following fractional-order financial crises contagions model:

    {du0.671(t)dt0.67=2u1(t)u22(tσ),du0.672(t)dt0.67=u2(t)(1.45+u1(tσ)u2(t)). (5.1)

    Apparently, system (5.1) owns the unique positive equilibrium point (1.0513,1.3793). Making use of Matlab software, we get γ0=0.5091 and σ0=0.0332. The conditions (K1)(K3) of Theorem 3.1 are fulfilled. In order to check the stability of the positive equilibrium point (1.0513,1.3793) and the appearance of Hopf bifurcation of system (5.1), we choose two unequal delay values. Let σ=0.025<σ0=0.0332, we obtain the computer simulation results that are presented in Figure 1. According to Figure 1, one can clearly see that the positive equilibrium point (1.0513,1.3793) keeps locally asymptotically stable situation. Figure 1 contains 4 subfigures. Subfigure 1 of Figure 1 shows that the state variable u11.0513 when the time increases. Subfigure 2 of Figure 1 implies that the state variable u21.3793 when the time increases. Subfigure 3 of Figure 1 manifests the numerical relation of u1 and u2. Subfigure 4 of Figure 1 display the numerical relation of t-u1-u2. Let σ=0.045>σ0=0.0332, we get the computer simulation results which are presented in Figure 2. According to Figure 2, we can clearly see that a Hopf bifurcation arises around the positive equilibrium point (1.0513,1.3793). Figure 2 contains 4 subfigures. Subfigure 1 of Figure 2 implies that the state variable u1 will keep a periodic oscillatory level around the value 1.0513 when the time increases. Subfigures 2 of Figure 2 implies that the state variable u2 will keep a periodic oscillatory state near the value 1.3793 when the time increases. Subfigures 3 of Figure 2 manifests the numerical relation of u1 and u2. Subfigures 4 of Figure 2 displays the numerical relation of t-u1-u2. The correlation for μ, γ0 and σ0 is listed in Table 1. Also, the bifurcation plots are presented to show that the bifurcation value is approximately equal to 0.0332 (see Figures 3 and 4).

    Figure 1.  The computer simulation results of system (5.1) when σ=0.025<σ0=0.0332.
    Figure 2.  The computer simulation results of system (5.1) when σ=0.045>σ0=0.0332.
    Table 1.  The correlation of μ,γ0 and σ0 in system (5.1).
    μ γ0 σ0
    0.18 0.9713 0.2109
    0.23 0.8107 0.2655
    0.36 0.7861 0.2872
    0.48 0.6723 0.0301
    0.67 0.5091 0.0332
    0.76 0.4728 0.0457
    0.81 0.4011 0.0504
    0.89 0.3856 0.0609
    0.94 0.2781 0.00821

     | Show Table
    DownLoad: CSV
    Figure 3.  Bifurcation plot of system (5.1): σ versus u1.
    Figure 4.  Bifurcation plot of system (5.1): σ versus u2.

    Example 5.2 Consider the following fractional-order controlled financial crises contagions model with delays:

    {du0.671(t)dt0.67=2u1(t)u22(tσ)+θ[u1(tσ)u1(t)],du0.672(t)dt0.67=u2(t)(1.45+u1(tσ)u2(t)), (5.2)

    Apparently, system (5.2) owns the unique positive equilibrium point (1.0513,1.3793). Let θ=3. Making use of Matlab software, one gets ϱ0=0.9012 and σ0=0.0583. The conditions (K4)-(K6) in Theorem 4.1 are satisfied. In order to check the stability of the positive equilibrium point (1.0513,1.3793) and the onset of Hopf bifurcation of system (5.2), we choose two unequal delay values. Let σ=0.045<σ0=0.0583, we get the computer simulation results that are presented in Figure 5. According to Figure 1, one can distinctly see that the equilibrium point (1.0513,1.3793) keeps locally asymptotically stable situation. Figure 5 includes 4 subfigures. Subfigure 1 of Figure 1 shows that the state variable u11.0513 when the time increases. Subfigure 2 of Figure 5 implies that the state variable u21.3793 when the time increases. Subfigure 3 of Figure 5 manifests the numerical relation of u1 and u2. Subfigure 4 of Figure 5 display the numerical relation of t-u1-u2. Let σ=0.075>σ0=0.0583, we obtain the computer simulation results which are presented in Figure 6. According to Figure 6, we can clearly see that a Hopf bifurcation arises around the positive equilibrium point (1.0513,1.3793). Figure 2 contains 4 subfigures. Subfigure 1 of Figure 6 implies that the state variable u1 will keep a periodic oscillatory level around the value 1.0513 when the time increases. Subfigures 2 of Figure 6 implies that the state variable u2 will keep a periodic oscillatory state near the value 1.3793 when the time increases. Subfigures 3 of Figure 6 manifests the numerical relation of u1 and u2. Subfigures 4 of Figure 6 displays the numerical relation of t-u1-u2. The correlation for μ, ϱ0 and σ0 is listed in Table 2. Also, the bifurcation plots are presented to show that the bifurcation value is approximately equal to 0.0583 (see Figures 7 and 8).

    Figure 5.  The computer simulation results of system (5.2) when σ=0.045<σ0=0.0583.
    Figure 6.  The computer simulation results of system (5.2) when σ=0.075>σ0=0.0583.
    Table 2.  The correlation of μ,ϱ0 and σ0 in system (5.2).
    μ ϱ0 σ0
    0.25 1.5209 0.0278
    0.38 1.4155 0.0357
    0.43 1.2376 0.0433
    0.55 0.9904 0.0502
    0.67 0.9012 0.0583
    0.73 0.8155 0.0624
    0.82 0.7466 0.0743
    0.90 0.6123 0.0829
    0.96 0.5842 0.0925

     | Show Table
    DownLoad: CSV
    Figure 7.  Bifurcation plot of system (5.2): σ versus u1.
    Figure 8.  Bifurcation plot of system (5.2): σ versus u2.

    Remark 5.1 From the computer numerical simulation results of Example 5.1 and Example 5.2, we know that the stability region of system (5.1) is [0,σ0=0.0332) and the Hopf bifurcation value is 0.0332, the stability region of system (5.2) is [0,σ0=0.0583) and the Hopf bifurcation value is 0.0583. Thus the stability region of system (5.1) is enlarged and the time of the onset Hopf bifurcation is postponed by designing a suitable delayed feedback controller.

    Fractional-order differential system has displayed underlying application prospect in the economic sphere. Based on the previous publications, we propose a new fractional-order delayed financial crises contagions model. By virtue of the stability theory and bifurcation knowledge of fractional-order differential equation, we derive a novel delay-independent stability and bifurcation condition to remain the stability and generate Hopf bifurcation for the involved fractional-order delayed financial crises contagions model. Through designing an appropriate delayed feedback controller, we can successfully control the stability region and the time of onset of Hopf bifurcation for the involved fractional-order delayed financial crises model. The investigated fruits are helpful for us to grasp the inherent law of economic operation and then serve mankind effectively. Also, the research approach can be applied to control the dynamical peculiarity of numerous other dynamical models in lots of disciplines.

    This work is supported by National Natural Science Foundation of China (No.61673008, No.62062018), Project of High-level Innovative Talents of Guizhou Province ([2016]5651), Key Project of Hunan Education Department (17A181), University Science and Technology Top Talents Project of Guizhou Province (KY[2018]047), Foundation of Science and Technology of Guizhou Province ([2019]1051), Guizhou University of Finance and Economics(2018XZD01). The authors would like to thank the referees and the editor for helpful suggestions incorporated into this paper.

    The authors declare that they have no conflict of interest.



    [1] G. S. Bloom, S. W. Glomb, Application of numbered undirected graphs, IEEE, 65 (1977), 562–570. doi: 10.1109/PROC.1977.10517. doi: 10.1109/PROC.1977.10517
    [2] J. Gross, J. Yellen, Graph theory and its applications, 3 Eds., London, UK: Chapman and Hall/CRC Press, 2018.
    [3] G. U. Maheswari, G. M. Jebarani, V. Balaji, Coding techniques through Fibonacci webs, difference cordial labeling and GMJ code method, J. Phys. Conf. Ser., 1139 (2018), 012077. doi: 10.1088/1742-6596/1139/1/012077. doi: 10.1088/1742-6596/1139/1/012077
    [4] G. Prasad, G. U. Maheswari, Matrix coding technique on sunflower graphs with edge product cordial labeling, IOP Conf. Ser. Mater. Sci. Eng., 872 (2020), 012004. doi: 10.1088/1757-899X/872/1/012004. doi: 10.1088/1757-899X/872/1/012004
    [5] J. A. Bondy, U. S. Murty, Graph theory, (Graduate texts in mathematics 244), Springer, New York, 2008.
    [6] B. D. Acharya, S. Arumugam, A. Rosa, Labeling of discrete structures and applications, New Delhi, India: Narosa Publishing House, 2008.
    [7] S. P. Lo, On edge-graceful labeling of graphs, Congr. Number., 50 (1985), 231–241.
    [8] M. A. Seoud, M. A. Salim, Further results on edge-odd graceful graphs, Turkish J. Math., 40 (2016), 647–656.
    [9] S. N. Daoud, Edge odd graceful labeling of some path and cycle-related graphs, AKCE Int. J. Graphs Comb., 14 (2017), 178–203. doi: 10.1016/j.akcej.2017.03.001. doi: 10.1016/j.akcej.2017.03.001
    [10] M. R. Zeen El Deen, Strong k-edge odd graceful labeling of graphs, Int. Math. Forum, 13 (2018), 393–405. doi: 10.12988/imf.2018.8634. doi: 10.12988/imf.2018.8634
    [11] M. R. Zeen El Deen, Edge-even graceful labeling of some graphs, J. Egypt. Math. Soc., 27 (2019), 20. doi: 10.1186/s42787-019-0025-x. doi: 10.1186/s42787-019-0025-x
    [12] M. R. Zeen El Deen, N. A. Omar, Further results on edge even graceful labeling of the join of two graphs, J. Egypt. Math. Soc., 28 (2020), 21. doi: 10.1186/s42787-020-00077-5. doi: 10.1186/s42787-020-00077-5
    [13] M. R. Zeen El Deen, N. A. Omar, Extending of edge even graceful labeling of graphs to strong r-edge even graceful labeling, J. Math., 2021 (2021), 1–19. doi: 10.1155/2021/6643173. doi: 10.1155/2021/6643173
    [14] J. A. Gallian, A dynamic survey of graph labeling, Electron. J. Comb., 2018.
    [15] M. R. Zeen El Deen, Edge δ graceful labeling for some cyclic-related graphs, Adv. Math. Phys., 2020 (2020), 6273245. doi: 10.1155/2020/6273245. doi: 10.1155/2020/6273245
  • This article has been cited by:

    1. Zitao Chen, Quanbin Deng, Kairui Chen, Adaptive Dynamic Programming and Its Application to Economic Dispatch in Microgrid: A Brief Overview, 2022, 9, 2409-5761, 13, 10.15377/2409-5761.2022.09.2
    2. Fouad Mohammad Salama, Nur Nadiah Abd Hamid, Umair Ali, Norhashidah Hj. Mohd Ali, Fast hybrid explicit group methods for solving 2D fractional advection-diffusion equation, 2022, 7, 2473-6988, 15854, 10.3934/math.2022868
    3. Wenni Lei, Zhe Li, Dongzhou Mei, Ricky Chee Jiun Chia, Financial crisis, labor market frictions, and economic volatility, 2023, 18, 1932-6203, e0291106, 10.1371/journal.pone.0291106
    4. Saleh Fahad Aljurbua, Extended existence results of solutions for FDEs of order 1<γ2, 2024, 9, 2473-6988, 13077, 10.3934/math.2024638
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3355) PDF downloads(194) Cited by(3)

Figures and Tables

Figures(28)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog