key: cord-0955999-clbut1ot authors: Liu, Qiong; Zhang, Meng; Chen, Lansun title: State feedback impulsive therapy to SIS model of animal infectious diseases() date: 2019-02-15 journal: Physica A DOI: 10.1016/j.physa.2018.09.161 sha: 39ca0b9d227bdaefed47859d6026007859774a57 doc_id: 955999 cord_uid: clbut1ot Controlling animal infectious diseases and its related infra microbes is of great significance to public health, since a lot of infectious diseases originate from animal epidemics and they often threaten human health. A state feedback impulsive model is constructed to depict the transmission and treatment of animal epidemics. Basing on the impulsive model, the existence of order-1 periodic solution and its stability are proved with a novel method. The theoretical results indicate that the impulsive treatment triggered by the number of infectious is an efficient approach to control animal infectious disease from breakout. Numerical simulation is presented to support the theoretical conclusion in the end. Livestock, poultry and wildlife help human in supplying food resource and keeping balance of ecosystem etc that people cannot subsist without getting support from animals. However, animal infectious throw mankind a deadly menace with infectious diseases initiating from animals [1] [2] [3] . Some diseases such as African swine fever and bovine leukemia only influence animals. But to some epidemic, those pathogenic bacteria can be spilled from domestic animals or wildlife to human beings, and the development and mutation during the transmission will make the pathogenic bacteria difficult to control [4] [5] [6] . A lot of diseases besetting human beings such as human immunodeficiency virus (HIV), SARS corona virus, Ebola virus and avian influenza viruses (AIV) all can trace their filiation from animals [7] . Taylor's achievement showed that about 75% of human epidemics stem from an animal source, and these zoo-noses are quite difficult to control after its prevalence among a crowd of human [8] . So preventing the formation and spread of zoonoses is quite important to public health of human beings. In recent achievements, some scholars have got some useful results in impulsive dynamic system. Among them, Jianjun Jiao et al. established an SIR model with impulsive vaccination and impulsive dispersal to analyze the effect of impulsive vaccination and restricting infected individuals boarding transports on disease spread [9] . Linfei Nie et al. proposed a state-dependent impulsive model with two control threshold values to study two-languages competition and they found the order-1 periodic solutions using qualitative analysis method [10] . Shengqiang Zhang et al. proposed a stochastic nonautonomous Lotka-Volterra predator-prey model with impulsive effects and investigate its stochastic dynamics [11] . Meng Zhang et al. considered the preventive control of computer worm's transmission through Internet [12] . Recently discussions regarding impulsive disturbance based on state feedback have dominated the research in field of control theory and optimization [13] [14] [15] [16] [17] [18] . Enlightened by those research about impulsive model, we construct a state feedback impulsive SIS model to analyze the therapy control of animal infectious diseases. In practice, animal epidemic is hardly to be detected during its latent period because no obvious symptoms can be observed. Veterinarians can find the disease only after the eruption and then can treatment measures be carried out. Nevertheless, a certain amount of animals have been infected by the virus at that time. Thus, it is more appropriate to represent the treatment of animal epidemics as a state feedback impulse which depends on the variable of infected animals. Without loss of generality, we assume that the disease can be detected when the amount of infected animals cumulate to a certain level I * , and the treatment will be carried out at the same time. We also assume that the infected animals recover after treatment and become susceptible again. Then the model is established as follows: where S(t) denotes the susceptible animals; I(t) represents the infected ones. According to the transmission feather, we make the following assumption: (1) animals increase in logistic pattern, and the environmental capacity is K ; (2) only the healthy animals can propagate and the newborn animals are susceptible ones; (3) the infected rate of susceptible is linear with both susceptible and infected animals at β; (4) natural death rate is d 0 and death rate caused by the epidemic is d ; (5) I * is the threshold that the infectious disease can be discovered, when the amount of infections is smaller than I * , the infectious disease develops according to the common way which is expressed as the first two equations of model (1.1); (6) once the threshold I * is reached, the infectious disease can be spotted and efficient treatment will be performed immediately. A portion of infections, here we assume the ratio is θ (0 < θ < 1), will be treated and turn into susceptible ones and can be infected again. Considering the practical significance, all the parameters in this paper are positive. The sketch of this paper is organized as the following: some preliminaries about the state feedback impulsive dynamic system were presented in Section 2; the existence and stability of the order-1 periodic solution to system (1.1) were analyzed in Sections 3 and 4 separately. Finally, we examine the theoretical results with numerical simulation in Section 5. In this section some definitions and lemmas about state feedback impulsive dynamic system are listed. These preliminaries will be used in Sections 3 and 4. (2.1) The dynamic system defined by (2.1) is called a state feedback impulsive dynamic system, and it is denoted by (Ω, f , ϕ, M). The initial point P of system (2.1) is defined anywhere in the plane R 2 but M{x, y}, i.e. P ∈ Ω = R 2 \M{x, y}. ϕ is a continuous mapping which satisfies ϕ(M) = N, and we called it impulse mapping. M{x, y} and N{x, y} are impulse set, and phase set respectively, and they are lines or curves in the plane R 2 . The order-1 periodic solution and successor function is defined as follows. For practice meaning, we only consider system (1. Without considering the impulse mapping, system (1.1) reduces into in which we only focus on the situation in the first quadrant. It is obvious that O(0, 0) and E 1 (K (r − d 0 ), 0) are two marginal equilibrium of system (3.1). Following the theory of stability, ) , then E 1 is a stable node (see Fig. 2 ) and the infectious disease will eliminate automatically without any human intervention. If R > 1, then E 1 is a saddle and system (3.1) has a unique positive equilibrium βK )), and the Jacobin matrix of point E 2 is The characteristic equation is λ 2 βK )) is a locally stable focus or node (see Fig. 3 ). This is the certain case we focus on in the following of this paper. To prove the existence of order-1 periodic solution of system (3.1), we need the following two lemmas. is a periodic trajectory of system (3.1) and its period is T , then the periodic trajectory is stable. is globally asymptotically stable. The proves of those lemmas are similar with the process in [12] and they are omitted here. βK ), then we can prove that system (1.1) has an order-1 periodic solution. In the following we will prove the existence of point C satisfying F (C ) < 0. First, we prove that in the first quadrant the line whose slope is −1 cannot be tangent with any trajectory of system (1.1). To system (1.1), on the point where the slope of direction vector is −1 the coordinate must satisfy Then we have Its graph is a parabola with the origin as vertex, downward opening and located under the x axis. So it is obvious that the line whose slope is −1 cannot be tangent with any trajectory of system (1.1) in the first quadrant. Second, there exists a line L whose slope is −1 and the trajectory of system (1.1) runs across line L from right side to left side at the points whose vertical coordinate satisfying I < I * . Assuming V = S + I, to system (1.1) we have If the parameter e satisfying e > (r − d 0 )K + I * , there must exist dV dt < 0. Then on the points of line L satisfying I < I * , the trajectories of system (1.1) pass line L from right side to left side. Assuming the line L passes the impulsive set, the phase set and x axis at D, C and E respectively (see Fig. 5 ). There exists a trajectory L C of system (1.1) which passes the line L from right side to left side at point C and runs upward-left until reaching the impulsive set at M C , and M C must be located on the left side of point D. Then the impulsive function ϕ maps M C to N C . For M C N C is parallel to L, the phase point N C must be located on the left side of point C . So the successor function of point C From the continuity of successor function, there must exist a point A between B and C in phase set which satisfies F (A) = 0. Then the trajectory passing point A forms the order-1 periodic solution of system (1.1) (see Fig. 6 ). In the following part we focus on the uniqueness of the order-1 periodic solution to system (1.1). If R > 1, for any θ(0 < θ < 1), the order-1 periodic solution of system (1.1) is unique. Case II: Select two different points I and J in phase set N satisfying x I < x J < q, and they are located on different sides of horizontal isoclinic line dI dt = 0 (see Fig. 8 ). Then, F (I) − F (J) => 0. The prove is similar with that in case I. Case III: Select two different points I and J in phase set N satisfying x I < x J < q, and both I and J are located on right side of horizontal isoclinic line dI dt = 0 (see Fig. 9 ). Then, there exist trajectories L I and L J initiating from I and J, and they cross The proofs of Case IV and Case V are similar to that of Case III and they are omitted here. By all accounts above, the successor function F (·) is monotonically decreasing in the phase set. Thus there exists only one point A satisfying F (A) = 0. Hence, the order-1 periodic solution to system (1.1) is unique. This completes the proof. In this section, we study the stability of the order-1 periodic solution to system (1.1). To the order-1 periodic solution of a differential dynamic system with state feedback impulsive, it is fantastic to prove the stability with the technique of ordinary differential dynamic system. To explain a novel method, we restate the concept of stability and establish an orthogonal coordinate before proving. ). Assume Γ is an order-1 periodic solution of a certain state feedback impulsive system. If there exists a sufficiently small neighborhood U(Γ ) such that the ω limit set of any trajectory whose starting point P ∈ U(Γ ) is Γ , then Γ is stable. Otherwise, it is defined unstable. Following the consequence of theorems in part 3, for any θ (0 < θ < 1), system (1.1) has an exclusive order-1 periodic solution denoted by Γ , and we assume A is the intersection of Γ and phase set N (see Fig. 10 ). Then to any point S 0 satisfying |AS 0 | → 0, we define a point-series where S i+1 is the subsequent point of S i , i = 0, 1, . . . , k, . . .. Denote their coordinates with s 0 , s 1 , . . . , s k , s k+1 , . . ., where In this expression, d k is the distance between S k and A, then we can easily find 0 is the coordinate of A. Without loss of generality, we assume that functions P(x, y) and Q (x, y) of system (2.1) have arbitrary order partial derivatives. We also assume the time period from A to C on the certain trajectory of (2.1) is T , then the period of the closed orbit Γ is T too. Then we can set up a new coordinate system (s, n) in adequately closed scope of orbit Γ . To any point in that field, there exists a certain point in the closed orbit Γ and its normal line passes the point. We define the first coordinate s as the length of the trajectory between the certain point and A on Γ , while the second coordinate n as the number of the normal vector (here we assume upside as the positive direction) (See Fig. 11 ). With the first coordinate s we have set up, we can establish the parametric equation of curveÂBC x = φ(s), y = ψ(s), (4.1) where s is the parameter. Then we can expressed the rectangular coordinate (x, y) with the orthogonal coordinate (s, n) as Following (2.1) and (4.2), we have It is obvious that n = 0 is a solution of (4.3) and the first-order partial derivative of F (s, n) with respect to n is continuous, then (4.3) can be expressed briefly as where P 0 and Q 0 represent the values of P and Q at point A, P x0 , P y0 , Q x0 and Q y0 are the partial derivatives of P and Q at n = 0 with respect to x, y respectively. It is obvious that P 0 φ ′′ + Q 0 ψ ′′ = 0. Then the order-1 linear approximate equation of (4.4) is and its solution is The proofs of Theorems 4.1 and 4.2 are similar with that in [12] , and omitted here. Then we can get the corollary directly. ≜ q(S, I). In this system, we can calculate where p S0 and q I0 are the partial derivatives of p and q at n = 0 with respect to S and I respectively. Inspired with Corollary 4.1, we can reach the verdict that the order-1 periodic solution of system (1.1) is stable. This completes the proof. System (1.1) describes the development of animal epidemics and its impulsive treatments, it is the most advanced model which can be regarded as an analog of the procedure. We have proved that system (1.1) has a unique stable periodic solution which means the animal epidemics will not break out with factitious treatment basing on the flip-flop of the infectious' number. And in this part, we examine the theoretical results with real parameter values. Let r = 0.2, K = 5, β = 0.7, d = 0.2, d 0 = 0.08, θ = 0.7, the impulse set and phase set are M = {(x, y)|x ≥ 0, y = 0.0457} and N = {(x, y)|x ≥ 0, y = 0.0137} respectively, and the condition R > 1 is satisfied. Here, we choose two groups of initial values, S 0 = 0.4, I 0 = 0.01 and S 0 = 0.65, I 0 = 0.01. The numerical simulation results verify the conclusion of Theorems 3.1, 3.2 and 4.3, that means system (1.1) has a unique stable periodic solution (see Fig. 12 ). From the simulation results, we can find that system (1.1) has a unique limited cycle, and it does not depend on the initial values. In practice, the stable periodic solution means the animal epidemic can be controlled by intermittent state feedback treatment. Apparently the number of infectious animals must be observed and recorded, since the treatment is activated by the value. Once the number of infectious animals accumulates to the threshold of impulsive set, human will carry out some treatment measures and the number of the infectious drop to the imagine set. That work costs huge manpower and material resources. In fact, since the periodic solution is unique and stable, the period of the periodic solution is a constant. If we observe the phenomenon that the interval between each time the amount of the infectious reaches the threshold is a constant, it forebodes the periodic solution has formed. In that case, we can carry out artificial intervention depending on the fix period instead of observing the infectious and record the amount. This theoretical result can improve the efficiency of animal epidemic's treatment and control. Spatio-Temporal Attributes of Pandemic and Epidemic Diseases Antigenic and molecular heterogeneity in recent swine influenza A(H1N1) virus isolates with possible implications for vaccination policy The evolutionary dynamics of influenza A virus adaptation to mammalian hosts Influenza virus non-structural protein NS1: Interferon antagonism and beyond Phosphorylation of Influenza A virus NS1 protein at threonine 49 suppresses its interferon antagonistic activity Intuitive representation of surface properties of biomolecules using BioBlender The role of wildlife in transboundary animal diseases Risk Factors for Human Disease Emergence Impulsive vaccination and dispersal on dynamics of an SIR epidemic model with restricting infected individuals boarding transports State impulsive control strategies for a two-languages competitive model with bilingualism and interlinguistic similarity Dynamics analysis and numerical simulations of a stochastic non-autonomous predator-prey system with impulsive effects A state feedback impulse model for computer worm control Dynamics of a nonautonomous stochastic SIS epidemic model with double epidemic hypothesis Stochastic inequalities and applications to dynamics analysis of a novel SIVS epidemic model with jumps Application of stochastic inequalities to global analysis of a nonlinear stochastic SIRS epidemic model with saturated treatment function Global properties of a cell mediated immunity in HIV infection model with two classes of target cells and distributed delays An impulsive state feedback control model for releasing white-headed langurs in captive to the wild Global properties of a class of virus infection models with multitarget cells, Nonlinear Dynam Optimal control strategy for an impulsive stochastic competition system with time delays and jumps Modeling the state dependent impulse control for computer virus propagation under media coverage A novel method for analyzing the stability of periodic solution of impulsive state feedback model Qualitative analysis of impulsive state feedback control to an algae-fish system with bistable property The periodic solutions of the impulsive state feedbck dynamical system Dynamics on a Holling II Predator-Prey Model With State-Dependent Impulsive Control Periodic Solution of Prey-Predator Model with Beddington-DeAngelis Functional Response and Impulsive State Feedback Control Homoclinic bifurcation of prey-predator model with impulsive state feedback control Periodic solution of the system with impulsive state feedback control Limit.cycle. theory, Limit cycle theory We would like to sincerely thank the reviewers for their careful reading and constructive opinions of the original manuscript. This work is supported by NSFC, China