
Citation: Qichao Ying, Jingzhi Lin, Zhenxing Qian, Haisheng Xu, Xinpeng Zhang. Robust digital watermarking for color images in combined DFT and DT-CWT domains[J]. Mathematical Biosciences and Engineering, 2019, 16(5): 4788-4801. doi: 10.3934/mbe.2019241
[1] | Luis Almeida, Michel Duprez, Yannick Privat, Nicolas Vauchelet . Mosquito population control strategies for fighting against arboviruses. Mathematical Biosciences and Engineering, 2019, 16(6): 6274-6297. doi: 10.3934/mbe.2019313 |
[2] | Rajivganthi Chinnathambi, Fathalla A. Rihan . Analysis and control of Aedes Aegypti mosquitoes using sterile-insect techniques with Wolbachia. Mathematical Biosciences and Engineering, 2022, 19(11): 11154-11171. doi: 10.3934/mbe.2022520 |
[3] | Chen Liang, Hai-Feng Huo, Hong Xiang . Modelling mosquito population suppression based on competition system with strong and weak Allee effect. Mathematical Biosciences and Engineering, 2024, 21(4): 5227-5249. doi: 10.3934/mbe.2024231 |
[4] | Yang Li, Jia Li . Stage-structured discrete-time models for interacting wild and sterile mosquitoes with beverton-holt survivability. Mathematical Biosciences and Engineering, 2019, 16(2): 572-602. doi: 10.3934/mbe.2019028 |
[5] | Arti Mishra, Benjamin Ambrosio, Sunita Gakkhar, M. A. Aziz-Alaoui . A network model for control of dengue epidemic using sterile insect technique. Mathematical Biosciences and Engineering, 2018, 15(2): 441-460. doi: 10.3934/mbe.2018020 |
[6] | Luis F. Gordillo . Optimal sterile insect release for area-wide integrated pest management in a density regulated pest population. Mathematical Biosciences and Engineering, 2014, 11(3): 511-521. doi: 10.3934/mbe.2014.11.511 |
[7] | Liming Cai, Shangbing Ai, Guihong Fan . Dynamics of delayed mosquitoes populations models with two different strategies of releasing sterile mosquitoes. Mathematical Biosciences and Engineering, 2018, 15(5): 1181-1202. doi: 10.3934/mbe.2018054 |
[8] | Mingzhan Huang, Shouzong Liu, Xinyu Song . Study of the sterile insect release technique for a two-sex mosquito population model. Mathematical Biosciences and Engineering, 2021, 18(2): 1314-1339. doi: 10.3934/mbe.2021069 |
[9] | Mugen Huang, Zifeng Wang, Zixin Nie . A stage structured model for mosquito suppression with immigration. Mathematical Biosciences and Engineering, 2024, 21(11): 7454-7479. doi: 10.3934/mbe.2024328 |
[10] | Zhenzhen Shi, Huidong Cheng, Yu Liu, Yanhui Wang . Optimization of an integrated feedback control for a pest management predator-prey model. Mathematical Biosciences and Engineering, 2019, 16(6): 7963-7981. doi: 10.3934/mbe.2019401 |
The sterile insect technique (SIT) is presently one of the most ecological methods for controlling insect pests responsible for disease transmission or crop destruction worldwide. This technique consists in releasing sterile males into the insect pest population [1,2,3]. This approach aims at reducing fertility and, consequently, reducing the target insect population after a few generations. Classical SIT has been modeled and studied theoretically in a large number of papers to derive results to study the success of these strategies using discrete, continuous, or hybrid modeling approaches (for instance, the recent papers[4,5,6,7,8,9,10]).
Despite this extensive research, little has been done concerning the stabilization of the target population near extinction after the decay caused by the massive initial SIT intervention, and there are still major difficulties due to the complexity of the dependency on climate, landscape, and many other parameters which would be difficult to be integrated into the mathematical models studied. Not being able to consider all these parameters in our mathematical models and knowing that these external factors strongly impact the evolution of the density of the target population, we focus our studies on releases that now depend on the target population density measurements since, as we will see below, this makes our control more robust. Indeed, several monitoring tools can provide information on the size of the wild population throughout the year. So, a control that considers this information to adapt the size of the releases is possible and useful. This was already the case of [9,11] in which a state feedback control law gives significant robustness qualities to the mathematical model of SIT. Although this approach provides evidence in terms of robustness because the control is directly adjusted according to the density of the population, its application requires to continuously measure the different states of the model. In practice, traps allow data to be collected to analyze the control's impact and technology is being developed that may allow us to obtain continuous data in the near future.
However, specific categories of data are still problematic or very expensive to obtain. For example, during an SIT intervention, it is difficult to measure the density of young females that have not yet been fecundated or of females that were fecundated by wild males. In this work we use another control theory tool, which consists of constructing a state estimator for a dynamical system and using this estimator to apply feedback control. A state observer or state estimator is a system that provides an estimate of the natural state using some partial measurements of the real system. In our case, using traps, wild males as well as sterile males, can be measured. Using the observer system technique, we have built a system that allows us to estimate all other states.
The problem of observer design for linear systems was established and solved by [12] and [13]. While Kalman's observer [12] was highly successful for linear systems, extending it to nonlinear systems took a lot of work. In several cases, the observer can be obtained from the extended Kalman filter by a particular choice of the matrix gain using linear matrix inequalities (LMIs). The development of the observer in this paper was motivated by its application to the SIT model. A model of this process can be written as
˙x=Ax+B(y)x+Du, | (1.1) |
y=Cx, | (1.2) |
where y∈Rm is the output, x∈Rn is the state vector, and u∈Rp is the input. The output matrix B(y) is such that the coefficients b(y)ij are bounded for all i,j.
Our paper has three parts. In the first part, thanks to the backstepping approach, we build a feedback control law that stabilizes the zero population state for the SIT model for the mosquito population, which considers only the compartments of young females and fertilized females presented in [14]. In the second part we construct a state estimator for the SIT model. Finally, in the third part we show that the application of this feedback, depending on the measured states and the ones estimated thanks to the state estimator, globally stabilizes the system.
The mosquito life cycle has several phases. The aquatic stage comprises eggs, larvae, and pupa, followed by the adult stage, where we consider both wild males and females. After emergence from the pupa, a female mosquito needs to mate and then to get a blood meal before it can start laying eggs. Then, every 4−5 days, it will take a blood meal and lay 100−150 eggs at different places (10−15 per place). For the mathematical description, we will consider the following compartments [14].
● E the density of population in aquatic stage,
● Y the density of young females, not yet laying eggs,
● F the density of fertilized and egg-laying females,
● M the density of males,
● Ms the density of sterile males,
● U the density of females that mate with sterile males.
The Y compartment represents the stage of the young females before the start of their gonotropic cycle, i.e., before they mate and take their first blood meal. It generally lasts for 3 to 4 days. The sterile insect technique introduces male mosquitoes to compete with wild males. We denote by Ms the density of sterile mosquitoes and by U the density of females that have mated with them. We assume that a female mating mosquito has probability MM+Ms to mate with a wild male and probability MsM+Ms to mate with a sterile one. Hence, the transfer rate η from the compartment Y splits into transfer rate of η1MM+Ms to compartment F and a transfer rate of η2MsM+Ms to compartment U of females that will be laying sterile (non-hatching) eggs. The mathematical model is the system of ordinary differential equations presented in [15]
˙E=βEF(1−EK)−(δE+νE)E, | (2.1) |
˙M=(1−ν)νEE−δMM, | (2.2) |
˙Y=ννEE−η1MM+MsY−η2MsM+MsY−δYY, | (2.3) |
˙F=η1MM+MsY−δFF, | (2.4) |
˙U=η2MsM+MsY−δUU, | (2.5) |
˙Ms=u−δsMs. | (2.6) |
The parameter δY is the mortality rate for young females (they can die without mating for diverse reasons like predators or other hostile environmental conditions). Male mosquitoes can mate for most of their lives. A female mosquito needs a successful mating to be able to reproduce for the rest of her life, βE>0 is the oviposition rate; δE,δM,δF,δY,δs>0 are the death rates, respectively, for eggs, wild adult males, fertilized females, young females, and sterile males; νE>0 is the hatching rate for eggs; ν∈(0,1), the probability that a pupa gives rise to a female, and (1−ν) is, therefore, the probability of giving rise to a male. K>0 is the environmental capacity for eggs. It can be interpreted as the maximum density of eggs that females can lay in breeding sites. Since here the larval and pupal compartments are not present, we consider that E represents all the aquatic compartments, in which case, this term K represents a logistic law's carrying capacity for the aquatic phase, which also includes the effects of competition between larvae. The control function u represents the number of mosquitoes released during the SIT intervention. It is interesting to follow the evolution of the state U because female mosquitoes, once fertilized by sterile males, will continue their gonotrophic cycle normally and, therefore, can still transmit disease. We will assume in this work that
δs≥δM. | (2.7) |
In [14,15], equilibria and their stability property were studied for the system without control.
˙E=βEF(1−EK)−(δE+νE)E, | (2.8) |
˙M=(1−ν)νEE−δMM, | (2.9) |
˙Y=ννEE−(η1+δY)Y, | (2.10) |
˙F=η1Y−δFF. | (2.11) |
Its basic offspring number is R0=η1βEννEδF(νE+δE)(η1+δY). For the rest of our work, we assume that
R0>1. | (2.12) |
We assume that wild males are more likely to fertilize young females because they are born in the same egg-laying site. We define
Δη=η1−η2≥0. | (3.1) |
Other authors, such as in [14], have already studied the stability of this type of model. The difference in our approach lies in the kind of control used initially for global stabilization. Indeed, in most of the prior studies, the controls u studied were independent of system states. Some previous works have considered certain simple applications of feedback control to SIT (see, for instance, [9,16,17]). In a previous paper [11], we used the backstepping method to build a feedback control system that simplifies the SIT model, which is presented in [4], assuming that all females are immediately fertilized. Here, we consider the system
˙E=βEF(1−EK)−(δE+νE)E, | (3.2) |
˙M=(1−ν)νEE−δMM, | (3.3) |
˙Y=ννEE−ΔηMM+MsY−(η2+δY)Y, | (3.4) |
˙F=η1MM+MsY−δFF, | (3.5) |
˙U=η2MsM+MsY−δUU, | (3.6) |
˙Ms=u−δsMs. | (3.7) |
Let N:=[0,+∞)6 and X:=(E,M,Y,F,U,Ms)T. When applying a feedback law u:N→[0,+∞), the closed-loop system is the system
˙X=H(X,u(X)), | (3.8) |
where H is the righthand side of Eqs (3.2)–(3.7). The construction method remains the same as in our previous paper [11]. In this work, we also consider solutions in the Filippov sense of our discontinuous closed-loop system (see, for instance, [18,19,20,21,22,23]). Let us define x:=(E,M,Y,F,U)T. We must rewrite the target system (3.2)–(3.6) in the following form to apply the backstepping method (see, for instance, [35]HY__HY, Theorem 12.24, page 334]):
{˙x=f(x,Ms),˙Ms=u−δsMs, | (3.9) |
where f:R6→R5 represents the righthand side of (3.2)–(3.6). We then consider the control system ˙x=f(x,Ms) with the state being x∈D:=[0,+∞)5 and the control being Ms∈[0,+∞). We assume that Ms is of the form Ms=θM for a constant θ>0. Then, we define and study the closed-loop system
˙x=f(x,θM). | (3.10) |
Its offspring number is
R(θ):=βEη1ννEδF(νE+δE)(Δη+(1+θ)(η2+δY)). | (3.11) |
Note that if R(θ)≤1, 0∈R5 is the only equilibrium point of the system in D. Our next proposition shows that the feedback law Ms=θM stabilizes our control system ˙x=f(x,Ms) if R(θ)<1.
Proposition 3.1. Assume that
R(θ)<1. | (3.12) |
Then, 00 is globally exponentially stable in D for system (3.10). The exponential convergence rate is bounded from above by the positive constant c defined by relation (3.16).
Proof. We apply Lyapunov's second theorem. To do so, we define V:[0,+∞)5→R+, x↦V(x),
V(x):=(1+2R(θ))ννE(νE+δE)(1−R(θ))E+νM+3R(θ)(1−R(θ))Y+(2+R(θ))βEννEδF(νE+δE)(1−R(θ))F+σU, | (3.13) |
where σ>0 is a constant that we will choose later.
As (3.12) holds, V is of class C1, V(x)>V((0,0,0,0,0)T)=0,∀x∈[0,+∞)5∖{(0,0,0,0,0)T}, V(x)→+∞ when ‖x‖→+∞ with x∈D and
˙V(x)=−βEννE(νE+δE)F−νδMM−(1+2R(θ))ννE(νE+δE)(1−R(θ))βEKFE−ν2νEE−η1βEννEδF(νE+δE)(1+θ)Y−ση21+θY+ση2Y−σδUU. |
By choosing
σ:=η1βEννER(θ)(1+θ)η2(νE+δE)δF | (3.14) |
we get
˙V(x)=−βEννE(νE+δE)F−νδMM−(1+2R(θ))ννE(νE+δE)(1−R(θ))βEKFE−ν2νEE−η1βEννE(1+θ(1−R(θ)))δF(νE+δE)(1+θ)2Y−σδUU. |
Using (3.12) once more, we get
˙V(x)≤−cV(x),∀x∈[0,+∞)5, | (3.15) |
with
c:=min{ν(νE+δE)(1−R(θ))(1+2R(θ)),δM,δF(1−R(θ))2+R(θ),η1βEννE(1+θ(1−R(θ)))δF(νE+δE)(1+θ)2(1−R(θ))3R(θ),δU}>0. | (3.16) |
This concludes the proof of Proposition 3.1.
Remark 3.1. When the Allee effect is included in the model (for instance, [4,Eq (2.5),Page 25]), the control Ms=θM can still be used, and the proof of the stability result can still be done using the same Lyapunov function (3.13).
We define
ϕ:=(2+R(θ))η1βEννE−3R(θ)ΔηδF(νE+δE)δF(νE+δE)(1−R(θ))(1+θ)−η1βEννER(θ)(1+θ)2(δE+νE)δF, | (3.17) |
Q:=3(η2+δY)(1+θ)(νE+δE)δF−(1−R(θ))η1βEννE, | (3.18) |
and for α>0, the map G:N:=[0,+∞)6→R, (xT,Ms)T↦G((xT,Ms)T) by
G((xT,Ms)T):=ϕY(θM+Ms)2α(M+Ms)(3θM+Ms)+((1−ν)νEθE−θδMM)(θM+3Ms)3θM+Ms+δsMs+1α(θM−Ms) if M+Ms≠0, | (3.19) |
G((xT,Ms)T):=0 if M+Ms=0. | (3.20) |
Finally, let us define the feedback law u:N→[0,+∞), (xT,Ms)T↦u((xT,Ms)T), by
u((xT,Ms)T):=max(0,G((xT,Ms)T)). | (3.21) |
The global stability result is the following.
Theorem 3.1. Assume that (3.12) holds. Then, 00∈N is globally exponentially stable in N for system (3.2)–(3.6) with the feedback law (3.21). The exponential convergence rate is bounded by the positive constant cp defined by
cp:=min{c,1α,δM,Q3(1+θ)δF(νE+δE),δU}. | (3.22) |
Lemma 3.1. Assume that (2.13) and (3.12) hold, then ϕ>0.
Proof.
Let us define ϕ1:=(2+R(θ))η1βEννE−3R(θ)ΔηδF(νE+δE)δF(νE+δE)(1−R(θ))(1+θ). We get from the relation (2.13) that η1βEννE>δF(νE+δE)(η1+δY). So,
ϕ1>2η1(1+θ)+(2+R(θ))δY+3R(θ)η2(1−R(θ))(1+θ). | (3.23) |
From relation (3.12), we get βEη1ννEδF(νE+δE)<Δη+(1+θ)(η2+δY). Thus,
ϕ>2η1(1+θ)+(2+R(θ))δY+3R(θ)η2(1−R(θ))(1+θ)−ΔηR(θ)+(1+θ)(η2+δY)R(θ)(1+θ)2,>2η1R(θ)(1+θ)+(2+R(θ))δY+3R(θ)η2(1−R(θ))(1+θ)+η2R(θ)(1+θ)2−η1R(θ)(1+θ)2−(η2+δY)R(θ)1+θ,>η1R(θ)(1+2θ)(1+θ)2+2R(θ)η2+2δY+R(θ)2(η2+δY)(1−R(θ))(1+θ)+η2R(θ)(1+θ)2,>0. |
Proof of Theorem 3.1. Let α>0 and define W:N→R by
W((xT,Ms)T):=V(x)+α(θM−Ms)2θM+Ms if M+Ms≠0, | (3.24) |
W((xT,Ms)T):=V(x) if M+Ms=0. | (3.25) |
We have
W is continuous, | (3.26) |
W is of classC1onN∖{(E,M,Y,F,U,Ms)T∈N;M+Ms=0}, | (3.27) |
W((xT,Ms)T)→+∞as‖x‖+Ms→+∞, withx∈DandMs∈[0,+∞), | (3.28) |
W((xT,Ms)T)>W(00)=0,∀(xT,Ms)T∈N∖{00}. | (3.29) |
From now on, and until the end of this proof, we assume that (xT,Ms)T is in N, and until (3.43) below, we further assume that
(M,Ms)≠(0,0). | (3.30) |
One has
˙W((xT,Ms)T)=∇V(x)T⋅f(x,Ms)+α(θM−Ms)2(θ˙M−˙Ms)(θM+Ms)−(θ˙M+˙Ms)(θM−Ms)(θM+Ms)2,=∇V(x)T⋅f(x,θM)+∇V(x)T⋅(f(x,Ms)−f(x,θM))+α(θM−Ms)θ˙M(θM+3Ms)−˙Ms(3θM+Ms)(θM+Ms)2. |
Since
∇V(x)T⋅(f(x,Ms)−f(x,θM))=((1+2R(θ))ννE(νE+δE)(1−R(θ))ν3R(θ)(1−R(θ))(2+R(θ))βEννEδF(νE+δE)(1−R(θ))η1βEννER(θ)(1+θ)η2(νE+δE)δF)⋅(00−Δη(θM−Ms)(M+Ms)(1+θ)Yη1(θM−Ms)(M+Ms)(1+θ)Y−η2(θM−Ms)(M+Ms)(1+θ)Y)=ϕY(θM−Ms)M+Ms, |
˙W((xT,Ms)T)=∇V(x)T⋅f(x,θM)+α(θM−Ms)(θM+Ms)2[(∇V(x)⋅(f((xT,Ms)T)−f(x,θM)))(θM+Ms)2α(θM−Ms)+θ˙M(θM+3Ms)−˙Ms(3θM+Ms)]=˙V(x)+α(θM−Ms)(θM+Ms)2[ϕY(θM+Ms)2α(M+Ms)+((1−ν)νEθE−θδMM)(θM+3Ms)−u(3θM+Ms)+δsMs(3θM+Ms)]. | (3.31) |
We take u as given by (3.21). Therefore, in the case where
ϕY(θM+Ms)2α(M+Ms)+((1−ν)νEθE−θδMM)(θM+3Ms)+δsMs(3θM+Ms)+1α(θM−Ms)(3θM+Ms)>0, | (3.32) |
u=13θM+Ms[ϕY(θM+Ms)2α(M+Ms)+((1−ν)νEθE−θδMM)(θM+3Ms)+δsMs(3θM+Ms)+1α(θM−Ms)(3θM+Ms)], |
which, together with (3.31), leads to
˙W((xT,Ms)T)=˙V(x)−(θM−Ms)2(3θM+Ms)(θM+Ms)2. | (3.33) |
Otherwise, i.e., if (3.32) does not hold,
ϕY(θM+Ms)2α(M+γsMs)+((1−ν)νEθE−θδMM)(θM+3Ms)+δsMs(3θM+Ms)+1α(θM−Ms)(3θM+Ms)≤0, | (3.34) |
so, by (3.21),
u=0. | (3.35) |
We consider two cases. If θM>Ms using (3.31), (3.34), and (3.35),
˙W((xT,Ms)T)≤˙V(x)−(θM−Ms)2(3θM+Ms)(θM+Ms)2,≤−cV(x)−(θM−Ms)2θM+Ms,≤−c1W((xT,Ms)T), | (3.36) |
with
c1:=min{c,1α}>0. | (3.37) |
Otherwise, if θM≤Ms, using once more (3.31) and (3.35),
˙W((xT,Ms)T)=˙V(x)+α(θM−Ms)(θM+Ms)2[ϕY(θM+Ms)2α(M+Ms)+θ((1−ν)νEE−δMM)(θM+3Ms)+δsMs(3θM+Ms)]. | (3.38) |
˙W((xT,Ms)T)=˙V(x)+α(θM−Ms)(θM+Ms)2[ϕY(θM+Ms)2α(M+Ms)+θ((1−ν)νEE]+α(θM−Ms)(θM+Ms)2[−δMM(θM+3Ms)+δsMs(3θM+Ms)]. | (3.39) |
From Lemma 3.1, we deduce that ϕ>0, and as (xT,Ms)T∈N, one has ϕY(θM+Ms)2α(M+Ms)+θ(1−ν)νEE≥0.
So,θM−Ms≤0⟹α(θM−Ms)(θM+Ms)2[ϕY(θM+Ms)2α(M+Ms)+θ(1−ν)νEE]≤0. | (3.40) |
Equation (3.39) becomes
˙W((xT,Ms)T)≤˙V(x)+α(θM−Ms)(θM+Ms)2[−θδMM(θM+3Ms)+δsMs(3θM+Ms)]. | (3.41) |
The inequality (2.7) gives δs≥δM, and one has
−θδMM(θM+3Ms)+δsMs(3θM+Ms)≥−θδMM(θM+3Ms)+δMMs(3θM+Ms)≥δM(Ms−θM)(Ms+θM). | (3.42) |
(3.42) together with θM−Ms≤0 implies that
˙W((xT,Ms)T)≤˙V(x)−αδM(θM−Ms)2(θM+Ms),≤−c2W((xT,Ms)T), | (3.43) |
with
c2:=min{c,δM}>0. | (3.44) |
Let us now deal with the case where (3.30) is not satisfied. Note that, for every τ≥0, M(τ)+Ms(τ)>0 implies that M(t)+Ms(t)>0 for all t≥τ. Thus, if M(0)+Ms(0)=0, there exists ts∈[0,+∞] such that M(t)+Ms(t)=0 if, and only if, t∈[0,ts]∖{+∞}. Let us study only the case ts∈(0,+∞) (the case ts=0 is obvious and the case ts=+∞ is a corollary of our study of the case ts∈(0,+∞)). Let us first point out that, for every (M,Ms)T∈[0,+∞)2 such that M+Ms>0, one has
MM+Ms≤1,(θM+Ms)2≤(3θM+Ms)2and(θM+Ms)2(M+Ms)(3θM+Ms)≤(3θM+Ms)M+Ms≤3θ+1,θM+3Ms3θM+Ms=θM3θM+Ms+3Ms3θM+Ms≤13+3≤4. |
So,
MM+Ms∈[0,1],(θM+Ms)2(M+Ms)(3θM+Ms)∈[0,3θ+1], and θM+3Ms3θM+Ms∈[0,4]. | (3.45) |
Let t↦X(t)=(E(t),M(t),Y(t),F(t),U(t),Ms(t))T be a solution (in the Filippov sense) of the closed-loop system (3.2)–(3.6) such that, for some ts∈(0,+∞),
M(t)+Ms(t)=0,∀t∈[0,ts]. | (3.46) |
Note that (3.46) implies that
M(t)=Ms(t)=0,∀t∈[0,ts]. | (3.47) |
From (3.45), (3.47), and the definition of a Filippov solution, one has on (0,ts)
(˙E˙M˙Y˙F˙U˙Ms)=(βEF(1−EK)−(νE+δE)E(1−ν)νEE−δMMννEE−κ(t)ΔηY−(η2+δY)Yη1Yκ(t)−δFFη2(1−κ(t))Y−δUUYg1(t)+Eg2(t)−δsMs) | (3.48) |
with
κ(t)∈[0,1],g1(t)∈ϕα[0,3θ+1] and g2(t)∈(1−ν)νEθ[0,4]. | (3.49) |
From (3.47) and the second line of (3.48), one has
E(t)=0,∀t∈[0,ts]. | (3.50) |
From the first line of (3.48) and (3.50), we get
F(t)=0,∀t∈[0,ts]. | (3.51) |
Let us first consider the case where Y(0)=0. Then, from the third line of (3.48) and (3.50), one has
Y(t)=0,∀t∈[0,ts]. | (3.52) |
To summarize, from (3.47), the fifth line of (3.48), (3.50), (3.51), and (3.52),
E(t)=M(t)=Y(t)=F(t)=Ms(t)=0 and ˙U(t)=−δUU(t),∀t∈[0,ts], | (3.53) |
which, with (3.13), (3.16), and (3.25) gives
˙W(t)=−σδUU(t)≤−δUW(t),∀t∈[0,ts]. | (3.54) |
Let us finally consider the case where Y(0)>0. Then, from the third line of (3.48),
Y(t)>0,∀t∈[0,ts], | (3.55) |
which, together with the fourth line of (3.48) and (3.51), implies
κ(t)=0,∀t∈[0,ts]. | (3.56) |
To summarize, from (3.47), the third and the fifth line of (3.48), (3.50), (3.51) and (3.56),
E(t)=M(t)=F(t)=Ms(t)=0,˙Y(t)=−(η2+δY)Y(t), and ˙U(t)=η2Y−δUU(t),∀t∈[0,ts], |
which, with (3.13), (3.14), (3.16), and (3.25) gives
˙W(t)=−(η2+δY)3R(θ)(1−R(θ))Y(t)+η2σY(t)−σδUU(t)=−R(θ)((η2+δY)3(1−R(θ))−η1βEννE(1+θ)(νE+δE)δF)Y(t)−σδUU(t)=−R(θ)(Q(1−R(θ))(1+θ)(νE+δE)δF)Y(t)−σδUU(t) | (3.57) |
where
Q:=3(η2+δY)(1+θ)(νE+δE)δF−(1−R(θ))η1βEννE. | (3.58) |
To end the proof, we have to prove that Q>0. Using the relation (3.11) and (3.12), we have
βEη1ννE<R(θ)δF(νE+δE)Δη+δF(νE+δE)(1+θ)(η2+δY). | (3.59) |
Recall that Δη=η1−η2. One has
Q=3(η2+δY)(1+θ)(νE+δE)δF−η1βEννE+R(θ)η1βEννE>2(η2+δY)(1+θ)(νE+δE)δF−R(θ)Δη(νE+δE)δF+R(θ)η1βEννE>2(η2+δY)(1+θ)(νE+δE)δF−R(θ)η1(νE+δE)δF+R(θ)η1βEννE+R(θ)η2(νE+δE)δF. |
From the relation (2.13), η1βEννE>δF(νE+δE)(η1+δY).
Q>2(η2+δY)(1+θ)(νE+δE)δF−R(θ)η1(νE+δE)δF+R(θ)δF(νE+δE)(η1+δY)+R(θ)η2(νE+δE)δF>2(η2+δY)(1+θ)(νE+δE)δF+R(θ)(η2+δY)(νE+δE)δF>0. |
We get
˙W(t)≤−c′W(t),∀t∈[0,ts], | (3.60) |
where
c′:=min{Q3(1+θ)δF(νE+δE),δU}. | (3.61) |
This proves Theorem 3.1 and gives the global exponential stability. From (3.37), (3.44), and (3.61), we obtain an estimate on the exponential decay rate
cp:=min{c,1α,δM,Q3(1+θ)δF(νE+δE),δU}. | (3.62) |
Note that η1 represents the natural fertility rate in the mosquito population. Wild males have a shorter maturity time in their life cycle than females. Thus, the fertilization phase is essentially around the hatching site. Sterile males are artificially released into the intervention region. We denote by p with (0≤p≤1) the proportion of sterile males that are released. Also, the effective fertilization during the mating could be diminished due to the sterilization, which leads us to assume that the effective mating rate of sterile insects is given by qη1 with 0≤q≤1. Putting together these assumptions, we get that the probability for a young female to mate with sterile males is η2MM+Ms with η2=pqη1. For the numerical simulation, we take η1=1 and η2=0.7. The numerical simulations of the dynamics when applying the feedback (3.21) is given in Figure 1. The parameters we use are given in the following table.
With the parameters given in Table 1, condition (3.12) is θ>102,06. We fix K=21000 and we consider the persistence equilibrium z0=(E0,M0,Y0,F0,U0,M0s) as initial condition. That gives E0=20700,M0=5300,Y0=1500,F0=13000, and U0=M0s=0. We take θ=290 and α=90.
Parameters | Description | Value |
βE | Effective fecundity | 10 |
νE | Hatching parameter | 0.05 |
δE | Mosquitoes in aquatic phase death rate | 0.03 |
δF | Fertilized female death rate | 0.04 |
δY | Young female death rate | 0.04 |
δM | Male death rate | 0.1 |
δs | Sterilized male death rate | 0.12 |
ν | Probability of emergence | 0.49 |
Remark 3.2. Note that the feedback satisfies
supε→0{|u(X)|:X∈N,‖X‖1∈B(0,ε)}⟶0. | (3.63) |
The advantage of applying feedback control is that when the density of the target population decreases, the control also decreases.
Remark 3.3. It is important to note that the backstepping feedback control (3.21) does not depend on the environmental capacity K, which is also an interesting feature for the field applications. In the case K=+∞, the Eq (2.8) becomes
˙E=βEF−(δE+νE)E, | (3.64) |
and we prove by the same process that the same feedback law (3.21) ensures the exponential stability of the SIT system (3.64), (3.3)–(3.7) with the same lower bound of the exponential convergence rate.
Our stabilization result is the following one.
Theorem 3.2. Assume that (3.12) holds and K=+∞. Then, 00∈N is globally exponentially stable in N for system (3.64), (3.3)–(3.7) with the feedback law (3.21). The exponential convergence rate is bounded by cp>0 and defined in (3.62).
Remark 3.4. Let us assume that the heterogeneity of the intervention zone strongly impacts the mating of female mosquitoes with sterile males more than we would have estimated. Suppose the estimated mating rate for the control (3.21) is ηe2=0.7, and let the mating rate be ηr2=0.4 for the dynamics. Keeping the other parameters and the same initial condition, we obtain the following figure. This parameter considerably impacts the convergence time of the states of the system. Note that with e=3×10−1 of error difference, we still have convergence. Estimation errors of the order of e=10−2 will have a negligible impact on the convergence time. This is because the backstepping control also depends on the states of the system. Thus, the states make a correction that can compensate for a certain margin of error. Unlike control, which only depends on the parameters, estimation errors have no correction from the dynamics. Therefore, this can be fatal to the success of the intervention. In practice, many external factors impact the life cycle of mosquitoes. These factors modify parameters such as birth, hatching, and fertilization rates. These factors are, for example, rainfall and the topography of the region. An SIT model that can integrate these factors is challenging to study (see [24]). The success of an SIT intervention depends strongly on the robustness of the control strategy. The results of our previous test that are reported in Figure 2 show us the advantage feedback control can provide in terms of robustness.
The application of feedback control requires measuring states such as eggs E and young females Y of the intervention zone over time. In practice, it is always important to estimate the density of adult mosquitoes to intervene in an area. This data is collected using mosquito traps distributed throughout the region. Despite various technological advances to improve these traps, it should be noted that some data is still easier to be measured than others. Measuring mosquito density in the aquatic phase E is difficult, especially in a heterogeneous area. It is also challenging to measure young females Y because females come in three categories, and we need to distinguish between unfertilized and fertilized females. Males are more easily measured because they are distinguishable. It can be also easy to distinguish wild males from laboratory males by marking processes applied to laboratory males. In this part of our paper, we will assume that the density of wild males and that of sterile males can be measured continuously. Our objective is to estimate the other densities. Observer design for nonlinear dynamic systems is a technique used in control theory to estimate the states of a system when only partial or indirect measurements are available. The difficulties in dealing with observer problems for general nonlinear systems is the proof of global convergence of the estimation error. Much literature exists on state observers and filters for nonlinear systems as they play crucial roles in control theory. To simplify the nonlinearity F(1−EK) of the SIT model, in this section we consider the simplified SIT model for environmental capacity K=+∞. On the one hand, the reason for studying such a model is that the simplified model can be considered relevant from a biological point of view within a large intervention domain or in areas where environmental capacity is difficult to estimate. On the other hand, based on the result presented in Theorem 3.2, the proposed feedback law (3.21) still stabilizes the simplified model around zero with the same convergence rate. We consider the following output control system:
˙E=βEF−(δE+νE)E, | (4.1) |
˙M=(1−ν)νEE−δMM, | (4.2) |
˙Y=ννEE−ΔηMM+MsY−(η1+δY)Y, | (4.3) |
˙F=η1MM+MsY−δFF, | (4.4) |
˙U=η2MsM+MsY−δUU, | (4.5) |
˙Ms=u−μsMs, | (4.6) |
y1=M, | (4.7) |
y2=Ms, | (4.8) |
where the states are X=(E,M,Y,F,U,Ms)T∈N, the control is u∈[0,+∞), and the output is y=(M,Ms)T∈R2+.
In particular, in this model we are confronted with a difficulty in which most observer construction theories are invalid because of the singularity at the origin. To go around this difficulty, we will use the fact that the main nonlinearity term MM+Ms is bounded and essentially the most accessible data to measure. This leads us to develop an observer for this type of system.
The usual observers for linear systems are the Luenberger observer and the Kalman observer. Observer design for a nonlinear system is a complex problem in control theory and has received much attention from many authors yielding a large literature of methods. Among than, the most famous are the change of coordinates to transform the nonlinear system into a linear system [25,26,27,28,29] and a second approach consists in using the extended Kalman filter (EKF) [30,31,32,33]. The state observer is called an exponential state observer if the observer error converges exponentially to zero. In this section we provide an explicit construction of a global observer for the following system.
{˙x(t)=Ax(t)+B(y(t))x(t)+Du(t),y(t)=Cx(t), | (4.9) |
where x(t)∈Rn is the state vector, u(t)∈Rp is the input vector, and y(t)∈Rm is the output vector. A∈Rn×n and C∈Rm×n are the appropriate matrices. The matrice B(y(t)) is in the form
B(y(t))=n,n∑i,j=1bij(y(t))en(i)eTn(j). | (4.10) |
We assume that for all y(t)∈Rm, the coefficients bij are bounded for all i=1,⋯,n and j=1,⋯,n, and denote
¯bij=maxt(bij(y(t)))andb_ij=mint(bij(y(t))). | (4.11) |
Then, the parameter vector b(t) remains in a bounded convex domain Sn,n of which 2(n2) vertices are defined by:
VSn,n={η=(η11,⋯,η1n,⋯,ηnn)|ηij∈{b_ij,¯bij}}. |
A state observer corresponding to (4.9) is given as follows:
{˙ˆx(t)=Aˆx(t)+B(y(t))ˆx+Du(t)−L(Cˆx−y(t)),ˆy(t)=Cˆx(t), | (4.12) |
where ˆx(t) denotes the estimate of the state x(t). The dynamics of the observer error e(t):=ˆx(t)−x(t) are ˙e(t)=(A−LC)e(t)+B(y(t))e(t)=(A+B(y(t))−LC)e(t). We define
A(b(t))=A+n,n∑i,j=1bij(y(t))eq(i)eTn(j). | (4.13) |
The dynamics of the observer error becomes
˙e(t)=(A(b(t))−LC)e(t). | (4.14) |
The observation problem consists in finding a gain L such that (4.14) converges exponentially toward zero. We use the following results in [34].
Theorem 4.1. The observer error converges exponentially toward zero if there exist matrices P=PT>0 and R of appropriate dimensions such that following LMIs are feasible:
AT(η)P−CTR+PA(η)−RTC+ξI<0, | (4.15) |
∀η∈VSn,n, | (4.16) |
for some constant ξ>0. When these LMIs are feasible, the observer gain L is given by L=P−1RT.
Proof. We follow [34] and consider the following quadratic Lyapunov function:
V(e)=eTPe, | (4.17) |
where P is the matrix in Theorem 4.1. We have ˙V(e)(t)=e(t)TF(b(t))e(t), where F(b(t))=(A(b(t))−LC)TP+P(A(b(t))−LC). For e(t)≠0, the condition V(e(t))>0 is satisfied because P>0, and the condition ˙V(e(t))<0 is satisfied if we have
F(b(t))<0for allb(t)∈Sn,n. | (4.18) |
Since the matrix function F is affine in b(t), using a convexity argument we deduce that ∀t≥0
˙V(e(t))<−ξ‖e(t)‖2P, | (4.19) |
if the following condition is satisfied F(η)<−ξI, ∀η∈Vn,n. Thus, if (4.15) holds, this inequality is also satisfied.
We rewrite the output SIT models (4.1)–(4.7) as
{˙X=AX+B(y)X+Du,y=CX, | (4.20) |
where X=(E,M,Y,F,U,Ms)T,
A=(−(δE+νE)000βE0(1−ν)νE−δM0000ννE0−(η2+δY)000000−δF000000−δU000000−δs),B(y)=(00000000000000−Δηy1y1+y200000η1y1y1+y200000η2y2y1+y2000000000)C=(010000000001),D=(0,0,0,0,0,1)T. |
As, N is an invariant set, one has 0≤y1y1+y2≤1. Solving the corresponding equation of (4.15) with ξ=1 in MATLAB, we get
P=104(0.0219−0.1567−0.1531−0.1703−0.03440−0.15678.9301−0.8472−0.8081−0.49290−0.1531−0.84724.57160.92771.08450−0.1703−0.80810.92774.3088−2.30120−0.0344−0.49291.0845−2.30124.74130000003.7267), | (4.21) |
R=103(0.23520.9704−0.4415−1.14010.06900000001.4162), | (4.22) |
L=(50.634201.415000.942602.654701.6023000.3800). | (4.23) |
With the parameters given in Table 1, the result for a simulation run with x0=(400,100,150,120,120,50)T, ˆx0=(120,70,70,50,60,0)T, and u=500000 is plotted in Figure 3. The asymptotic behavior of the different estimates ˆE,ˆF,ˆY, and ˆU (dashed) illustrates the exponential convergence of the estimation error shown in Theorem 4.1.
The feedback control (3.21) depends on the states E, M, Y, and Ms. From the measurement of states M and Ms, an observer system has been built in the previous section. This state observer is used to estimate both eggs E and young females Y. In this section, we show that u(ˆX,y) stabilizes the dynamics at the origin. We consider the coupled system
{˙X=f(X,ˆu(ˆX,y)),˙ˆX=f(ˆX,ˆu(ˆX,y))−L(Cˆx−y), | (5.1) |
with
ˆu(ˆX,y)=max(0,S(ˆX,y)). | (5.2) |
where S:R4×R2+→R, (ˆX,y)T↦S(ˆX,y) is defined by
S(ˆX,y):=G(ˆE,M,ˆY,Ms) | (5.3) |
The main result of this section is the following theorem.
Theorem 5.1. Assume that (3.12) holds. Then, 00∈E=N×R6 is globally exponentially stable in E for system (5.1) with the feedback law (5.2). The convergence rate is bounded by the positive constant ce defined by
ce:=min{c1,c2,c′,ξ4}. | (5.4) |
Proof. Let λ>0, and we define H:E→R by
H(X,ˆX)=W(X)+λ√V(e) | (5.5) |
with e=ˆX−X.
H is continuous onEandC1onE∖{(X,ˆX)∈E;M+Ms=0}, | (5.6) |
H(X,ˆX)→+∞as‖(X,ˆX)‖→+∞, | (5.7) |
H(X,ˆX)>H(0)=0,∀(X,ˆX)∈E∖{0}. | (5.8) |
In this proof, from now on we assume that (X,ˆX)T is in E. Until (5.16) is included, we also assume that
(M,Ms)≠(0,0). | (5.9) |
One has
˙H(X,ˆX)=˙V(X)+α(θM−Ms)(θM+Ms)2[ϕY(θM+Ms)2α(M+Ms)+((1−ν)νEθE−θδMM)(θM+3Ms)−ˆu(ˆX,y)(3θM+Ms)+δsMs(3θM+Ms)]+λ˙V(e)2√V(e). |
Replacing the term ˆu(ˆX,y) by ˆu(ˆX,y)−ˆu(X,y)+ˆu(X,y), we get
˙H(X,ˆX)=˙W(X)+α(θM−Ms)(3θM+Ms)(θM+Ms)2(u(ˆX,y)−u(X,y))+λ˙V(e)2√V(e). | (5.10) |
Lemma 5.1. There exists C>0 such that, for all (X,ˆX)∈E and for all y∈R2+,
‖ˆu(ˆX,y)−ˆu(X,y)‖≤C‖ˆX−X‖. | (5.11) |
Note that ˙V(e)≤−ξ‖e‖2P. Thanks to this lemma, there exists C′>0 independent of y such that
˙H(X,ˆX)≤˙W(X)+C′‖e‖−ξλ‖e‖2P2√V(e). | (5.12) |
Note that there exists a constant β>0 such that ‖e‖≤β‖e‖P. So,
˙H(X,ˆX)≤˙W(X)−(λξ2−βC′)‖e‖P. | (5.13) |
Hence, for λ=4C′β/ξ, and using the relation (3.37) and (3.44),
˙H(X,ˆX)≤−min{c1,c2}W(X)−λξ4‖e‖P. | (5.14) |
We conclude that there exists a constant
cs:=min{c1,c2,ξ4} | (5.15) |
such that
˙H(X,ˆX)<−csH(X,ˆX),ifM+Ms≠0. | (5.16) |
Let us now deal with the case where (3.30) is not satisfied. As we explained previously in the proof of the Theorem 3.1, it is sufficient to study only the case ts∈(0,+∞). Let t↦(E(t),M(t),Y(t),F(t),U(t),Ms(t),ˆE(t),ˆM(t),ˆY(t),ˆF(t),ˆU(t),ˆMs(t))T be a solution (in the Filippov sense) of the closed-loop system (5.1) such that, for some ts∈(0,+∞),
M(t)+Ms(t)=0∀t∈[0,ts] | (5.17) |
Note that (5.17) implies that
M(t)=Ms(t)=0,∀t∈[0,ts] | (5.18) |
From (3.45), (3.47), and the definition of a Filippov solution, one has on (0,ts),
(˙E˙M˙Y˙F˙U˙Ms)=(βEF(1−EK)−(νE+δE)E(1−ν)νEE−δMMννEE−κ(t)ΔηY−(η2+δY)Yη1Yκ(t)−δFFη2(1−κ(t))Y−δUUmax(0,ˆYg1+ˆEg2)−δsMs) | (5.19) |
(˙ˆE˙ˆM˙ˆY˙ˆF˙ˆU˙ˆMs)=(βEˆF−(νE+δE)ˆE(1−ν)νEˆE−δMˆMννEˆE−κ(t)ΔηˆY−(η2+δY)ˆYη1ˆYκ(t)−δFˆFη2(1−κ(t))ˆY−δUˆUmax(0,ˆYg1+ˆEg2)−δsˆMs)−LCˆX, | (5.20) |
with
κ(t)∈[0,1],g1(t)∈ϕα[0,3θ+1] and g2(t)∈(1−ν)νEθ[0,4]. | (5.21) |
From (5.18) and the second line of (5.19), one has
E(t)=0,∀t∈[0,ts] | (5.22) |
From the first line of (5.19) and (5.22), we get
F(t)=0,∀t∈[0,ts]. | (5.23) |
In the case where Y(0)=0, from the third line of (5.19) and (5.22), one has
Y(t)=0,∀t∈[0,ts]. | (5.24) |
To summarize, from (5.18), the fifth line of (5.19), (5.22), (5.23), and (5.24),
E(t)=M(t)=Y(t)=F(t)=Ms(t)=0 and ˙U(t)=−δUU(t),∀t∈[0,ts], | (5.25) |
which, with (3.13), (3.16), and (3.25) gives
˙W(t)=−σδUU(t)≤−δUW(t),∀t∈[0,ts]. | (5.26) |
In the case where Y(0)>0, from the third line of (5.19),
Y(t)>0,∀t∈[0,ts], | (5.27) |
which, together with the fourth line of (5.19) and (5.23) implies
κ(t)=0,∀t∈[0,ts]. | (5.28) |
Referring to this case already studied in the proof of Theorem 3.1, we get
˙W(t)≤−c′W(t),∀t∈[0,ts]. | (5.29) |
κ(t)∈[0,1],g1(t)∈ϕα[0,3θ+1] and g2(t)∈(1−ν)νEθ[0,4], | (5.30) |
˙Ms(t)=max(0,ˆYg1+ˆEg2)−δsMs | (5.31) |
Since Ms(t)=0∀t∈[0,ts], max(0,ˆYg1+ˆEg2)=0. For all κ(t)∈[0,1], in these two cases, the dynamics of the observation error remains
˙e=(A(κ(t))−LC)e, | (5.32) |
and one has
˙H(X,ˆX)=−c′W(X)−λξ2‖e‖P. | (5.33) |
We conclude that there exists a constant
cw:=min{c′,ξ2}, | (5.34) |
such that
˙H(X,ˆX)≤−cwH(X,ˆX). | (5.35) |
This proves Theorem 5.1 and gives the global exponential stability with the exponential decay rate ce given by relation (5.4).
We apply the backstepping control u function of the measured states y and the estimated states ˆE and ˆY given by the relation (5.2) with the following initial condition x0=(20000,5000,1500,12000,500) and ˆx0=(2000,500,150,1200,0).
The response of system (5.1) to the backstepping control (5.2) is illustrated in the Figure 4. The asymptotic convergence in large t of the different estimated ˆE,ˆF,ˆY, and ˆU (dashed) to their corresponding state variables E,F,Y, and U, respectively, illustrates the exponential convergence of the estimation error stated in Theorem 4.1. The convergence of the states and their estimates, toward zero proves the efficiency of the control (5.2) as shown in the Theorem 5.1. Figure 5 shows that the applied control function decreases when the density of the target population decreases.
In this work, we have built a feedback control law to stabilize the SIT model presented in [14,15] at extinction. Control by state feedback is a type of control rarely proposed in the literature for the overall stabilization of the SIT model. The feedback control (3.21) developed in this work has many advantages, including robustness to changing parameters. We have shown in Remark 3.4 that despite the margin of error that can be made in the estimation of the parameters, this feedback control still makes the system converge to extinction. Moreover, it does not depend on environmental capacity and this control law ensures exponential stability with the same convergence rate for the SIT system even in the high environmental capacity limit (see Theorem 3.2). Remark 3.2 shows that when the density of the target population decreases, the control also decreases
In Section 4 of our work, we built an observer for the SIT model where, using the measurement of male mosquitoes, our state estimator gives us an estimate of the other states of the system. This aspect is rarely studied for this type of dynamics. An accurate estimate of the mosquito population enables resources to be allocated more efficiently. If the intervention is effective in some areas but not in others, resources can be reallocated to maximize impact. On the other hand, the data collected during the SIT intervention provides essential information on the impact of the control in the conditions of the intervention area. This will enable informed decisions on future control strategies to be adopted according to conditions in the intervention zone by adding complementary methods or by adapting existing approaches.
One of the applications we made was to show in Section 5 that by using the data estimated via our observer to adjust the feedback control, we globally stabilize the system upon extinction. Figure 4 shows that the difficulty of estimating eggs and young females during an intervention can be compensated by the application of the observer system. Data collected on the mosquito population is also used in epidemic prevention programs. They help to adapt public health programs for better control of mosquito-borne diseases.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The author wishes to thank Luis Almeida and Jean-Michel Coron for having drawn his attention to this problem and for the many enlightening discussions during this work.
The authors declare there is no conflict of interest.
[1] | H. Larijani and G. Rad, A new spatial domain algorithm for gray scale images watermarking, International Conference on Computer and Communication Engineering, (2008), 157–161. |
[2] | H. Nemade and V. Kelkar, Reversible watermarking for colored medical images using histogram shifting method, 3rd International Conference on Computing for Sustainable Global Development, (2016), 2664–2668. |
[3] | A. Al-Nu'aimi and R. Qahwaji, Adaptive watermarking for digital colored images based on the energey of edges, IEEE International Conference on Signal Processing and Communications, (2007), 1371–1374. |
[4] | J. O'Ruanaidh and T. Pun, Rotation, scale, and translation invariant digital image watermarking, Sign. Process., 66 (1998), 303–317. |
[5] | F. N. Thakkar and V. K. Srivastava, A blind medical image watermarking: DWT-SVD based robust and secure approach for telemedicine applications. Multimed. Tools Appl., 76 (2017), 3669–3697. |
[6] | H. Hu and H. Lyu, Collective blind image watermarking in DWT-DCT domain with adaptive embedding strength governed by quality metrics, Multimed. Tools Appl., 76 (2017), 6575–6594. |
[7] | S. Horng, A blind image copyright protection scheme for e-government. J. Vis. Communic. Image, 24 (2013), 1099–1105. |
[8] | S. A. Parah, J. A. Sheikh, N. A. Loan, et al., Robust and blind watermarking technique in DCT domain using inter-block coefficient differencing, Dig. Signal Process., 53 (2016), 11–24. |
[9] | X. Kang, J. Huang and W. Zeng, Efficient general print-scanning resilient data hiding based on uniform log-polar mapping, IEEE T. Inf. Foren. Sec., 5 (2010), 1–12. |
[10] | D. Zheng, J. Zhao and A. El. Saddik, RST invariant digital image watermarking based on log-polar mapping and phase correlation, IEEE T. Circ. Syst. Vid., 13 (2003), 753–765. |
[11] | P. Niu, X. Wang, H. Jin, et al., A feature-based robust digital image watermarking scheme using Bandelet transform, Optics Laser Tech., 43 (2011), 437–450. |
[12] | M. Amini, M. O. Ahmad and M. N. S. Swamy, A robust multibit multiplicative watermark decoder using vector-based hidden Markov model in wavelet domain, IEEE T. Circ. Syst. Vid., 28 (2018), 402–413. |
[13] | M. Amini, M. O. Ahmad and M. N. S. Swamy, Digital watermark extraction in wavelet domain using hidden Markov model, Multimed. Tools Appl., 76 (2017), 3731–3749. |
[14] | E. Nezhadarya, Z. J. Wang and R. K. Ward, Robust image watermarking based on multiscale gradient direction quantization, IEEE T. Inf. Foren. Sec., 6 (2011), 1200–1213. |
[15] | V Ananthaneni, U. R.Nelakuditi, Hybrid digital image watermarking using contourlet transform (CT), DCT and SVD. Int. J. Image Process., 11 (2017), 85–93. |
[16] | A. K. Singh, M. Dave and A. Mohan, Hybrid technique for robust and imperceptible image watermarking in DWT–DCT–SVD Domain. Natl. Acad. Sci. Lett., 37 (2014) 351–358. |
[17] | Z. Wang, A. C. Bovik, H. R. Sheik, et al., Image quality assessment: from error visibility to structural similarity, IEEE T. Image Process., 13 (2004), 600–612. |
[18] | UCID-an Uncompressed Colour Image Database–a benchmark database for image retrieval with predefined ground truth, Accessed: 2019. (available: http://imagedatabase.cs.washington.edu/groundtruth/) |
[19] | T. Lin, M. Maire, S. Belongie, et al., Microsoft COCO: Common objects in context, European Conference on Computer Vision, (2014), 740–755, Accessed: 2019. (available: http://cocodataset.org/#home) |
[20] | T. Pevny, P. Bas, and J. Fridrich, Steganalysis by subtractive pixel adjacency matrix, IEEE T. Inf. Foren. Sec., 5 (2010), 215–224. |
[21] | J. Kodovsky, J. Fridrich and V. Holub, Ensemble classifiers for steganalysis of digital media, IEEE T. Inf. Foren. Sec., 7 (2012), 432–444. |
1. | Kala Agbo bidi, Luís Almeida, Jean-Michel Coron, Global Stabilization of a Sterile Insect Technique Model by feedback Laws, 2025, 204, 0022-3239, 10.1007/s10957-024-02566-4 | |
2. | Kala Agbo bidi, Jean-Michel Coron, Amaury Hayat, Nathan Lichtlé, A novel approach to feedback control with deep reinforcement learning, 2025, 202, 01676911, 106102, 10.1016/j.sysconle.2025.106102 |
Parameters | Description | Value |
βE | Effective fecundity | 10 |
νE | Hatching parameter | 0.05 |
δE | Mosquitoes in aquatic phase death rate | 0.03 |
δF | Fertilized female death rate | 0.04 |
δY | Young female death rate | 0.04 |
δM | Male death rate | 0.1 |
δs | Sterilized male death rate | 0.12 |
ν | Probability of emergence | 0.49 |