Processing math: 100%
Research article

Existence and global exponential stability of compact almost automorphic solutions for Clifford-valued high-order Hopfield neutral neural networks with D operator

  • Received: 13 November 2021 Revised: 31 December 2021 Accepted: 12 January 2022 Published: 18 January 2022
  • MSC : 34K14, 34K20, 34K40, 92B20

  • In this paper, a class of Clifford-valued higher-order Hopfield neural networks with D operator is studied by non-decomposition method. Except for time delays, all parameters, activation functions and external inputs of this class of neural networks are Clifford-valued functions. Based on Banach fixed point theorem and differential inequality technique, we obtain the existence, uniqueness and global exponential stability of compact almost automorphic solutions for this class of neural networks. Our results of this paper are new. In addition, two examples and their numerical simulations are given to illustrate our results.

    Citation: Yuwei Cao, Bing Li. Existence and global exponential stability of compact almost automorphic solutions for Clifford-valued high-order Hopfield neutral neural networks with D operator[J]. AIMS Mathematics, 2022, 7(4): 6182-6203. doi: 10.3934/math.2022344

    Related Papers:

    [1] Mingzhan Huang, Shouzong Liu, Ying Zhang . Mathematical modeling and analysis of biological control strategy of aphid population. AIMS Mathematics, 2022, 7(4): 6876-6897. doi: 10.3934/math.2022382
    [2] Mingzhan Huang, Xiaohuan Yu . Dynamic analysis of a mosquito population model with a stage structure and periodic releases of sterile males. AIMS Mathematics, 2023, 8(8): 18546-18565. doi: 10.3934/math.2023943
    [3] Zhongcai Zhu, Xue He . Rich and complex dynamics of a time-switched differential equation model for wild mosquito population suppression with Ricker-type density-dependent survival probability. AIMS Mathematics, 2023, 8(12): 28670-28689. doi: 10.3934/math.20231467
    [4] Weili Kong, Yuanfu Shao . The effects of fear and delay on a predator-prey model with Crowley-Martin functional response and stage structure for predator. AIMS Mathematics, 2023, 8(12): 29260-29289. doi: 10.3934/math.20231498
    [5] Eric M. Takyi, Charles Ohanian, Margaret Cathcart, Nihal Kumar . Sex-biased predation and predator intraspecific competition effects in a prey mating system. AIMS Mathematics, 2024, 9(1): 2435-2453. doi: 10.3934/math.2024120
    [6] Junjie He, Di Li, Shouzong Liu . Global dynamics of a mosquito population suppression model with stage and sex structure. AIMS Mathematics, 2023, 8(6): 14027-14046. doi: 10.3934/math.2023717
    [7] Hui Huang, Kaihong Zhao, Xiuduo Liu . On solvability of BVP for a coupled Hadamard fractional systems involving fractional derivative impulses. AIMS Mathematics, 2022, 7(10): 19221-19236. doi: 10.3934/math.20221055
    [8] Folashade Agusto, Daniel Bond, Adira Cohen, Wandi Ding, Rachel Leander, Allis Royer . Optimal impulse control of West Nile virus. AIMS Mathematics, 2022, 7(10): 19597-19628. doi: 10.3934/math.20221075
    [9] Rongjie Yu, Hengguo Yu, Min Zhao . Steady states and spatiotemporal dynamics of a diffusive predator-prey system with predator harvesting. AIMS Mathematics, 2024, 9(9): 24058-24088. doi: 10.3934/math.20241170
    [10] Saiwan Fatah, Arkan Mustafa, Shilan Amin . Predator and n-classes-of-prey model incorporating extended Holling type Ⅱ functional response for n different prey species. AIMS Mathematics, 2023, 8(3): 5779-5788. doi: 10.3934/math.2023291
  • In this paper, a class of Clifford-valued higher-order Hopfield neural networks with D operator is studied by non-decomposition method. Except for time delays, all parameters, activation functions and external inputs of this class of neural networks are Clifford-valued functions. Based on Banach fixed point theorem and differential inequality technique, we obtain the existence, uniqueness and global exponential stability of compact almost automorphic solutions for this class of neural networks. Our results of this paper are new. In addition, two examples and their numerical simulations are given to illustrate our results.



    In reality to show fractional order demeanor which can change with time and space in case of large number of physical processes. The operations of differentiation and integration of fractional order are authorized by fractional calculus. The fractional order may be taken on imaginary and real values. The theory of fuzzy sets is continuously drawing the attention of researchers towards itself due to its rich application in many fields including mechanics, electrical, engineering, processing signals, thermal system, robotics and control, signal processing and many other fields. Therefore, it has been notice that it is the centre of increasing interest of researcher during past few years.

    Differential equations are effective tools for describing a wide range of phenomena in modern-world problems. There has been meaningful progress in the study of different classes of differential equations. Traditional integer-order derivatives have recently lost popularity in recent decades in favor of fractional-order derivatives. This is because a variety of mathematical models for current issues involving fractional-order derivatives have been investigated, and their findings have been considerable. In contrast to integer-order derivatives, which are local operators, noninteger-order derivatives have the advantage of being global operators that yield precise and consistent results. Numerous classes of differential equations have been reorganized and constructed in terms of fractional-order derivatives as a result of these great benefits. We discuss some issues where fractional-order derivatives are excellent resources, such as models of global population growth, issues involving blood alcohol content, video tape models, etc. Similarly, fractional-order derivatives [3,4] are used in fields such as electrodynamics, fluid dynamics, fluid mechanics, and so on, see [1,2].

    One of the major classes of differential equations is the class of implicit differential equations. These equations are useful in the management and economic sciences. The differential equations in the equilibrium state are typical of the implicit type in economic difficulties. Related to this, we can use implicit functions to explore important aspects of most real-world graphs or surface geometry.

    Differential equations with impulsive conditions, on the other hand, are important in almost every branch of science. Dynamical systems with impulsive phenomena are used in physics, biology, economics, engineering, and other fields. See [5,6,7] for more information on how to model procedures with discontinuous jumps and disruptions using differential equations with impulsive conditions. Although impulsive differential equations have received considerable attention, it is worth noting that many aspects of these equations remain to be studied and explored. Delay differential equations can take many different forms. Pantograph differential equations, also known as proportional delay differential equations, are one type. These equations are significant because they can be used to simulate a wide range of issues in fields like population studies, physics, chemistry, economics, infectious diseases, biology, medicine, physiological and pharmaceutical kinetics, chemical kinetics, light absorption by interstellar matter, navigational control of ships and airplanes, electronic systems, electrodynamics, quantum mechanics, and more.

    The pantograph equation is the most effective form of delay differential equations; these equations have gotten a lot of attention because of the numerous applications in which they appear [8,9,10]. Recently many of the existence-uniqueness of solutions for different classes of fractional pantograph equations, see, for instance, [11,12,13]. Also, we refer to the articles [14,15] and the references therein for a rigorous inspection of the Ulam-stability of fractional pantograph differential equations [16,17]. Balachandran et al. [18], Hashemi et al. [19] and Alzabut et al. [20] worked on the fractional pantograph equations. Abbas et al. [21,22] solve ordinary differential equations. Niazi et al.[23], Iqbal et al. [24], Shafqat et al. [25,26], Alnahdi [27], Khan et al. [28], Boulares et al. [29] and Abuasbeh et al. [30,31,32] existence-uniqueness of the fuzzy fractional evolution equations were investigated. M. Houas [33] worked on the existence and Ulam stability of solution for FPDEs with two Caputo-hadamard type derivatives:

    cHDα(cHDβ+γ)u(ω))=f(ω,u(ω),u(λω)),ω[1,],γRm,0<λ<1,xu(1)=θ,u()=ϑ,θ,ϑRm,

    where 0<α,β1,cHDα,cHDβ are Caputo-Hadamard fractional derivatives and f:[1,]×Rm×RmRm, is given continuous function.

    We discuss the existence-uniqueness and Ulam-stability of solutions for the FOPDEs with Caputo fractional derivatives in light of the mentioned research:

    cDϖ(cDϱx(ω))=Aψ(ω,x(ω),x(μω))+Iϖ1f(ω,x(ω),x(λω)),[0,1],x(0)=0,Dϖ1x(1)=bυ0x(s)ds,ϱRm, (1.1)

    where 1<ϖ,ϱ2,cDϖ,cDϱ denote the Caputo fractional derivatives, ϖ1>0,A is real constant, ψ,f:[1,2]×Rm×Rm×RmRm, x:C(0,1)L1(0,1) and x=supω|x(ω)| are given continuous functions. The goal of this work is to investigate the existence, uniqueness and Ulam-Hyers-Rassias stability of results to fractional integral pantograph differential equations. Some scholars discovered FDE results in the literature, although the vast majority of them were differential equations of the first order. We discovered the results for Caputo derivatives of order (1,2) in our research. Stability, as a part of differential equation theory, is vital in both theory and application. As a result, stability is a key subject of study for researchers, and research papers on stability for FDE have been published in the last two decades, for example, essential conditions for solution stability and asymptotic stability of FDEs. We use fractional integral pantograph differential equations and apply them on various kinds of fractional derivatives and studied the existence and stability of Ulam-Hyers. Pantograph equations also play a pivotal role in pure and applied mathematics and physics. Motivated by their significance, a ton of scientists generalized these equations into different types and presented the solvability aspect of such problems both numerically and theoretically; for additional subtleties. The remaining of this paper is as below. In Section 2, we discuss and outline that are the most important features. The existence-uniqueness of the solution is demonstrated in Section 3. Section 4 investigates the Ulam-Hyers-Rassias stability (UHRS) of the solution. In addition, an example is presented in Section 5 and applications in Section 6. Finally, in Section 6, a conclusion is given.

    The operator cDϖ is the fractional derivative in the sense of Caputo [34,35], defined by

    cDϖφ(ω)=1Γ(nϖ)ω0(ωs)nϖ1φ(n)(s)ds,ϖ>0,n=[ϖ]+1,

    and the RL fractional integral [34,35] of order ϖ>0, denoted by

    Iϖφ(ω)=1Γ(α)ω0(ωs)ϖ1φ(s)ds,ω>0.

    The following lemmas are required.

    Lemma 2.1. [35,36] Assume ϖ>ϱ>0 and φL1([a,b]). Then DςIϖφ(ω=Iϖϱφ(ω),ω[,].

    Lemma 2.2. [35,36] For ϖ>0 and σ>1, we have

    Iϖ[(ωx)ϱ]=Γ(ϱ+1)Γ(ϖ+ϱ+1)(ωx)ϖ+ϱ.

    In particular, for x=0 and σ>1, we have

    Iϖ[1]=1Γ(ϖ+1)ωϖ.

    We also require the below lemmas.

    Lemma 2.3. [36] Suppose ϱ>0 and xC(0,1)L1(0,1). Then the fractional differential equation RLDϱu(ω)=0 has a unique solution

    x(ω)=nȷ=1ciωϱȷ,

    where cȷRm,ȷ=1,2,...,n,n=[ϱ]+1.

    Lemma 2.4. [36] Let ϱ>0. Then for xC(0,1)L1(0,1) and RLDϱxC(0,1)L1(0,1), we have

    Iϱ[RLDϱx(ω)]=x(ω)+nȷ=1ciωϱi,

    where cȷRm,ȷ=1,2,...,n and n1<ϱ<n.

    Lemma 2.5. [36] For ϑ>0, the general solution of the fractional differential equation cDϑu(ω)=0 is given by

    u(ω)=n1ȷ=0cȷωȷ,

    where cȷRm,ȷ=0,1,2,...,n1,n=[ϑ]+1.

    Lemma 2.6. [36] Assume ϖ>0. Then

    Iϖ[cDϖu(ω)]=u(ω)+n1ȷ=0cȷωȷ,

    for some cȷRm,ȷ=0,1,2,...,n1,n1<ϖ<n.

    Lemma 2.7. [37] Let O:EE be a completely continuous operator (that is, a compact map that is constrained to any bounded set in E). For some 0<ρ<1}, let Θ(O)={xE:x=ρO(x). The set Θ(O) is then either unbounded or has at least one fixed point.

    Assume W=C([0,1],Rm) denote the Banach space of continuous functions with the norm ||u|| from [0,1] to Rm, where ||x||=supω[0,]|x(ω)|.

    In this section, we will look at the UH and UHRS for the fractional problem 1.1. We examine the following inequalities for ω[0,1]:

    |cDϖ[cDϱu(ω)]Aψ(ω,u(ω),u(μω))Iϖ1f(ω,u(ω),u(λω))|ς, (2.1)

    and

    |cDϖ[cDϱu(ω)]Aψ(ω,u(ω),u(μω))Iϖ1f(ω,u(ω),u(λω))|ςg(ω), (2.2)

    Definition 2.1. The problem 1.1 is UHS if a real number ηφ,ϕȷ>0(ȷ=1,...,m) exists such that for each ς>0 and for each solution vW of the inequality 2.1, there exists a solution uW of the problem 1.1:

    |v(ω)u(ω)|ηφ,ϕiς,ω[0,1],ȷ=1,...,m.

    Definition 2.2. The problem 1.1 is UHRS with respect to gC([0,1],Rm+) if there exists a real number ηφ,ϕȷ>0(ȷ=1,...,m) such that for each ς>0 and for each solution vW of the inequality 2.2, there exists a solution uW of problem 1.1 with

    |v(ω)u(ω)|ηφ,ϕȷςg(ω),ω[0,1],ȷ=1,...,m.

    Remark 2.8. If and only if there exists a function F:[0,1]Rm, a function vW is a solution of the inequality 2.1 as

    |F(ω)|η,ω[0,1]

    and

    cDϖ[cDϱu(ω)]A(ω,u(ω),u(λω))Iϖ+ϱf(ω,u(ω),u(λω))=F(ω),ω[0,1].

    Lemma 3.1. Consider the fractional problem with h(ω)C([0,1],Rm),

    Dϖ[Dϱx(ω)]=h(ω),ω[0,1],1<ϖ,ϱ<2, (3.1)

    with the condition

    x(0)=0,Dϖ1x(1)=bυ0x(s)ds,ϖRm. (3.2)

    Then,

    x(ω)=1Γ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1h(s)ds+c1Γ(ϱ)Γ(ϖ+ϱ)ωϖ+ϱ1+c2. (3.3)

    Proof. Using Lemma 2.4,

    cDϱx(ω)=1Γ(ϖ)ω0(ωs)ϖ1h(s)ds+c1ωϖ1, (3.4)

    where c1Rm. Now, by Lemma 2.6, we have

    x(ω)=1Γ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1h(s)ds+c1Γ(ϱ1+1)Γ(ϖ+ϱ1+1)ωϖ+ϱ1+c2=1Γ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1h(s)ds+c1Γ(ϱ)Γ(ϖ+ϱ)ωϖ+ϱ1+c2. (3.5)

    where c2Rm. By using 3.2, we get

    c2=x(0)=0,

    and

    c1=b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0(ωs)ϖ+ϱ1h(s)ds.

    Substituting the value of c1 and c2 in 3.6 yields the solution 3.4. We define

    Qx(ω)=AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1φ(s,x(s),x(μs)ds+AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1f(s,x(s),x(λs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0φ(s,x(s),x(μs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0f(s,x(s),x(λs))ds. (3.6)

    We require the following conditions for the detail:

    (C1):φ,ϕȷ:[0,1]×Rm×RmRm,ȷ=1,...,m, are continuous functions and there exist a constants κȷ>0,(ȷ=1,...,m+1) such that for all ω[0,1] and xı,yıRm,ı=1,2.

    |φ(ω,x1,x2)φ(ω,y1,y2)|κ1(|x1y1|+|x2y2|),|ϕȷ(ω,x1,x2)φ(ω,y1,y2)|κȷ+1(|x1y1|+|x2y2|),ȷ=1,...,m.

    (C2):ψ:C([0,1],Rm)Rm is continuous function with ψ(0)=0 and there exists a constant ω>0,

    |ψ(x)ψ(y)|ω|xy|,x,yC([0,1],Rm).

    Using BFPT, we will develop a special remedy for the FPE (1).

    Theorem 3.2. Suppose that (C1) and (C2) hold. If

    14ω8k((|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]         )+(|A|Γ(ϖ+ϱ+1)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+2)bυϖ+ϱ+1   Γ(ϱ+1)]         )32+(2|A|Γ(ϖ+ϱ)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ)]         )+(2|A|Γ(ϖ+ϱ+1)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+2)bυϖ+ϱ+1   Γ(ϱ+1)]+      ω   )         ), (3.7)

    where κ=max{κȷ,ȷ=1,2,...,m+1}, then fractional problem 1.1 has a unique solution on [1,2].

    Proof. Assume that L1=max{Lȷ:ȷ=1,2,...,m+1}, where Lȷ are finite numbers given by supω[0,1]|φ(ω,0,0)| and Lȷ+1=supω[0,1]|φȷ(ω,0,0)|. Define

    r2L[|A|Γ(ϖ+ϱ+1)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+2)bυϖ+ϱ+1   Γ(ϱ+1)]         ]32(2|A|Γ(ϖ+ϱ+1)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+2)bυϖ+ϱ+1   Γ(ϱ+1)]+      ω         ),

    we show that QBrBr, where Br={uW:ur}. For uBr and by (C1) and (C2),

    |φ(ω,x(ω),x(λω)||φ(ω,x(ω),x(λω)φ(s,0,0)|+|φ(s,0,0)|2κ1x+L12κ1r+L1,|φi(ω,x(ω),x(λω)||φi(ω,x(ω),x(λω)φȷ(s,0,0)|+|φȷ(s,0,0)|2κȷ+1u+Lȷ+12κȷ+1r+Lȷ+1,ȷ=1,2,...,m, (3.8)

    and

    |ψ(x)|ωxωr. (3.9)

    Using 3.8 and 3.9,

    |Qx(ω)|AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1φ(s,x(s),x(μs))ds+AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1f(s,x(s),x(λs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0φ(s,x(s),x(μs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0f(s,x(s),x(λs))ds.

    By (C1) and (C2), we can write

    |Qx(ω)|(2kr+L)[2|A|Γ(ϖ+ϱ)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]+2|A|Γ(ϖ+ϱ+1)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ+1Γ(ϱ+1)]]=2[2|A|kΓ(ϖ+ϱ)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]]+2L[|A|Γ(ϖ+ϱ+1)+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+2)bυϖ+ϱ+1Γ(ϱ+1)]]r.

    Thus

    Q(x)r.

    This implies that QBrBr. Now we have u,vBr and all ω[0,],

    |Qx(ω)Qy(ω)|=|A|Γ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1|φ(s,x(s),x(μs)φ(s,y(s),y(μs)|ds+AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1|f(s,x(s),x(λs))f(s,y(s),y(λs))|ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0|φ(s,x(s),x(μs)φ(s,y(s),y(μs)|ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0|f(s,x(s),x(λs))f(s,y(s),y(λs))|ds.

    So, according to (C1) and (C2), we get

    QxQy(4k[|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]])xy.

    We can see from 3.7 that Q is a contractive operator. As a result the BFPT, has a fixed point that is a solution of 1.1.

    By Lemma 3.1, we demonstrate the existence of solutions to the FPE 1.1 in the below result.

    For the forthcoming result, we need to provide the below conditions:

    (C3):φ,ϕȷ:[0,1]×Rm×RmRm,ȷ=1,...,m, are continuous and there exist real constants bȷ,dȷ0 and aȷ>0,ȷ=1,...,m+1 such that for any x1,y1Rm, we have

    |φ(ω,x1,y2)|a1+b1|x1|+d1|x2|,

    (C4):ψ:C([0,1],Rm)Rm is continuous function with ψ(0)=0 and there exists constant ς>0 such that

    |ψ(x)|ςuxC([0,1],Rm).

    Theorem 3.3. Suppose that the conditions (C3) and (C4) are satisfied. If the inequality

    |A|b1+d1Γ(α+β)+b(α+β)Γ(β+1)(bȷ+1+dȷ+1   )Γ(β)[Γ(α+β+1)bυα+β   Γ(β+1)]<3ϖ2, (3.10)

    at least one solution to the 1.1 problem can be found in [1,2].

    Proof. In the first step, we demonstrate the complete continuity of the operator Q:WW. It follows that the operator Q is continuous because the functions φ,ϕȷ(ȷ=1,...,m) and ψȷ are continuous.

    Assume ΩW be bounded. Then there exist positive constants Mȷ,(ȷ=1,...,m+1) such that |φ(ω,x,y)|M1,|ϕi(ω,x,y)|Mȷ+1, for each x,yΩ and constants Π such that |ψ(x)|Π for all xC([0,1],Rm). Then we have for any xΩ, we have

    |Qx(ω)|AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1φ(s,x(s),x(μs)ds+AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1f(s,x(s),x(λs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0φ(s,x(s),x(μs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0f(s,x(s),x(λs))ds2|A|M1Γ(ϖ+ϱ+1)+2b(ϖ+ϱ)Γ(ϱ+1)Mi+1Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)],

    which assumes that

    Q(u)2(|A|M1Γ(ϖ+ϱ+1)+b(ϖ+ϱ)Γ(ϱ+1)Mȷ+1Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]).

    As a result of the preceding inequality, the operator Q is uniformly bounded.

    Following that, Q is an equicontinuous set of W. Assume, ω1,ω2(0,1] be replaced by ω1<ω2. Then,

    |Qx(ω2)Qx(ω1)||A|Γ(ϖ+ϱ)ω10[(ω2s)ϖ+ϱ1(ω1s)ϖ+ϱ1]|ψ(s,x(s),x(μs)|ds+|A|Γ(ϖ+ϱ)ω2ω1[(ω2s)ϖ+ϱ1]|ψ(s,x(s),x(μs)|ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω10[(ω2s)ϖ+ϱ1(ω1s)ϖ+ϱ1]f(s,x(s),x(λs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω2ω1[(ω2s)ϖ+ϱ1]f(s,x(s),x(λs))ds+[ωϖ+ϱ12ωϖ+ϱ11](|A|Γ(ϖ+ϱ)ω0(1s)ϖ+ϱ1|ψ(s,x(s),x(μs))|dsb(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]10(1s)ϖ+ϱ1|f(s,x(s),x(λs))|ds).

    Hence, we obtain

    |Qx(ω2)Qx(ω1)||A|M1Γ(ϖ+ϱ)ω10[(ω2s)ϖ+ϱ1(ωϖ+ϱ2ωϖ+ϱ1)]+b(ϖ+ϱ)Γ(ϱ+1)Mȷ+1Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)][(ω2s)ϖ+ϱ1+(ωϖ+ϱ2ωϖ+ϱ1)][ωϖ+ϱ12ωϖ+ϱ11](|A|M1Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Mȷ+1Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]).

    According to the above inequality, we can state that Qx(ω2)Qx(ω1)0 as ω2ω10. By applying the Arzela-Ascoli theorem, Q:WW is completely continuous.

    Finally, we demonstrate that the Θ set defined by

    Θ={xW:x=ρQ(x),1<ρ<2},

    is bounded. Suppose xΘ, then x=ρQ(x). For each ω[0,1],

    x(ω)=ρQx(ω).

    Then

    |x(ω)|=|Qx(ω)|AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1φ(s,x(s),x(μs))ds+AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1f(s,x(s),x(λs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0φ(s,x(s),x(μs))ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]υ0f(s,x(s),x(λs))ds.

    Hence

    x2|A|Γ(ϖ+ϱ)[|a1+(b1+d1)x]+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)][aȷ+1+(bȷ+1+dȷ+1)x]=(2|A|(b1+d1)Γ(ϖ+ϱ)+2b(ϖ+ϱ)Γ(ϱ+1)(bȷ+1+dȷ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)])u+2|A|Γ(ϖ+ϱ)a1+2b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]aȷ+1,

    which implies that

    x2|A|a1Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)aȷ+1Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]32(2|A|(bȷ+dȷ)Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)(bȷ+1+dȷ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]         ).

    This demonstrates that Θ is bounded. As a result, according to Lemma 2.7, operator O has at least one fixed point. As a result, there is at least one solution to the fractional pantograph problem 1.1 on [1,2]. The proof is finally completed.

    In this portion, we will look at the problem's UH and UHRS 1.1.

    Theorem 4.1. If conditions (C1) and (C2) are satisfied and also

    b(ϖ+ϱ)Γ(ϱ+1)Γ(ϖ+ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]<Γ(ϖ+ϱ+1)2k|A|. (4.1)

    Then the problem 1.1 is UHS, and thus generalized UHS.

    Proof. Presume that xW is the only solution to the issue and that yW is a solution to the inequality 2.1:

    Dϖ(Dϱx(ω))=Aψ(ω,x(ω),x(μω))+Iϖ1f(ω,x(ω),x(λω)),ω[0,1],1<ϖ,ϱ<2,x(0)=φ(0),Dϖ1x(1)=bυ0x(s)ds,ϱRm. (4.2)

    Using Eq 4.2, we can write

    x(ω)=Iϖ+ϱhx(ω)+c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]Iϖ+ϱhx(ω).

    The inequality 1.1 has now been integrated, and we have

    |y(ω)Iϖ+ϱhx(ω)+c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]Iϖ+ϱhx(ω)|ϖΓ(ϖ+ϱ)ωϖ+ϱϖΓ(ϖ+ϱ+1).

    For any ω[0,1], we have

    |y(ω)x(ω)|=|y(ω)Iϖ+ϱhx(ω)c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]Iϖ+ϱhx(ω)||y(ω)Iϖ+ϱhy(ω)+c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]Iϖ+ϱhy(ω)|,

    where

    hx(ω)=Aψ(ω,x(ω),x(μω))+Iϖ1f(ω,x(ω),x(λω)),

    and

    hy(ω)=Aψ(ω,y(ω),y(μω))+Iϖ1f(ω,y(ω),y(λω)),

    then

    Iϖ+ϱ[hy(ω)hx(ω)]=Iϖ+ϱ[Aψ(ω,y(ω),y(μω))Aψ(ω,x(ω),x(μω))]+Iϖ+ϱ[f(ω,y(ω),y(λt))f(ω,x(ω),x(λω))]=AΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1[ψ(s,y(s),y(μω))ψ(s,x(s),x(μω))]ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω0(ωs)ϖ+ϱ1[f(ω,y(ω),y(λω))f(ω,x(ω),x(λω))]ds.

    According to (C1), we get

    |Iϖ+ϱ[hy(ω)hx(ω)]|2κΓ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1|y(ω)x(ω)|ds+2kb(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω0(ωs)ϖ+ϱ1|y(ω)x(ω)|ds.

    This leads to

    |y(ω)x(ω)||x(ω)Iϖ+ϱhy(ω)c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]Iϖ+ϱhy(ω)|+2κ[|A|Γ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1|y(ω)x(ω)|ds+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω0(ωs)ϖ+ϱ1|y(ω)x(ω)|ds],

    which indicates

    |y(ω)x(ω)|ςΓ(ϖ+ϱ+1)+2k[|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]]||y(s)x(s)||.

    Also, one can see that

    ||y(s)x(s)||(12κ[|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)   Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]])ςΓ(ϖ+ϱ+1).

    Thus, we obtain the inequality

    |y(ω)x(ω)|1Γ(ϖ+ϱ)2κ[|A|+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]]ς=ε.

    Problem 1.1 is UHS, we conclude.

    Theorem 4.2. If conditions (C1),(C2) and 3.10 are satisfied. Suppose there exists Λg>0,

    Iϖ+ϱ[g(ω)]Λgg(ω), (4.3)

    for all ω[0,1], where gC([0,1],Rm+) is increasing. The problem 1.1 then becomes UHRS concerning g.

    Proof. Assume yW be a solution to the inequality 2.2, and also presume xW be the unique solution to the problem:

    Dϖ(Dϱx(ω))=Aψ(ω,x(ω),x(μω))+Iϖ1f(ω,x(ω),x(λω)),ω[0,1],1<ϖ,ϱ<2,x(0)=φ(0),Dϖ1x(1)=bυ0x(s)ds,ϱRm. (4.4)

    We obtain by integrating the inequality 2.2:

    |y(ω)Iϖ+ϱhy(ω)+c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]Iϖ+ϱhy(ω)|ςIϖ+ϱ[g(ω)]ςΛgg(ω).

    In addition, we have

    |y(ω)x(ω)||y(ω)Iϖ+ϱhy(ω)+c1b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]Iϖ+ϱhy(ω)|+2k|A|Γ(α+β)ω0(ωs)ϖ+ϱ1y(ω)x(ω)ds+2bk(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω0(ωs)ϖ+ϱ1y(ω)x(ω)ds,

    which implies that

    |y(ω)x(ω)|ςΛgg(ω)+2κ|A|Γ(ϖ+ϱ)ω0(ωs)ϖ+ϱ1y(ω)x(ω)ds+2bk(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱΓ(ϱ+1)]ω0(ωs)ϖ+ϱ1y(ω)x(ω)ds.

    Hence

    |y(ω)x(ω)|ςΛgg(ω)+2κ[|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]]y(ω)x(ω).

    This implies

    y(ω)x(ω)(12κ[|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]])ςΛgg(ω).

    Then, for each ω[0,1], there is

    |y(ω)x(ω)|Λg12k[|A|Γ(ϖ+ϱ)+b(ϖ+ϱ)Γ(ϱ+1)Γ(ϱ)[Γ(ϖ+ϱ+1)bυϖ+ϱ   Γ(ϱ+1)]        ]ςg(ω):=ϕ(ω)g(ω).

    As a result, we conclude that problem 1.1 has UHRS.

    We will use the below example to demonstrate our main points.

    Example 5.1. Assume the Caputo-type fractional pantograph equation

    cD52[CD32u(ω)]=35e3t+5u(ω)+35e3t+5u(32ω)+278,ω[0,e],u(0)=75,u(e)=7. (5.1)

    We have α52,β=32,λ=12 and =e. On the other hand,

    f(ω,u,v)=35e3t+5u(ω)+35e3ω+5u(32ω)+278.

    For ω[0,] and (u1,v1) and (u2,v2)Rm,

    |f(ω,u1,v1)f(ω,u2,v2)|35e5(|u1u2|+|v1v2|).

    Hence, the condition (H1) holds with K=217e5. Thus conditions

    4KΓ(ϖ+ϱ)ωα+β8.1035×102

    are satisfied. According to Theorem 3.2, the problem 5.1 has a unique solution on [1,e], and according to Theorem 4.1, fractional problem 5.1 is UH stable.

    Pantograph is a linkage constituting of five link connected with pin joints to form revolute pairs. It is connected in a manner based on parallelograms so that the movement of one point, in tracing an image, produces identical movements by second point. A pantograph is used to reproduce to an enlarged or a reduced scale and as exactly as possible the path described by a given point. Pantograph equation is one of the most prominent delay differential equation that play a significant role in mathematical physics. This equation exists in several branches of pure and applied mathematics including dynamical systems, control system, probability, number theory, quantum mechanics, electrodynamics. In particular, Taylor and Ockendon formulated such type of equation to describe how to receive electric current from the pantograph of an electric locomotive. Figure 1 displays the pantograph model.

    Figure 1.  Pantograph model.

    In this work, we have explained the existence-uniqueness and Ulam-type stability of solutions for FOPDEs with fractional Caputo derivative. We have established the existence-uniqueness results applying the BFPT and Leray-Schauder's alternative. Moreover, the UHS and the UHRS have been discussed. To illustrate the theoretical results we have given an example. Future work may also involve generalizing other tasks, adding observability, and developing the concept that was introduced in this mission. This is a fertile field with numerous research initiatives that have the potential to produce a wide range of theories and applications.

    The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University of funding this work through Small Groups (Project under grant number (1/350/43)).

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper



    [1] X. Z. Liu, K. L. Teo, B. J. Xu, Exponential stability of impulsive high-order Hopfield-type neural networks with time-varying delays, IEEE T. Neural Networ., 16 (2005), 1329–1339. https://doi.org/10.1109/TNN.2005.857949 doi: 10.1109/TNN.2005.857949
    [2] A. M. Alimi, C. Aouiti, F. Chérif, F. Dridi, M. S. M'hamdi, Dynamics and oscillations of generalized high-order Hopfield neural networks with mixed delays, Neurocomputing, 321 (2018), 274–295. https://doi.org/10.1016/j.neucom.2018.01.061 doi: 10.1016/j.neucom.2018.01.061
    [3] C. Aouiti, E. A. Assali, Stability analysis for a class of impulsive high-order Hopfield neural networks with leakage time-varying delays, Neural Comput. Appl., 31 (2019), 7781–7803. https://doi.org/10.1007/s00521-018-3585-z doi: 10.1007/s00521-018-3585-z
    [4] B. Li, Y. K. Li, Existence and global exponential stability of almost automorphic solution for Clifford-valued high-order Hopfield neural networks with leakage delays, Complexity, 2019 (2019), 6751806. https://doi.org/10.1155/2019/6751806 doi: 10.1155/2019/6751806
    [5] Y. K. Li, J. L. Xiang, B. Li, Almost periodic solutions of quaternion-valued neutral type high-order Hopfield neural networks with state-dependent delays and leakage delays, Appl. Intell., 50 (2020), 2067–2078. https://doi.org/10.1007/s10489-020-01634-2 doi: 10.1007/s10489-020-01634-2
    [6] Y. K. Li, X. F. Meng, Almost automorphic solutions for quaternion-valued Hopfield neural networks with mixed time-varying delays and leakage delays, J. Syst. Sci. Complex., 33 (2020), 100–121. https://doi.org/10.1007/s11424-019-8051-1 doi: 10.1007/s11424-019-8051-1
    [7] J. K. Hale, S. M. V. Lunel, Introduction to functional differential equations, New York: Springer, 2013.
    [8] Y. L. Xu, Exponential stability of pseudo almost periodic solutions for neutral type cellular neural networks with D operator, Neural Process. Lett., 46 (2017), 329–342. https://doi.org/10.1007/s11063-017-9584-8 doi: 10.1007/s11063-017-9584-8
    [9] L. G. Yao, Global convergence of CNNs with neutral type delays and D operator, Neural Comput. Appl., 29 (2018), 105–109. https://doi.org/10.1007/s00521-016-2403-8 doi: 10.1007/s00521-016-2403-8
    [10] L. G. Yao, Global exponential convergence of neutral type shunting inhibitory cellular neural networks with D operator, Neural Process. Lett., 45 (2017), 401–409. https://doi.org/10.1007/s11063-016-9529-7 doi: 10.1007/s11063-016-9529-7
    [11] X. J. Guo, C. X. Huang, J. D. Cao, Nonnegative periodicity on high-order proportional delayed cellular neural networks involving D operator, AIMS Mathematics, 6 (2021), 2228–2243. https://doi.org/10.3934/math.2021135 doi: 10.3934/math.2021135
    [12] A. P. Zhang, Pseudo almost periodic solutions for neutral type SICNNs with D operator, J. Exp. Theor. Artif. Intell., 29 (2017), 795–807. https://doi.org/10.1080/0952813X.2016.1259268 doi: 10.1080/0952813X.2016.1259268
    [13] A. P. Zhang, Almost periodic solutions for SICNNs with neutral type proportional delays and D operators, Neural Process. Lett., 47 (2018), 57–70. https://doi.org/10.1007/s11063-017-9631-5 doi: 10.1007/s11063-017-9631-5
    [14] Y. K. Li, L. Zhao, X. R. Chen, Existence of periodic solutions for neutral type cellular neural networks with delays, Appl. Math. Model., 36 (2012), 1173–1183. https://doi.org/10.1016/j.apm.2011.07.090 doi: 10.1016/j.apm.2011.07.090
    [15] H. D. Yang, Weighted pseudo almost periodicity on neutral type CNNs involving multi-proportional delays and D operator, AIMS Mathematics, 6 (2021), 1865–1879. https://doi.org/10.3934/math.2021113 doi: 10.3934/math.2021113
    [16] Y. K. Li, N. N. Huo, B. Li, On μ-pseudo almost periodic solutions for Clifford-valued neutral type neural networks with delays in the leakage term, IEEE T. Neur. Net. Lear., 32 (2021), 1365–1374. https://doi.org/10.1109/TNNLS.2020.2984655 doi: 10.1109/TNNLS.2020.2984655
    [17] S. Buchholz, G. Sommer, On Clifford neurons and Cliford multi-layer perceptrons, Neural Networks, 21 (2008), 925–935. https://doi.org/10.1016/j.neunet.2008.03.004 doi: 10.1016/j.neunet.2008.03.004
    [18] Y. Liu, P. Xu, J. Q. Lu, J. L. Liang, Global stability of Clifford valued recurrent neural networks with time delays, Nonlinear Dyn., 84 (2016), 767–777. https://doi.org/10.1007/s11071-015-2526-y doi: 10.1007/s11071-015-2526-y
    [19] J. W. Zhu, J. T. Sun, Global exponential stability of Clifford-valued recurrent neural networks, Neurocomputing, 173 (2016), 685–689. https://doi.org/10.1016/j.neucom.2015.08.016 doi: 10.1016/j.neucom.2015.08.016
    [20] Y. K. Li, J. L. Xiang, Existence and global exponential stability of anti-periodic solution for Clifford-valued inertial Cohen–Grossberg neural networks with delays, Neurocomputing, 332 (2019), 259–269. https://doi.org/10.1016/j.neucom.2018.12.064 doi: 10.1016/j.neucom.2018.12.064
    [21] N. Huo, B. Li, Y. K. Li, Anti-periodic solutions for Clifford-valued high-order Hopfield neural networks with state-dependent and leakage delays, Int. J. Appl. Math. Comput. Sci., 30 (2020), 83–98. https://doi.org/10.34768/amcs-2020-0007 doi: 10.34768/amcs-2020-0007
    [22] Y. K. Li, S. P. Shen, Compact almost automorphic function on time scales and its application, Qual. Theory Dyn. Syst., 20 (2021), 86. https://doi.org/10.1007/s12346-021-00522-5 doi: 10.1007/s12346-021-00522-5
    [23] Y. K. Li, B. Li, Pseudo compact almost automorphy of neutral type Clifford-valued neural networks with mixed delays, DCDS-B, 2021, https://doi.org/10.3934/dcdsb.2021248
    [24] Y. K. Li, S. P. Shen, Pseudo almost periodic synchronization of Clifford-valued fuzzy cellular neural networks with time-varying delays on time scales, Adv. Differ. Equ., 2020 (2020), 593. https://doi.org/10.1186/s13662-020-03041-w doi: 10.1186/s13662-020-03041-w
    [25] F. Brackx, R. Delanghe, F. Sommen, Clifford analysis, Boston: Pitman, 1982.
    [26] S. Bochner, A new approach to almost periodicity, Proc. Natl. Acad. Sci. USA, 48 (1962), 2039–2043. https://doi.org/10.1073/pnas.48.12.2039 doi: 10.1073/pnas.48.12.2039
    [27] S. Bochner, Continuous mappings of almost automorphic and almost periodic functions, Proc. Natl. Acad. Sci. USA, 52 (1964), 907–910. https://doi.org/10.1073/pnas.52.4.907 doi: 10.1073/pnas.52.4.907
    [28] T. Diagana, Almost automorphic type and almost periodic type functions in abstract spaces, New York: Springer, 2013.
    [29] B. Es-sebbar, Almost automorphic evolution equations with compact almost automorphic solutions, C. R. Math., 354 (2016), 1071–1077. https://doi.org/10.1016/j.crma.2016.10.001 doi: 10.1016/j.crma.2016.10.001
  • 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(2129) PDF downloads(91) Cited by(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog