
In this paper, we investigate the complex dynamics of a classical discrete-time prey-predator system with the cost of anti-predator behaviors. We first give the existence and stability of fixed points of this system. And by using the central manifold theorem and bifurcation theory, we prove that the system will experience flip bifurcation and Neimark-Sacker bifurcation at the equilibrium points. Furthermore, we illustrate the bifurcation phenomenon and chaos characteristics via numerical simulations. The results may enrich the dynamics of the prey-predator systems.
Citation: Ceyu Lei, Xiaoling Han, Weiming Wang. Bifurcation analysis and chaos control of a discrete-time prey-predator model with fear factor[J]. Mathematical Biosciences and Engineering, 2022, 19(7): 6659-6679. doi: 10.3934/mbe.2022313
[1] | Minjie Zheng, Yulai Su, Guoquan Chen . An improved sampled-data control for a nonlinear dynamic positioning ship with Takagi-Sugeno fuzzy model. Mathematical Biosciences and Engineering, 2024, 21(5): 6019-6041. doi: 10.3934/mbe.2024265 |
[2] | Luyao Li, Licheng Fang, Huan Liang, Tengda Wei . Observer-based event-triggered impulsive control of delayed reaction-diffusion neural networks. Mathematical Biosciences and Engineering, 2025, 22(7): 1634-1652. doi: 10.3934/mbe.2025060 |
[3] | Yilin Tu, Jin-E Zhang . Event-triggered impulsive control for input-to-state stability of nonlinear time-delay system with delayed impulse. Mathematical Biosciences and Engineering, 2025, 22(4): 876-896. doi: 10.3934/mbe.2025031 |
[4] | Xinggui Zhao, Bo Meng, Zhen Wang . Event-triggered integral sliding mode control for uncertain networked linear control systems with quantization. Mathematical Biosciences and Engineering, 2023, 20(9): 16705-16724. doi: 10.3934/mbe.2023744 |
[5] | Jingya Wang, Ye Zhu . L2−L∞ control for memristive NNs with non-necessarily differentiable time-varying delay. Mathematical Biosciences and Engineering, 2023, 20(7): 13182-13199. doi: 10.3934/mbe.2023588 |
[6] | Wei Zou, Yanxia Shen, Lei Wang . Design of robust fuzzy iterative learning control for nonlinear batch processes. Mathematical Biosciences and Engineering, 2023, 20(11): 20274-20294. doi: 10.3934/mbe.2023897 |
[7] | Gavin D'Souza, Jenna Osborn, Shayna Berman, Matthew Myers . Comparison of effectiveness of enhanced infection countermeasures in different scenarios, using a dynamic-spread-function model. Mathematical Biosciences and Engineering, 2022, 19(9): 9571-9589. doi: 10.3934/mbe.2022445 |
[8] | Mingxia Gu, Zhiyong Yu, Haijun Jiang, Da Huang . Distributed consensus of discrete time-varying linear multi-agent systems with event-triggered intermittent control. Mathematical Biosciences and Engineering, 2024, 21(1): 415-443. doi: 10.3934/mbe.2024019 |
[9] | Wenjing Wang, Jingjing Dong, Dong Xu, Zhilian Yan, Jianping Zhou . Synchronization control of time-delay neural networks via event-triggered non-fragile cost-guaranteed control. Mathematical Biosciences and Engineering, 2023, 20(1): 52-75. doi: 10.3934/mbe.2023004 |
[10] | Ramesh Garimella, Uma Garimella, Weijiu Liu . A theoretic control approach in signal-controlled metabolic pathways. Mathematical Biosciences and Engineering, 2007, 4(3): 471-488. doi: 10.3934/mbe.2007.4.471 |
In this paper, we investigate the complex dynamics of a classical discrete-time prey-predator system with the cost of anti-predator behaviors. We first give the existence and stability of fixed points of this system. And by using the central manifold theorem and bifurcation theory, we prove that the system will experience flip bifurcation and Neimark-Sacker bifurcation at the equilibrium points. Furthermore, we illustrate the bifurcation phenomenon and chaos characteristics via numerical simulations. The results may enrich the dynamics of the prey-predator systems.
Dynamic positioning system (DPS) is a system that can automatically control the position and course of a ship [1]. It only depends on its own propulsion system. The position and course of the ship's motion are measured by various sensors equipped by the ship itself, and the computer is used to control the ship's thrust device to generate propulsion force and torque to resist the interference force caused by the external environment. Then the ship can maintain the target's positions and heading. The DPS can fix the ship in a certain position, and can also automatically adjust the heading of the ship to the optimal desired position according to the real-time direction of wind, wave and current. Additionally, it has been widely used in marine investigation ship, drilling ship, salvage ship, mining ship, cable laying ship and so on. By far, severable important control methods have been reported for DPS ([2,3,4,5,6,7]). Liang et al. [2] studies the state space model of DPS by using Cummins equation and convolution integral term. Based on this model, a static output feedback controller based on L∞ is designed. Based neural fuzzy algorithm, Fang and Lee [3] use a self-tuning controller to control the rotation speed of the outboard thruster to achieve the goal of controlling the ship's position, heading and path tracking. In Zhang et al. [4], a robust neural control algorithm for DPS based on the combination of dynamic surface control, neural network and robust damping technology is proposed. This algorithm solves the strong coupling between state variables and the uncertainty of actuator gain, reduces the parameter requirements of the DP model. In Ngongi et al. [5], a robust fuzzy controller using optimal control technology is proposed. The controller uses T-S fuzzy model to estimate the nonlinear part of DPS, and uses the generalized eigenvalue method to solve the controller gain. To deal with the issue of the model uncertainties, Li et al. [6] investigated a new collaborative design for path tracking of hybrid USVs and UAVs in the presence of structural uncertainty and external disturbances. In order to control the convergence of USV-UAV to the desired path, an adaptive fuzzy control algorithm was designed by integrating DSC and backstepping technology. In Zhang et al. [7], by combining robust neural damping with DSC technology, the gain dependent adaptive laws for the ships have been developed to address constraints of gain uncertainty and environmental disturbances.
Due to the wide application of digital circuits, microelectronic devices and modern computers, the sampled-data system has developed rapidly. The sampled-data control system is a kind of hybrid system, which combines the characteristics of discrete system and continuous system. Therefore, it can effectively reduce the consumption of computing resources, reduce the burden of communication transmission and is easier to implement in engineering. Recently, the sampled-data control theory has become one of the research hotspots in multiagent systems [8,9], fuzzy systems [10,11,12], chaotic systems [13,14,15], neural networks systems [16,17,18] and so on ([19,20,21,22,23,24]), and has wide applications in train systems [25], airship systems [26] and ship systems [27]. Scholars have made corresponding achievements in sampled-data control ([28,29,30,31,32,33]). Qian et al. [29] discussed the global stabilization of a class of nonlinear systems with higher order powers based on output feedback sampling control. A linear output feedback sampled-data controller is constructed in Zhang et al. [30], which can semi globally asymptotically stabilize the system with high order and linear growth nonlinearity. Zou et al. [31] designed a sampling control protocol using the backstep-ping method to ensure the leader follower consistency of second-order multi-agent with external interference and unpredictable speed information. In Zhang et al. [32], the global output feedback decentralized control problem for a class of uncertain interconnected systems based on sampled-data control is studied. By constructing sampled high gain observers and controllers, the corresponding system is guaranteed to be globally asymptotically stable.
DPS is a system controller by computer. Recently, sampled-data theory has been used to solve the control issue of DPS. In Katayama [34], the nonlinear sampled-data control issue of DPS is discussed by using the integrator backstepping technique and a semiglobal stability controller is designed. In Katayama and Aoki [35], the linear trajectory tracking issue for sampled-data DPS is discussed. Additionally, a state-feedback controller and observer are designed based on the Euler approximation model. In Zheng et al. [36], robust sampled-data control of neutral system for DPS is discussed. the state-derivative control law is designed. By combining the delayed decomposition technique and Wirtinger based integral inequality, less conservative results are obtained. In Yang and Zheng [37], the fault-tolerant control for sampled-data DPS with actuator failures is discussed. Zheng et al. [38] discusses the nonlinear fuzzy sampled-data control issue for DPS, and a free weighting matrix is used to ensure the stability of the system. In Chen et al. [39], the tracking control of DPS with actuator-failure and aperiodic measurement information is studied. Additionally, aperiodic sampling data controller is designed to guarantee that the DPS has excellent tracking performance.
Although some researches have been reported on the sampled-data control DPS, there is still much room for improvement. One reason is that the term in LKF share common quadratic function, which may lead to conservative result. The other reason is that due to the neglect of some important and useful terms in LKF, the available features of the actual sampling mode are not fully captured, which may lead to conservatism to a certain extent.
Motivated by the discussion, the design and stabilization of sampled-data for LC with DPS is discussed in this paper. First, the design method for LC is given. Second, the motion model of LC with DPS is established. Third, a novel LKF is introduced to capture the sampling pattern fully. Then, the stability conditions of the system are derived by LMIs and corresponding design algorithm of the controller is given. Finally, a specific example verifies the effectiveness of the proposed method.
To solve the current problem such as ship function, aesthetic and humanization, we should base on the related theories of product system design and integrate design art and ship engineering knowledge in a scientific, systematic way to obtain the method of ship design innovation.
Due to the highly humanization, highly esthetic, highly functional and highly technology-integrated characteristics, it is of typical significance that we take cruise ship as representatives to do research and for educational means.
Then, a 150 m five-star luxury cruise is designed as follow. Figure 1 shows the three-dimensional modelling, Figure 2 is the panoramic effect picture and Figure 3 is nightscape effect picture of the luxury cruise.
This article focuses on the fixed point control of ships on the sea surface, so the influence on heave, sway and pitch is ignored, and only the three degrees of freedom of heave, sway and yaw which have significant effects on ship are considered. Hence, the following mathematical model of the LC with DPS is considered.
M˙υ(t)+Dυ(t)=u(t)+w(t),˙η(t)=J(ψ(t))υ(t), | (1) |
where
M=[m−X˙u000m−Y˙vmxG−Y˙r0mxG−Y˙rIz−N˙r],D=[−Xu000−Yv−Yr0−Nv−Nr],J(ψ(t))=[cos(ψ(t))−sin(ψ(t))0sin(ψ(t))cos(ψ(t))0001]. |
where η(t)=[xa(t)ya(t)ψ(t)]T shows the position and heading angle of the ship in the northeast coordinate system, υ(t)=[p(t)v(t)r(t)]T describes the speed and heading rotation rate of the ship in the hull coordinate system (see Figure 4). J(ψ) is the rotation matrix, u(t)=[u1(t),u2(t),u3(t)] are the control force vectors provided by the thruster, where u1(t) is the forward control force, u2(t) is the lateral drift control force, and u3(t) is the yaw control torque. M is the inertial matrix, and D is damping matrix. m is the ship's mass, and Iz is the rotational inertia; Xu,Yv,Nr are the linear damping coefficients of the ship in three degrees of freedom, and X˙u,Y˙v,N˙r are the additional mass generated by the ship in the surge, sway and yaw. w(t)=[w1(t),w2(t),w3(t)] is an unknown disturbance term that represents the environmental disturbances such as wind, waves and currents.
Assume that the environment disturbance is unknown but energy bounded, that is, |wi(t)|⩽h,h is upper bound of the disturbances.
Assume that yaw angle ψ(t) is small enough, which means that cos(ψ(t))≈1,sin(ψ(t))≈0, then
J(ψ(t))≈I. | (2) |
Define
x(t)=[η(t)υ(t)]T=[xp(t)yp(t)ψ(t)p(t)q(t)r(t)]T | (3) |
Substitute (3) into (1) that
˙x(t)=Ax(t)+Bu(t)+Bww(t), | (4) |
where
A=[03×3I3×303×3−M−1D],B=[03×3M−1],Bw=[03×3M−1]. |
Assume that the state variable of the DPS is sampled by the time 0=t0<t1<⋯<tk<⋯<limk→∞tk=+∞. Assume the sampling period is
tk+1−tk⩽d,∀k⩾0,d>0, |
where d is the upper bound of the sampling period. The sampled-data controller is designed
u(t)=Kx(tk),tk⩽t<tk+1, | (5) |
where K is controller gain matrix.
Substitute (5) into (4), it can be obtained
˙x(t)=Ax(t)+BKx(tk)+Bwx(t) | (6) |
The paper's aim is to propose a sampled-data controller to guarantee the stability of the system. To derive our main theorems, the following lemma is introduced.
Lemma 1([40]): For scalars τ2>τ1, any constant matrix Z∈Rn,Z>0, vector function w:[τ1,τ2]∈Rn, the following inequalities are established
−∫t−τ1t−τ2wT(α)Zw(α)dα⩽−1τ2−τ1(∫t−τ1t−τ2wT(α)dα)TZ(∫t−τ1t−τ2w(α)dα) | (7) |
−∫−τ1−τ2∫tt+αwT(s)Zw(s)dsdα⩽−2τ22−τ21(∫−τ1−τ2∫tt+αwT(s)dsdα)TZ(∫−τ1−τ2∫tt+αw(s)dsdα) | (8) |
In the section, the stability conditions for system (6) are proposed. Consequently, a sampled-data controller is developed. First, the notations are defined as follow:
τ(t)=t−tk,d(t)=d−τ(t),ς(t)=[xT(t)xT(tk)∫ttkxT(s)ds]Tζ(t)=[xT(t)˙x(t)TxT(tk)∫ttkxT(s)ds] |
Theorem 1: The system (6) with w(t)=0 is asymptotically stable, if there exist symmetric positive matrices P>0,Q>0,Z>0,R>0,U>0, X11,X13,X22,X23,X33,G,Mi,Ni(i=1,2,3,4) and scales d>0,εi(i=1,2,3,4), such that the following LMIs hold:
Ψ1=[P+d(X11+XT11+Q)d(X11+X22−Q)dX13∗d(X11+XT22+Q)dX23∗∗dX33]>0 | (9) |
Ψ2=[Ξ11+Ω11Ξ12+Ω12Ξ13+Ω13Ξ14+dXT33∗Ξ22+Ω22Ξ23+Ω23Ξ24+dX13∗∗Ξ33+Ω33Ξ34∗∗∗Ξ44]<0 | (10) |
Ψ3==[Π11Π12∗Π22]<0. | (11) |
where
Π11=[Ξ11+Ω11Ξ12+dNT2Ξ13+dNT3Ξ14+dNT4∗Ξ22+d24ZΞ23Ξ24∗∗Ξ33+Ω33Ξ34∗∗∗Ξ44],Π12=[−dM1dN1−dM2dN2−dRT12−dM3dN3−dM4dN4],Π22=[−dQ−dR110∗−2Z],Ξ11=−(X11+XT11)+MT1+M1+ε1AGT+ε1GATΞ12=P+MT2−ε1G−ε2ATGTΞ13=(X11+X22)+MT3−M1+ε1GBK+ε3ATGTΞ14=−X13+MT4−N1+ε4ATGTΞ22=−ε2G−ε2GTΞ23=−M2+ε2GBK−ε3GTΞ24=−ε4GT−N2Ξ33=−(X22+XT22)−MT3−M3+ε3GBK+ε3KTBTGTΞ34=−X23−U12−MT4−N3+ε4KTBTGTΞ44=−X33−1dU22−NT4−N4Ω11=d(X13+XT13)+dU22−ZΩ12=dX11+dXT11Ω13=dXT23+dUT12+ZΩ22=dR11+dQ+d24ZΩ23=−d(X11+X22)+dR12Ω33=dR22+dU11−Z. |
Proof. Consider the following LKF
V(t)=4∑i=1Vi(t),t∈[tk,tk+1)V1(t)=x(t)TPx(t)+d(t)ς(t)TXς(t)+d(t)∫ttk˙x(s)TQ˙x(s)dsV2(t)=d(t)∫ttk[˙x(s)x(tk)]TR[˙x(s)x(tk)]dsV3(t)=d(t)∫ttk[x(tk)x(s)]TU[x(tk)x(s)]dsV4(t)=d(t)∫0−d(t)∫tt+θ˙x(s)TZ˙x(s)dsdθ | (12) |
where
X=[X11+XT11−X11−X22X13∗X22+XT22X23∗∗X33],R=[R11R12∗R22],U=[U11U12∗U22] |
By Lemma 1, it can be verified that
V1(t)⩾x(t)TPx(t)+d(t)ς(t)TXς(t)+d(t)d[x(t)−x(tk)]TQ[x(t)−x(tk)]=d(t)dςT(t)Ψij1ς(t) | (13) |
From LMI (9), V1(t)⩾0 can be concluded, which means that V(t)⩾0.
Calculate the derivative of V(t), it yields that
˙V1(t)=2˙x(t)TPx(t)+2d(t)ςT(t)X[˙xT(t)0xT(t)]T−ςT(t)Xς(t)+d(t)˙xT(t)Q˙x(t)−∫ttk˙xT(s)Q˙x(s)ds,˙V2(t)=d(t)[˙x(t)x(tk)]TR[˙x(t)x(tk)]−∫tk[˙x(s)x(tk)]TR[˙x(s)x(tk)]ds,˙V3(t)=d(t)[x(tk)x(t)]TU[x(tk)x(t)]−∫ttk[x(tk)x(s)]TU[x(tk)x(s)]ds=d(t)[x(tk)x(t)]TU[x(tk)x(t)]−τ(t)xT(tk)U11x(tk)=−2xT(tk)U12∫ttkx(s)ds−∫ttkxT(s)U22x(s)ds˙V4(t)=d(t)τ(t)˙xT(t)Z˙x(t)−d(t)∫ttk˙xT(s)Z˙x(s)ds−∫0−τ(t)∫tt+α˙xT(s)Z˙x(s)dsdα⩽(d(t)+τ(t))24˙xT(t)Z˙x(t)−d(t)∫ttk˙xT(s)Z˙x(s)ds−∫0−τ(t)∫tt+α˙xT(s)Z˙x(s)dsdα⩽d4d(t)˙xT(t)Z˙x(t)+d4τ(t)˙xT(t)Z˙x(t)−d(t)d[x(tk)x(t)]T[Z−Z−ZZ][x(tk)x(t)]−∫0−τ(t)∫tt+α˙xT(s)Z˙x(s)dsdα | (14) |
By Lemma 1, one can obtain
0=2ξT(t)N×[τ(t)x(t)−∫ttkx(s)ds−∫0−τ(t)∫tt+α˙x(s)dsdα]⩽2τ(t)ξT(t)Nx(t)−2ξT(t)N∫ttkx(s)ds+τ2(t)2ξT(t)NZ−1NTξ(t)+2τ2(t)[∫0−τ(t)∫tt+α˙xT(s)dsdα]Z[∫0−τ(t)∫tt+α˙x(s)dsdα]⩽2τ(t)ξT(t)Nx(t)−2ξT(t)N∫ttkx(s)ds+τ2(t)2ξT(t)NZ−1NTξ(t)+∫0−τ(t)∫tt+α˙xT(s)Z˙x(s)dsdα | (15) |
which implies
∫0−τ(t)∫tt+α˙x(s)Z˙x(s)dsdα⩽2τ(t)ξT(t)Nx(t)−2ξT(t)N∫ttkx(s)ds+dτ(t)2ξT(t)NZ−1NTξ(t) | (16) |
Combining (14) with (16) that
˙V4(t)⩽d4d(t)˙xT(t)Z˙x(t)+d4τ(t)˙xT(t)Z˙x(t)−d(t)d[x(tk)x(t)]T[Z−Z−ZZ][x(tk)x(t)]+2τ(t)ξT(t)Nx(t)−2ξT(t)N∫ttkx(s)ds+dτ(t)2ξT(t)NZ−1NTξ(t) | (17) |
For any matrixes G,M, and scalars εi,i=1,2,3,4, we have
0=ξT(t)M×[x(t)−x(tk)−∫ttk˙x(s)ds] | (18) |
0=2[ε1xT(t)G+ε2˙xT(t)G+ε3x(tk)G+ε4∫ttkxT(s)Gds]×[−˙x(t)+Ax(t)+BKx(tk)] | (19) |
From (13), (14) and (17)–(19), we obtain that
˙V(t)⩽d(t)dξT(t)Ψ3ξ(t)+1d∫ttk[ξ(t)˙x(s)]TˆΨ4[ξ(t)˙x(s)]ds | (20) |
where
ˆΨ4=[Π11Θ∗−dQ−R11]+d22[N0]Z−1[N0]TΘij=[−dMT1−dMT2−dR12−dMT3−dMT4]T | (21) |
Following the Schur complement, (21) implies ˆΨij4<0. which implies that ˙V(t)<−σ‖ when x(t) \neq 0, \sigma>0. Thus, the system (6) is asymptotically stable. This completed the proof.
Remark 1: Note that some novel terms such as d(t)\int_{ - d(t)}^0 {\int_{t + \theta }^t {\dot x{{(s)}^T}Z\dot x(s)dsd\theta } } are added in the constructed LKF, which means that the characteristics about sampling patterns are fully captured.
Remark 2: If {\mathit{\boldsymbol{Q}}} = {{\mathit{\boldsymbol{R}}}_{12}} = {{\mathit{\boldsymbol{R}}}_{22}} = {\mathit{\boldsymbol{Z}}} = 0 and {{\mathit{\boldsymbol{X}}}_{13}} = {{\mathit{\boldsymbol{X}}}_{23}} = {{\mathit{\boldsymbol{X}}}_{33}} = {{\mathit{\boldsymbol{U}}}_{13}} = {{\mathit{\boldsymbol{U}}}_{12}} = {{\mathit{\boldsymbol{U}}}_{22}} = 0 , the LKF is simplified to that in [38] and [39]. Therefore, the proposed LKF are general and have wider application scopes. Moreover, according to Newton-Leibniz formulas, free matrices M and N have been considered to handle the derivation of LKF to avoid using the bounding techniques. Therefore, the conservatism is reduced.
Remark 3: Compared with [38] and [39], the more relaxed constraint conditions are introduced by adding X-dependent term in {V_1}(t) , which means that the involved symmetric matrices are not required to be positive definite. Besides, the term d(t)\int_{{t_k}}^t {{{\dot x}^T}(s){\mathit{\boldsymbol{Q}}}\dot x(s)ds} is added in the LKF to obtain more relaxed condition and longer sampling period. Hence, the conservativeness can be further reduced.
Remark 4: The number of decision variables in the paper is 12n2 + 13n, which depends on the system order n. When n increase, it will consume longer central processing unit time when solving stability conditions. Furthermore, it is worth noting that the number of decision variables in [41] and [42] is 70n2 + 12n and 142n2 + 18n, which is bigger than that in the paper. This means that the solution of the condition in [41] and [42] will waste more time to obtain less conservative results. It illustrates that the proposed methodology has lower computational complexity than [41] and [42].
Furthermore, the controller design algorithm is introduced via the following theorem.
Theorem 2: The system (6) with w(t) = 0 is asymptotically stable, if there exist symmetric positive matrices \tilde P > 0, \; \tilde Q > 0, \; \tilde Z > 0, \tilde R > 0, \tilde U > 0, {\tilde X_{11}}, {\tilde X_{13}}, {\tilde X_{22}}, {\tilde X_{23}}, {\tilde X_{33}}, \tilde G, {\tilde M_i}, {\tilde N_i}\; (i = 1, 2, 3, 4) and scales d > 0, {\varepsilon _i}\; (i = 1, 2, 3, 4) , such that:
\Psi _1^{ij} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{\tilde P}}} + d\left( {{{{\mathit{\boldsymbol{\tilde X}}}}_{11}} + {\mathit{\boldsymbol{\tilde X}}}_{11}^T + {\mathit{\boldsymbol{\tilde Q}}}} \right)}&{d\left( {{{{\mathit{\boldsymbol{\tilde X}}}}_{11}} + {{{\mathit{\boldsymbol{\tilde X}}}}_{22}} - {\mathit{\boldsymbol{\tilde Q}}}} \right)}&{d{{{\mathit{\boldsymbol{\tilde X}}}}_{13}}} \\ *&{d\left( {{{{\mathit{\boldsymbol{\tilde X}}}}_{11}} + {\mathit{\boldsymbol{\tilde X}}}_{22}^T + {\mathit{\boldsymbol{\tilde Q}}}} \right)}&{d{{{\mathit{\boldsymbol{\tilde X}}}}_{23}}} \\ *&*&{d{{{\mathit{\boldsymbol{\tilde X}}}}_{33}}} \end{array}} \right] > 0 | (22) |
{\tilde \Psi _2} = \left[ {\begin{array}{*{20}{c}} {{{\tilde \Xi }_{11}} + {{\tilde \Omega }_{11}}}&{{{\tilde \Xi }_{12}} + {{\tilde \Omega }_{12}}}&{{{\tilde \Xi }_{13}} + {{\tilde \Omega }_{13}}}&{{{\tilde \Xi }_{14}} + d\tilde {\mathit{\boldsymbol{X}}}_{33}^T} \\ *&{{{\tilde \Xi }_{22}} + {{\tilde \Omega }_{22}}}&{{{\tilde \Xi }_{23}} + {{\tilde \Omega }_{23}}}&{{{\tilde \Xi }_{24}} + d{{\tilde {\mathit{\boldsymbol{X}}}}_{13}}} \\ *&*&{{{\tilde \Xi }_{33}} + {{\tilde \Omega }_{33}}}&{{{\tilde \Xi }_{34}}} \\ *&*&*&{{{\tilde \Xi }_{44}}} \end{array}} \right] < 0 | (23) |
{\Psi _3} = \left[ {\begin{array}{*{20}{c}} {{{\tilde \Pi }_{11}}}&{{{\tilde \Pi }_{12}}} \\ *&{{{\tilde \Pi }_{22}}} \end{array}} \right] < 0, | (24) |
where
\begin{array}{l} {{\tilde \Pi }_{11}} = \left[ {\begin{array}{*{20}{c}} {{{\tilde \Xi }_{11}} + {{\tilde \Omega }_{11}}}&{{{\tilde \Xi }_{12}} + d{\mathit{\boldsymbol{\tilde N}}}_2^T}&{{{\tilde \Xi }_{13}} + d{\mathit{\boldsymbol{\tilde N}}}_3^T}&{{{\tilde \Xi }_{14}} + d{\mathit{\boldsymbol{\tilde N}}}_4^T} \\ *&{{{\tilde \Xi }_{22}} + \frac{{{d^2}}}{4}{\mathit{\boldsymbol{\tilde Z}}}}&{{{\tilde \Xi }_{23}}}&{{{\tilde \Xi }_{24}}} \\ *&*&{{{\tilde \Xi }_{33}} + {{\tilde \Omega }_{33}}}&{{{\tilde \Xi }_{34}}} \\ *&*&*&{{{\tilde \Xi }_{44}}} \end{array}} \right], \hfill \\ {{\tilde \Pi }_{12}} = \left[ {\begin{array}{*{20}{c}} { - d{{{\mathit{\boldsymbol{\tilde M}}}}_2}}&{d{{{\mathit{\boldsymbol{\tilde N}}}}_2}} \\ { - d{{{\mathit{\boldsymbol{\tilde M}}}}_3}}&{d{{{\mathit{\boldsymbol{\tilde N}}}}_3}} \\ { - d\tilde R_{12}^T - d{{{\mathit{\boldsymbol{\tilde M}}}}_3}}&{d{{{\mathit{\boldsymbol{\tilde N}}}}_3}} \\ { - d{{{\mathit{\boldsymbol{\tilde M}}}}_4}}&{d{{{\mathit{\boldsymbol{\tilde N}}}}_4}} \end{array}} \right], \hfill \\ {{\tilde \Pi }_{22}} = \left[ {\begin{array}{*{20}{c}} { - d{\mathit{\boldsymbol{\tilde Q}}} - d{{{\mathit{\boldsymbol{\tilde R}}}}_{11}}}&0 \\ *&{ - 2{\mathit{\boldsymbol{\tilde Z}}}} \end{array}} \right]. \hfill \\ {{\tilde \Xi }_{11}} = - {{{\mathit{\boldsymbol{\tilde X}}}}_{11}} - {\mathit{\boldsymbol{\tilde X}}}_{11}^T + {\mathit{\boldsymbol{\tilde M}}}_1^T + {{{\mathit{\boldsymbol{\tilde M}}}}_1} - {\varepsilon _1}{\mathit{\boldsymbol{A}}}{{{\mathit{\boldsymbol{\tilde G}}}}^T} - {\varepsilon _1}{\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{A}}}^T} \hfill \\ {{\tilde \Xi }_{12}} = {\mathit{\boldsymbol{\tilde P}}} + {\mathit{\boldsymbol{\tilde M}}}_2^T + {\varepsilon _1}{\mathit{\boldsymbol{\tilde G}}} - {\varepsilon _2}{{\mathit{\boldsymbol{A}}}^T}{{{\mathit{\boldsymbol{\tilde G}}}}^T} \hfill \\ {{\tilde \Xi }_{13}} = {{{\mathit{\boldsymbol{\tilde X}}}}_{11}} + {{{\mathit{\boldsymbol{\tilde X}}}}_{22}} + {\mathit{\boldsymbol{\tilde M}}}_3^T - {{{\mathit{\boldsymbol{\tilde M}}}}_1} - {\varepsilon _1}{\mathit{\boldsymbol{B\tilde K}}} - {\varepsilon _3}{\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{A}}}^T} \hfill \\ {{\tilde \Xi }_{14}} = - {{{\mathit{\boldsymbol{\tilde X}}}}_{13}} + {\mathit{\boldsymbol{\tilde M}}}_4^T - {{{\mathit{\boldsymbol{\tilde N}}}}_1} - {\varepsilon _4}{\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{A}}}^T} \hfill \\ {{\tilde \Xi }_{22}} = {\varepsilon _2}{\mathit{\boldsymbol{\tilde G}}} + {\varepsilon _2}{{{\mathit{\boldsymbol{\tilde G}}}}^T} \hfill \\ {{\tilde \Xi }_{23}} = - {{{\mathit{\boldsymbol{\tilde M}}}}_2} - {\mathit{\boldsymbol{B\tilde K}}} + {\varepsilon _3}{\mathit{\boldsymbol{\tilde G}}} \hfill \\ {{\tilde \Xi }_{24}} = {\varepsilon _4}{{{\mathit{\boldsymbol{\tilde G}}}}^T} - {{{\mathit{\boldsymbol{\tilde N}}}}_2} \hfill \\ {{\tilde \Xi }_{33}} = - {{{\mathit{\boldsymbol{\tilde X}}}}_{22}} - {\mathit{\boldsymbol{\tilde X}}}_{22}^T - {\mathit{\boldsymbol{\tilde M}}}_3^T - {{{\mathit{\boldsymbol{\tilde M}}}}_3} - {\varepsilon _3}{\mathit{\boldsymbol{B\tilde K}}} - {\varepsilon _3}{{{\mathit{\boldsymbol{\tilde K}}}}^T}{{\mathit{\boldsymbol{B}}}^T} \hfill \\ {{\tilde \Xi }_{34}} = - {{{\mathit{\boldsymbol{\tilde X}}}}_{23}} - {{{\mathit{\boldsymbol{\tilde U}}}}_{12}} - {\mathit{\boldsymbol{\tilde M}}}_4^T - {{{\mathit{\boldsymbol{\tilde N}}}}_3} - {\varepsilon _4}{{{\mathit{\boldsymbol{\tilde K}}}}^T}{{\mathit{\boldsymbol{B}}}^T} \hfill \\ {{\tilde \Xi }_{44}} = - {{{\mathit{\boldsymbol{\tilde X}}}}_{33}} - \frac{1}{d}{{{\mathit{\boldsymbol{\tilde U}}}}_{22}} - {\mathit{\boldsymbol{\tilde N}}}_4^T - {{{\mathit{\boldsymbol{\tilde N}}}}_4} \hfill \\ {{\tilde \Omega }_{11}} = d\left( {{{{\mathit{\boldsymbol{\tilde X}}}}_{13}} + {\mathit{\boldsymbol{\tilde X}}}_{13}^T} \right) + d{{{\mathit{\boldsymbol{\tilde U}}}}_{22}} - {\mathit{\boldsymbol{\tilde Z}}} \hfill \\ {{\tilde \Omega }_{12}} = d{{{\mathit{\boldsymbol{\tilde X}}}}_{11}} + d{\mathit{\boldsymbol{\tilde X}}}_{11}^T \hfill \\ {{\tilde \Omega }_{13}} = d{\mathit{\boldsymbol{\tilde X}}}_{23}^T + d{\mathit{\boldsymbol{\tilde U}}}_{12}^T + {\mathit{\boldsymbol{\tilde Z}}} \hfill \\ {{\tilde \Omega }_{22}} = d{{{\mathit{\boldsymbol{\tilde R}}}}_{11}} + d{\mathit{\boldsymbol{\tilde Q}}} + \frac{{{d^2}}}{4}{\mathit{\boldsymbol{\tilde Z}}} \hfill \\ {{\tilde \Omega }_{23}} = - d\left( {{{{\mathit{\boldsymbol{\tilde X}}}}_{11}} + {{{\mathit{\boldsymbol{\tilde X}}}}_{22}}} \right) + d{{{\mathit{\boldsymbol{\tilde R}}}}_{12}} \hfill \end{array} |
{\tilde \Omega _{33}} = d{{\mathit{\boldsymbol{\tilde R}}}_{22}} + d{{\mathit{\boldsymbol{\tilde U}}}_{11}} - {\mathit{\boldsymbol{\tilde Z}}} |
The controller is derived that
K = \tilde K{\tilde G^{ - T}} | (25) |
Proof: Denoting
\begin{array}{l} {\mathit{\boldsymbol{\tilde G}}} = {{\mathit{\boldsymbol{G}}}^{ - 1}}, \;{\mathit{\boldsymbol{\tilde K}}} = {\mathit{\boldsymbol{K}}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, \;{\mathit{\boldsymbol{\tilde P}}} = {\mathit{\boldsymbol{\tilde GP}}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, \;{{{\mathit{\boldsymbol{\tilde X}}}}_{11}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{X}}}_{11}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde X}}}}_{22}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{X}}}_{22}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde X}}}}_{13}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{X}}}_{13}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, \hfill \\ {{{\mathit{\boldsymbol{\tilde X}}}}_{23}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{X}}}_{23}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde X}}}}_{33}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{X}}}_{33}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde R}}}}_{11}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{R}}}_{11}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde R}}}}_{12}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{R}}}_{12}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde R}}}}_{22}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{R}}}_{22}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, \hfill \\ {\mathit{\boldsymbol{\tilde Q}}} = {\mathit{\boldsymbol{\tilde GQ}}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde U}}}}_{11}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{U}}}_{11}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde U}}}}_{12}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{U}}}_{12}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {{{\mathit{\boldsymbol{\tilde U}}}}_{22}} = {\mathit{\boldsymbol{\tilde G}}}{{\mathit{\boldsymbol{U}}}_{22}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, {\mathit{\boldsymbol{\tilde Z}}} = {\mathit{\boldsymbol{\tilde GZ}}}{{{\mathit{\boldsymbol{\tilde G}}}}^T}, \hfill \\ {\Upsilon _1} = diag\left\{ {{\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}} \right\}, {\Upsilon _2} = diag\left\{ {{\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}} \right\}, {\Upsilon _3} = diag\left\{ {{\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}, {\mathit{\boldsymbol{\tilde G}}}} \right\}, \hfill \\ {\mathit{\boldsymbol{\tilde M}}} = {\Upsilon _2}{\mathit{\boldsymbol{M}}}\Upsilon _2^T, {\mathit{\boldsymbol{\tilde N}}} = {\Upsilon _2}{\mathit{\boldsymbol{N}}}\Upsilon _2^T. \hfill \end{array} |
Pre- and post-multiplying (9)–(11) by {\Upsilon _1}, {\Upsilon _2}, {\Upsilon _3} and \Upsilon _1^T, \Upsilon _2^T, \Upsilon _3^T respectively, (22)–(24) are obtained. This completed the proof.
In this section, to vertify the performance of the given method, a simulation example for a LC is given [43]. The LC model parameters are considered in Table 1 as follow.
Items | Symbol | Values |
Length between perpendiculars | {L_{pp}} | 220.3 m |
Breadth at waterline | {B_{wl}} | 32.2 m |
Draft amidships | {T_m} | 7.2 m |
Displacement | \nabla | 33,229 t |
Vert. location of canter of gravity | KG | 15.1 m |
Metacentric height | GM | 2.75 m |
Transverse gyradius | {k_{xx}} | 12.9 m |
Longitudinal gyradius | {k_{yy}} | 55.1 m |
Vertical gyradius | {k_{zz}} | 55.1 m |
And the M and D in model (1) are
{\mathit{\boldsymbol{M}}} = \left[ {\begin{array}{*{20}{c}} {0.754}&0&0 \\ 0&{1.199}&{0.211} \\ 0&{0.029}&{0.524} \end{array}} \right], {\mathit{\boldsymbol{D}}} = \left[ {\begin{array}{*{20}{c}} {0.014}&0&0 \\ 0&{0.102}&{ - 0.0024} \\ 0&{0.192}&{0.095} \end{array}} \right] |
Noted that {\mathit{\boldsymbol{A}}} = \left[{\begin{array}{*{20}{c}} 0 & {\mathit{\boldsymbol{I}}} \\ 0 & { - {{\mathit{\boldsymbol{M}}}^{ - 1}}{\mathit{\boldsymbol{D}}}} \end{array}} \right], \quad {\mathit{\boldsymbol{B}}} = \left[{\begin{array}{*{20}{c}} 0 \\ {{{\mathit{\boldsymbol{M}}}^{ - 1}}} \end{array}} \right], \quad {{\mathit{\boldsymbol{B}}}_w} = \left[{\begin{array}{*{20}{c}} 0 \\ {{{\mathit{\boldsymbol{M}}}^{ - 1}}} \end{array}} \right], then
\begin{array}{l} {\mathit{\boldsymbol{A}}} = \left[ {\begin{array}{*{20}{c}} 0&0&0&1&0&0 \\ 0&0&0&0&1&0 \\ 0&0&0&0&0&1 \\ 0&0&0&{ - 0.0186}&0&0 \\ 0&0&0&0&{ - 0.0208}&{0.0342} \\ 0&0&0&0&{ - 0.3653}&{ - 0.1832} \end{array}} \right], \hfill \\ {\mathit{\boldsymbol{B}}} = {{\mathit{\boldsymbol{B}}}_w} = \left[ {\begin{array}{*{20}{c}} 0&0&0 \\ 0&0&0 \\ 0&0&0 \\ {1.3263}&0&0 \\ 0&{0.8422}&{ - 0.3391} \\ 0&{ - 0.0466}&{1.9272} \end{array}} \right]. \hfill \end{array} |
First, from Table 2, the maximum sampling period obtained by theorem 2 is d = 0.652 , which improves [37,38,39] above 160.8,146.97 and 22.56% respectively. One can find that the controller in this paper can obtain longer sampling period than that in [37,38,39].
The initial value of system state is assumed as {x_s}(t) = \left[{\begin{array}{*{20}{c}} {10} & {10} & {0.1} & 0 & 0 & 0 \end{array}} \right] , and d = 3.5513, {\varepsilon _1} = {\varepsilon _2} = {\varepsilon _3} = {\varepsilon _4} = 1, The desired value of system state is {x_r}(t) = \left[{\begin{array}{*{20}{c}} 0 & 0 & 0 & 0 & 0 & 0 \end{array}} \right]. Then, the gain can be computed that
{\mathit{\boldsymbol{K}}} = \left[ {\begin{array}{*{20}{c}} { - 9.6426}&0&0&{ - 20.1006}&0&0 \\ 0&{ - 15.3203}&{ - 0.5965}&0&{ - 32.0584}&{ - 1.4204} \\ 0&{1.1998}&{ - 5.7224}&{ - 0.0022}&{0.8754}&{ - 12.7550} \end{array}} \right]. |
The external environment disturbance is considered that
w(t) = \left[ {\sin (0.3t), \sin (0.2t), \sin (0.1t)} \right] |
Compared the results with the reference [39] which used the normal LKF to deal with the sampled-data issue of the DPS, the simulation results are shown in Figures 5–10 respectively. Figures 5–8 are the responses of the ship's positions and yaw angle. It can be seen that the stable time required to achieve the expected positions and yaw angle using normal LKF in [39] is 20 seconds, while the stable time based on improved LKF in the paper is 12 seconds, which improves the [39] about 40%.
Figures 8–10 are the responses of the velocities. It is shown that the stable time required to achieve the expected velocities using normal LKF in [39] is 22 seconds, while the stable time based on improved LKF in the paper is 15 seconds, which improves the [39] about 31.81%.
From Figures 5–10, we can see that in the presence of external environmental disturbances, the sampled-data controller designed in the paper can track the desired goal of the ship more smoothly and quickly. Besides, the fluctuations of the system are fewer. Therefore, the proposed sampled-data controller is more effective and robustness.
This article focused on the sampled-data issue for the LC with DPS. By constructing an improved LKF, the information about sampling modes is fully captured. Then, the conditions of asymptotical stability are present by means of LMI. Additionally, the sampled-data controller for LC with DPS is designed by analyzing the stabilization conditions. A numerical example is provided to illustrate that the proposed methods are effective and less conservativeness can be achieved. It is worth mentioning that the ship motion model only considers the linear problems, which may result in approximate errors in the actual object and reduce the robustness performance. In future research, the nonlinear systems and the modeling errors will be considered in the control synthesis process. In addition, other practical situations such as input saturation and incomplete state information will also be studied.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This study was funded by Undergraduate Education and Teaching Reform Research Project of Fujian Province in 2022 (FBJG20220194); Graduate Education and Teaching Reform Research project of Jimei University in 2021 (YJG2105); National Natural Science Foundation of China (51879119); Youth Innovation Foundation of Xiamen (3502Z20206019); The Natural Science Foundation of Fujian Province (2021J01822).
The authors declare there is no conflict of interest.
[1] | H. I. Freedman, Deterministic Mathematical Models in Population Ecology, Marcel Dekker Inc., New York, 1980. |
[2] | M. Kot, Elements of Mathematical Ecology, Cambridge University Press, 2001. https://doi.org/10.1017/CBO9780511608520 |
[3] |
M. A. Aziz-Alaoui, O. M. Daher, Boundedness and global stability for a predator-prey model with modified Leslie-Gower and Holling-Type Ⅱ schemes, Appl. Math. Lett., 16 (2003), 1069–1075. https://doi.org/10.1016/S0893-9659(03)90096-6 doi: 10.1016/S0893-9659(03)90096-6
![]() |
[4] |
M. Sen, M. Banerjee, A. Morozov, Bifurcation analysis of a ratio-dependent prey-predator model with the Allee effect, Ecol. Complexity, 11 (2012), 12–27. https://doi.org/10.1016/j.ecocom.2012.01.002 doi: 10.1016/j.ecocom.2012.01.002
![]() |
[5] |
J. C. Huang, D. M. Xiao, Analyses of bifurcations and stability in a predator-prey system with Holling Type-Ⅳ functional response, Acta Math. Appl. Sinica, 20 (2004), 167–178. https://doi.org/10.1007/s10255-004-0159-x doi: 10.1007/s10255-004-0159-x
![]() |
[6] |
J. Wang, Y. Cai, S. Fu, W. Wang, The effect of the fear factor on the dynamics of a predator-prey model incorporating the prey refuge, Chaos, 29 (2019), 83109. https://doi.org/10.1063/1.5111121 doi: 10.1063/1.5111121
![]() |
[7] |
T. Qiao, Y. Cai, S. Fu, W. Wang, Stability and Hopf bifurcation in a predator-prey model with the cost of anti-predator behaviors, Int. J. Bifurcation Chaos, 29 (2019), 1950185. https://doi.org/10.1142/S0218127419501852 doi: 10.1142/S0218127419501852
![]() |
[8] |
H. Zhang, Y. Cai, S. Fu, W. Wang, Impact of the fear effect in a prey-predator model incorporating a prey refuge, Appl. Math. Comput., 356 (2019), 328–337. https://doi.org/10.1016/j.amc.2019.03.034 doi: 10.1016/j.amc.2019.03.034
![]() |
[9] |
R. M. May, Simple mathematical models with very complicated dynamics, Nature, 261 (1976), 459–467. https://doi.org/10.1038/261459a0 doi: 10.1038/261459a0
![]() |
[10] |
L. Yuan, Q. Yang, Bifurcation, invariant curve and hybrid control in a discrete-time predator-prey system, Appl. Math. Modell., 39 (2015), 2345–2362. https://doi.org/10.1016/j.apm.2014.10.040 doi: 10.1016/j.apm.2014.10.040
![]() |
[11] |
M. Zhao, C. Li, J. Wang, Complex dynamic behaviors of a discrete-time predator-prey system, J. Appl. Anal. Comput., 7 (2017), 478–500. https://doi.org/10.11948/2017030 doi: 10.11948/2017030
![]() |
[12] | S. M. Rana, Bifurcation and complex dynamics of a discrete-time predator-prey system, Comput. Ecol. Software, 5 (2015), 187–200. |
[13] |
A. Q. Khan, K. Tanzeela, Neimark-Sacker bifurcation and hybrid control in a discrete-time Lotka-Volterra model, Math. Methods Appl. Sci., 43 (2020), 5887–5904. https://doi.org/10.1002/mma.6331 doi: 10.1002/mma.6331
![]() |
[14] |
H. N. Agiza, E. M. ELabbasy, H. EL-Metwally, A. A. Elsadany, Chaotic dynamics of a discrete prey-predator model with Holling-Type Ⅱ, Nonlinear Anal. Real World Appl., 10 (2019), 116–129. https://doi.org/10.1016/j.nonrwa.2007.08.029 doi: 10.1016/j.nonrwa.2007.08.029
![]() |
[15] |
G. Feng, X. Song, Bifurcation and chaos of a discrete-time population model, Discrete Dyn. Nat. Soc., 2020 (2020), 1–7. https://doi.org/10.1155/2020/8474715 doi: 10.1155/2020/8474715
![]() |
[16] |
A. G. M. Selvam, S. B. Jacob, R Dhineshbabu, Bifurcation and chaos control for prey predator model with step size in discrete time, J. Phys. Conf. Ser., 1543 (2020), 12010. https://doi.org/10.1088/1742-6596/1543/1/012010 doi: 10.1088/1742-6596/1543/1/012010
![]() |
[17] |
E. M. Elabbasy, A. A. Elsadany, Y. Zhang, Bifurcation analysis and chaos in a discrete reduced Lorenz system, Appl. Math. Comput., 228 (2014), 184–194. https://doi.org/10.1016/j.amc.2013.11.088 doi: 10.1016/j.amc.2013.11.088
![]() |
[18] |
L. F. Cheng, H. J. Cao, Bifurcation analysis of a discrete-time ratio-dependent predator-prey Model with Allee effect, Commun. Nonlinear Sci. Numer. Simul., 38 (2016), 288–302. https://doi.org/10.1016/j.cnsns.2016.02.038 doi: 10.1016/j.cnsns.2016.02.038
![]() |
[19] |
Z. M. He, X. Lai, Bifurcation and chaotic behavior of a discrete-time predator-prey system, Nonlinear Anal. Real World Appl., 12 (2011), 403–417. https://doi.org/10.1016/j.nonrwa.2010.06.026 doi: 10.1016/j.nonrwa.2010.06.026
![]() |
[20] |
X. Wang, L. Zanette, X. Zou, Modelling the fear effect in predator-prey interactions, J. Math. Biol., 73 (2016), 1179–1204. https://doi.org/10.1007/s00285-016-0989-1 doi: 10.1007/s00285-016-0989-1
![]() |
[21] | R.M. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, 1973. |
[22] |
X. Li, W. Jiang, J. Shi, Hopf bifurcation and turing instability in the reaction-diffusion Holling-Tanner predator-prey model, IMA J. Appl. Math., 78 (2013), 287–306. https://doi.org/10.1093/imamat/hxr050 doi: 10.1093/imamat/hxr050
![]() |
[23] | J. Guckenheimer, P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Springer, New York, 1983. https://doi.org/10.1007/978-1-4612-1140-2 |
[24] | Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, Springer-Verlag, New York, 1998. https://doi.org/10.1007/b98848 |
[25] | S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos, Springer-Verlag, New York, 1990. https://doi.org/10.1007/978-1-4757-4067-7 |
[26] | G. Chen, X. Dong, From Chaos to Order: Perspectives, Methodologies, and Applications, World Scientific, Singapore, 1998. https://doi.org/10.1142/3033 |
[27] | S. Elaydi, An Introduction to Difference Equations, 3rd edition, Springer-Verlag, New York, 2005. https://doi.org/10.1007/0-387-27602-5 |
[28] | S. Lynch, Dynamical Systems with Applications Using Mathematica, 1st edition, Birkhauser, Boston, 2007. https://doi.org/10.1007/978-0-8176-4586-1 |
1. | Minjie Zheng, Yulai Su, Changjian Yan, Further Stability Criteria for Sampled-Data-Based Dynamic Positioning Ships Using Takagi–Sugeno Fuzzy Models, 2024, 16, 2073-8994, 108, 10.3390/sym16010108 | |
2. | Minjie Zheng, Yulai Su, Guoquan Chen, An improved sampled-data control for a nonlinear dynamic positioning ship with Takagi-Sugeno fuzzy model, 2024, 21, 1551-0018, 6019, 10.3934/mbe.2024265 | |
3. | Zihe Qin, Feng Zhang, Wenlin Xu, Yu Chen, Jinyu Lei, A Model-Free Adaptive Positioning Control Method for Underactuated Unmanned Surface Vessels in Unknown Ocean Currents, 2024, 12, 2077-1312, 1801, 10.3390/jmse12101801 |
Items | Symbol | Values |
Length between perpendiculars | {L_{pp}} | 220.3 m |
Breadth at waterline | {B_{wl}} | 32.2 m |
Draft amidships | {T_m} | 7.2 m |
Displacement | \nabla | 33,229 t |
Vert. location of canter of gravity | KG | 15.1 m |
Metacentric height | GM | 2.75 m |
Transverse gyradius | {k_{xx}} | 12.9 m |
Longitudinal gyradius | {k_{yy}} | 55.1 m |
Vertical gyradius | {k_{zz}} | 55.1 m |
Items | Symbol | Values |
Length between perpendiculars | {L_{pp}} | 220.3 m |
Breadth at waterline | {B_{wl}} | 32.2 m |
Draft amidships | {T_m} | 7.2 m |
Displacement | \nabla | 33,229 t |
Vert. location of canter of gravity | KG | 15.1 m |
Metacentric height | GM | 2.75 m |
Transverse gyradius | {k_{xx}} | 12.9 m |
Longitudinal gyradius | {k_{yy}} | 55.1 m |
Vertical gyradius | {k_{zz}} | 55.1 m |
Method | [37] | [38] | [39] | Theorem1 |
d_2 | 0.25 | 0.264 | 0.532 | 0.652 |