# Chemostat Model of Competition between Plasmid-Bearing and Plasmid-Free Organism with the Impulsive State Feedback Control

Authors: Fengmei Tao; Zhong Zhao; Lansun Chen
Journal: Discrete Dynamics in Nature and Society
 (2018)

---

## Abstract

In this paper, we propose a chemostat model of competition between plasmid-bearing and plasmid-free organism with the impulsive state feedback control. The sufficient condition for existence of the positive period-1 solution is obtained by means of successor function and the qualitative properties of the corresponding continuous system. We show that the impulsive control system is more effective than the corresponding continuous system if we choose a suitable threshold value of the state feedback control in the process of manufacturing the desired products through genetically modified techniques. Furthermore, a new method of proving the stability of the order-1 periodic solution is given based on the theory of the limit cycle of the continuous dynamical system. Finally, mathematical results are justified by some numerical simulations.


---

## Body

## 1. Introduction

With the rapid development of biotechnology, manufacturing the desired products through genetically modified techniques has been widely applied in many fields, such as agriculture, industrial biotechnology, and medicine. In general, genetical alteration is carried out by the introduction of a recombinant DNA into the cell in the form of a plasmid [1]. However, the plasmid may be lost in the reproductive process, which will introduce the plasmid-free organism into the process [2]. The genetically altered organism (plasmid-bearing organism) is less competitor than the plasmid-free organism because of the burden imposed by production [3]. At the same time, plasmid-free organism not only consumes the nutrient but also reduces the efficiency of the desired products. Therefore, the key problem is how to improve the efficiency of the desired products and to reduce the negative effect of the plasmid-free organism on the plasmid-bearing organism. Researchers [2–8] have done a lot of work to reduce the negative effect of the plasmid-free organism on the production efficiency of the desired products. In [2], a selective media was added from an external source and authors investigated the global stability of the equilibria by using the Dulac criterion and the Poincaré-Bendixson Theorem. Ai [5] introduced an inhibitor into the model of competition between plasmid-bearing and plasmid-free organisms to reduce the effect of the plasmid-free organism on the plasmid-bearing organism. Yuan et al. [6] proposed a model for competition between plasmid-bearing and plasmid-free organisms in the chemostat with an external inhibitor and obtained some sufficient conditions of global attractivity to the extinction equilibria by constructing appropriate Lyapunov-like functionals. Based on a feedback control, Dimitrova [9] formulated a dynamical model of plasmid-bearing, plasmid-free competition in the chemostat with general specific growth rates, and global stabilization of the dynamics towards a practically important coexistence equilibrium point is obtained. Shi et al. [10] investigated a chemostat model with plasmid-bearing, plasmid-free competition and impulsive effect, and they obtained the invasion threshold of the plasmid-free organism and plasmid-bearing organism.Recently, impulsive state feedback control has been widely used in chemostat model [11–14]. Tian et al. [12] studied the dynamics of the bioprocess with the impulsive state feedback control and obtained the existence and stability of period-1 solution of the bioprocess. In [13], authors proposed a turbidostat model with the feedback control, and the sufficient conditions of existence of positive order-1 periodic solution were obtained by using the existence criteria of periodic solution of a general planar impulsive autonomous system. As far as the authors know, little information is given about the introduction of the impulsive state feedback control into chemostat model of competition between plasmid-bearing and plasmid-free organism. In this paper, we will formulate a mathematical model of competition between plasmid-bearing and plasmid-free organism with the impulsive state feedback control so as to further improve the efficiency of the desired products in the process of the genetic alteration.The paper is organized as follows: a mathematical model of competition between plasmid-bearing and plasmid-free organism with the impulsive state feedback control is proposed in Section2. In Section 3, the qualitative analysis of system without impulsive control is given. Furthermore, the existence and stability of order-1 periodic solution are investigated in Section 4. Finally, we give some numerical simulations and a brief discussion.


## 2. Model Description and Preliminaries

In the process of manufacturing the desired products through genetically modified techniques, an important way to control the negative effect is to reduce the competitive advantage of the plasmid-free organism.The process of the genetic alteration is complex because it includes a series of biological and chemical reactions. Therefore, understanding the dynamic mechanism of genetic alteration is the most important factor for improving the efficiency of the desired products. Luo et al. [15] formulated a model of plasmid-bearing and plasmid-free competition in a chemostat as follows: (1)dSdt=DS0-S-σ1Sx1-σ2Sx2,dx1dt=x1f1S1-q-D,dx2dt=x2f2S-D+qx1f1S,where S(t) is the nutrient concentration at time t. x1 denotes the concentration of plasmid-bearing organisms at time t, and x2 is the concentration of plasmid-free organisms at time t. The consumption rates and the specific growth rates of plasmid-bearing and plasmid-free organisms are σi and fi (i=1,2), respectively. The probability that a plasmid is lost in reproduction is represented by q(0<q<1).Although the reason why plasmid is lost in the reproductive process is not completely understood by theoretical biologists even at present time, it has been experimentally validated that the plasmid-bearing cells have at least for some substrate concentrations a lower maximum specific growth rate than the plasmid-free counterpart [9, 16, 17]. To enhance the efficiency of the desired products and prevent the negative effect of the plasmid-free organism on the plasmid-bearing organism, we try to find the critical threshold value (where plasmid-bearing cells reach the maximum specific growth rate), which can be obtained by the experiment. When the growth rate of the plasmid-free organism reaches the critical threshold value, we begin to extract the plasmid-free organism so as to control the concentration of the plasmid-free organism to be lower than the critical threshold value. Based on [11–14], we introduce the impulsive state feedback control into system (1): (2)dSdt=DS0-S-μ1Sx1δ-μ2Sx2δ,dx1dt=x1μ1S1-q-D,dx2dt=x2μ2S-D+μ1qSx1,x2<h,Δx1=-E1x1,Δx2=-E2x2,x2=h,where μ1 and μ2 denote the maximal growth rate of the plasmid-bearing and plasmid-free organisms. δ is called a yield constant. h is a critical threshold value which can be obtained by experiment. E1 and E2(0<Ei<1,i=1,2) are the ratios of the plasmid-bearing and plasmid-free organisms extracted from the chemostat.We firstly consider the qualitative property of (2) without the impulsive effect. The variables in the above system may be rescaled by measuring S=S0S~,x1=δS0x,x2=δS0y,t=τ/D, and then system (2) becomes (3)dS~dτ=1-S~-μ1S0S~xD-μ2S0S~yD,dxdτ=xμ11-qS0DS~-1,dydτ=yμ2S0S~D-1+μ1qS0S~xD,y<H,Δx=-E1x,Δy=-E2y,y=H,where H=h/δS0.Lemma 1.
Supposeω(τ)=(S~(τ),x(τ),y(τ)) is a solution of (3) subject to ω(0+)≥0, and then ω(τ)≥0 for all τ≥0, and further ω(τ)>0,τ≥0 if ω(0+)>0.The proof is obvious. Hence we omit it.From the first three equations, we have(4)dS~dτ+dxdτ+dydτ=1-S~-x-y,which shows S~+x+y=(S~(0)+x(0)+y(0))e-τ+(1-e-τ) and we obtain (5)limt→∞S~+x+y=1.From system (3) and (5), we obtain that the dynamical behavior of system (2) can be determined by the following system: (6)dxdτ=xA1-q1-x-y-1,dydτ=yB1-x-y-1+Aq1-x-yx,y<H,Δx=-E1x,Δy=-E2y,y=H,where A=μ1S0/D, B=μ2S0/D.Definition 2 (see [18]).
LetM denote the impulsive set and N be the phase set. Suppose that g:N→N is a mapping. For any point P∈N, there exists a t1>0 such that F(P)=f(P,t1)=P1∈M, P1+=φ(P1)∈N; then g(P)=l(P1+)-l(P) is called the successor function of point P and the point P1+ is called the successor point of P.


## 3. Qualitative Analysis

We firstly consider the qualitative property of (6) without the impulsive effect. (7)dxdτ=xA1-q1-x-y-1=Px,y,dydτ=yB1-x-y-1+Aq1-x-yx=Qx,y.System (7) has three equilibria: E0(0,0), E1(0,B-1/B), and a positive equilibrium E2(x∗,y∗), (8)x∗=A1-q-1A1-q-BA1-qA-B,y∗=Aqx∗A1-q-B,where A(1-q)>max⁡B,1.Lemma 3.
System (7) has no periodic trajectories in the positive quadrant.Proof.
We construct the Dulac functionF(x,y)=1/xy, and then we have (9)∂PFx,y∂x+∂FQx,y∂y=A1-qy-Bx-Aq1-xy2<0.According to the Bendixson-Dulac Theorem [19], there is no periodic orbit in the positive quadrant.
Next, we discuss the stability of the equilibria. The Jacobian matrix evaluated at the pointE0(0,0) is (10)A1-q0AqB-1and the eigenvalues are λ1=A(1-q),λ2=B-1. Hence, the equilibrium E0(0,0) is saddle for B>1.Theorem 4.
The equilibriumE1(0,B-1/B) is globally asymptotically stable if B>1.Proof.
We can compute Jacobian matrix of the equilibriumE1(0,B-1/B) as follows: (11)1-qB-101-B+AqB1-Band the eigenvalues are λ3=(1-q)/B-1, λ4=1-B. Obviously, the equilibrium E1(0,B-1/B) is asymptotically stable if B>1. Again from Lemma 3, the equilibrium E1(0,B-1/B) is globally asymptotically stable for B>1.Theorem 5.
The positive equilibriumE2(x∗,y∗) is globally asymptotically stable if A(1-q)>max⁡B,1.Proof.
Firstly, we show the positive equilibriumE2(x∗,y∗) is locally asymptotically stable. Jacobian matrix of the equilibrium E2(x∗,y∗) is given as follows: (12)JE2x∗,y∗=-A1-qx∗-A1-qx∗-By∗+Aq1-x∗-y∗-Aqx∗B1-x∗-y∗-By∗-Aqx∗.We perform a basic elementary column transformation for the matrix J(E2(x∗,y∗)) and obtain (13)J′E2x∗,y∗=-A1-qx∗0-By∗+Aq1-x∗-y∗-Aqx∗B-Aq1-x∗-y∗-1.Hence, the eigenvalues of matrix J(E2(x∗,y∗)) have the same eigenvalues as the matrix J′(E2(x∗,y∗)), and the eigenvalues are λ5=-A(1-q)x∗<0, λ6=(B-Aq)(1-x∗-y∗)-1=B-A/A(1-q)<0. In view of Lemma 3, the equilibrium E2(x∗,y∗) is globally asymptotically stable if (14)A1-q>max⁡B,1.


## 4. Existence and Stability of the Order-1 Periodic Solution

Case 1 (existence and stability of semitrivial periodic solution).
Supposex=0, and we have the following subsystem: (15)dydτ=yB-1-By,y<H,Δy=-E2y,y=H.System (15) possesses a periodic solution η(t) defined on t∈(nτ,(n+1)τ] for y(0+)=(1-E2)H, where (16)η2t=B-1H1-E2eB-1t-nTB-1-B1-E2H+B1-E2HeB-1t-nτ.If we let (17)T=1B-1ln⁡B-1-B1-E2HB-1-BH1-E2,(18)yT=hand y(T+)=(1-E2)h, which shows that system (6) has the following semitrivial periodic solution (ξ(t),η(t)), where (19)ξt=0,ηt=B-1H1-E2eB-1t-nτB-1-B1-E2H+B1-E2HeB-1t-nτfor t∈(nτ,(n+1)τ]. Therefore, we obtain that (0,η(t)) is the semitrivial periodic solution of system (6).Next we give one lemma firstly to discuss the stability of this semitrivial periodic solution of system (6).Lemma 6.
Theτ-periodic solution x(t)=ξ(t),y(t)=η(t) of the system (20)dxdt=Px,y,dydt=Qx,y,ϕx,y≠0,Δx=αx,y,Δy=βx,y,ϕx,y=0is orbitally asymptotically stable if the Floquet multiplier μ2 satisfies the condition |μ2|<1, where (21)μ2=Πk=1qΔkexp⁡∫0τ∂P∂xξt,ηt+∂Q∂yξt,ηtdt,with (22)Δk=P+∂β/∂y∂ϕ/∂x-∂β/∂x∂ϕ/∂y+∂ϕ/∂x+Q+∂α/∂x∂ϕ/∂y-∂α/∂y∂ϕ/∂x+∂ϕ/∂yP∂ϕ/∂x+Q∂ϕ/∂yand P,Q,∂α/∂x, ∂α/∂y,∂β/∂x, ∂β/∂y,∂ϕ/∂x, and ∂ϕ/∂y are calculated at the point (ξ(τk),η(τk)), P+=P(ξ(τk+),η(τk+)), Q+=Q(ξ(τk+),η(τk+)).ϕ(x,y) is a sufficiently smooth function with grad ϕ(x,y)≠0 and τk(k∈N) is the time of kth jump.The proof of this lemma is referred to in [20].Theorem 7.
The semitrivial periodic solution(0,η(t)) of system (6) is orbitally asymptotically stable if 1-Φ-1<E1, where (23)Φ=1-E2B-1-BH1-E2B-1-BH×1-B-BH1-E2A1-q-2B/BB-1-BH1-E2B-1-BH1-E2A1-q+B-2/B-1.Proof.
SinceP(x,y)=x(A(1-q)(1-x-y)-1) and Q(x,y)=y(B(1-x-y)-1)+Aq(1-x-y)x, α(x,y)=-E1x, β(x,y)=-E2y, ϕ(x,y)=y-H. Using Lemma 6, we obtain (24)∂P∂x=A1-q1-x-y-1-A1-qx,∂Q∂x2=B1-x-y-1-By,∂α∂x=-E1,∂α∂y=0,∂β∂x=0,∂β∂y=-E2,∂ϕ∂x=0,∂ϕ∂y=1.Δ1=P+∂β/∂y∂ϕ/∂x-∂β/∂x∂ϕ/∂y+∂ϕ/∂x+Q+∂α/∂x∂ϕ/∂y-∂α/∂y∂ϕ/∂x+∂ϕ/∂yP∂ϕ/∂x+Q∂ϕ/∂y=Q+ξT+,ηT+1-E1QξT,ηT=1-E11-E2B-1-BH1-E2B-1-BHand (25)exp∫0T∂P∂xξt,ηt+∂Q∂yξt,ηtdt=exp∫0TA1-q+B-2-A1-q-2Bηtdt=1-B-BH1-E2A1-q-2B/BB-1-BH1-E2B-1-BH1-E2A1-q+B-2/B-1.Hence (26)μ=Πk=1nΔk×exp∫0T∂P∂xξt,ηt+∂Q∂yξt,ηtdt=1-E11-E2B-1-BH1-E2B-1-BH×1-B-BH1-E2A1-q-2B/BB-1-BH1-E2B-1-BH1-E2A1-q+B-2/B-1=1-E1Φ.Therefore, if |μ|<1 holds, the semitrivial periodic solution (0,η(t)) is orbitally asymptotically stable.Remark. Defining E1∗=Φ-1-1, we obtain that a fold bifurcation occurs at E1=E1∗ due to |μ|=1, and a positive periodic solution may emerge for E1<E1∗.Case 2 (existence and stability of the positive periodic solution).
In this section, we denote the impulsive setM=x,y∈R+2∣x≥0,y=H, the impulsive function φ(x,H)=((1-E1)x,(1-E2)H), and the phase set N=x,y∈R+2∣x≥0,y=1-E2H. In the following, we will prove the existence and stability of the order-1 periodic solution of system (6).
Obviously, the trajectories starting from the regiony>y∗ will tend to the positive equilibrium E(x∗,y∗) after impulsive effect of at most finite times.Theorem 8.
IfH<y∗, then system (6) has an order-1 periodic solution.Proof.
Suppose the impulsive setM intersects the y-axis at the point E. The phase set N intersects the x-axis at the point F and the isoclinal line dy/dt=0 at the point C. Suppose the point A is close enough to the point F and the trajectory starting from the point A intersects the impulsive set M at the point B and reaches the point A1 due to the impulsive effect Δx=-E1x, Δy=-E2y. The point A1 is the successor point of A. Thus the successor function of A satisfies f(A)=xA1-xA>0. Similarly, the trajectory from the point C inevitably intersects the impulsive set at the point D. The point D is mapped into the point C1 after pulses. Furthermore, point C1 is surely on the left of point C from the property of the vector field; hence the successor function of the point C becomes f(C)=xC1-xC<0. According to the continuity of the successor function, we obtain that there exists a point G such that fG=0; that is, system (6) presents an order-1 periodic solution (see Figure 1).
In the following, we will investigate the stability of the order-1 periodic solution by means of stability of the limit cycle of continuous dynamic system. We consider the general continuous system as follows:(27)dxdt=Px,y,dydt=Qx,y,where functions P(x,y) and Q(x,y) are continuous and derivative about x,y.Figure 1
The existence of order-1 periodic solution of system (6) for H<y∗.The stability of the limit cycle is investigated by Poincaré map [21] and we give some theorems as follows.Lemma 9 (see [21]).
Suppose the functions¯=f(s) is a continuous map from the line segment N to itself and s=0 is a fixed point of the continuous map f(s). Then the fixed point s=0 is stable(unstable) if the part near origin of curve s¯=f(s) on the plane (s,s¯) lies in the interior of the domain s¯/s≤1-ε(≥1+ε),ε>0.Lemma 10 (see [21]).
Supposes¯=f(s) is continuous and derivative at the neighborhood of the point s=0; then the fixed point s=0 is stable (unstable) for ds¯/ds<1(>1).The detailed proofs of Lemmas9 and 10 can be seen in [21].Since the order-1 periodic solution is not a closed curve, we will make good use of the method of the square approximation to investigate the stability of the order-1 periodic solution.Theorem 11.
Supposel is the length of the closed orbit Γ¯ of system (6) and the closed order-1 periodic limit Γ¯ is stable for ∫0lH(s)ds<1.Proof.
LetΓABC denote the order-1 periodic solution of system (6) (see Figure 2(a)). The impulsive orbit ΓABC intersects the impulsive set M at point A and then jumps into the phase set N, and intersection point is B. The closed orbit including the curve ΓBCA and line segment AB is a closed order-1 periodic limit, which is denoted as Γ¯(see Figure 2(a)). Draw normal line n passing through B∈Γ and establish coordinate system (s,n) on point B. Choose any point D∈U(A,ε),ε>0 and the trajectory starting from the point D intersects vertically n-axis at point Bk and then intersects impulsive set M at the point F. The point F is mapped into the point E due to impulsive effect. The trajectory passing through point E intersects vertically n-axis at Bk+1 as t increases again.
Suppose the rectangular coordinate ofB is (φ(s),ψ(s)); we can obtain the relation of the point Bk between the rectangular coordinates (x,y) and curvilinear coordinates (s,n) by using the geometric method. (28)x=φs-nψ′s,y=ψs+nφ′s,where φ′(s)=dx/dsB=P0/P02+Q02, ψ′(s)=dy/dsB=Q0/P02+Q02. P0 and Q0 are values of P,Q at the point B; that is, P0=P(φ(s),ψ(s)), Q0=Q(φ(s),ψ(s)).
From (28), we deduce that dy/dx=ψ′(s)+φ′sdn/ds+nφ′′(s)/φ′s-ψ′sdn/ds-nψ′′s = Qφs-nψ′s,ψs+nφ′s/Pφs-nψ′s,ψs+nφ′s,
and we obtain(29)dnds=Qφ′-Pψ′-nPφ′′+Qψ′′Pφ′+Qψ′≜Fs,n.It is easy to obtain that n=0 is a special solution of system (29). Since functions P and Q are continuous and derivative, F(s,n) presents continuous first-order partial derivative about the parameter n. Thus, the Taylor expansion of F(s,n) is given as follows: (30)dnds=Fn′s,nn=0·n+on,where Fn′(s,n)=d/dn(F(s,n)). By computing, we have (31)Fn′s,nn=0=P02Qy0-P0Q0Py0+Qx0+Q02Px0P02+Q023/2≜Hs,where Py0,Px0,Qy0, and Qx0 denote partial derivatives of P,Q when n=0.
Therefore, (29) becomes dn/ds=H(s)n, and we obtain n(s)=n0exp∫0lHsds, n(0)=n0, where l is the length of the curve.
In Figure2(a), n0 represents the ordinate of Bk, and n denotes the ordinate of Bk+1. The Poincaré map will map Bk into Bk+1 and the successor function f(n0)=n. According to Lemma 9 and Lemma 10, if n/n0<1, the curve BkEFBK+1 starting from the neighborhood of point B will tend to the closed order-1 periodic limit Γ¯ABC, as k→∞. From Lemma 9 and Lemma 10, we can obtain that the closed order-1 periodic limit is stable. The proof is completed.Figure 2
The stability of the order-1 periodic solution of system (3) for x2<x2∗. (b) Γn is the order-1 periodic solution of system (32).
(a) (b)According to the arc length formulads=P02+Q02dt we obtain H(s)=∫0τPx0+Qy0dt-1/2∫0τd/dt[lnP02+Q02)]dt.Corollary 12 (see [21]).
Suppose the integral along the order-1 periodic solutionΥ satisfies H(s)<0; then the order-1 periodic solution Υ is stable.
SupposeΥ is the τ-periodic solution of the continuous dynamic system (27); we obtain that the value of the integral equals zero along the periodic solution of the continuous dynamic system (27); that is, JΥ=∫0τ1/2d/dt[ln⁡P02+Q02]dt=0. Next we will prove that the value of integral along the closed order-1 periodic limit Γ¯ has similar result. Denote F(x,y)=1/2d/dtln⁡P02+Q02.Lemma 13.
Assume the functionF(x,y) is continuous and differentiable; then the integral along the closed order-1 periodic limit Γ of system (6) satisfies ∫0τdF(x,y)/dt=0, where τ is the period of the closed order-1 periodic limit Γ¯.Proof.
SupposeΓ¯ is a closed order-1 periodic limit of system (6). ΓBCA^ denotes the curve of system (6) from the point B(xB,yB) to the point A(xA,yA), where ΓBCA^(0)=B, SBCA^(T)=A. ΓAB denotes the line segment AB¯. Make a time transformation τ~=n-1/nτ, and system (6) becomes(32)dxdτ~=xA1-q1-x-y-1,dydτ~=yB1-x-y-1+Aq1-x-yx,x2<H,Δx=-E1x,Δy=-E1y,x2=H.Let Γn denote the closed order-1 periodic limit of system (32). Since system (32) is topologically equivalent to system (6), the stability of the closed order-1 periodic limit Γn is similar to Γ¯.
SupposeΓn¯BCA^ denotes the curve of system (32) from the point B(xB,yB) to the point A(xA,yA), where Γn¯BCA^(0)=B, Γn¯BCA^(n-1/nτ)=A. Γn¯AB denotes the line segment AB¯. Equation of the line segment AB is y=yA-yB/xA-xB(x-xA)+yA, xB≤x≤xA. Obviously, system (32) is the square approximation of system (6), which implies Γn→Γ as t→∞. Thus (33)∮0τdFx,ydtdt=∫0n-1/nτdFx,ydtdt=12∫BCAdln⁡P02+Q02dtdt+12∫ABdln⁡P02+Q02dtdt→0,n→∞, which shows ∫0τdF(x,y)/dtdt=0.Therefore, we have the following result.Theorem 14.
Assume the integral along the order-1 periodic solutionΓ of system (6) has (34)∫0TPx0+Qy0dt<0;then the order-1 periodic solution Γ is stable.Theorem 15.
The order-1 periodic solution of system (6) is stable.Proof.
From the first two equations, we have(35)∂P∂x+∂Q∂y=A1-q1-x-y-1-A1-qx+B1-x-y-1-Aqx.Sign of system (35) is uncertainty and we make an equivalent transformation v(x,y)=1/xy using the Dulac’s theorem [19]; system (6) is topologically equivalent to the following system: (36)dxdτ~=A1-q1-x-y-1y=P1x,y,dydτ~=B1-x-y-1x+Aq1-x-yy=Q1x,y.We can compute ∂P1/∂x+∂Q1/∂y=-(A(1-q)/y+B/x+Aq(1-x)/y2)<0. Using Theorem 14, we obtain that the order-1 periodic solution is stable.


## 5. Discussion

The ultimate goal of the genetic alteration produces good desired products by introducing a recombinant DNA into the cell in the form of a plasmid. However, the plasmid is easily lost in the reproductive process, which will bring some negative effects on the desired products. Therefore, how to effectively control the concentration of the plasmid-free organism is important to improve the efficiency of the desired products. In this paper, we formulate a mathematical model of competition between plasmid-bearing and plasmid-free organism in order to answer the question of how to reduce the negative effect of the plasmid-free organism on the plasmid-bearing organism by means of the impulsive state feedback control. The positive equilibriumE(x∗,y∗) of system (7) is a globally asymptotically stable, which is simulated in Figure 3 with the parameters A=12,q=0.4,B=3.5, E(x∗,y∗)=(0.3748366014,0.4862745100). From the practical view of production, plasmid-free organism will bring a loss of production by consuming the nutrient. Figure 3 shows that the concentration of the plasmid-bearing organism reaches the maximum value x∗=0.3748366014, while the concentration of the plasmid-free organism still increases and reaches the maximum value y∗=0.4862745100. Hence, the impulsive state feedback control must be carried out to prohibit the rapid growth of the plasmid-free organism. In Theorem 7, the semitrivial periodic solution (0,η(t)) of system (6) is orbitally asymptotically stable if 1-Φ-1<E1 (see Figure 4) and the parameters are given A=5,q=0.3,B=3.5,H=0.4,E1=0.3,E2=0.6. Theorem 7 and Figure 4 show that the plasmid-bearing organism will become extinct if the extraction ratio of the plasmid-bearing organism is higher than E∗=1-Φ-1. Hence, we should try to reduce the effect on the plasmid-bearing organism when we implement the impulsive state feedback control on the plasmid-free organism.Figure 3
Time series and phase portrait of the globally asymptotical stability of the positive equilibrium(x∗,y∗) with the parameters A=12,q=0.4,B=3.5.
(a) (b) (c)Figure 4
Time series and phase portrait of the semitrivial periodic solution of system (6) with the parameters A=5, q=0.3, B=3.5, H=0.4, E1=0.3, E2=0.6.
(a) (b) (c)To effectively inhibit the negative effect of the plasmid-free organisms on the plasmid-bearing organism, the existence and stability of the order-1 periodic solution are given in Theorem8 and Theorem 15. Let the parameters A=12,q=0.4,B=3.5 be the same as Figure 3 and E1=0.01,E3=0.6. In view of Figure 3, we choose the critical threshold value H=0.3748366014 and simulate existence of the order-1 periodic solution in Figure 5. We show that the impulsive state feedback control is carried out once the growth rate of the plasmid-bearing organism reaches the maximum value and the negative effect of the plasmid-free organism is minimal. From Figure 5, we see that the value of the plasmid-bearing organism is 0.45≤x≤0.58 and the value of the plasmid-bearing organism is 0.15≤y≤0.42 under the condition of the impulsive state feedback control. Comparing Figure 3 with Figure 5, we find that the impulsive control is more effective in maximizing the concentration of plasmid-bearing organism.Figure 5
Time series and phase portrait of the order-1 periodic solution of system (6) with the parameters A=12,q=0.4,B=3.5, H=0.3748366014, E1=0.01,E2=0.6.
(a) (b) (c)



---

*Source: 6401059-2018-10-23.xml*