Processing math: 100%
Research article Special Issues

Energy supplier selection using Einstein aggregation operators in an interval-valued q-rung orthopair fuzzy hypersoft structure

  • The selection of energy suppliers is important for sustainable energy management, as selecting the most appropriate suppliers reduces the environmental impact and improves resource optimization through sustainable practices. Our primary objective of this work was to develop a system for identifying energy suppliers by assessing various characteristics and their associated sub-attributes. Interval-valued q-rung orthopair fuzzy hypersoft sets (IVq-ROFHSS) originate by developing an association among interval-valued q-rung orthopair fuzzy sets and hypersoft sets. It is a crucial resource to handle unpredictable situations, mainly when presenting a component in a real-life scenario. IVq-ROFHSS is a new structure developed to manage the sub-parametric values of the alternatives. We developed the Einstein operational laws for IVq-ROFHSS and extended the Interval-valued q-rung ortho-pair fuzzy hypersoft Einstein weighted average (IVq-ROFHSEWA) and interval-valued q-rung ortho-pair fuzzy hypersoft Einstein weighted geometric (IVq-ROFHSEWG) operators. Moreover, we used the developed operators to formulate a multi-attribute group decision-making strategy to choose the ideal provider in sustainable energy management. The presented fuzzy robust approach reliably reiterated the challenged energy supplier selection in supply chain management to regular activities while alleviating overall expenses and promising stable reliability.

    Citation: Muhammad Saqlain, Xiao Long Xin, Rana Muhammad Zulqarnain, Imran Siddique, Sameh Askar, Ahmad M. Alshamrani. Energy supplier selection using Einstein aggregation operators in an interval-valued q-rung orthopair fuzzy hypersoft structure[J]. AIMS Mathematics, 2024, 9(11): 31317-31365. doi: 10.3934/math.20241510

    Related Papers:

    [1] A. Q. Khan, Ibraheem M. Alsulami . Discrete Leslie's model with bifurcations and control. AIMS Mathematics, 2023, 8(10): 22483-22506. doi: 10.3934/math.20231146
    [2] Xiongxiong Du, Xiaoling Han, Ceyu Lei . Dynamics of a nonlinear discrete predator-prey system with fear effect. AIMS Mathematics, 2023, 8(10): 23953-23973. doi: 10.3934/math.20231221
    [3] Fatao Wang, Ruizhi Yang, Yining Xie, Jing Zhao . Hopf bifurcation in a delayed reaction diffusion predator-prey model with weak Allee effect on prey and fear effect on predator. AIMS Mathematics, 2023, 8(8): 17719-17743. doi: 10.3934/math.2023905
    [4] Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Neimark-Sacker bifurcation, chaos, and local stability of a discrete Hepatitis C virus model. AIMS Mathematics, 2024, 9(11): 31985-32013. doi: 10.3934/math.20241537
    [5] Xiaoming Su, Jiahui Wang, Adiya Bao . Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491. doi: 10.3934/math.2024656
    [6] Ahmad Suleman, Rizwan Ahmed, Fehaid Salem Alshammari, Nehad Ali Shah . Dynamic complexity of a slow-fast predator-prey model with herd behavior. AIMS Mathematics, 2023, 8(10): 24446-24472. doi: 10.3934/math.20231247
    [7] Weili Kong, Yuanfu Shao . Bifurcations of a Leslie-Gower predator-prey model with fear, strong Allee effect and hunting cooperation. AIMS Mathematics, 2024, 9(11): 31607-31635. doi: 10.3934/math.20241520
    [8] Ali Al Khabyah, Rizwan Ahmed, Muhammad Saeed Akram, Shehraz Akhtar . Stability, bifurcation, and chaos control in a discrete predator-prey model with strong Allee effect. AIMS Mathematics, 2023, 8(4): 8060-8081. doi: 10.3934/math.2023408
    [9] Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Discrete Hepatitis C virus model with local dynamics, chaos and bifurcations. AIMS Mathematics, 2024, 9(10): 28643-28670. doi: 10.3934/math.20241390
    [10] Fethi Souna, Salih Djilali, Sultan Alyobi, Anwar Zeb, Nadia Gul, Suliman Alsaeed, Kottakkaran Sooppy Nisar . Spatiotemporal dynamics of a diffusive predator-prey system incorporating social behavior. AIMS Mathematics, 2023, 8(7): 15723-15748. doi: 10.3934/math.2023803
  • The selection of energy suppliers is important for sustainable energy management, as selecting the most appropriate suppliers reduces the environmental impact and improves resource optimization through sustainable practices. Our primary objective of this work was to develop a system for identifying energy suppliers by assessing various characteristics and their associated sub-attributes. Interval-valued q-rung orthopair fuzzy hypersoft sets (IVq-ROFHSS) originate by developing an association among interval-valued q-rung orthopair fuzzy sets and hypersoft sets. It is a crucial resource to handle unpredictable situations, mainly when presenting a component in a real-life scenario. IVq-ROFHSS is a new structure developed to manage the sub-parametric values of the alternatives. We developed the Einstein operational laws for IVq-ROFHSS and extended the Interval-valued q-rung ortho-pair fuzzy hypersoft Einstein weighted average (IVq-ROFHSEWA) and interval-valued q-rung ortho-pair fuzzy hypersoft Einstein weighted geometric (IVq-ROFHSEWG) operators. Moreover, we used the developed operators to formulate a multi-attribute group decision-making strategy to choose the ideal provider in sustainable energy management. The presented fuzzy robust approach reliably reiterated the challenged energy supplier selection in supply chain management to regular activities while alleviating overall expenses and promising stable reliability.



    The well-known Leslie-Gower model with a Holling Ⅱ functional response is expressed as follows:

    {dxdt=rx(1xK)bxya+x,dydt=sy(1yx), (1.1)

    where the biological interpretation of each parameter is provided in Table 1. It is noteworthy that the Leslie-Gower model incorporates a response function indicating the predator's carrying capacity, which is proportional to the prey population, a feature absent in the Lotka-Volterra model. This distinction has been extensively investigated in many references, such as [1,2,3,4,5].

    Table 1.  Biological meaning of parameters in model (1.1).
    Parameter Interpretation
    x prey population density
    y predator population density
    r intrinsic growth rate of the prey
    K maximum prey carrying capacity of the environment
    b maximum predation rate when prey is abundant
    a prey population at which predator attack capability saturates
    s intrinsic growth rate of the predator

     | Show Table
    DownLoad: CSV

    Previous studies predominantly focused on direct predation effects. In 2016, Wang et al. [6] demonstrated through experiments that predator-induced fear (indirect effects) in prey leads to a reduction in prey birth rates. They introduced the fear factor F(k,y)=11+ky, where k represents the intensity of fear driving antipredator behavior. To further investigate the influence of fear on population dynamics, we enhance the system (1.1) by incorporating the fear factor F(k,y). The modified model is formulated as follows:

    {dxdt=rx1+ky(1xK)bxya+x,dydt=sy(1yx). (1.2)

    For additional significant findings on the impact of fear in predator-prey models, refer to [7,8,9].

    There has been increasing recognition that traditional integer-order differential equations may not sufficiently capture the complexity of biological systems [10,11,12,13,14]. In reality, the behaviors of most organisms in nature are influenced by their historical context. Fractional-order derivatives, which extend the concept of integer-order differentiation, offer a more flexible and accurate framework for modeling memory and hereditary properties in population dynamics. Therefore, the authors of this paper aim to apply the Caputo fractional derivative to (1.2), thereby extending it into a fractional model. The Caputo fractional derivative of a function u(t) of order α(0,1] is given in [15] as follows

    Dαtu(t)=1Γ(1α)t0u(τ)(tτ)αdτ,

    where Γ is the Gamma function. By replacing the integer-order derivative with the Caputo fractional derivative, the following model is obtained:

    {Dαtx=rx1+ky(1xK)bxya+x,Dαty=sy(1yx). (1.3)

    Specifically, when α=1, the fractional-order model (1.3) reduces to the integer-order model (1.2), demonstrating that the fractional-order model serves as a generalization of the integer-order model. Moreover, in fractional calculus, the rate of change at any given moment, expressed by the fractional-order derivative, depends on the population density over a specified time interval. This feature gives the fractional-order model (1.3) a distinct advantage in capturing memory effects within populations.

    When studying populations with nonoverlapping generations or small sizes, mathematical models are often expressed in discrete terms. Although some variables may change continuously in real-world scenarios, the recording of these changes typically happens at specific time intervals during data collection. Thus, employing discrete systems to examine the dynamic behavior of biological populations is highly practical and significant. In references [16,17], the authors employed the piecewise constant approximation method to discretize continuous fractional predator-prey models and explored their dynamical properties. Singh and Sharma [18] examine a discrete prey-predator model with Holling type Ⅱ functional response and prey refuge, identifying bifurcations and controlling chaos through state feedback, pole placement, and hybrid techniques. Berkal and Almatrafi [19] used the exponential piecewise constant argument to discretize continuous fractional activator-inhibitor system. Their analysis includes stability assessments, investigations of Neimark-Sacker and period-doubling bifurcations, and numerical simulations that validate the theoretical findings on the system's dynamics. For a more comprehensive exploration of discrete model studies, readers are advised to consult references [20,21,22,23,24,25] and the related literature cited therein. However, current literature lacks studies on fractional discrete-time predator-prey Leslie-Gower systems that incorporate the fear effect in prey population. Using the same method as in references [16,17] to discretize model (1.3), we can obtain the following discrete model:

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn],yn+1=yn+hαΓ(α+1)[syn(1ynxn)]. (1.4)

    Here h>0 represents the time interval of production.

    Furthermore, the key contributions and findings of this study are summarized as follows:

    ● The Caputo fractional derivative of order (0,1] is utilized to incorporate the memory effect into the dynamical behavior of the proposed model.

    ● The existence and stability of fixed points are investigated.

    ● Conditions for the occurrence and direction of period-doubling bifurcation and Neimark–Sacker bifurcation at the positive fixed point are established.

    ● State feedback and hybrid control strategies are employed to manage bifurcations and chaotic behavior in the model.

    ● To validate the accuracy of our theoretical findings, numerical examples for the fractional-order discrete-time Leslie-Gower model with a fear factor are provided.

    The remainder of this paper is structured as follows: In Section 2, we investigate the existence and stability of the equilibrium points of model (1.4). Section 3 analytically demonstrates that model (1.4), under specific parametric conditions, undergoes period-doubling or Neimark-Sacker bifurcation. Section 4 explores the control of chaos toward an unstable equilibrium point using feedback control or hybrid control approaches. Section 5 presents a quantitative analysis of the dynamics of model (1.4) to validate our analytical findings. Finally, Section 6 offers brief conclusions.

    This section explores the existence of equilibrium points in model (1.4) and assesses their stability by evaluating the eigenvalues of the Jacobian matrix at these points. The following definition and lemma are introduced to assist in this stability analysis of equilibrium points.

    Definition 1. [26] Let λ1 and λ2 denote the two roots of the characteristic equation F(λ)=λ2+pλ+q=0 associated with the Jacobian matrix J(x,y). The equilibrium point (x,y) is termed

    (1) sink if |λ1|<1 and |λ2|<1, and the sink is locally asymptotically stable;

    (2) source if |λ1|>1 and |λ2|>1, and the source is locally unstable;

    (3) saddle if |λ1|>1 and |λ2|<1 (or |λ1|<1 and |λ2|>1);

    (4) non-hyperbolic if either |λ1|=1 or |λ2|=1.

    Lemma 1. [26] Let F(1)>0 in F(λ)=λ2Mλ+N, where λ1,λ2 are the two roots of F(λ)=0. Then, the following results hold true:

    (1) |λ1|<1 and |λ2|<1 if, and only if, F(1)>0 and N<1;

    (2) |λ1|<1 and |λ2|>1 (or |λ1|>1 and |λ2|<1) if, and only if, F(1)<0;

    (3) |λ1|>1 and |λ2|>1 if, and only if, F(1)>0 and N>1;

    (4) λ1=1 and |λ2|1 if, and only if, F(1)=0 and M0,2;

    (5) λ1 and λ2 are conjugate complex and |λ1,2|=1 if, and only if, M2<4N and N=1.

    It is evident that the equilibrium points of model (1.4) satisfy the following equations:

    {x=x+hαΓ(α+1)[rx1+ky(1xK)bxya+x],y=y+hαΓ(α+1)[sy(1yx)].

    This algebraic system is satisfied if x=K and y=0, indicating that the model (1.4) has a boundary equilibrium point E0(K,0) for all model parameters. To find the interior equilibrium point E1, we will solve the following system simultaneously:

    r1+ky(1xK)bya+x=0,s(1yx)=0.

    From the second equation, we find y=x. Substituting y=x into the first equation yields:

    (rK+bk)x2+(arK+br)xar=0,

    Given that all model parameters are positive, we obtain the following results through direct calculations.

    Theorem 1. The model (1.4) always has a boundary equilibrium point E0(K,0) and a positive equilibrium point E1=(x,y), where

    x=y=(arK+br)+(arK+br)2+4ar(rK+bk)2(rK+bk).

    The Jacobian matrix of model (1.4) evaluated at any point (x,y) is as follows:

    J(x,y)=[1+hαΓ(α+1)[r1+ky(12xK)aby(a+x)2]hαΓ(α+1)[rkx(1+ky)2(1xK)+bxa+x]hαΓ(α+1)sy2x21+hαΓ(α+1)(s2syx)]. (2.1)

    Theorem 2. The boundary equilibrium point E0(K,0) exhibits the following behaviors:

    (1) It is a source point if r>2Γ(α+1)hα;

    (2) It is a saddle point if 0<r<2Γ(α+1)hα;

    (3) It is a non-hyperbolic point if r=2Γ(α+1)hα.

    Proof. The Jacobian matrix (2.1) evaluated at the boundary equilibrium point E0(K,0) is given by:

    J(E0)=[1hαrΓ(α+1)hαΓ(α+1)bKa+K01+hαsΓ(α+1)].

    The eigenvalues of J(E0) are λ1=1hαrΓ(α+1) and λ2=1+hαsΓ(α+1). Clearly, |λ2|>1 and:

    |λ1|={<1if r>2Γ(α+1)hα,>1if 0<r<2Γ(α+1)hα,=1if r=2Γ(α+1)hα.

    By applying Definition 1, the proof is complete.

    Next, we analyze the local dynamics of model (1.4) at the positive equilibrium point E1(x,y). The Jacobian matrix (2.1) evaluated at E1(x,y) is expressed as:

    J(E1)=(1+Aa11Aa12As1As), (2.2)

    where

    A=hαΓ(α+1),a11=r1+kx(12xK)abx(a+x)2,a12=rkx(1+kx)2(1xK)bxa+x.

    The characteristic equation of J(E1) is given by:

    λ2Mλ+N=0, (2.3)

    where

    M=2+(a11s)A,N=1+(a11s)A(a11+a12)sA2.

    Let F(λ)=λ2Mλ+N, then:

    F(0)=N,F(1)=1+M+N,F(1)=1M+N.

    As a result, we establish the following theorem.

    Theorem 3. Let E1 be the unique positive equilibrium point of model (1.4). Then:

    (1) E1 is a sink point if |M|<1+N<2;

    (2) E1 is a source point if |M|<|1+N| and |N|>1;

    (3) E1 is a saddle point if M2>4N and |M|>|1+N|;

    (4) E1 is non-hyperbolic if |M|=|1+N| or N=1 and |M|<2.

    Proof. (1) According to Lemma 1, E1 is a sink point if and only if F(1)>0, F(1)>0, and N<1. This condition is satisfied when |M|<1+N<2. Similarly, the proofs for Theorem 3 (2)–(4) follow straightforwardly from the definitions and properties of M and N.

    In this section, we will analyze the period-doubling and Neimark-Sacker bifurcation behaviors of model (1.4) at the positive equilibrium point E1(x,y). The reason for not analyzing the boundary equilibrium point is that, at E0(K,0), the predators have become extinct, leaving only the prey.

    Assume M2>4N and

    s=s1:=4+2Aa112A+A2(a11+a12),

    and we can ascertain that the characteristic equation (2.3) satisfies

    F(1)=1+M+N=1+[2+(a11s1)A]+[1+(a11s1)A(a11+a12)s1A2]=4+2a11A[2A+(a11+a12)A2]s1=0.

    Therefore, the eigenvalues of the characteristic equation (2.3) are

    λ1=1,λ2=1+M=3+(a11s1)A.

    Let

    s2=a11+2A,s3=a11+4A.

    The condition |λ2|1 implies that s1s2,s3.

    Based on the above analysis, we can conclude that when the parameters vary within the set

    P.D={(h,α,r,k,K,b,a,s):M2>4N,s=s1andss2,s3}, (3.1)

    model (1.4) will undergo period-doubling bifurcation at E1(x,y).

    Select parameters (h,α,r,k,K,b,a,s)P.D and consider s as a small perturbation of s, i.e., s=ss1, where |s|1. After this perturbation, model (1.4) can be represented as follows:

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn],yn+1=yn+hαΓ(α+1)[(s1+s)yn(1ynxn)]. (3.2)

    Let un=xnx and vn=yny, transforming the equilibrium point E1(x,y) into the origin O(0,0). The model (3.2) can thus be expressed as:

    {un+1=un+hαΓ(α+1)[r(un+x)1+k(vn+y)(1un+xK)b(un+x)(vn+y)a+(un+x)],vn+1=vn+hαΓ(α+1)[(s1+s)(vn+y)(1vn+yun+x)]. (3.3)

    The Taylor expansion of model (3.3) around (un,vn)=(0,0) yields the following form:

    [un+1vn+1]=[1+Aa11Aa12As11As1][unvn]+[f(un,vn,s)g(un,vn,s)], (3.4)

    where

    f(un,vn,s)=c13u2n+c14unvn+c15v2n+c16u3n+c17u2nvn+c18unv2n+c19v3n+O((|un|+|vn|+|s|)4),g(un,vn,s)=c23u2n+c24unvn+c25v2n+c26u3n+c27u2nvn+c28unv2n+c29v3n+d1uns+d2vns+d3u2ns+d4unvns+d5v2ns+O((|un|+|vn|+|s|)4),

    and

    c13=ArK(1+kx)+abAx(a+x)3,c14=Ark(1+kx)2(12xK)abA(a+x)2,c15=Ark2x(1+kx)3(1xK),c16=Aabx(a+x)4,c17=ArkK(1+kx)2+abA(a+x)3,c18=Ark2(1+kx)3(12xK),c19=Ark3x(1+kx)4(1xK),c23=c25=As1x,c24=2As1x,c26=c28=As1x2,c27=2As1x2,c29=0,d1=A,d2=A,d3=d5=Ax,d4=2Ax.

    Define

    T=[Aa12Aa122Aa11λ21Aa11],

    then

    T1=1Aa12(1+λ2)[λ21Aa11Aa122+Aa11Aa12].

    Using the transformation:

    [unvn]=T[˜un˜vn],

    model (3.4) is transformed into the following form:

    [˜un+1˜vn+1]=[100λ2][˜un˜vn]+[˜f(˜un,˜vn,s)˜g(˜un,˜vn,s)], (3.5)

    where

    ˜f(˜un,˜vn,s)=(λ21Aa11)c13Aa12c23Aa12(1+λ2)u2n+(λ21Aa11)c14Aa12c24Aa12(1+λ2)unvn+(λ21Aa11)c15Aa12c25Aa12(1+λ2)v2n+(λ21Aa11)c16Aa12c26Aa12(1+λ2)u3n+(λ21Aa11)c17Aa12c27Aa12(1+λ2)u2nvn+(λ21Aa11)c18Aa12c28Aa12(1+λ2)unv2n+(λ21Aa11)c19Aa12c29Aa12(1+λ2)v3nd11+λ2unsd21+λ2vnsd31+λ2u2nsd41+λ2unvnsd51+λ2v2ns+O((|un|+|vn|+|s|)4),
    ˜g(˜un,˜vn,s)=(2+Aa11)c13+Aa12c23Aa12(1+λ2)u2n+(2+Aa11)c14+Aa12c24Aa12(1+λ2)unvn+(2+Aa11)c15+Aa12c25Aa12(1+λ2)v2n+(2+Aa11)c16+Aa12c26Aa12(1+λ2)u3n+(2+Aa11)c17+Aa12c27Aa12(1+λ2)u2nvn+(2+Aa11)c18+Aa12c28Aa12(1+λ2)unv2n+(2+Aa11)c19+Aa12c29Aa12(1+λ2)v3n+d11+λ2uns+d21+λ2vns+d31+λ2u2ns+d41+λ2unvns+d51+λ2v2ns+O((|un|+|vn|+|s|)4),

    and

    un=Aa12˜un+Aa12˜vn,vn=(2+Aa11)˜un+(λ21Aa11)˜vn.

    Next, we apply the center manifold theorem [26] to analyze the dynamics around the equilibrium point (˜un,˜vn)=(0,0) at s=0. According to the theorem, the model (3.5) has a center manifold, which can be represented as:

    Wc(0,0,0)={(˜un,˜vn,s)R3+:˜vn=w(˜un,s),w(0,0)=0,Dw(0,0)=0}.

    Assume that

    w(˜un,s)=η1˜u2n+η2˜uns+η3(s)2+O((|˜un|+|s|)3).

    Then, the center manifold must satisfy

    w(˜un+˜f(˜un,w(˜un,s),s),s)λ2w(˜un,s)˜g(˜un,w(˜un,s),s)=0.

    By comparing the coefficients, it can be obtained that

    η1=Aa12[(2+Aa11)c13+Aa12c23](2+Aa11)[(2+Aa11)c14+Aa12c24]1λ22+(2+Aa11)2[(2+Aa11)c15+Aa12c25]Aa12(1λ22),η2=(2+Aa11)d2Aa12d1(1+λ2)2,η3=0.

    Thus, the model (3.5), when restricted to the center manifold Wc(0,0,0), is given by:

    ˜W:˜un+1=˜un+m1˜u2n+m2˜uns+m3˜u3n+m4˜u2ns+m5˜un(s)2+O((|˜un|+|s|)4) (3.6)

    where

    m1=Aa12[(λ21Aa11)c13Aa12c23]1+λ2(2+Aa11)[(λ21Aa11)c14Aa12c24]1+λ2+(2+Aa11)2[(λ21Aa11)c15Aa12c25]Aa12(1+λ2)m2=d1Aa121+λ2+d2(2+Aa11)1+λ2m3=2Aa12η1[(λ21Aa11)c13Aa12c23]1+λ2+Aa12η1(λ232Aa11)[(λ21Aa11)c14Aa12c24]1+λ22η1(2+Aa11)(λ21Aa11)[(λ21Aa11)c15Aa12c25]Aa12(1+λ2)+(Aa12)2[(λ21Aa11)c16Aa12c26]1+λ2Aa12(2+Aa11)[(λ21Aa11)c17Aa12c27]1+λ2+(2+Aa11)2[(λ21Aa11)c18Aa12c28]1+λ2(2+Aa11)3[(λ21Aa11)c19Aa12c29]Aa12(1+λ2)m4=2Aa12η2[(λ21Aa11)c13Aa12c23]1+λ2+η2(λ232Aa11)[(λ21Aa11)c14Aa12c24]1+λ22η2(2+Aa11)(λ21Aa11)[(λ21Aa11)c15Aa12c25]Aa12(1+λ2)Aa12η1d11+λ2(λ21Aa11)η1d21+λ2(Aa12)2d31+λ2+Aa12(2+Aa11)d41+λ2(2+Aa11)2d51+λ2m5=Aa12η2d11+λ2(λ21Aa11)η2d21+λ2.

    In order for Eq (3.6) to undergo period-doubling bifurcation, it is necessary that the following two quantities possess nonzero values.

    β1=(2˜W˜uns+12˜Ws2˜W˜u2n)|(0,0)=m2β2=[163˜W˜u3n+(122˜W˜u2n)2]|(0,0)=m3+m21.

    Summarize the above analysis into the following theorem.

    Theorem 4. If β1β20, then model (1.4) undergoes a period-doubling bifurcation at the positive equilibrium point E1(x,y) when parameters vary in a small neighborhood of P.D. Additionally, when β2>0 (respectively, β2<0), model (1.4) bifurcates from the equilibrium point E1(x,y) to a stable (respectively, unstable) 2-periodic orbit.

    Assume M2<4N and

    s=s4:=a111+A(a11+a12),

    and we can determine that the eigenvalues of the characteristic equation (2.3) satisfy

    |λ1,2|2=|M2±4NM22i|2=1+(a11s4)A(a11+a12)s4A2=1+a11A[1+(a11+a12)A]As4=1.

    Namely, Eq (2.3) has two complex conjugate roots with unit modulus. It is clear from the above discussion that

    d|λ1,2|ds|s=s4=A2[1+(a11+a12)A]0.

    In addition, it is crucial that when s=s4, λθ1,2(s4)1(θ=1,2,3,4), which is equivalent to M(s4)2,1,0,2. Since M2<4N, we deduce M(s4)2,2. Additionally, we necessitate that M(s4)0,1, which leads to

    s4s5:=a11+3A,s4s6:=a11+2A.

    Based on the preceding analysis, we conclude that when the parameters vary within a small neighborhood of the set

    N.S={(h,α,r,k,K,b,a,s):M2<4N,s=s4,ss5,ss6}, (3.7)

    model (1.4) undergoes a Neimark-Sacker bifurcation at E1(x,y).

    Next, assuming (h,α,r,k,K,b,a,s)N.S. and |s|1 represents a small perturbation of s4, model (1.4) can be described as follows:

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn],yn+1=yn+hαΓ(α+1)[(s4+s)yn(1ynxn)]. (3.8)

    Let un=xnx and vn=yny, then we have

    [un+1vn+1]=[1+Aa11Aa12As41As4][unvn]+[F(un,vn)G(un,vn)] (3.9)

    where

    F(un,vn)=c13u2n+c14unvn+c15v2n+c16u3n+c17u2nvn+c18unv2n+c19v3n+O((|un|+|vn|)4),G(un,vn)=c23u2n+c24unvn+c25v2n+c26u3n+c27u2nvn+c28unv2n+c29v3n+O((|un|+|vn|)4),

    and c13, c14, c15, c16, c17, c18, c19, c23, c24, c25, c26, c27, c28, c29 are given in (2.3) by substituting s1 for s4+s.

    In order to obtain the normal form of model (3.9) at s=0, we use the following transformation:

    [unvn]=[ω1+Aa11ρ0As4][ˆunˆvn],

    with

    ρ=M2,ω=4NM22.

    Using this transformation, model (3.9) will transform as follows:

    [ˆun+1ˆvn+1]=[ρωωρ][ˆunˆvn]+[ˆF(ˆun,ˆvn)ˆG(ˆun,ˆvn)], (3.10)

    where

    ˆF(ˆun,ˆvn)=As4c13+(ρ1Aa11)c23As4ωu2n+As4c14+(ρ1Aa11)c24As4ωunvn+As4c15+(ρ1Aa11)c25As4ωv2n+As4c16+(ρ1Aa11)c26As4ωu3n+As4c17+(ρ1Aa11)c27As4ωu2nvn+As4c18+(ρ1Aa11)c28As4ωunv2n+As4c19+(ρ1Aa11)c29As4ωv3n+O((|un|+|vn|)4),ˆG(ˆun,ˆvn)=1As4(c23u2n+c24unvn+c25v2n+c26u3n+c27u2nvn+c28unv2n+c29v3n)+O((|un|+|vn|)4),

    and

    un=ωˆun+(1+Aa11ρ)ˆvn,vn=As4ˆvn.

    Next, a nonzero real number is defined as follows:

    Ω=Re[(12λ1)λ221λ1ξ11ξ12]12|ξ11|2|ξ21|2+Re(λ2ξ22) (3.11)

    where

    ξ11=14[2ˆFˆu2n+2ˆFˆv2n+i(2ˆGˆu2n+2ˆGˆv2n)]|s=0ξ12=18[2ˆFˆu2n2ˆFˆv2n+22ˆGˆunˆvn+i(2ˆGˆu2n2ˆGˆv2n22ˆFˆunˆvn)]|s=0ξ21=18[2ˆFˆu2n2ˆFˆv2n22ˆGˆunˆvn+i(2ˆGˆu2n2ˆGˆv2n+22ˆFˆunˆvn)]|s=0ξ22=116[3ˆFˆu3n+3ˆFˆunˆv2n+3ˆGˆu2nˆvn+3ˆGˆv3n+i(3ˆGˆu3n+3ˆGˆunˆv2n3ˆFˆu2nˆvn3ˆFˆv3n)]|s=0.

    Through some complicated calculations, we get

    2ˆFˆu2n|s=0=2ωAs4[c13As4+c23(ρ1Aa11)],2ˆFˆunˆvn|s=0=1As4[c14(As4)2As4(ρ1Aa11)(2c13c24)3c23(ρ1Aa11)2],2ˆFˆv2n|s=0=1ωAs4[2c15(As4)22(As4)2(ρ1Aa11)(c14c25)2As4(c24c13)(ρ1Aa11)2+2c23(ρ1Aa11)3],3ˆFˆu3n|s=0=6ω2As4[c16As4+c26(ρ1Aa11)],3ˆFˆu2nˆvn|s=0=ωAs4[2c17(As4)2As4(ρ1Aa11)(6c162c27)6c26(ρ1Aa11)2],3ˆFˆunˆv2n|s=0=1As4[2c18(As4)2+2(As4)2(c282c17)(ρ1Aa11)+6As4(c16c27)(ρ1Aa11)2+6c26(ρ1Aa11)3],3ˆFˆv3n|s=0=6(1+Aa11ρ)ωAs4[c18(As4)3+(As4)2(c28c17)(ρ1Aa11)+As4(c16c27)(ρ1Aa11)2+c26(ρ1Aa11)3],

    and

    2ˆGˆu2n|s=0=2c23ω2As4,2ˆGˆunˆvn|s=0=ωAs4[c24As42c23(ρ1Aa11)],2ˆGˆv2n|s=0=2As4[c25(As4)2As4c24(ρ1Aa11)+c23(ρ1Aa11)2],3ˆGˆu3n|s=0=6c26ω3As4,3ˆGˆu2nˆvn|s=0=ω2As4[2c27As46c26(ρ1Aa11)],3ˆGˆunˆv2n|s=0=ωAs4[2c28(As4)24c27As4(ρ1Aa11)+6c26(ρ1Aa11)2],3ˆGˆv3n|s=0=6(1+Aa11ρ)As4[c26(1+Aa11)2+c27As4(1+Aa11)+c28(As4)2(2(1+Aa11)c26+As4c27)ρ+c26ρ2].

    Based on the aforementioned calculations, we establish the theorem regarding the existence and direction of Neimark-Sacker bifurcation.

    Theorem 5. If (h,α,r,k,K,b,a,s)N.S. and Ω0, then model (1.4) undergoes a Neimark-Sacker bifurcation at the equilibrium point E1(x,y) when the parameter s varies in the vicinity of s4. Moreover, if Ω<0 (respectively, Ω>0), then an attracting (respectively, repelling) closed invariant curve bifurcates from the equilibrium point for s>0 (respectively, s<0).

    Chaos often has detrimental effects on biological systems, disrupting the ecological balance of populations and directly influencing long-term population growth projections. Implementing effective control policies not only safeguards the size of ecological populations but also establishes a strong foundation for the sustainable exploitation of ecological resources [27,28]. This section explores two control methods aimed at effectively managing the chaos generated by model (1.4).

    In this subsection, the state feedback control method [22] will be employed to regulate the chaos exhibited by model (1.4). To achieve this, we introduce the following controlled model.

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn]+Sn,yn+1=yn+hαΓ(α+1)[syn(1ynxn)], (4.1)

    which corresponds to model (1.4). The feedback controlling force is defined as

    Sn=p1(xnx)p2(yny), (4.2)

    where (x,y) represents the positive equilibrium point of the model (1.4), and p1, p2 stand for the feedback gains. The Jacobian matrix of the controlled model (4.1) evaluated at the positive equilibrium point E1(x,y) is given by

    J1(x,y)=(1+Aa11p1Aa12p2As1As), (4.3)

    where the variables A, a11, and a12 are defined in Eq (2.2). The corresponding characteristic equation of the Jacobian matrix J1(x,y) is

    λ2(2+Aa11Asp1)λ+(1+Aa11p1)(1As)(Aa12p2)As=0. (4.4)

    Let λ1 and λ2 be the roots of the Eq (4.3), then

    λ1λ2=(1+Aa11p1)(1As)(Aa12p2)As. (4.5)

    The lines of marginal stability l1, l2, and l3 are derived by solving λ1λ2=1, λ1=1 and λ2=±1, respectively. These conditions ensure that |λ1,2|=1. Then, we derive the marginal stability lines as follows:

    l1:(1As)p1Asp2=A2s(a11+a12)+A(a11s), (4.6)
    l2:p1+p2=A(a11+a12), (4.7)
    l3:(2As)p1Asp2=A2s(a11+a12)+2A(a11s)+4. (4.8)

    Therefore, l1, l2, and l3 in the (p1,p2)-plane form a triangular region which leads to |λ1,2|<1.

    Theorem 6. If p1 and p2 lie within a triangular region bounded by the lines l1, l2, and l3, it can be concluded that the model (4.1) is stable.

    Next, we apply the hybrid control approach proposed by [23] to control chaos. The controlled model of (1.4) with the hybrid control approach is depicted below:

    {xn+1=ρxn+ρhαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn]+(1ρ)xn,yn+1=ρyn+ρhαΓ(α+1)[syn(1ynxn)]+(1ρ)yn. (4.9)

    where 0<ρ<1, and the controlled strategy in (4.9) combines feedback control and parameter perturbation. By appropriately choosing the controlled parameter ρ, the chaotic behaviors of the equilibrium point (x,y) of the controlled model (4.9) can be accelerated (delayed) or even entirely eliminated. The Jacobian matrix of the controlled model (4.9), evaluated at the positive equilibrium point (x,y), is given by

    J2(x,y)=(1+Aρa11Aρa12Aρs1Aρs), (4.10)

    where the variables A, a11, and a12 are defined in Eq (2.2). Then, the positive equilibrium point (x,y) of the controlled model (4.9) is locally asymptotically stable if the roots of the characteristic polynomial of (4.10) lie within the open unit disk. According to the Jury condition, the equilibrium point of the model remains stable if, and only if, the following conditions are met:

    |2+Aρa11Aρs|<1+(1+Aρa11)(1Aρs)(Aρa12)(Aρs)<2.

    Theorem 7. If |2+Aρa11Aρs|<1+(1+Aρa11)(1Aρs)(Aρa12)(Aρs)<2 can hold, it can be concluded that the model (4.9) is stable.

    In this section, by exploring specific cases of model (1.4), we confirm the theoretical analysis above and discover new intriguing complex dynamic behaviors. Additionally, we validate the effectiveness of linear feedback techniques and hybrid control strategies for chaos control through numerical simulations.

    Take

    h=0.69,α=0.85,r=2.27,k=1.64,K=2.16,b=0.24,a=2.02. (5.1)

    After straightforward calculations, we determine the positive equilibrium point E1=(1.7498,1.7498) and the critical bifurcation value s1=2.8472. The Jacobian matrix of model (1.4) evaluated at E1 with s=s1 is given by

    J(E1)=(0.67330.14972.19651.1965),

    whose characteristic equation is

    F(λ)=λ2+0.5232λ0.4768=0. (5.2)

    The roots of (5.2) are λ1=1 and λ2=0.4768. Moreover, we have β1=0.79590 and β2=1.0080>0. According to Theorem 4, model (1.4) undergoes period-doubling bifurcation at E1 as s passes through s1. This behavior, verified through corresponding bifurcation diagrams shown in Figure 1, utilizes initial conditions (x0,y0)=(1.74,1.74) and varies s in the range [2.7,3.7].

    Figure 1.  Bifurcation diagrams of model (1.4) with parameter values as given in (5.1) and initial conditions (1.74,1.74).

    Take

    h=0.75,α=0.8,r=4.2,k=0.5,K=3.5,b=3.6,a=8.1. (5.3)

    After straightforward calculations, we determine the positive equilibrium point E1=(2.1746,2.1746) and the critical bifurcation value s4=1.1872. The Jacobian matrix of model (1.4) evaluated at E1 with s=s4 is given by

    J(E1)=(0.07120.98841.01260.0126),

    whose characteristic equation is

    F(λ)=λ20.0586λ+1=0. (5.4)

    The roots of (5.4) are λ1,2=0.02930±0.9996i. Moreover, we have d=0.39120 and Ω=1392.0464<0. According to Theorem 5, model (1.4) undergoes Neimark-Sacker bifurcation at E1 as s passes through s4. This behavior, verified through corresponding bifurcation diagrams shown in Figure 2, utilizes initial conditions (x0,y0)=(2.2,2.2) and varies s in the range [1.1,1.6].

    Figure 2.  Bifurcation diagrams of model (1.4) with parameter values as given in (5.3) and initial conditions (2.2,2.2).

    We chose the parameter values as follows:

    h=0.75,α=0.8,r=4.2,k=0.5,K=3.5,b=3.6,a=8.1,s=1.58. (5.5)

    By Figure 3, we can get that the variables xn and yn in the model (1.4) are in a chaotic state.

    Figure 3.  Plots for model (1.4) with parameter values as given in (5.5) and initial conditions (2.2,2.2).

    We use the same parameter values as in (5.5). In Figure 4, the triangular region defined by Theorem 6 bounds the parameters p1 and p2. Inside this region, the chaotic behavior produced by model (1.4) is effectively managed, resulting in asymptotic convergence toward the equilibrium point E1=(2.1746,2.1746).

    Figure 4.  Stability region for the controlled model (4.1).

    In this case, with feedback gains set to p1=0.57 and p2=0.38 and the controller activated at the 3000th iteration of the model, Figure 5 demonstrates the control effect. A chaotic trajectory is successfully stabilized at the equilibrium point E1=(2.1746,2.1746). This indicates that the feedback control approach is effective in mitigating bifurcation and chaos.

    Figure 5.  Plots for controlled model (4.1) with parameter values as given in (5.5), p1=0.57, p2=0.38, and initial conditions (2.2,2.2).

    Finally, using the same parameter values as in (5.5), the Jacobian matrix of the controlled model (4.9), evaluated at E1, is given by:

    (10.9288ρ0.988439ρ1.34765ρ11.34765ρ). (5.6)

    The characteristic polynomial of (5.6) is given by

    λ2(22.27645ρ)λ+2.58377ρ22.27645ρ+1=0. (5.7)

    The roots of (5.7) lie within the open unit disk if, and only if, 0<ρ<0.881058. Additionally, the plots for xn and yn of the controlled model (4.9) are shown in Figure 6 with ρ=0.87. From Figure 6, it is clear that the positive equilibrium point E1 is stable. Therefore, it can be concluded that employing the hybrid control approach is effective in mitigating bifurcation and chaos.

    Figure 6.  Plots for controlled model (4.9) with parameter values as given in (5.5), ρ=0.87, and initial conditions (2.2,2.2).

    In this study, we proposed a fractional Leslie-Gower model with a Holling type Ⅱ functional response and antipredator behavior. By employing the piecewise constant approximation method, we derived the discrete model (1.4) and analyzed its dynamical behavior, including the existence and stability of equilibrium points and the possibility of local bifurcations. The following conclusions can be drawn from our research:

    (1) The discrete model (1.4) has two equilibrium points: E0(K,0) and E1(x,y). The only positive coexistence equilibrium point, E1(x,y), reflects the coexistence of predators and prey.

    (2) Our theoretical analysis and numerical simulations of the positive equilibrium point E1(x,y) indicate that the model undergoes period-doubling bifurcation and Neimark-Sacker bifurcation under specific parameter conditions. Figures 1 and 2 illustrate how these bifurcations can lead to chaotic behavior at E1.

    (3) By applying state feedback control and hybrid control methods, we effectively managed the chaotic behavior generated by the model (1.4), as shown in Figures 36. These interventions mitigated the adverse effects of chaos and bifurcations, consequently enhancing ecosystem resilience.

    This study advances our understanding of the complex dynamics in ecological models influenced by fear effects and provides practical techniques for controlling chaotic behavior in such models. Future work could explore different discretization methods for the model, and new parameters could be chosen to study the influence of various ecological effects on population dynamics.

    Yao Shi: Formal analysis, validation, writing-original draft & editing, visualization, methodology; Zhenyu Wang: Conceptualization, investigation, writing-review & editing, software, supervision. All authors have read and approved the final version of the manuscript for publication.

    This work was supported by the Innovation Foundation of Hebei University of Engineering (Grant Nos. SJ2401002097), the National Natural Science Foundation of China (Grant Nos. 12401519), and the Central Guidance on Local Science and Technology Development Fund of Hebei Province (Grant Nos. 246Z1825G).

    The authors declare no conflict of interest.



    [1] L. A. Zadeh, Fuzzy sets, Inform. Control, 8 (1965), 338–353. https://doi.org/10.1016/S0019-9958(65)90241-X doi: 10.1016/S0019-9958(65)90241-X
    [2] I. B. Turksen, Interval valued fuzzy sets based on normal forms, Fuzzy Set. Syst., 20 (1986), 191–210. https://doi.org/10.1016/0165-0114(86)90077-1 doi: 10.1016/0165-0114(86)90077-1
    [3] K. T. Atanassov, Intuitionistic fuzzy sets, In: Intuitionistic Fuzzy Sets, Studies in Fuzziness and Soft Computing, Physica, Heidelberg, 35 (1999), 1–137. https://doi.org/10.1007/978-3-7908-1870-3_1
    [4] K.T. Atanassov, Interval valued intuitionistic fuzzy sets, In: Intuitionistic Fuzzy Sets, Studies in Fuzziness and Soft Computing, Physica, Heidelberg, 35 (1999). https://doi.org/10.1007/978-3-7908-1870-3_2
    [5] R. R. Yager, Pythagorean membership grades in multi-criteria decision making, IEEE T. Fuzzy Syst., 22 (2013), 958–965. https://doi.org/10.1109/TFUZZ.2013.2278989 doi: 10.1109/TFUZZ.2013.2278989
    [6] X. Peng, Y. Yang, Fundamental properties of interval-valued Pythagorean fuzzy aggregation operators, Int. J. Intell. Syst., 31 (2016), 444–487. https://doi.org/10.1002/int.21790 doi: 10.1002/int.21790
    [7] R. R. Yager, Generalized orthopair fuzzy sets, IEEE T. Fuzzy Syst., 25 (2016), 1222–1230. https://doi.org/10.1109/TFUZZ.2016.2604005 doi: 10.1109/TFUZZ.2016.2604005
    [8] B. P. Joshi, A. Singh, P. K. Bhatt, K. S. Vaisla, Interval valued q-rung orthopair fuzzy sets and their properties, J. Intell. Fuzzy Syst., 35 (2018), 5225–5230. https://doi.org/10.3233/JIFS-169806 doi: 10.3233/JIFS-169806
    [9] D. Molodtsov, Soft set theory—first results, Comput. Math. Appl., 37 (1999), 19–31. https://doi.org/10.1016/S0898-1221(99)00056-5 doi: 10.1016/S0898-1221(99)00056-5
    [10] K. Hayat, M. I. Ali, B. Y. Cao, F. Karaaslan, New results on type-2 soft sets, Hacet. J. Math. Stat., 47 (2018), 855–876. https://doi.org/10.15672/HJMS.2017.484 doi: 10.15672/HJMS.2017.484
    [11] K. Hayat, T. Mahmood, Some applications of bipolar soft set: Characterizations of two isomorphic hemi-rings via BSI-h-ideals, British J. Math. Comput. Sci., 13 (2016), 1–21. https://doi.org/10.9734/BJMCS/2016/22028 doi: 10.9734/BJMCS/2016/22028
    [12] Y. Jiang, Y. Tang, Q. Chen, H. Liu, J. Tang, Interval-valued intuitionistic fuzzy soft sets and their properties, Comput. Math. Appl., 60 (2010), 906–918. https://doi.org/10.1016/j.camwa.2010.05.036 doi: 10.1016/j.camwa.2010.05.036
    [13] X. Peng, Y. Yang, J. Song, Pythagoren fuzzy soft set and its application, Comput. Eng., 41 (2015), 224–229.
    [14] R. M. Zulqarnain, I. Siddique, A. Iampan, D. Baleanu, Aggregation operators for Interval-valued Pythagorean fuzzy soft set with their application to solve multi-attribute group decision making problem, CMES-Comput. Model. Eng., 131 (2022), 1717–1750. https://doi.org/10.32604/cmes.2022.019408 doi: 10.32604/cmes.2022.019408
    [15] A. Hussain, M. I. Ali, T. Mahmood, M. Munir, q-Rung orthopair fuzzy soft average aggregation operators and their application in multicriteria decision-making, Int. J. Intell. Syst., 35 (2020), 571–599. https://doi.org/10.1002/int.22217 doi: 10.1002/int.22217
    [16] X. Yang, K. Hayat, M. S. Raja, N. Yaqoob, C. Jana, Aggregation and interaction aggregation soft operators on interval-valued q-rung orthopair fuzzy soft environment and application in automation company evaluation, IEEE Access, 10 (2022), 91424–91444. https://doi.org/10.1109/ACCESS.2022.3202211 doi: 10.1109/ACCESS.2022.3202211
    [17] F. Smarandache, Extension of soft set to hypersoft set, and then to plithogenic hypersoft set, Neutrosophic Sets Sy., 22 (2018), 168–170.
    [18] S. Debnath, Interval-valued intuitionistic hypersoft sets and their algorithmic approach in multi-criteria decision making, Neutrosophic Sets Sy., 48 (2022), 226–250.
    [19] R. M. Zulqarnain, X. L. Xin, M. Saeed, A development of Pythagorean fuzzy hypersoft set with basic operations and decision-making approach based on the correlation coefficient, Theory Appl. Hypersoft Set, 2021 (2021), 85–106.
    [20] R. Zulqarnain, I. Siddique, R. Ali, F. Jarad, A. Iampan, Aggregation operators for interval-valued pythagorean fuzzy hypersoft set with their application to solve MCDM problem, CMES-Comput. Model. Eng., 135 (2023), 619–651. https://doi.org/10.32604/cmes.2022.022767 doi: 10.32604/cmes.2022.022767
    [21] S. Khan, M. Gulistan, N. Kausar, S. Kousar, D. Pamucar, G. M. Addis, Analysis of cryptocurrency market by using q-rung orthopair fuzzy hypersoft set algorithm based on aggregation operator, Complexity, 2022 (2022). https://doi.org/10.1155/2022/7257449 doi: 10.1155/2022/7257449
    [22] R. M. Zulqarnain, I. Siddique, A. Mahboob, H. Ahmad, S. Askar, S. H. Gurmani, Optimizing construction company selection using einstein weighted aggregation operators for q-rung orthopair fuzzy hypersoft set, Sci. Rep., 13 (2023), 6511. https://doi.org/10.1038/s41598-023-32818-8 doi: 10.1038/s41598-023-32818-8
    [23] C. Samantra, S. Datta, S. S. Mahapatra, Application of fuzzy based VIKOR approach for multi-attribute group decision making (MAGDM): A case study in supplier selection, Decis. Mak. Manuf. Serv., 6 (2012), 25–39. https://doi.org/10.7494/dmms.2012.6.1.25 doi: 10.7494/dmms.2012.6.1.25
    [24] S. Zhang, Y. Hou, S. Zhang, M. Zhang, Fuzzy control model and simulation for nonlinear supply chain system with lead times, Complexity, 2017 (2017), 2017634. https://doi.org/10.1155/2017/2017634 doi: 10.1155/2017/2017634
    [25] S. Zhang, C. Zhang, S. Zhang, M. Zhang, Discrete switched model and fuzzy robust control of dynamic supply chain network, Complexity, 2018 (2018), 3495096. https://doi.org/10.1155/2018/3495096 doi: 10.1155/2018/3495096
    [26] C. Yu, Y. Shao, K. Wang, L. Zhang, A group decision making sustainable supplier selection approach using extended TOPSIS under interval-valued Pythagorean fuzzy environment, Expert Syst. Appl., 121 (2019), 1–17. https://doi.org/10.1016/j.eswa.2018.12.010 doi: 10.1016/j.eswa.2018.12.010
    [27] S. Zhang, P. Zhang, M. Zhang, Fuzzy emergency model and robust emergency strategy of supply chain system under random supply disruptions, Complexity, 2019 (2019), 3092514. https://doi.org/10.1155/2019/3092514 doi: 10.1155/2019/3092514
    [28] S. Zhang, S. Li, S. Zhang, M. Zhang, Decision of lead-time compression and stable operation of supply chain, Complexity, 2017 (2017), 7436764. https://doi.org/10.1155/2017/7436764 doi: 10.1155/2017/7436764
    [29] S. H. Gurmani, Z. Zhang, R. M. Zulqarnain, S. Askar, An interaction and feedback mechanism-based group decision-making for emergency medical supplies supplier selection using T-spherical fuzzy information, Sci. Rep., 13 (2023), 8726. https://doi.org/10.1038/s41598-023-35909-8 doi: 10.1038/s41598-023-35909-8
    [30] H. Wang, Sustainable circular supplier selection in the power battery industry using a linguistic T-spherical fuzzy MAGDM model based on the improved ARAS method, Sustainability, 14 (2022), 7816. https://doi.org/10.3390/su14137816 doi: 10.3390/su14137816
    [31] J. Wu, Y. Jin, M. Zhou, M. Cao, Y. Liu, A group consensus decision making based sustainable supplier selection method by combing DEMATEL and VIKOR, J. Intell. Fuzzy Syst., 42 (2022), 2595–2613. https://doi.org/10.3233/JIFS-211929 doi: 10.3233/JIFS-211929
    [32] Y. Diao, Q. Zhang, Optimization of management mode of small- and medium- sized enterprises based on decision tree model, J. Math., 2021 (2021), 2815086. https://doi.org/10.1155/2021/2815086 doi: 10.1155/2021/2815086
    [33] J. Wang, H. Gao, G. Wei, Y. Wei, Methods for multiple-attribute group decision making with q-rung interval-valued orthopair fuzzy information and their applications to the selection of green suppliers, Symmetry, 11 (2019), 56. https://doi.org/10.3390/sym11010056 doi: 10.3390/sym11010056
    [34] A. Ali, K. Ullah, A. Hussain, An approach to multi-attribute decision-making based on intuitionistic fuzzy soft information and Aczel-Alsina operational laws, J. Decis. Anal. Intell. Comput., 3 (2023), 80–89. https://doi.org/10.31181/jdaic10006062023a doi: 10.31181/jdaic10006062023a
    [35] R. M. Zulqarnain, X. L. Xin, H. Garg, W. A. Khan, Aggregation operators of pythagorean fuzzy soft sets with their application for green supplier chain management, J. Intell. Fuzzy Syst., 40 (2021), 5545–5563. https://doi.org/10.3233/JIFS-202781 doi: 10.3233/JIFS-202781
    [36] R. M. Zulqarnain, I. Siddique, S. Ahmad, A. Iampan, G. Jovanov, Đ. Vranješ, et al., Pythagorean fuzzy soft Einstein ordered weighted average operator in sustainable supplier selection problem, Math. Probl. Eng., 2021 (2021), 1–16. https://doi.org/10.1155/2021/2559979 doi: 10.1155/2021/2559979
    [37] S. Zhang, M. Zhang, Mitigation of bullwhip effect in closed‐loop supply chain based on fuzzy robust control approach, Complexity, 2020 (2020), 1085870. https://doi.org/10.1155/2020/1085870 doi: 10.1155/2020/1085870
    [38] A. Hussain, K. Ullah, T. Senapati, S. Moslem, Energy supplier selection by TOPSIS method based on multi-attribute decision-making by using novel idea of complex fuzzy rough information, Energy Strateg. Rev., 54 (2024), 101442. https://doi.org/10.1016/j.esr.2024.101442 doi: 10.1016/j.esr.2024.101442
    [39] W. Wang, X. Liu, Intuitionistic fuzzy information aggregation using Einstein operations, IEEE T. Fuzzy Syst., 20 (2012), 923–938. https://doi.org/10.1109/TFUZZ.2012.2189405 doi: 10.1109/TFUZZ.2012.2189405
    [40] W. Wang, X. Liu, Intuitionistic fuzzy geometric aggregation operators based on Einstein operations, Int. J. Intell. Syst., 26 (2011), 1049–1075. https://doi.org/10.1002/int.20498 doi: 10.1002/int.20498
    [41] P. Liu, Y. Li, Y. Chen, Some generalized Einstein aggregation operators based on the interval-valued intuitionistic fuzzy numbers and their application to group decision making, Sci. Iran., 22 (2015), 2684–2701. https://doi.org/10.1109/TFUZZ.2013.2248736 doi: 10.1109/TFUZZ.2013.2248736
    [42] W. Wang, X. Liu, Some interval-valued intuitionistic fuzzy geometric aggregation operators based on einstein operations, In 2012 9th International Conference on Fuzzy Systems and Knowledge Discovery, IEEE, 2012,604–608. https://doi.org/10.1109/FSKD.2012.6234364
    [43] H. Garg, A new generalized Pythagorean fuzzy information aggregation using Einstein operations and its application to decision making, Int. J. Intell. Syst., 31 (2016), 886–920. https://doi.org/10.1002/int.21809 doi: 10.1002/int.21809
    [44] M. Asif, U. Ishtiaq, I. K. Argyros, Hamacher aggregation operators for Pythagorean fuzzy set and its application in multi-attribute decision-making problem, Spectrum Oper. Res., 2 (2025), 27–40. https://doi.org/10.31181/sor2120258 doi: 10.31181/sor2120258
    [45] K. Rahman, S. Abdullah, M. S. A. Khan, Some interval-valued Pythagorean fuzzy Einstein weighted averaging aggregation operators and their application to group decision making, J. Intell. Syst., 29 (2018), 393–408. https://doi.org/10.1515/jisys-2017-0212 doi: 10.1515/jisys-2017-0212
    [46] K. Rahman, A. Ali, S. Abdullah, Multiattribute group decision making based on interval-valued Pythagorean fuzzy Einstein geometric aggregation operators, Granular Comput., 5 (2020), 361–372. https://doi.org/10.1007/s41066-019-00154-w doi: 10.1007/s41066-019-00154-w
    [47] M. Deveci, I. Gokasar, D. Pamucar, D. M. Coffman, E. Papadonikolaki, Safe E-scooter operation alternative prioritization using a q-rung orthopair fuzzy Einstein based WASPAS approach, J. Clean. Prod., 347 (2022), 131239. https://doi.org/10.1016/j.jclepro.2022.131239 doi: 10.1016/j.jclepro.2022.131239
    [48] Y. Xu, A two-stage multi-criteria decision-making method with interval-valued q-Rung Orthopair fuzzy technology for selecting bike-sharing recycling supplier, Eng. Appl. Artif. Intell., 119 (2023), 105827. https://doi.org/10.1016/j.engappai.2023.105827 doi: 10.1016/j.engappai.2023.105827
    [49] R. M. Zulqarnain, I. Siddique, F. Jarad, Y. S. Hamed, K. M. Abualnaja, A. Iampan, Einstein aggregation operators for Pythagorean fuzzy soft sets with their application in multiattribute group decision-making, J. Funct. Space., 2022 (2022), 1–21. https://doi.org/10.1155/2022/1358675 doi: 10.1155/2022/1358675
    [50] R. M. Zulqarnain, H. K. U. Rehman, J. Awrejcewicz, R. Ali, I. Siddique, F. Jarad, et al., Extension of Einstein average aggregation operators to medical diagnostic approach under q-rung orthopair fuzzy soft set, IEEE Access, 10 (2022), 87923–87949. https://doi.org/10.1109/ACCESS.2022.3199069 doi: 10.1109/ACCESS.2022.3199069
    [51] R. M. Zulqarnain, R. Ali, J. Awrejcewicz, I. Siddique, F. Jarad, A. Iampan, Some Einstein geometric aggregation operators for Q-rung orthopair fuzzy soft set with their application in MCDM, IEEE Access, 10 (2022), 88469–88494. https://doi.org/10.1109/ACCESS.2022.3199071 doi: 10.1109/ACCESS.2022.3199071
    [52] R. M. Zulqarnain, H. Naveed, I. Siddique, J. C. R. Alcantud, Transportation decisions in supply chain management using interval-valued q-rung orthopair fuzzy soft information, Eng. Appl. Artif. Intell., 133 (2024), 108410. https://doi.org/10.1016/j.engappai.2024.108410 doi: 10.1016/j.engappai.2024.108410
    [53] P. Sunthrayuth, F. Jarad, J. Majdoubi, R. M. Zulqarnain, A. Iampan, I. Siddique, A novel multicriteria decision-making approach for einstein weighted average operator under Pythagorean fuzzy hypersoft environment, J. Math., 2022 (2022), 1951389. https://doi.org/10.1155/2022/1951389 doi: 10.1155/2022/1951389
    [54] R. M. Zulqarnain, I. Siddique, R. Ali, F. Jarad, A. Iampan, Einstein weighted geometric operator for Pythagorean fuzzy hypersoft with its application in material selection, Comput. Model. Eng. Sci., 135 (2023), 2557–2583. https://doi.org/10.32604/cmes.2023.023040 doi: 10.32604/cmes.2023.023040
    [55] M. Sajid, K. A. Khan, A. U. Rahman, S. A. Bajri, A. Alburaikan, H. A. E. W. Khalifa, A novel algorithmic multi-attribute decision-making framework for solar panel selection using modified aggregations of cubic intuitionistic fuzzy hypersoft set, Heliyon, 10 (2024). https://doi.org/10.1016/j.heliyon.2024.e36508 doi: 10.1016/j.heliyon.2024.e36508
    [56] B. Huang, J. Miao, Q. Li, A vetoed multi-objective grey target decision model with application in supplier choice, J. Grey Syst., 34 (2022).
    [57] M. Saqlain, Sustainable hydrogen production: a decision-making approach using VIKOR and intuitionistic hypersoft sets, J. Intell. Manag. Decis., 2 (2023), 130–138. https://doi.org/10.56578/jimd020303 doi: 10.56578/jimd020303
    [58] M. T. Hamid, M. Abid, Decision support system for mobile phone selection utilizing fuzzy hypersoft sets and machine learning, J. Intell. Manag. Decis., 3 (2024), 104–115. https://doi.org/10.56578/jimd030204 doi: 10.56578/jimd030204
    [59] M. Sarwar, T. Li, Fuzzy fixed point results and applications to ordinary fuzzy differential equations in complex valued metric spaces, Hacet. J. Math. Stat., 48 (2019), 1712–1728. https://doi.org/10.15672/HJMS.2018.633 doi: 10.15672/HJMS.2018.633
    [60] Y. Xia, J. Wang, B. Meng, X. Chen, Further results on fuzzy sampled-data stabilization of chaotic nonlinear systems, Appl. Math. Comput., 379 (2020), 125225. https://doi.org/10.1016/j.amc.2020.125225 doi: 10.1016/j.amc.2020.125225
    [61] M. Gao, L. Zhang, W. Qi, J. Cao, J. Cheng, Y. Kao, et al., SMC for semi-Markov jump TS fuzzy systems with time delay, Appl. Math. Comput., 374 (2020), 125001. https://doi.org/10.1016/j.amc.2019.125001 doi: 10.1016/j.amc.2019.125001
    [62] J. Ge, S. Zhang, Adaptive inventory control based on fuzzy neural network under uncertain environment, Complexity, 2020 (2020), 6190936. https://doi.org/10.1155/2020/6190936 doi: 10.1155/2020/6190936
    [63] N. Zhang, W. Qi, G. Pang, J. Cheng, K. Shi, Observer-based sliding mode control for fuzzy stochastic switching systems with deception attacks, Appl. Math. Comput., 427 (2022), 127153. https://doi.org/10.1016/j.amc.2022.127153 doi: 10.1016/j.amc.2022.127153
    [64] Q. Sun, J. Ren, F. Zhao, Sliding mode control of discrete-time interval type-2 fuzzy Markov jump systems with the preview target signal, Appl. Math. Comput., 435 (2022), 127479. https://doi.org/10.1016/j.amc.2022.127479 doi: 10.1016/j.amc.2022.127479
    [65] Z. X. Duan, J. L. Liang, Z. R. Xiang, H∞ control for continuous-discrete systems in TS fuzzy model with finite frequency specifications, Discrete Cont. Dyn.-S, 64 (2022), 1–18.
  • This article has been cited by:

    1. Yao Shi, Xiaozhen Liu, Zhenyu Wang, Bifurcation Analysis and Chaos Control of a Discrete Fractional-Order Modified Leslie–Gower Model with Nonlinear Harvesting Effects, 2024, 8, 2504-3110, 744, 10.3390/fractalfract8120744
    2. Rutaba Saleem, Ahmad Javid, Interaction solutions and qualitative analysis of the Painlevé integrable (2+1)-dimensional extended Kadomtsev–Petviashvili equation, 2025, 128, 11100168, 611, 10.1016/j.aej.2025.05.013
  • Reader Comments
  • © 2024 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(813) PDF downloads(60) Cited by(2)

Figures and Tables

Figures(3)  /  Tables(6)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog