Research article

Extinction and stationary distribution of stochastic predator-prey model with group defense behavior


  • Received: 10 June 2022 Revised: 17 August 2022 Accepted: 22 August 2022 Published: 06 September 2022
  • Considering that many prey populations in nature have group defense behavior, and the relationship between predator and prey is usually affected by environmental noise, a stochastic predator-prey model with group defense behavior is established in this paper. Some dynamical properties of the model, including the existence and uniqueness of global positive solution, sufficient conditions for extinction and unique ergodic stationary distribution, are investigated by using qualitative theory of stochastic differential equations, Lyapunov function analysis method, Itô formula, etc. Furthermore, the effects of group defense behavior and environmental noise on population stability are also discussed. Finally, numerical simulations are carried out to illustrate that the effects of environmental noise on both populations are negative, the appropriate group defense level of prey can maintain the stability of the relationship between two populations, and the survival threshold is strongly influenced by the intrinsic growth rate of prey population and the intensity of environmental noise.

    Citation: Yansong Pei, Bing Liu, Haokun Qi. Extinction and stationary distribution of stochastic predator-prey model with group defense behavior[J]. Mathematical Biosciences and Engineering, 2022, 19(12): 13062-13078. doi: 10.3934/mbe.2022610

    Related Papers:

    [1] Shuangjie Yuan, Jun Zhang, Yujia Lin, Lu Yang . Hybrid self-supervised monocular visual odometry system based on spatio-temporal features. Electronic Research Archive, 2024, 32(5): 3543-3568. doi: 10.3934/era.2024163
    [2] Yanling Dong, Xiaolan Zhou . Advancements in AI-driven multilingual comprehension for social robot interactions: An extensive review. Electronic Research Archive, 2023, 31(11): 6600-6633. doi: 10.3934/era.2023334
    [3] Yu Lei, Zhi Su, Chao Cheng . Virtual reality in human-robot interaction: Challenges and benefits. Electronic Research Archive, 2023, 31(5): 2374-2408. doi: 10.3934/era.2023121
    [4] Jiarui Chen, Aimin Tang, Guanfeng Zhou, Ling Lin, Guirong Jiang . Walking dynamics for an ascending stair biped robot with telescopic legs and impulse thrust. Electronic Research Archive, 2022, 30(11): 4108-4135. doi: 10.3934/era.2022208
    [5] Mingxing Xu, Hongyi Lin, Yang Liu . A deep learning approach for vehicle velocity prediction considering the influence factors of multiple lanes. Electronic Research Archive, 2023, 31(1): 401-420. doi: 10.3934/era.2023020
    [6] Ismail Ben Abdallah, Yassine Bouteraa, Saleh Mobayen, Omar Kahouli, Ali Aloui, Mouldi Ben Amara, Maher JEBALI . Fuzzy logic-based vehicle safety estimation using V2V communications and on-board embedded ROS-based architecture for safe traffic management system in hail city. Electronic Research Archive, 2023, 31(8): 5083-5103. doi: 10.3934/era.2023260
    [7] Yong Zhao, Shanshan Ren . Synchronization for a class of complex-valued memristor-based competitive neural networks(CMCNNs) with different time scales. Electronic Research Archive, 2021, 29(5): 3323-3340. doi: 10.3934/era.2021041
    [8] Lijing Ma, Shiru Qu, Lijun Song, Junxi Zhang, Jie Ren . Human-like car-following modeling based on online driving style recognition. Electronic Research Archive, 2023, 31(6): 3264-3290. doi: 10.3934/era.2023165
    [9] Jiping Xing, Xiaohong Jiang, Yu Yuan, Wei Liu . Incorporating mobile phone data-based travel mobility analysis of metro ridership in aboveground and underground layers. Electronic Research Archive, 2024, 32(7): 4472-4494. doi: 10.3934/era.2024202
    [10] Longkui Jiang, Yuru Wang, Weijia Li . Regress 3D human pose from 2D skeleton with kinematics knowledge. Electronic Research Archive, 2023, 31(3): 1485-1497. doi: 10.3934/era.2023075
  • Considering that many prey populations in nature have group defense behavior, and the relationship between predator and prey is usually affected by environmental noise, a stochastic predator-prey model with group defense behavior is established in this paper. Some dynamical properties of the model, including the existence and uniqueness of global positive solution, sufficient conditions for extinction and unique ergodic stationary distribution, are investigated by using qualitative theory of stochastic differential equations, Lyapunov function analysis method, Itô formula, etc. Furthermore, the effects of group defense behavior and environmental noise on population stability are also discussed. Finally, numerical simulations are carried out to illustrate that the effects of environmental noise on both populations are negative, the appropriate group defense level of prey can maintain the stability of the relationship between two populations, and the survival threshold is strongly influenced by the intrinsic growth rate of prey population and the intensity of environmental noise.



    Robotics is a critical component of industrial development and automation. It enables significant growth in manufacturing efficiency, productivity, precision and quality, while also increasing safety when performing tasks that would be dangerous for a human operator to perform.

    Automation solutions incorporating mobile robots in industry have grown in popularity in the last decades. This is due to an increasing demand for flexible solutions. Autonomous mobile robots are used in multiple areas, such as industrial automation, logistics and warehouse management, health care and services and many other industrial and non-industrial applications [1,2].

    The basics of mobile robotics can be divided into locomotion, perception, localization and planning and navigation. Locomotion is required for the robot to move in its working environment. It relies on suitable motion control schemes that generate the necessary control actions to the motors' drivers, based on set points and feedback data [1,2,3,4].

    There are many works focused on the analysis of kinematic models of differential drive vehicles [5,6,7,8,9]. Additionally, multiple studies show methods to test the accuracy of these models and how to modify them to reduce errors in odometry [10,11,12,13]. On the other hand, robot path planning and optimization either in static or dynamic environments represent key issues. Path planning strategies are often classified as (ⅰ) analytical, (ⅱ) enumerative, (ⅲ) meta-heuristic and (ⅳ) evolutionary [14]. In particular, the evolutionary methods can be regarded as optimization procedures, and many algorithms have been proposed [15,16,17,18,19,20,21,22].

    A quality control system is necessary to ensure that the robot can move with high accuracy, which is very important for industrial environments. Motion control systems for mobile robots with various drive types utilizing Beckhoff's TwinCAT automation software have been addressed in a small number of works. For instance, in [23] a new remote operated vehicle (ROV) control system based on TwinCAT was developed. The solution allowed the connection of functional modules and hardware and/or software modifications, while reducing the time needed for maintenance. A graphical user interface (GUI) was also created for real-time data monitoring and controlling the ROV. In [24] a heavy-duty omni-directional Mecanum robot was presented, together with its control system and simulation design. The control system combined a Beckhoff module and the Robot Operating System (ROS) for low-level motion execution and high-level intelligent navigation tasks, respectively. In addition, a virtual simulation environment was validated. In [25] real-time sensor fusion for smooth position feedback of a heavy-duty field robot was addressed. The robot autonomous driving system was based on a Beckhoff hardware platform running TwinCAT. The programming was performed in a MATLAB/Simulink environment, while the functional modules were built on a target TwinCAT system. The robot used wireless communication between the real-time hardware and the host PC. In [26] a solution for a network of real-time-capable modules was derived and illustrated on a prototype of a mobile robot using TwinCAT CNC. In [27] a novel wheel-legged robotic system was developed. The robot control software was implemented using TwinCAT. In [28] a path planning and navigation control system was developed for a 12 m length driverless electric bus. The TwinCAT software system was utilized for path planning and trajectory tracking controller, while ensuring real-time update and synchronization rate. However, to the best of the authors' knowledge, none of these TwinCAT implementations has addressed motion accuracy testing and odometry of the robot, nor motion correction in order to obtain minimum errors. On the other hand, the proposed approach is flexible, completely modular and extendable to accommodate different sensors and actuators, as well as other robotic differential drive platforms utilizing TwinCAT 3.

    The present paper proposes a motion control scheme for a low-cost differential drive mobile robot. The robot locomotion platform consists of a differential drive setup using two wheels powered by Beckhoff motors and a free caster wheel. The motion control is developed using Beckhoff's TwinCAT 3 automation software in an industrial PC (IPC). The system is tested and experimentally tuned to achieve excellent performance. Field tests of the robot performing linear and angular trajectories reveal errors below 0.02 and 0.03%, respectively, after tuning its motion parameters. The proposed approach can be extended and applied to other robotic differential drive platforms utilizing TwinCAT 3, being a cost-effective solution for real-world applications.

    The paper is organized as follows: Section 2 describes the robot and its hardware; Section 3 presents the development of the motion control system, the kinematics of the robot and the testing and correction of the odometry. Finally, Section 4 outlines the main conclusions.

    The mobile robot uses a differential drive platform with two independently driven wheels and one free caster wheel. Each of the driven wheels is equipped with an AS1050-0120 Beckhoff motor with an incorporated incremental encoder. The resolution of the encoders is 1.8° or 200 increments per revolution. The motors are coupled with a gearbox with a gear ratio of 5. The robotic platform is also equipped with an anthropomorphic robot arm, which is not addressed in this paper. A picture of the robot and the two independently driven wheels can be seen in Figure 1, along with their respective Beckhoff stepper motors and gearboxes [29].

    Figure 1.  Front picture of the robot.

    The electrical cabinet can be seen in Figure 2. It contains three different power supplies, one with 48 V, one with 24 V and one with 5 V:

    Figure 2.  Electrical cabinet of the robot.

    1) The 48 V power supply is used to power 5 EtherCAT terminals (3 EL7047 and 2 EL7041) that connect to 3 stepper motors used in the robot arm and 2 stepper motors used for the robot's locomotion.

    2) The 24 V power supply is used to power the EtherCAT terminal EL9800, which supplies 2 EtherCAT EL7037 modules that connect to 2 other stepper motors used in the robot arm. It also powers their encoders, the robot's IPC, the EtherCAT coupler EK1100 and the input (EL1008) and output terminals (EL2008).

    3) The 5 V power supply is used for three of the encoders in the robot arm and for the two encoders used for locomotion.

    Figure 3 presents a schema of the electrical power used in the robot. The IPC used in this robot is the C6015-0010 model from Beckhoff, running with Windows 10 as its operating system and with TwinCAT 3 installed. The robot is also equipped with an emergency button that, when activated, turns off all circuits that send power to the motors, while keeping the command circuits on, including the IPC and the EtherCAT communication.

    Figure 3.  Electrical power schema of the robot.

    The software used for the control and command of the robotic system is TwinCAT 3. TwinCAT 3 is an automation software solution for PC-based control that can turn a PC with EtherCAT communication into a real time controller with multiple subsystems running simultaneously. It consists of a XAR (eXtended Automation Runtime) that runs the already programmed control modules and the XAE (eXtended Automation Engineering) that, using Microsoft Visual Studio, allows the user to setup, modify and program the control system to be used.

    TwinCAT 3 uses a number of modules that communicate with one another to form a complete controller. There are modules for the system configuration (SYSTEM), motion control (MOTION), PLC programming (PLC), safety (SAFETY), C++ programming (C++), analytics (ANALYTICS) and I/O configuration (I/O). This paper does not make use of the SAFETY, C++ or ANALYTICS modules.

    When configuring the system, the software in the I/O module recognized all connected EtherCAT devices. Furthermore, the motor drivers were configured, which included entering multiple motor properties, such as maximum current, nominal voltage, number of steps per revolution and encoder resolution. The different movement axes were configured in the MOTION module. This module is responsible for the numerical control of the movement axes. It can use point-to-point control to generate the signals needed to send to the appropriate I/O device. It has several tabs where different fields must be configured correctly. The PLC module implements the system logic and is responsible for the overall control of the robot's motion system. This module is programmable according to the IEC 61131-3 international standard for programmable PLCs, which supports the following programming languages: ladder diagram (LD), function block diagram (FBD), sequential function chart (SFC), structured text (ST) and instruction list (IL). In the case of this paper, every program in the PLC module was developed in ST.

    The robot used in this paper was also simultaneously used for the development of an autonomous navigation system using ROS. As such, the objective was for TwinCAT to receive velocity commands from the ROS system while simultaneously sending odometry information to the ROS environment. This information flow is shown in Figure 4.

    Figure 4.  Flow of information between TwinCAT 3 and ROS.

    The ROS Navigation Stack uses the move_base package. This package generates the following commands: linear velocity in x (Vx), linear velocity in y (Vy) and angular velocity in θ (Vθ). These directions are relative to the robot's current position, in accordance with Figure 5. This specific robot can only have velocities in the x and θ directions since it is a differential drive robot. Because there are no other possible velocities that a differential drive robot can take, Vx and Vθ are also respectively known as linear and angular velocity.

    Figure 5.  Robot's movement axes.

    Because TwinCAT 3 controls the movement of both wheel axes independently, the robot's linear and angular velocities must be converted into velocities for the left and right wheels. Furthermore, odometry information about the robot must be obtained using information from the encoders.

    TwinCAT's MOTION module was configured with information about the motors and transmission systems used. As a result, it can use wheel velocities specified in mm/s and automatically read information from encoders in mm. As such, it is only necessary to transform Vx and Vθ into VR (velocity of right wheel) and VL (velocity of left wheel), and 𝛿R and 𝛿L into 𝛿x, 𝛿y and 𝛿θ, where 𝛿 means the variation of position between each control cycle.

    In differential drives, linear and angular velocities are given as functions of the velocities of the wheels, according to [5]:

    Vx=VR+VL2 (1)
    Vθ=VRVLB (2)

    where B is the axial distance between the driven wheels. Adding and rearranging Eqs (1) and (2) results in

    VR=Vx+VθB2 (3)
    VL=(2Vx)VR (4)

    Since the increments in position of each wheel can be obtained in mm from the encoder data in TwinCAT, linear and angular displacements of the robot are given by

    δx=δR+δL2 (5)
    δθ=δRδLB (6)

    The length of the red arc in Figure 6 is given by 𝛿x (the size is exaggerated for purposes of easier understanding). The linear distance, d, between the two points (green line) can be obtained with the trigonometric relations:

    d=2Rsin(δθ2) (7)
    R=δxδθ (8)
    Figure 6.  Arc traveled by the robot between two points.

    By adding the previous two equations, Eq (9) is obtained:

    d=2δxδθsin(δθ2) (9)

    This equation, however, can only be used if the robot did not move in a straight line between two cycles of the program. If 𝛿θ is 0, then Eq (10) is used instead:

    d=δx (10)

    The position of the robot in relation to its starting frame can be obtained over time by adding the linear and angular displacements between two cycles of the program (k + 1 and k) to the last calculated position of the robot, using

    θk+1=θk+δθ (11)
    xk+1=xk+δxcos(θk+θk+12) (12)
    yk+1=yk+δxsin(θk+θk+12) (13)

    The angle used to calculate x and y is the average value of the angular position in each of the two cycles to reduce errors. Initial position is initialized to (0, 0, 0). Since usually these cycles are calculated at very high frequencies, the distance between two points in an arc given by Eq (9) can be approximated by the length of the arc itself with negligible error. Lower cycle frequency and higher linear and angular velocities result in a higher error when using this approximation. Herein, the cycle time of the program is 10 ms. At the robot's maximum linear velocity (Vx) of 0.3 m/s and maximum angular velocity (Vx) of 0.7 rad/s, the length of the arc traveled between each cycle (𝛿x) is 3 mm, and the angular displacement (𝛿θ) is 0.007 rad. Using these values in Eq (9) results in a linear distance traveled (d) of approximately 2.999994 mm. This means that, in this case, the maximum error obtainable from the use of this approximation is around 0.0002%.

    The kinematics was implemented in a program developed in ST language in TwinCAT. It was necessary to take the velocity commands and make the wheels move at the respective wheel velocities until a new velocity command is obtained. To achieve this, the MC_MoveVelocity function block was used. This function block starts a continuous movement with specified velocity and direction. However, this function block cannot receive a null velocity as an input, which means that an additional function block must be used for when the velocity of any of the wheels is zero. As such, the MC_Halt function block was used for these cases, as it stops an axis's movement with a defined breaking ramp.

    Since the MC_MoveVelocity function block only takes velocity inputs larger than 0 and a forward or negative direction value, IF conditions are written for all nine possible combinations of positive, null, and negative velocity values of each wheel. The velocity command variables Vx and Vθ and the odometry variables x, y and θ are all set as global variables to allow them to be written and read by ROS. The program also reads the current velocity of each wheel and calculates the linear and angular velocities the robot is moving at currently, as these are necessary odometry data for ROS. The odometry is calculated using the difference between the absolute values of the encoders in the current and last cycles. During the first cycle of the program, odometry values are set to (0, 0, 0), the motors are powered, and their axes are reset.

    Because the MC_MoveVelocity and MC_Halt function blocks require a rising edge signal of the Execute Boolean input, it is also necessary to add the following lines present in Figure. These lines feed the mentioned function blocks an input of "Execute = FALSE". This allows the robot's movement to be flexible and able to change velocities at any moment.

    Figure 7.  Necessary lines to allow new movement commands.

    In order to minimize the trajectory errors, the values of the robot parameters need to be properly determined. This problem is challenging, but a deep analysis is out of the scope of this paper. As such, herein, an empirical tuning approach was adopted based on experimental testing and error analysis, by comparing the actual odometry trajectories with those recorded by an external measurement system. Naturally, this method is simple, but it does not guarantee optimal robot behavior. Different methods, namely, parameters' optimization based on evolutionary algorithms, can be an effective means to address the issue [30,31].

    The quality, and magnitude of the errors, of the implemented odometry were tested by comparing the odometry values with the real trajectory of the robot. As such, a pen was attached to the robotic platform at the robot's rotation axis. A component to attach the pen to the robot in the position of its frame was modeled and 3D printed (Figure 8). Its geometry only allows the pen to move in the vertical axis, ensuring that it can stay in contact with the ground while the robot moves. A spring was inserted inside to guarantee that there is tension between the pen's ballpoint and the ground.

    Figure 8.  Pen gripper installed on the robot.

    To make the robot perform specific movements that can be repeated, another TwinCAT program was also written in ST. This program includes the odometry functionalities already implemented, but the robot movements are specified for a set distance and relative to the position at the start of the movement. To achieve this, the program utilizes the MC_MoveRelative function block. It can be used to perform linear, circular or rotational trajectories.

    The quality of the linear odometry was tested by measuring the length of the line drawn when instructing the robot to move 1 m forward (Figure 9). This test was done at a linear velocity of 0.025 m/s and repeated three times. The results of the linear odometry tests are shown in Table 1.

    Figure 9.  Picture of the lines produced by the robot and ruler used to measure them.
    Table 1.  Results of the initial linear odometry tests.
    x given by the odometry Length of line Error %Error
    Test 1 1000.0 mm 1009.5 mm -9.5 mm -0.94%
    Test 2 1000.0 mm 1011.0 mm -11.0 mm -1.09%
    Test 3 1000.0 mm 1010.5 mm -10.5 mm -1.04%
    Average 1000.0 mm 1010.3 mm -10.3 mm -1.02%

     | Show Table
    DownLoad: CSV

    The results of the linear odometry test show that the robot's odometry was giving values lower than the measured values by around 1.03%. To counteract this error, the scaling factor parameter of the motor's encoder was changed in TwinCAT. This scaling factor specifies how much each wheel moves in mm per increment of the encoder. The scaling factor is given by:

    SF=πD4ei (14)

    where SF is the scaling factor, D is the diameter of the wheel, i is the gear ratio of the transmission system, and e is the number of encoder increments per rotation.

    To correct the linear odometry of the robot, it is necessary to divide the scaling factor by (1 + %Error):

    SFcorrected=SF1+%Error (15)

    Since e and i are constant, any changes to the SF can also be interpreted as changes to the wheel diameter considered. The original SF used in this robot was 0.015340 mm/Inc (mm per increment of the encoder). This was obtained using a wheel diameter of 100 mm in Eq (14). In this case, using Eq (15) a corrected SF of 0.015498 mm/Inc is obtained. This equates to considering a wheel diameter of around 101.0 mm. The previous linear odometry tests were repeated using the corrected scaling factor. The results are shown in Table. We verify now an average error of around 0.02%, or 0.2 mm, which is within the resolution of the ruler used to measure the line length.

    Table 2.  Results of the linear odometry tests after correction.
    x given by the odometry Length of line Error %Error
    Test 1 1000.0 mm 1000.0 mm 0.0 mm 0.0%
    Test 2 1000.0 mm 1000.0 mm 0.0 mm 0.0%
    Test 3 1000.0 mm 1000.5 mm -0.5 mm -0.05%
    Average 1000.0 mm 1000.2 mm -0.2 mm -0.02%

     | Show Table
    DownLoad: CSV

    Angular odometry tests were performed after the correction of the linear odometry with the new scaling factor. The robot was programmed to perform an in-place rotation of 360°. A line parallel to the wheel's initial and final positions was drawn on a piece of paper on the floor. The angle between these two lines was measured and compared to the angle given by the robot's odometry. This test was done at an angular velocity of 0.182 rad/s and repeated three times. The results are summarized in Table.

    Table 3.  Results of the first series of angular odometry tests.
    θ given by the odometry Angle measured Error %Error
    Test 1 360.0° 363.1° -3.1° -0.85%
    Test 2 360.0° 362.9° -2.9° -0.80%
    Test 3 360.0° 361.8° -1.8° -0.50%
    Average 360.0° 362.6° -2.6° -0.72%

     | Show Table
    DownLoad: CSV

    Angular displacement is given by Eq (6). Considering that the linear odometry is already corrected, angular odometry can only be corrected by changing the value of the distance between wheels (B), which is inversely proportional to the angular displacement. It can be corrected by:

    Bcorrected=B(1+%Error) (16)

    In this case, the measuring error was considered too large in comparison with the odometry error to take any conclusions. Because of this, a new series of tests was conducted by commanding the robot to perform a 1080° in-place rotation with the same angular velocity of 0.182 rad/s. This larger movement results in an increased odometry error, while maintaining similar measurement errors, which allows for a more accurate percentual error to be calculated. Results of these tests are shown in Table 4.

    Table 4.  Results of the second series of angular odometry tests.
    θ given by the odometry Angle measured Error %Error
    Test 1 1080.0° 1089.1° -9.1° -0.84%
    Test 2 1080.0° 1090.0° -10.0° -0.92%
    Test 3 1080.0° 1088.5° -8.5° -0.78%
    Average 1080.0° 1089.2° -9.2° -0.84%

     | Show Table
    DownLoad: CSV

    By applying the average error obtained from these tests to Eq (16) using the original value of B of 550 mm, a corrected distance between wheels of 545.38 mm is determined. The same angular odometry tests were repeated using the corrected value of B. The results are shown in Table 5.

    Table 5.  Results of the angular odometry tests after correction.
    θ given by the odometry Angle measured Error %Error
    Test 1 1080.0° 1079.2° +0.8° +0.07%
    Test 2 1080.0° 1079.8° +0.2° +0.02%
    Test 3 1080.0° 1080.1° -0.1° -0.01%
    Average 1080.0° 1079.7° +0.3° +0.03%

     | Show Table
    DownLoad: CSV

    In a 1080° trajectory, the results now indicate an average error of 0.03°, or approximately 0.03%, which is within the repeatability errors of the tests. As a result, the value of 545.38 mm for the corrected distance between wheels is accepted.

    When using different values for the velocities, the odometry errors appeared to be similar due to low wheel slippage occurring in the testing environment.

    A motion control system for a low-cost differential drive mobile robot was proposed, based on Beckhoff's TwinCAT 3 automation software running on an IPC. An empirical tuning approach based on experimental testing and error analysis, by comparing the odometry trajectories with those measured with an external measurement system, was implemented. The approach led to good robot performance, with errors below 0.02 and 0.03% while performing linear and angular trajectories, respectively. The method is simple, but it does not guarantee optimal robot behavior. Further work will address parameter optimization, namely, the adoption of evolutionary algorithms, and additional tests to further assess the accuracy of the system, including verifying the accuracy of each wheel individually and testing differences between clockwise and counterclockwise movement. The proposed approach is flexible and modular, and it can be extended, adjusted and used in other robotic differential drive based solutions utilizing TwinCAT 3.

    The authors declare there is no conflict of interest.



    [1] S. L. Lima, L. M. Dill, Behavioral decisions made under the risk of predation: a review and prospectus, Can. J. Zool., 68 (1990), 619–640. https://doi.org/10.1139/z90-092 doi: 10.1139/z90-092
    [2] P. Cong, M. Fan, X. Zou, Dynamics of a three-species food chain model with fear effect, Commun. Nonlinear Sci. Numer. Simul., 99 (2021), 105809. https://doi.org/10.1016/j.cnsns.2021.105809 doi: 10.1016/j.cnsns.2021.105809
    [3] H. Qi, X. Meng, T. Hayat, A. Hobiny, Stationary distribution of a stochastic predator–prey model with hunting cooperation, Appl. Math. Lett., 124 (2022), 107662. https://doi.org/10.1016/j.aml.2021.107662 doi: 10.1016/j.aml.2021.107662
    [4] H. Qi, X. Meng, Threshold behavior of a stochastic predator-prey system with prey refuge and fear effect, Appl. Math. Lett., 113 (2021), 106846. https://doi.org/10.1016/j.aml.2020.106846 doi: 10.1016/j.aml.2020.106846
    [5] Y. Wang, X. Zou, On a predator–prey system with digestion delay and anti-predation strategy, J. Nonlinear Sci., 30 (2020), 1579–1605. https://doi.org/10.1007/s00332-020-09618-9 doi: 10.1007/s00332-020-09618-9
    [6] G. Tang, S. Tang, R. A. Cheke, Global analysis of a Holling type Ⅱ predator-prey model with a constant prey refuge, Nonlinear Dyn., 76 (2014), 635–647. https://doi.org/10.1007/s11071-013-1157-4 doi: 10.1007/s11071-013-1157-4
    [7] W. Cresswell, J. L. Quinn, Faced with a choice, sparrowhawks more often attack the more vulnerable prey group, Oikos, 104 (2004), 71–76. https://doi.org/10.1111/j.0030-1299.2004.12814.x doi: 10.1111/j.0030-1299.2004.12814.x
    [8] J. Wei, W. Shao, M. Cao, J. Ge, P. Yang, L. Chen, et al., Phenylacetonitrile in locusts facilitates an antipredator defense by acting as an olfactory aposematic signal and cyanide precursor, Sci. Adv., 5 (2019), eaav5495. https://doi.org/10.1126/sciadv.aav5495 doi: 10.1126/sciadv.aav5495
    [9] A. A. Salih, M. Baraibar, K. K. Mwangi, G. Artan, Climate change and locust outbreak in East Africa, Nat. Clim. Change, 10 (2020), 584–585. https://doi.org/10.1038/s41558-020-0835-8 doi: 10.1038/s41558-020-0835-8
    [10] C. N. Meynard, M. Lecoq, M. P. Chapuis, C. Piou, On the relative role of climate change and management in the current desert locust outbreak in East Africa, Global Change Biol., 26 (2020), 3753–3755. https://doi.org/10.1111/gcb.15137 doi: 10.1111/gcb.15137
    [11] H. I. Freedman, G. S. Wolkowicz, Predator-prey systems with group defence: the paradox of enrichment revisited, Bull. Math. Biol., 48 (1986), 493–508. https://doi.org/10.1016/S0092-8240(86)90004-2 doi: 10.1016/S0092-8240(86)90004-2
    [12] H. I. Freedman, S. Ruan, Hopf bifurcation in three-species food chain models with group defense, Math. Biosci., 111 (1992), 73–87. https://doi.org/10.1016/0025-5564(92)90079-C doi: 10.1016/0025-5564(92)90079-C
    [13] G. Gimmelli, B. W. Kooi, E. Venturino, Ecoepidemic models with prey group defense and feeding saturation, Ecol. Complexity, 22 (2015), 50–58. https://doi.org/10.1016/j.ecocom.2015.02.004 doi: 10.1016/j.ecocom.2015.02.004
    [14] C. Xu, S. Yuan, T. Zhang, Global dynamics of a predator-prey model with defense mechanism for prey, Appl. Math. Lett., 62 (2016), 42–48. https://doi.org/10.1016/j.aml.2016.06.013 doi: 10.1016/j.aml.2016.06.013
    [15] X. Cheng, J. Luo, Y. Zhao, Dynamic analysis of a population competition model with disease in one species and group defense in another species, Int. J. Bifurcation Chaos, 30 (2020), 2050181. https://doi.org/10.1142/S0218127420501813 doi: 10.1142/S0218127420501813
    [16] M. Das, G. P. Samanta, A prey-predator fractional order model with fear effect and group defense, Int. J. Dyn. Control, 9 (2021), 334–349. https://doi.org/10.1007/s40435-020-00626-x doi: 10.1007/s40435-020-00626-x
    [17] Y. Du, B. Niu, J. Wei, A predator-prey model with cooperative hunting in the predator and group defense in the prey, Discrete Contin. Dyn. Syst. B, (2021). https://doi.org/10.3934/dcdsb.2021298 doi: 10.3934/dcdsb.2021298
    [18] D. Xiao, S. Ruan, Codimension two bifurcations in a predator-prey system with group defense, Int. J. Bifurcation Chaos, 11 (2001), 2123–2131. https://doi.org/10.1142/S021812740100336X doi: 10.1142/S021812740100336X
    [19] P. H. Leslie, A stochastic model for studying the properties of certain biological systems by numerical methods, Biometrika, 45 (1958), 16–31. https://doi.org/10.2307/2333042 doi: 10.2307/2333042
    [20] X. Mao, G. Marion, E. Renshaw, Environmental Brownian noise suppresses explosions in population dynamics, Stochastic Processes Appl., 97 (2002), 95–110. https://doi.org/10.1016/s0304-4149(01)00126-0 doi: 10.1016/s0304-4149(01)00126-0
    [21] A. Gray, D. Greenhalgh, L. Hu, X. Mao, J. Pan, A stochastic differential equation SIS epidemic model, SIAM J. Appl. Math., 71 (2011), 876–902. https://doi.org/10.1137/10081856x doi: 10.1137/10081856x
    [22] D. Jiang, N. Shi, X. Li, Global stability and stochastic permanence of a non-autonomous logistic equation with random perturbation, J. Math. Anal. Appl., 340 (2008), 588–597. https://doi.org/10.1016/j.jmaa.2007.08.014 doi: 10.1016/j.jmaa.2007.08.014
    [23] C. Ji, D. Jiang, N. Shi, Analysis of a predator-prey model with modified Leslie-Gower and Holling-type Ⅱ schemes with stochastic perturbation, J. Math. Anal. Appl., 359 (2009), 482–498. https://doi.org/10.1016/j.jmaa.2009.05.039 doi: 10.1016/j.jmaa.2009.05.039
    [24] B. Wen, Z. Teng, Z. Li, The threshold of a periodic stochastic SIVS epidemic model with nonlinear incidence, Phys. A Stat. Mech. Appl., 508 (2018), 532–549. https://doi.org/10.1016/j.physa.2018.05.056 doi: 10.1016/j.physa.2018.05.056
    [25] Y. Chen, B. Wen, Z. Teng, The global dynamics for a stochastic SIS epidemic model with isolation, Phys. A, 492 (2018), 1604–1624. https://doi.org/10.1016/j.physa.2017.11.085 doi: 10.1016/j.physa.2017.11.085
    [26] X. Meng, S. Zhao, T. Feng, T. Zhang, Dynamics of a novel nonlinear stochastic SIS epidemic model with double epidemic hypothesis, J. Math. Anal. Appl., 433 (2016), 227–242. https://doi.org/10.1016/j.jmaa.2015.07.056 doi: 10.1016/j.jmaa.2015.07.056
    [27] S. He, S. Tang, L. Rong, A discrete stochastic model of the COVID-19 outbreak: forecast and control, Math. Biosci. Eng., 17 (2020), 2792–2804. https://doi.org/10.3934/mbe.2020153 doi: 10.3934/mbe.2020153
    [28] S. Li, S. Zhang, A research of pest management SI stochastic model with effect of pesticides function, J. Syst. Sci. Math. Sci., 37 (2017), 1379. https://doi.org/10.12341/jssms13165 doi: 10.12341/jssms13165
    [29] S. Zhang, S. Yuan, T. Zhang, A predator-prey model with different response functions to juvenile and adult prey in deterministic and stochastic environments, Appl. Math. Comput., 413 (2022), 126598. https://doi.org/10.1016/j.amc.2021.126598 doi: 10.1016/j.amc.2021.126598
    [30] M. Liu, H. Qiu, K. Wang, A remark on a stochastic predator-prey system with time delays, Appl. Math. Lett., 26 (2013), 318–323. https://doi.org/10.1016/j.aml.2012.08.015 doi: 10.1016/j.aml.2012.08.015
    [31] A. Skvortsov, B. Ristic, A. Kamenev, Predicting population extinction from early observations of the Lotka-Volterra system, Appl. Math. Comput., 320 (2018), 371–379. https://doi.org/10.1016/j.amc.2017.09.029 doi: 10.1016/j.amc.2017.09.029
    [32] F. Vadillo, Comparing stochastic Lotka-Volterra predator-prey models, Appl. Math. Comput., 360 (2019), 181–189. https://doi.org/10.1016/j.amc.2019.05.002 doi: 10.1016/j.amc.2019.05.002
    [33] G. Cai, Y. Lin, Stochastic analysis of the Lotka-Volterra model for ecosystems, Phys. Rev. E, 70 (2004), 041910. https://doi.org/10.1103/PhysRevE.70.041910 doi: 10.1103/PhysRevE.70.041910
    [34] Q. Liu, D. Jiang, Influence of the fear factor on the dynamics of a stochastic predator-prey model, Appl. Math. Lett., 112 (2021), 106756. https://doi.org/10.1016/j.aml.2020.106756 doi: 10.1016/j.aml.2020.106756
    [35] M. Liu, K. Wang, Global stability of a nonlinear stochastic predator-prey system with Beddington-DeAngelis functional response, Commun. Nonlinear Sci. Numer. Simul., 16 (2011), 1114–1121. https://doi.org/10.1016/j.cnsns.2010.06.015 doi: 10.1016/j.cnsns.2010.06.015
    [36] Q. Liu, L. Zu, D. Jiang, Dynamics of stochastic predator-prey models with Holling Ⅱ functional response, Commun. Nonlinear Sci. Numer. Simul., 37 (2016), 62–76. https://doi.org/10.1016/j.cnsns.2016.01.005 doi: 10.1016/j.cnsns.2016.01.005
    [37] M. Liu, M. Deng, Analysis of a stochastic hybrid population model with Allee effect, Appl. Math. Comput., 364 (2020), 124582. https://doi.org/10.1016/j.amc.2019.124582 doi: 10.1016/j.amc.2019.124582
    [38] S. Zhang, T. Zhang, S. Yuan, Dynamics of a stochastic predator-prey model with habitat complexity and prey aggregation, Ecol. Complexity, 45 (2021), 100889. https://doi.org/10.1016/j.ecocom.2020.100889 doi: 10.1016/j.ecocom.2020.100889
    [39] X. Mao, C. Yuan, Stochastic Differential Equations With Markovian Switching, Imperial College Press, 2006. http://doi.org/10.1142/p473
    [40] R. Khasminskii, Stochastic stability of differential equations, Springer Berlin, 2011. https://doi.org/10.1007/978-3-642-23280-0
    [41] S. Marino, I.B. Hogue, C. J. Ray, D. Kirschner, A methodology for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol., 254 (2008), 178–196. https://doi.org/10.1016/j.jtbi.2008.04.011 doi: 10.1016/j.jtbi.2008.04.011
  • This article has been cited by:

    1. Alexander Sterk-Hansen, Bjørn H.H. Saghaug, Daniel Hagen, Muhammad Faisal Aftab, 2023, A ROS 2 and TwinCAT Based Digital Twin Framework for Mechatronics Systems, 979-8-3503-1568-4, 485, 10.1109/ICCMA59762.2023.10374978
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

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

Metrics

Article views(2431) PDF downloads(114) Cited by(1)

Figures and Tables

Figures(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog