Research article Special Issues

Dynamics and asymptotic profiles of steady states of an SIRS epidemic model in spatially heterogenous environment

  • Received: 11 August 2019 Accepted: 24 October 2019 Published: 06 November 2019
  • This paper performs qualitative analysis on a reactionɃdiffusion SIRS epidemic system with ratioɃdependent incidence rate in spatially heterogeneous environment. The threshold dynamics in the term of the basic reproduction number R0 is established. And the asymptotic profile of endemic equilibrium is determined if the diffusion rate of the susceptible individuals is small. The results show that restricting the movement of susceptible individuals can effectively control the number of infectious individuals.

    Citation: Baoxiang Zhang, Yongli Cai, Bingxian Wang, Weiming Wang. Dynamics and asymptotic profiles of steady states of an SIRS epidemic model in spatially heterogenous environment[J]. Mathematical Biosciences and Engineering, 2020, 17(1): 893-909. doi: 10.3934/mbe.2020047

    Related Papers:

    [1] Santosh Kumar Sharma, D. K. Palwalia, Vivek Shrivastava . Distributed generation integration optimization using fuzzy logic controller. AIMS Energy, 2019, 7(3): 337-348. doi: 10.3934/energy.2019.3.337
    [2] Nagaraj C, K Manjunatha Sharma . Fuzzy PI controller for bidirectional power flow applications with harmonic current mitigation under unbalanced scenario. AIMS Energy, 2018, 6(5): 695-709. doi: 10.3934/energy.2018.5.695
    [3] Hari Charan Nannam, Atanu Banerjee . A novel control technique for a single-phase grid-tied inverter to extract peak power from PV-Based home energy systems. AIMS Energy, 2021, 9(3): 414-445. doi: 10.3934/energy.2021021
    [4] Rashid Al Badwawi, Mohammad Abusara, Tapas Mallick . Speed control of synchronous machine by changing duty cycle of DC/DC buck converter. AIMS Energy, 2015, 3(4): 728-739. doi: 10.3934/energy.2015.4.728
    [5] Katherin Indriawati, Bambang L. Widjiantoro, Ali Musyafa . Design of sensor and actuator fault tolerant control system on wind turbine benchmark for Region II. AIMS Energy, 2019, 7(2): 111-126. doi: 10.3934/energy.2019.2.111
    [6] Akash Talwariya, Pushpendra Singh, Mohan Lal Kolhe, Jalpa H. Jobanputra . Fuzzy logic controller and game theory based distributed energy resources allocation. AIMS Energy, 2020, 8(3): 474-492. doi: 10.3934/energy.2020.3.474
    [7] Imam Abadi, Chairul Imron, Mohammad Musa Bachrowi, Dwi Nur Fitriyanah . Design and implementation of battery charging system on solar tracker based stand alone PV using fuzzy modified particle swarm optimization. AIMS Energy, 2020, 8(1): 142-155. doi: 10.3934/energy.2020.1.142
    [8] Ahmed Neffati, Amira Marzouki . Local energy management in hybrid electrical vehicle via Fuzzy rules system. AIMS Energy, 2020, 8(3): 421-437. doi: 10.3934/energy.2020.3.421
    [9] Rasool M. Imran, Kadhim Hamzah Chalok . Innovative mode selective control and parameterization for charging Li-ion batteries in a PV system. AIMS Energy, 2024, 12(4): 822-839. doi: 10.3934/energy.2024039
    [10] Ahmed Eldessouky, Hossam Gabbar . SVC control enhancement applying self-learning fuzzy algorithm for islanded microgrid. AIMS Energy, 2016, 4(2): 363-378. doi: 10.3934/energy.2016.2.363
  • This paper performs qualitative analysis on a reactionɃdiffusion SIRS epidemic system with ratioɃdependent incidence rate in spatially heterogeneous environment. The threshold dynamics in the term of the basic reproduction number R0 is established. And the asymptotic profile of endemic equilibrium is determined if the diffusion rate of the susceptible individuals is small. The results show that restricting the movement of susceptible individuals can effectively control the number of infectious individuals.


    Over the years, the switched reluctance generator (SRG) has been a focal point for many research studies [1]. The SRG has a basic, rigid and simple structure. This device includes neither a permanent magnet nor twisting of the rotor. This construction decreases thus the expense of the SRG and its maintenance. The SRG can operate during high-velocity activities without the worry of mechanical problems [2].

    Accordingly, the SRG offers several advantages compared to other machines, e.g., the rotor does not require any winding and it is constituted only of ferromagnetic material. The major windings (losses) are concentrated in the stator. There are numerous sorts of structures that can drive the SRG [3]. The power electronic converter usually used for controlling the SRG is the asymmetric power electronic converter. It has the advantage of being able to independently control each phase of the SRG [4]. The stator windings are associated with the arrangement of the upper and lower switches of the inverter [5]. Due to the complex nonlinear characteristics of SRGs, there exist very few previous works dedicated to the identification and modeling of SRGs [6,7]. This was later described as nonlinear systems structured in blocks [8,9,10]. The most popular models of nonlinear systems in blocks can be found in [11,12,13]. The system nonlinearities can be static or dynamic [14,15,16]. Note that the problem of SRG speed control has been addressed by using several techniques, e.g., fuzzy logic control [17], artificial neural networks (ANNs) [18] and genetic algorithms [19,20]. Furthermore, the optimization techniques like particle swarm optimization and bacteria foraging can be used to regulate the SRG speed. Then, an ANN-proportional integral (ANN-PI) controller was introduced to regulate the power injected into the grid [21]. A proportional resonant controller that has been used to control the power produced by the SRG is given in [22]. To minimize the torque ripple of SRGs, an ANN control method has been developed, as presented in [23]. To control the position of the linear switched reluctance machine, a flower pollination algorithm was developed [24].

    The fuzzy logic technique was introduced by Zadeh in 1965 [25]. The fuzzy logic controller (FLC) has been widely investigated and used to regulate complex nonlinear systems. The FLC is an artificial decision-maker that relies on human decision-making behavior through the use of language rules rather than mathematical models [26]. An exact model is not necessarily required with the FLC technique and it is created using language information. Mathematical equations to describe the system to be controlled are not necessary for FLCs. The FLC has been used in many fields, e.g., in solar photovoltaic systems [27], wind energy [28], power systems [29] and smart agriculture systems [30]. Other previous works have yielded combinations of fuzzy logic with other methods. In [31], the authors used an adaptive neuro-fuzzy inference system to assist the power monitoring for a smart grid. Adaptive fuzzy control using a genetic algorithm has been used to overcome the influence of nonlinearities in the vessel dynamic positioning system [32]. In [33], a self-tuning FLC-based speed is proposed for the control of an SRG for wind energy applications. The ACO is a heuristic approach that can be used to find the global optimum of a given problem. In this respect, note that the ACO has already been employed in wind energy control [34] and for planning the wind farm layouts [35]. In [36], the authors discussed a design method entailing the use of the ACO based on a proportional-integral-derivative (PID) controller for a zeta converter.

    The aim of this study was to control the output voltage of the SRG. Techniques using the ACO-based PI controller and FLC have been developed. In this respect, the SRG output voltage is adjusted by using the turn-off angle θoff. To obtain the best value of θoff, two methods are proposed. The first one consists of a PI controller tuned by an ACO algorithm. In the second approach, the tun-off angle θoff is determined by using an FLC.

    For convenience, the main contributions of this work are summarized as follows:

    ● Unlike many previous methods dedicated to controlling the SRG voltage, the proposed methods are easy to implement.

    ● In this paper, very interesting tools and concepts are proposed, such as ACO for tuning the PI controller, as well as the FLC approach.

    ● The established algorithm, which applies an ACO technique for SRG voltage control, is new.

    ● To improve the results obtained via the ACO algorithm, four cost functions were applied and compared.

    ● In this study, a second solution to control the SRG output voltage by implementing a FLC approach was also developed. The best solution was chosen by comparing these two methods.

    The paper is organized as follows. The mathematical model and operation principle of SRGs are introduced in Section 2. In Section 3, the control methods for the SRG based on the ACO and FLC are described. An overview and description of the ACO are presented. Then, the cost functions used to optimize the ACO method are also proposed in this section. Then, the FLC is discussed. To show the effectiveness of these methods, simulation examples are presented in Section 4. Finally, comparisons between the results obtained using the two approaches are also discussed in this section. In Section 5, the concluding remarks of this paper are summarized.

    First, note that the SRG winding flux λ(.) depends on the phase current i and the rotor position θ. Then, for each stator phase, the voltage v can be expressed as the sum of the voltage drop on the resistor and the derivative of the flux:

    v=Rsi+dλ(θ,i)dt (1)

    where Rs is the phase resistor. The phase inductance L(θ,i) is a nonlinear function depending on the rotor position θ and the phase current i. Furthermore, the linkage flux λ(θ,i) can be expressed according to L and i as follows:

    λ(θ,i)=L(θ,i)i (2)

    Then, replacing the flux λ(θ,i) in Eq (2) with its expression according to L and i, the expression of phase voltage becomes

    v=Rsi+L(θ,i)didt+idθdtL(θ,i)θ (3)

    The last term on the right side of Eq (3) corresponds to the back electromotive force e. Specifically, one has

    e=idθdtL(θ,i)θ=iωL(θ,i)θ (4)

    where ω is the SRG's rotor speed. When the SRG operates as a generator, the current changes its direction. The back electromotive force thus becomes negative.

    The SRG stator is characterized by winding resistors with very small values [6]. The resistive voltage drop can be neglected [6]. Accordingly, the phase current i, the phase voltage, v and the back electromotive force e are related by the following equation:

    ve=L(θ)didt (5)

    Here, an asymmetric half-bridge (AHB) converter is proposed to drive the SRG (Figure 1). The SRG operates as a generator when the phase inductance L(θ) decreases; this means that dL(θ)/dt<0 (Figure 2).

    Figure 1.  SRG connected to an AHB converter.
    Figure 2.  Excitation and generation stages for one SRG phase.

    The SRG is interconnected via the DC link to a voltage source converter (VSC), which is linked to the electrical grid. Then, the DC link voltage VDC is regulated by the VSC control. The latter allows for the injection of the generated power from the SRG into the electrical grid.

    In this respect, note that the characteristics L(θ) and i(θ) are smooth curves and are often linearized in the study of generator and motor modes of SRGs, as shown in Figure 2 [37]. It can be easily seen in Figure 2 that the curve of the electrical current i can be subdivided into three main intervals, or the following three situations:

    ● The first situation is characterized by e<VDC. In this case, the phase current decreases after the excitation stage. This situation happens when the operating speed is decreased due to the electromotive force diminishing.

    ● The second situation can be obtained when e=VDC. Accordingly, after the duration of the excitation stage, the current tends to remain constant.

    ● The last situation can be seen when the SRG is working at high speeds; one thus has e>VDC. After the excitation stage, the phase current tends to rise.

    To determine the SRG's operation mode, the base speed should be evaluated. To improve the SRG performance, two techniques for SRG control can be used. For low operating speeds, the control of the SRG is performed by using the existing hysteresis method. The power control of the SRG is related to the adjustment of the reference current Iref. In this respect, the firing angle θon and turn-off angle θoff are the driving parameters (Figure 2). For high speeds, the control of the SRG is achieved by using a single-pulse technique. Accordingly, the power control is based on changing the turn-off angle θoff while the firing angle θon is kept constant.

    ACO is a meta-heuristic algorithm and probabilistic method. It is commonly used for solving hard combinatorial optimization problems based on graph representations. This method aims to find the best paths based on several possible graphs. This algorithm was initially proposed by Marco Dorigo in 1992 [38]. The artificial ants used in this technique were inspired by real ant colonies and designed to search for the best path. The shorter path is obtained by analyzing and combining the results of the path of each ant in the colony. The main issue is related to finding the shortest path among all of the trajectories traveled by different ants to reach the food. When the ants are moving, they leave a chemical pheromone trail on the ground. Depending on the distance of the path, the pheromone quality will change. Every ant chooses the path according to the intensity of the deposited pheromone. If there is no more deposit of pheromone, its intensity decreases according to the time. Specifically, the other ants are attracted by the pheromone; they thus choose the path where the pheromone is of high intensity. This path will be the best solution and has a great probability of being chosen.

    With this approach, the lifestyle behavior of real ant colonies is used to solve the optimization problem. Here, the artificial ant technique is proposed to seek the best solution by moving from one node to another. With this algorithm, the artificial ants move according to their previous positions, which have been stored in a specific data structure. Once the ants have accomplished their tour between the initial node and the last one, the pheromone consistency of each path is updated. The concentration of pheromone will be of high quality if the artificial ants finished their tour by taking the shortest path, and vice versa. The steps of the proposed algorithm are summarized in Figure 3.

    Figure 3.  ACO algorithm.

    For problem optimization, there are several error criteria commonly used in the literature, including the integral square error (ISE), integral of the absolute error (IAE), integral time of the absolute error (ITAE), integral of multiplied absolute error, quadratic error and total square variation. All of these error criteria have zero as the lower bound and can be considered as cost functions of the proposed ACO algorithm. To improve the overall performance obtained through the use of the PI controller, four cost functions were selected. The considered cost function criteria are ISE, mean squared error (MSE), IAE and ITAE. Specifically, the objective function of the ACO algorithm is based on minimizing the considered error criteria. The parameters of the PI controller can be obtained by using the optimization results.

    Here, the ACO-based control method for the SRG voltage is introduced. Figure 4 shows the closed-loop control corresponding to this technique, where Vref denotes the reference voltage and Vm stands for the SRG output voltage. With this method, the SRG is driven at high speeds. The firing angle θon is kept constant at 40°, while the turn-off angle θoff is adjusted by the algorithm. Accordingly, the magnetization cycle of the SRG phase is controlled by changing the value of θoff through the use of a PI controller.

    Figure 4.  SRG voltage control via PI controller implementation in ACO.

    Even if the PID controller would improve the system stability, the derivative action of the PID controller would likely increase the amplitude of disturbances. Furthermore, a couple of ripples is often produced in SRG applications. In order to avoid more disturbances, a PI controller was selected.

    In the proposed approach, the PI controller parameters are obtained by using an ACO algorithm. The transfer function of the proposed PI controller is as follows:

    Gc(S)=KP+KIS (6)

    where KP and KI are parameters of the PI controller.

    Generally, the FLC has three main stages, namely, fuzzification, fuzzy rule extraction and defuzzification. The inputs are turned into fuzzy sets by using linguistic elements and membership functions throughout the fuzzification process. The two most well-known fuzzy systems are the Mamdani and the Takagi-Sugeno-Kang models. Here, the Mamdani fuzzy systems have two inputs and one output is used. The voltage error signal ε (ε=VrefVout) and its derivative ˙ε are the inputs of the FLC. The output θFLC of the FLC corresponds to the turn-off angle θoff of the SRG.

    The design of the FLC was determined by assigning seven fuzzy sets for the inputs (ε, ˙ε) and the output (θFLC), namely, {NB (negative big), NM (negative middle), NS (negative small), ZZ (zero), PS (positive small), PM (positive middle), PB (positive big)}. In this study, the membership functions related to (ε, ˙ε) and θFLC are the chosen type of triangular shapes; they are given by Figures 5 and 6, respectively.

    Figure 5.  Membership functions corresponding to the inputs ε and ˙ε.
    Figure 6.  Membership functions corresponding to the output θFLC.

    The second method developed in this work consists of using an FLC to determine the best values of θoff while the angle θon is fixed at 40°. Figure 7 shows the closed-loop control corresponding to this technique.

    Figure 7.  SRG voltage control through the use of fuzzy logic.

    In this section, examples of simulations are presented to show the effectiveness of the proposed methods. The first method aims to control the output voltage of the SRG by using an ACO algorithm. This method has two steps. In the first step, the optimum parameters for each cost function (i.e., ISE, MSE, IAE and ITAE) of the PI controller are determined. The optimum PI controller parameters are chosen by taking the parameter values that yields the smallest error between the reference voltage and the output voltage of the SRG. In the second step, the obtained values are compared based on the time domain specifications. Specifically, the steady-state time, peak time and overshoot.

    In the second method, an FLC replaces the PI controller.

    Remark 1:

    The applied SRG is characterized by the parameters given in the Appendix (Part a). The initialization parameters of the ACO algorithm and PI controller are also given in the Appendix (Parts b and c).

    The approach used to initialize the ACO algorithm is based on results reported in the literature and simulations. First, we applied many combinations by changing the number of ants (10, 20, 30, 40, 50, 80) and the number of iterations (20, 50,100,200) in simulations. Second, in consideration of the obtained results and time constraints, we have chosen the optimal number of ants as 30 and optimal number of iterations as 100.

    Based on the ACO algorithm, the parameters for the optimal PI controllers corresponding to each cost function were obtained; they are presented in Table 1.

    Table 1.  Optimal PI parameters for each cost function.
    ACO KP KI
    ITAE 0.36977 2.497
    MSE 2.7879 5.8634
    IAE 1.802 0.82332
    ISE 1.2821 4.7698

     | Show Table
    DownLoad: CSV

    The responses of the SRG output voltage corresponding to the use of the optimal PI controllers have been plotted as shown in Figure 8 and Figure 9, where the reference voltage is constant (Vref=400V). The error variations for the ISE, IAE, MSE and IATE have been plotted as shown in Figure 10, Figure 11, Figures 12 and 13, respectively.

    Figure 8.  Output voltage and voltage reference based on ISE (a) and IAE (b).
    Figure 9.  Output voltage and voltage reference based on MSE (a) and ITAE (b).
    Figure 10.  Error variation for ISE.
    Figure 11.  Error variation for IAE.
    Figure 12.  Error variation for MSE.
    Figure 13.  Error variation for ITAE.

    To compare the performances of the PI controllers, their characteristic time-domain specifications (steady-state time, peak time, overshoot) were compared; the results are given in Table 2. Finally, the best PI controller was chosen. For convenience, the voltage ripple of the output voltage signal resulting from the use of the ACO-PI controller method and the cost error ITAE was estimated. The value of this voltage ripple was 0.76%.

    Table 2.  Time-domain characteristics of the output voltage according to the control method.
    Method Cost function type Peak time Peak Undershoot Overshoot Settling max Settling min Settling time Rise time
    PI ACO-ITAE 5.5277 400.1517 0 0.0015 400.1517 399.7504 9.9993 0.4553
    PI ACO-MSE 8.2291 400.1502 0 2.4386e-04 400.1502 399.7400 9.9997 1.2306
    PI ACO-IAE 9.9650 400.1138 0 1.4018e-04 400.1138 399.7126 9.9997 2.4563
    PI ACO-IAE 9.9650 400.1138 0 1.4018e-04 400.1138 399.7126 9.9997 2.4563
    FLC -------- ***** 400.8519 0 0.0596 400.8519 399.7932 ***** ******

     | Show Table
    DownLoad: CSV

    In the second method, the PI controller is replaced with an FLC. In this respect, two cases can be distinguished. The first one consists of taking the wind speed as constant; in the second case, the wind speed is variable. When the wind speed is taken as a constant (here, ω=157rad/s), the response of the output voltage of the SRG as a result of using the FLC was determined; it is shown in Figure 14; the error variation is shown in Figure 15. The results obtained for the time domain are summarized in Table 2. It can be seen that the output voltage signal of the FLC yields a voltage ripple with a value of 0.54%.

    Figure 14.  Output voltage obtained by using an FLC with ω = 157 rad/s and Vref = 400 V.
    Figure 15.  Error variation resulting from operating the FLC with ω = 157 rad/s and Vref = 400 V.

    Here, the wind speed was varied from 110rad/s to 170 rad/s. The wind speed curve is given in Figure 16. For this case, the SRG output voltage corresponding to the use of the FLC controller can be seen in Figure 17; the error variation is shown in Figure 18.

    Figure 16.  Profile of wind speed.
    Figure 17.  Output voltage resulting from implementation of the FLC with variable wind speed.
    Figure 18.  Error variation associated with the FLC and the condition of variable wind speed.

    The characteristics of the SRG output voltage signal corresponding to the ACO-PI control scheme with IATE, as well as that corresponding to the FLC scheme are summarized in Table 3.

    Table 3.  Characteristics of the SRG output voltage signal based on ACO-PI-IATE and FLC implementation.
    Output voltage ACO-PI-ITAE FLC
    Max 4.002e + 02 4.006e + 02
    Min 3.960e + 02 4.000e + 02
    Peak to peak 4.190e + 00 6.123e ­ 01
    Mean 3.999e + 02 4.003e + 02
    Median 4.000e + 02 4.003e + 02
    Root-mean-square 3.999e + 02 4.003e + 02

     | Show Table
    DownLoad: CSV

    ● Case of constant wind speed

    To show the effectiveness of this study, examples of simulations have been provided. The reference voltage was kept constant at Vref=400V. These simulations were established by applying a constant wind speed of 157rad/s and choosing a load of 1600 Ω. The results of the SRG output voltage that were obtained by using a conventional PI controller with different cost functions (i.e., ACO-ISE, ACO-MSE, ACO-IAE and ACO-ITAE), and as based on fuzzy control, are shown in Figures 8, 9 and 14. These results show that the SRG output voltage resulting from the use of the ACO-based PI controller and the FLC techniques converges to its reference value. It can be seen by using this method that a constant torque allows a constant voltage. On the other hand, the results summarized in Table 2 show that the ACO-based PI controller that utilizes ITAE provides the best results. Additionally, the time-domain characteristics show that the FLC improves the dynamic and steady-state features. Furthermore, the FLC yielded a smaller voltage ripple (0.54%) than the PI controller (0.76%).

    ● Case of variable wind speed

    The wind speed was varied in the range of 110–170 rad/s. The proposed simulations were established by taking the reference voltage Vref=400V and applying a load of 1600 Ω. Generally, the tuning of the ACO-based PI controller parameters requires significant computational time. Thus, to overcome this issue, an FLC controller was developed. It is shown in Figure 17 that the SRG output voltage converges fast to the reference voltage even if the wind speed changes.

    This article discusses the problem of SRG output voltage control. This SRG was implemented in a wind turbine operating at 75 KW that was connected to an electrical grid. Here, a PI controller based on the ACO technique has been proposed to regulate the SRG output voltage. To determine the PI controller parameters, several cost functions were used. By comparing the different results obtained by using the given cost functions, the best PI controller parameters were chosen.

    When the wind speed was taken as constant, the simulation results confirmed that the output voltage associated with the ACO-PI controller technique and a cost function (i.e., ITAE, MSE, IAE or ISE) converges to the reference voltage. Furthermore, it has been shown that the solution based on the ITAE gives the best results for the control solution based on the ACO-PI controller. It is interesting to note that this approach requires a short simulation time, unlike several other methods. On the other hand, the FLC yielded better results than the ACO-PI controller for variable wind speeds.

    The system parameters applied in the simulations are as shown below.

    (a) SRG parameters:

    Number of stator and rotor poles, 12/8 (respectively); Frequency, [F]=50Hz; DC supply voltage, [Vdc]=240V; Reference current, 200A; Hysteresis band, [+10,10]; Mechanical load torque, Tm = 10 N·m; DC link capacitor, C 3e3F; Speed of wind, ω=157rad/s; Voltage of excitation, Vewc=50V; Initial voltage, VCCinit=400V; Vdc reference, VCCref=400V; Load, R=1600Ω; Upper limit of turn-off angle, θoffmax=21; Firing angle, θon=40; Maximum current, Imax = 8.5 A.

    (b) PI controller parameters:

    The boundaries of Kp and Ki are [0.15] and [0.510], respectively.

    (c) ACO parameters:

    Number of iterations, nn = 100; Number of ants = 30; Pheromone decay parameter, α = 0.8; Relative importance of pheromone with respect to distance, β = 0.2; Evaporation rate = 0.7; Number of parameters = 2; Number of nodes = 1000.

    All authors declare no conflict of interest regarding this paper.



    [1] R. May and R. Anderson, Spatial heterogeneity and the design of immunization programs, Math. Biosci., 72 (1984), 83-111.
    [2] H. Hethcote and J. W. Van Ark, Epidemiological models for heterogeneous populations: proportionate mixing, parameter estimation, and immunization programs, Math. Biosci., 84 (1987), 85-118.
    [3] V. Capasso, Mathematical structures of epidemic systems, volume 97. Springer, 1993.
    [4] M. E. Alexander and S. M. Moghadas. Periodicity in an epidemic model with a generalized nonlinear incidence, Math. Biosci., 189 (2004), 75-96.
    [5] W. D. Wang, Epidemic models with nonlinear infection forces, Math. Biosci. Eng., 3 (2006), 267-279.
    [6] D. Xiao and S. Ruan, Global analysis of an epidemic model with nonmonotone incidence rate, Math. Biosci., 208 (2007), 419-429.
    [7] Y. Cai, Y. Kang and W. M. Wang, A stochastic SIRS epidemic model with nonlinear incidence rate, Appl. Math. Comp., 305 (2017), 221-240.
    [8] W. Liu, S. A. Levin and Y. Iwasa, Influence of nonlinear incidence rates upon the behavior of sirs epidemiological models, J. Math. Biol., 23 (1986), 187-204.
    [9] W. Liu, H. W. Hethcote and S. A. Levin, Dynamical behavior of epidemiological models with nonlinear incidence rates, J. Math. Biol., 25 (1987), 359-380.
    [10] H. W. Hethcote, The mathematics of infectious diseases, SIAM Rev., 42 (2000), 599-653.
    [11] B. Fred and C.-C. Carlos, Mathematical models in population biology and epidemiology (Second Edition). Springer, 2012.
    [12] Y. Cai, Y. Kang, M. Banerjee, et al., A stochastic SIRS epidemic model with infectious force under intervention strategies, J. Differ. Equations, 259 (2015), 7463-7502.
    [13] S. Yuan and B. Li, Global dynamics of an epidemic model with a ratio-dependent nonlinear incidence rate, Discrete Dyn. Nat. Soc., 2009 (2009), 609306.
    [14] C. Neuhauser, Mathematical Challenges in Spatial Ecology, Notices AMS, 48 (2001), 1304-1314.
    [15] S. Ruan, Spatial-Temporal Dynamics in Nonlocal Epidemiological Models, In: Takeuchi Y., Iwasa Y., Sato K. (eds) Mathematics for Life Science and Medicine. Springer, Berlin, Heidelberg, 2007.
    [16] W. E. Fitzgibbon, M. Langlais and J. J. Morgan, A reaction-diffusion system modeling direct and indirect transmission of diseases, Discrete Cont. Dyn.-B, 4 (2004), 893-910.
    [17] Z. Wang and J. Wu, Travelling waves of a diffusive Kermack-McKendrick epidemic model with non-local delayed transmission, P. Roy. Soc. A-Math. Phy., 466 (2010), 237-261.
    [18] Y. Cai and W. M. Wang, Dynamics of a parasite-host epidemiological model in spatial heterogeneous environment, Discrete Cont. Dyn.-B, 20 (2015), 989-1013.
    [19] Y. Cai and W. M. Wang, Fish-hook bifurcation branch in a spatial heterogeneous epidemic model with cross-diffusion, Nonl. Anal. Real World Appl., 30 (2016), 99-125.
    [20] W. M. Wang, X. Gao, Y. Cai, et al., Turing patterns in a diffusive epidemic model with saturated infection force, J. Franklin Inst., 355 (2018), 7226-7245. doi: 10.1016/j.jfranklin.2018.07.014
    [21] Y. Cai, Y. Kang, M. Banerjee, et al., Complex dynamics of a host-parasite model with both horizontal and vertical transmissions in a spatial heterogeneous environment, Nonl. Anal. Real World Appl., 40 (2018), 444-465.
    [22] P. Magal, G. Webb and Y. Wu, On a vector-host epidemic model with spatial structure, Nonlinearity, 31 (2018), 5589-5614.
    [23] Y. Cai, Z. Ding, B. Yang, et al., Transmission dynamics of Zika virus with spatial structure-A case study in Rio de Janeiro, Brazil, Phys. A, 514 (2019), 729-740.
    [24] J. Ge, K. Kim, Z. Lin, et al., An SIS reaction-diffusion-advection model in a low-risk and high-risk domain, J. Differ. Equations, 259 (2015), 5486-5509.
    [25] Y. Cai, X. Lian, Z. Peng, et al., Spatiotemporal transmission dynamics for influenza disease in a heterogenous environment, Nonl. Anal. Real, 46 (2019), 178-194.
    [26] E. E. Holmes, M. A. Lewis, J. E. Banks, et al., Partial differential equations in ecology: Spatial interactions and population dynamics, Ecology, 75 (1994), 17-29.
    [27] T. Caraco, M. Duryea, G. Gardner, et al., Host spatial heterogeneity and extinction of an SIS epidemic, J. Theor. Biol., 192 (1998), 351-361.
    [28] A. L. LLoyd and R. M. May, Spatial heterogeneity in epidemic models, J. Theor. Biol., 179 (1996), 1-11.
    [29] J. Dushoff and S. Levin, The effects of population heterogeneity on disease invasion, Math. Biosci., 128 (1995), 25-40.
    [30] S. Merler and M. Ajelli, The role of population heterogeneity and human mobility in the spread of pandemic influenza, Proc. R. Soc. B, 277 (2010), 557-565.
    [31] B. T. Grenfell, O. N. Bjornstad and J. Kappey, Travelling waves and spatial hierarchies in measles epidemics, Nature, 414 (2001), 716-723.
    [32] M. J. Keeling, M. E. Woolhouse; D. J. Shaw, et al., Dynamics of the 2001 UK foot and mouth epidemic: stochastic dispersal in a heterogeneous landscape, Science, 294 (2001), 813-817. doi: 10.1126/science.1065973
    [33] V. Colizza, A. Barrat, M. Barthelemy, et al., The role of the airline transportation network in the prediction and predictability of global epidemics, Proc. Natl. Acad. Sci. USA, 103 (2006), 2015-2020.
    [34] L. Hufnagel, D. Brockmann and T. Geisel, Forecast and control of epidemics in a globalized world, Proc. Natl. Acad. Sci. USA, 101 (2004), 15124-15129.
    [35] D. Henry and D. B. Henry, Geometric theory of semilinear parabolic equations, Springer-Verlag Berlin, 1981.
    [36] M. H. Protter and H. F. Weinberger, Maximum principles in differential equations, Prentice-Hall, New Jersey, 1967.
    [37] K. Yamazaki and X. Wang, Global well-posedness and asymptotic behavior of solutions to a reaction-convection-diffusion cholera epidemic model, Discrete Cont. Dyn-B, 21 (2017), 1297-1316.
    [38] Z. Du and R. Peng, A priori L estimates for solutions of a class of reaction-diffusion systems, J. Math. Biol., 72 (2016), 1429-1439.
    [39] N. K. Vaidya, F.-B. Wang and X. Zou, Avian influenza dynamics in wild birds with bird mobility and spatial heterogeneous environment, Discrete Cont. Dyn.-B, 17 (2013), 2829-2848.
    [40] O. Diekmann, J. A. P. Heesterbeek and J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365-382.
    [41] P. Van den Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29-48.
    [42] L. J. S. Allen, B. M. Bolker, Y. Lou, et al., Asymptotic profiles of the steady states for an SIS epidemic reaction-diffusion model, Discrete Cont. Dyn-A, 21 (2008), 1-20.
    [43] R. Peng and X. Zhao, A reaction-diffusion SIS epidemic model in a time-periodic environment, Nonlinearity, 25 (2012), 1451-1471.
    [44] W. D. Wang and X.-Q. Zhao, Basic reproduction numbers for reaction-diffusion epidemic models, SIAM J. Appl. Dyn. Sys., 11 (2012), 1652-1673.
    [45] Y. Lou and T. Nagylaki, Evolution of a semilinear parabolic system for migration and selection without dominance, J. Differ. Equations, 225 (2006), 624-665.
    [46] X. Zhao, J. Borwein and P. Borwein, Dynamical systems in population biology, volume 16. Springer, 2003.
    [47] P. Magal and X.-Q. Zhao, Global attractors and steady states for uniformly persistent dynamical systems, SIAM J. Math. Anal., 37 (2005), 251-275.
    [48] M. A. Pino, A priori estimates and applications to existence-nonexistence for a semilinear elliptic system, Indiana U. Math. J., 43 (1994), 77-129.
  • This article has been cited by:

    1. Hafid Oubouaddi, Fatima Ezzahra El Mansouri, Ali Bouklata, Ramzi Larhouti, Abdelmalek Ouannou, Adil Brouri, Parameter Estimation of Electrical Vehicle Motor, 2023, 18, 2224-2856, 430, 10.37394/23203.2023.18.46
  • Reader Comments
  • © 2020 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(4488) PDF downloads(535) Cited by(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog