
In response to the coronavirus disease 2019 (COVID-19) pandemic, Japan conducted mass vaccination. Seventy-two million doses of vaccine (i.e., for 36 million people if a double dose is planned per person) were obtained, with initial vaccination of the older population (≡ 65 years). Because of the limited number of vaccines, the government discussed shifting the plan to administering only a single dose so that younger individuals (<65 years) could also be vaccinated with one shot. This study aimed to determine the optimal vaccine distribution strategy using a simple mathematical method. After accounting for age-dependent relative susceptibility after single- and double-dose vaccination (vs and vd, respectively, compared with unvaccinated), we used the age-dependent transmission model to compute the final size for various patterns of vaccine distributions. Depending on the values of vs, the cumulative risk of death would be lower if all 72 million doses were used as a double dose for older people than if a single-dose program was conducted in which half is administered to older people and the other half is administered to adults (i.e., 1,856,000 deaths in the former program and 1,833,000-2,355,000 deaths [depending on the values of vs] in the latter). Even if 90% of older people were vaccinated twice and 100% of adults were vaccinated once, the effective reproduction number would be reduced from 2.50 to1.14. Additionally, the cumulative risk of infection would range from 12.0% to 54.6% and there would be 421,000-1,588,000deaths (depending on the values of vs). If an epidemic appears only after completing vaccination, vaccination coverage using a single-dose program with widespread vaccination among adults will not outperform a double-dose strategy.
Citation: Tetsuro Kobayashi, Hiroshi Nishiura. Prioritizing COVID-19 vaccination. Part 1: Final size comparison between a single dose and double dose[J]. Mathematical Biosciences and Engineering, 2022, 19(7): 7374-7387. doi: 10.3934/mbe.2022348
[1] | Guodong Li, Ying Zhang, Yajuan Guan, Wenjie Li . Stability analysis of multi-point boundary conditions for fractional differential equation with non-instantaneous integral impulse. Mathematical Biosciences and Engineering, 2023, 20(4): 7020-7041. doi: 10.3934/mbe.2023303 |
[2] | Luyao Li, Licheng Fang, Huan Liang, Tengda Wei . Observer-based event-triggered impulsive control of delayed reaction-diffusion neural networks. Mathematical Biosciences and Engineering, 2025, 22(7): 1634-1652. doi: 10.3934/mbe.2025060 |
[3] | Liping Wu, Zhongyi Xiang . A study of integrated pest management models with instantaneous and non-instantaneous impulse effects. Mathematical Biosciences and Engineering, 2024, 21(2): 3063-3094. doi: 10.3934/mbe.2024136 |
[4] | Trayan Stamov, Gani Stamov, Ivanka Stamova, Ekaterina Gospodinova . Lyapunov approach to manifolds stability for impulsive Cohen–Grossberg-type conformable neural network models. Mathematical Biosciences and Engineering, 2023, 20(8): 15431-15455. doi: 10.3934/mbe.2023689 |
[5] | Yang Wang, Yongyang Liu, Yansheng Liu . Total controllability of non-autonomous second-order measure evolution systems with state-dependent delay and non-instantaneous impulses. Mathematical Biosciences and Engineering, 2023, 20(2): 2061-2080. doi: 10.3934/mbe.2023095 |
[6] | Jian Meng, Bin Zhang, Tengda Wei, Xinyi He, Xiaodi Li . Robust finite-time stability of nonlinear systems involving hybrid impulses with application to sliding-mode control. Mathematical Biosciences and Engineering, 2023, 20(2): 4198-4218. doi: 10.3934/mbe.2023196 |
[7] | Ganesh Kumar Thakur, Sudesh Kumar Garg, Tej Singh, M. Syed Ali, Tarun Kumar Arora . Non-fragile synchronization of BAM neural networks with randomly occurring controller gain fluctuation. Mathematical Biosciences and Engineering, 2023, 20(4): 7302-7315. doi: 10.3934/mbe.2023317 |
[8] | Haitao Zhu, Xinrui Ji, Jianquan Lu . Impulsive strategies in nonlinear dynamical systems: A brief overview. Mathematical Biosciences and Engineering, 2023, 20(2): 4274-4321. doi: 10.3934/mbe.2023200 |
[9] | Yanjie Hong, Wanbiao Ma . Sufficient and necessary conditions for global attractivity and stability of a class of discrete Hopfield-type neural networks with time delays. Mathematical Biosciences and Engineering, 2019, 16(5): 4936-4946. doi: 10.3934/mbe.2019249 |
[10] | Biwen Li, Yujie Liu . Average-delay impulsive control for synchronization of uncertain chaotic neural networks with variable delay impulses. Mathematical Biosciences and Engineering, 2025, 22(6): 1382-1398. doi: 10.3934/mbe.2025052 |
In response to the coronavirus disease 2019 (COVID-19) pandemic, Japan conducted mass vaccination. Seventy-two million doses of vaccine (i.e., for 36 million people if a double dose is planned per person) were obtained, with initial vaccination of the older population (≡ 65 years). Because of the limited number of vaccines, the government discussed shifting the plan to administering only a single dose so that younger individuals (<65 years) could also be vaccinated with one shot. This study aimed to determine the optimal vaccine distribution strategy using a simple mathematical method. After accounting for age-dependent relative susceptibility after single- and double-dose vaccination (vs and vd, respectively, compared with unvaccinated), we used the age-dependent transmission model to compute the final size for various patterns of vaccine distributions. Depending on the values of vs, the cumulative risk of death would be lower if all 72 million doses were used as a double dose for older people than if a single-dose program was conducted in which half is administered to older people and the other half is administered to adults (i.e., 1,856,000 deaths in the former program and 1,833,000-2,355,000 deaths [depending on the values of vs] in the latter). Even if 90% of older people were vaccinated twice and 100% of adults were vaccinated once, the effective reproduction number would be reduced from 2.50 to1.14. Additionally, the cumulative risk of infection would range from 12.0% to 54.6% and there would be 421,000-1,588,000deaths (depending on the values of vs). If an epidemic appears only after completing vaccination, vaccination coverage using a single-dose program with widespread vaccination among adults will not outperform a double-dose strategy.
Modeling of activity represents, manipulates and communicates real-world objects of daily life [11]. One of the basic tools in mathematical models for representing the biological processes of the human brain are neural networks. These are also commonly referred to as parallel distributed processing models or simply neural nets. The first mathematical model of a neural network is attributed to McCulloch and Pitts [17] who in 1943 showed that any logical function could be configured by a parallel distributed processing system of interconnected digital neurons. A very concise summary of the foundations of neural networks was provided in the book [4]. Various models of neural networks are used to consider specific tasks and in order to specify a particular neural network model it suffices to define its architectural, computational, and adaptive dynamics. The computational dynamics specifies the network initial state and a rule for its updates in time, providing that the network topology and configuration are fixed. At the beginning of the computational mode, the states of input neurons are assigned to the network input and the remaining neurons find themselves in the initial state. All potential network inputs and states form the input and state space of neural network, respectively. In dynamic (recurrent, feedback) systems, the current output signals depend, in general, on current and past input signals. There are two equivalent classes of dynamic systems: continuous-time and discrete-time systems. Continuous-time dynamic systems operate with signals which are functions of a continuous variable, that is continuous functions of time and they are described by means of differential equations.
The main objective of this paper is to study stability of continuous-time neural networks without assuming the smoothness, monotonicity and unboundedness of the activation function. The model introduced is a generalization of some additive neural networks. These type of models can simulate biological neural networks and artificial intelligence systems, from the view point of reality, and it should also naturally take into account evolutionary processes of the biological systems as well as disturbances of external influence.
In controlling nonlinear systems, the Lyapunov second method provides a way to analyze the stability of the system without explicitly solving the differential equations. Stability results concerning neural networks can be found in [10,15]. One type of stability, very useful in real world problems, is the so called Lipschitz stability. Dannan and Elaydi [5] introduced the notion of Lipschitz stability for ordinary differential equations and as noted in [5] this type of stability is important only for nonlinear problems, since it coincides with uniform stability in linear systems.
Sometimes neural networks are subjected to certain impulsive state displacements at fixed moments of time. Stability problems of neural networks under the presence of instantaneous impulses was considered in [2,3,6]. Finite-time stability of systems with instantaneous impulses was studied in [16] and for delays see [13,14]. In the case when the time acting of these state displacements is not negligibly small, this phenomenon is modeled by so called non-instantaneous impulses.
In [19] the speed of the neural impulse is discussed and the reason why modeling with instantaneous impulses is not adequate is considered. As a result, we will consider non-instantaneous impulses in modeling neural impulses.
In this paper, we study stability properties of neural networks with non-instantaneous impulses, i.e. perturbations which start their actions abruptly at some points and then continue on given finite intervals. We study the general case when
− the self-regulating parameters of all units are variable in time;
− the functions of the connection between two neurons in the network are variable in time bounded functions;
− the activation functions are non-Lipschitz (usually in the literature only the case of Lipschitz activation functions is studied).
In this paper, the investigation is based on the Lyapunov method. Some Lipschitz stability sufficient criteria are provided by using Lyapunov functions (not only quadratic ones). The results are illustrated on several particular neural networks.
We will consider the model proposed by Hopfield [7] and known as Hopfield's graded response neural network in the case when the neurons are subject to a certain impulsive state displacement at fixed moments and the duration of this displacement is not negligible small (they are known as non-instantaneous impulses). The model is based on an analog circuit consisting of capacitors, resistors and amplifiers. The input-output characteristics of the amplifiers in the circuit are modeled by some functions known as activation functions. In the literature these functions are assumed to satisfy a wide variety of assumptions, mainly they are Lipschitz functions. A neuron in the circuit is modeled by an elementary sub-circuit composed of an amplifier, a reverse amplifier, a capacitor and a resistor. The neural network is itself constructed from such elementary neurons interconnected by resistors which are analogous to biological synapses. The biological counterpart of inhibitory synapses is realized by reverse amplifiers. Usually there are two sets of inputs to the neurons. The first set of inputs are the currents I1,I2,…,In which are external inputs and often in the literature are considered as constants. The second set consists of those arising from feedback connections from other op-amps including possible self-connections. Resistances are always positive and a negative resistance to a neuron can be implemented by choosing an inverted output from certain amplifiers.
In this paper, we will assume two increasing sequences of points {ti}∞i=1 and {si}∞i=1 are given such that 0<si≤ti<si+1, i=1,2,…, and limk→∞sk=∞.
Let t0∈R+=[0,∞) be a given arbitrary point. Without loss of generality we will assume that 0≤t0<s1.
Remark 1. In the case tk=sk, k=1,2,… we obtain the case of instantaneous impulses and the so called fixed points of impulses (see, for example [8,12]).
If for the natural number k: tk≠sk then the interval (sk,tk] is called the interval of non-instantaneous impulses (see [1]).
Consider the general model of Hopfield's graded response neural network with non-instantaneous impulses (NNNI)
x′i(t)=−ci(t)xi(t)+n∑j=1aij(t)fj(xj(t))+Ii(t) for t∈∪∞k=0(tk,sk+1], i=1,2,…nxi(t)=Φk,i(t,xi(sk−0)) for t∈(sk,tk], k=1,2,…xi(t0)=x0i, | (2.1) |
where n represents the number of neurons in the network, xi(t) corresponds to the membrane potential of i-th neuron at time t, x(t)=(x1(t),x2(t),…,xn(t))T∈Rn, xi(sk−0)=limt↑skxi(t), k=1,2,…,i=1,2,…,n, ci(t)>0,i=1,2,…,n, is the rate with which the unit self-regulates or resets its potential when isolated from other units and inputs of the i-th unit at time t, aij(t), i,j=1,2,…,n, correspond to the synaptic connection weights of the i-th neuron to the j-th neuron at time t, fj(xj) denotes a measure of response or activation of j-th neuron to its membrane potential, f(x)=(f1(x1),f2(x2),…,fn(xn))T and I=(I1,I2,…,In)T is an external bias vector, the intervals (sk,tk),k=1,2,…, are intervals of acting the state displacements which are impulsive activated at fixed moments sk, the functions Φk,i(t,xi), k=1,2,… are the impulsive functions giving the impulsive perturbation of the i-th neuron on the interval (sk,tk].
In applications, when Hopfield-type neural networks are used for associative memories, the equilibrium states of the network serve as stored patterns and their stability implies that the stored patterns can be retrieved in the presence of perturbations of the network which can be non-instantaneous impulsive state displacements. When Hopfield networks are used for the solution of optimization problems, the equilibrium states denote the possible optimal solutions and their stability will guarantee the convergence to the solutions beginning from an appropriate initial state.
Definition 1. A vector x∗∈Rn, x∗=(x∗1,x∗2,…,x∗n)T is an equilibrium of NNNI (2.1), iff the equalities
0=−ci(t)x∗i+n∑j=1aij(t)fj(x∗j)+Ii(t) for t∈∪∞k=0(tk,sk+1], i=1,2,…,n | (2.2) |
and
x∗i=Φk,i(t,x∗i) for t∈(sk,tk],k=1,2,…, i=1,2,…,n | (2.3) |
hold.
Remark 2. The problem (2.1) is studied in [6] in the special case of constant self-regulating parameters ci(t)≡ci, constant connection of the i-th neuron to the j-th neuron, i.e. aij(t)≡aij, i,j=1,2,…,n, constant external bias vectors I(t)≡I, instantaneous impulses tk=sk and Lipschitz activation functions. Several types of sufficient conditions for existence of a constant vector satisfying equality (2.2) for t≥t0 are given in [6].
We assume the following:
Assumption A1. The NNNI (2.1) has an equilibrium x∗∈Rn.
If assumption A1 is satisfied, then we can shift the equilibrium point x∗ of system (2.1) to the origin. The transformation y(t)=x(t)−x∗ is used to put system (2.1) in the following form:
y′i(t)=−ci(t)yi(t)+n∑j=1aij(t)Fj(yj(t)) for t∈∪∞k=0(tk,sk+1], i=1,2,…nyi(t)=ϕk,i(t,yi(sk−0)) for t∈(sk,tk], k=1,2,…yi(t0)=y0i, | (2.4) |
where Fj(u)=fj(u+x∗j)−fj(x∗j),j=1,2,…,n and ϕk,i(t,u)=Φk,i(t,u+x∗i)−Φk,i(t,x∗i), i=1,2,…,n, k=1,2,…, y0i=x0i−x∗i with u∈R.
Note that Fi(0)≡0 and Φj,i(t,0)=0 for j=1,2,…, i=1,2,…,n, and t∈[sj,tj] and if the point x∗∈Rn is an equilibrium of NNNI (2.1), then the point y∗=0 is an equilibrium of NNNI (2.4). This allows us to study the stability properties of the zero equilibrium of NNNI (2.4). Also, without loss of generality we can assume t0=0.
We will discuss the equilibrium of some NNNI with various activation functions. It will be useful for stability analysis.
EXAMPLE 1. Let n=1, Φk∈C([sk,tk]×R,R), k=1,2,…, f∈C(R,R) and consider the scalar NNNI
x′(t)=−c(t)x(t)+a(t)f(x(t))+I(t) for t∈∪∞k=0(tk,sk+1],x(t)=Φk(t,x(sk−0)) for t∈(sk,tk], k=1,2,…x(0)=x0, | (2.5) |
where c∈C(∪∞k=0(tk,sk+1],(0,∞)), I∈C(∪∞k=0(tk,sk+1],R).
Case 1.. Let I(t)=0.5πc(t), Φk(t,u)=tku+(1−tk)π2,t∈[sk,tk], and the activation function be the cosine function f(u)=cos(u) (see [18]).
The point x∗=0.5π is an equilibrium point of NNNI(2.5) because
−c(t)0.5π+a(t)cos(0.5π)+0.5πc(t)=0, t∈∪∞k=0(tk,sk+1] |
and 0.5π=Φk(t,0.5π),t∈[sk,tk],k=1,2,….
The transformation y(t)=x(t)−0.5π will shift the model (2.5) to the following model with zero equilibrium:
y′(t)=−c(t)y(t)+a(t)cos(y(t)+π2) for t∈∪∞k=0(tk,sk+1],y(t)=tky(sk−0) for t∈(sk,tk], k=1,2,…y(0)=y0, | (2.6) |
with y0=x0−π2.
Case 2. Let I(t)=0.5c(t), Φk(t,u)=tku+0.5−t2k,t∈[sk,tk], and the activation function be the Probit function f(u)=√2 erf−1(2u−1) with erf(u)=2√π∫u0e−t2dt is the error function or the activation function be the Logit function f(u)=log(u1−u).
The point x∗=0.5 is an equilibrium point of NNNI(2.5) because −c(t)0.5+a(t)f(0.5)+0.5c(t)=0 for all t∈∪∞k=0[tk,sk+1] and Φk(t,0.5)=0.5 for t∈[sk,tk].
The transformation y(t)=x(t)−0.5 will shift the model (2.5) to the following model with zero equilibrium:
y′(t)=−c(t)y(t)+a(t)√2 erf−1(2y(t)) for t∈∪∞k=0(tk,sk+1],y(t)=tky(sk−0) for t∈(sk,tk], k=1,2,…y(0)=y0, | (2.7) |
with y0=x0−0.5.
Case 3. If I(t)=0, Φk(t,u)=ξk(t)u for t∈[sk,tk],u∈R with ξk∈C([sk,tk],R) and any activation function such that f(0)=0 then the point x∗=0 is an equilibrium point of NNNI(2.5).
First we give some results for Lipschitz stability of the general type of nonlinear non-instantaneous impulsive differential equations (NIDE).
Consider the general type of initial value problem (IVP) for NIDE:
x′(t)=H(t,x(t)) for t∈∪∞k=0(tk,sk+1],x(t)=Ψk(t,x(sk−0)) for t∈(sk,tk], k=1,2,…x(t0)=x0, | (3.1) |
where x,x0∈Rn,H:∪∞k=0[tk,sk+1]×Rn→Rn,Ψk:[sk,tk]×Rn→Rn,(k=1,2,3,…).
Remark 3. The model MNNI (2.1) is a special case of NIDE (3.1).
A detailed description of the solutions of IVP for NIDE (3.1) and some qualitative properties of the solutions of (3.1) are given in the book [1].
Let J⊂R+ be a given interval. We introduce the following classes of functions
NPC(J)={u:J→Rn:u∈C(J/∪∞k=1{sk},Rn):aaaaaaaaau(sk)=u(sk−0)=limt↑sku(t)<∞, u(sk+0)=limt↓sku(t)<∞, k:sk∈J},NPC1(J)={u:J→Rn:u∈NPC(J), u∈C1(J⋂([0,s1]∪∞k=1[tk,sk+1],Rn):aaaaaaaaau′(sk)=u′(sk−0)=limt↑sku′(t)<∞, k:sk∈J}. |
Remark 4. According to the above description any solution of (3.1) might have a discontinuity at any point sk,k=1,2,….
Let J⊂R+, 0∈J, ρ>0. We introduce the following sets:
M(J)={a∈C[J,R+]:a(0)=0, a(r) is strictly increasing in J, and a−1(αr)≤rqa(α) for some function qa: qa(α)≥1, if α≥1},K(J)={a∈C[J,R+]:a(0)=0, a(r) is strictly increasing in J, and a(r)≤Kar for some constant Ka>0},Sρ={x∈Rn: ||x||≤ρ}. |
Remark 5. The function a(u)=K1u, K1>0 is from the class K(R+) with Ka=K1.
The function a(u)=K2√u, K2≥1 is from the class M([0,1]) with qa(α)=α2.
The function a(u)=K3u, K3≥1 is from the class M(R+) with qa(α)=K3α.
We will use the class Λ of Lyapunov-like functions, defined and used for impulsive differential equations in [12].
Definition 2. Let J⊂R+ be a given interval, and Δ⊂Rn be a given set. We will say that the function V(t,x):J×Δ→R+, belongs to the class Λ(J,Δ) if
− The function V(t,x) is a continuous on J/{sk∈J}×Δ and it is locally Lipschitz with respect to its second argument;
− For each sk∈J and x∈Δ there exist finite limits
V(sk,x)=V(sk−0,x)=limt↑skV(t,x) and V(sk+0,x)=limt↓skV(t,x). |
Lipschitz stability of systems of differential equations with non-instantaneous impulses is defined and studied in [9]:
Definition 3. (Lipschitz Stability) [9] The equilibrium x∗ of NNNI (2.1) is said to be
− uniformly Lipschitz stable if there exists M≥1 and δ>0 such that for any t0≥0 and x0∈Rn the inequality |x0−x∗|<δ implies |x(t;t0,x0)−x∗|≤M|x0−x∗| for t≥t0;
− globally uniformly Lipschitz stable if there exists M≥1 such that for any t0≥0 and x0∈Rn the inequality |x0−x∗|<∞ implies |x(t;t0,x0)−x∗|≤M|x0−x∗| for t≥t0.
Remark 6. In Definitions 1.1-1.4 [5] for Lipschitz stability the constants M>0. Note we use the constant M≥1 in Definition 3 since our inequality has to hold at t=t0.
Remark 7. Note the Lipschitz stability property of the equilibrium of NNNI (2.1) implies the corresponding stability property of the zero solution of (2.4) and vice versa.
We will use the following scalar comparison differential equation with non-instantaneous impulses:
u′=g(t,u) for t∈∪∞k=0(tk,sk+1],u(t)=ψk(t,u(sk−0)) for t∈(sk,tk],k=0,1,2,…,u(t0)=u0 | (3.2) |
where u,u0∈R,g:∪∞k=0[tk,sk]×R→R,ψk:[sk,tk+1]×R→R (k=0,1,2,3,…).
Note the Lipschitz stability of the scalar non-instantaneous impulsive equation (3.2) is similarly defined as in Definition 3.
We introduce the following conditions:
(H1). The function H∈C([0,s1]⋃∞k=1[tk,sk+1]×Rn,Rn), H(t,0)=0, and for any k=1,2,… the functions Ψk∈C([sk,tk]×Rn,Rn), Ψk(t,0)=0.
(H2). The function g∈C([0,s1]⋃∞k=1[tk,sk+1]×R+,R), g(t,0)=0, and for any k=1,2,… the functions ψk∈C([sk,tk]×R+,R+) are nondecreasing with respect to their second argument and ψk(t,0)=0.
EXAMPLE 2. Consider the scalar ordinary differential equation
u′=cos(t)u(t) for t≥t0, u(t0)=u0, | (3.3) |
with a solution u(t)=u0esin(t)−sin(t0) for t≥t0. Therefore, |u(t)|≤|u0|e2 for t≥t0, i.e. the solution is globally uniformly Lipschitz stable. (see Figure 1 for the graph of the solution with initial value u0=−0.2 and the graph of the solution with initial value u0=2 and M=1.1).
Now consider the scalar differential equation with non-instantaneous impulses:
u′=cos(t)u(t) for t∈∪∞k=0(tk,sk+1],u(t)=ψk(t,u(sk−0)) for t∈(sk,tk],k=1,2,…,u(t0)=u0, | (3.4) |
where tk=2k+1,k=0,1,2,…, sk=2k,k=1,2,….
Case 1. Let ψk(t,u)=tu. The solution is a nondecreasing function approaching infinity and the zero solution is not globally uniformly Lipschitz stable (see Figure 2 for the graph of the solution with initial value u0=1 and the graph of the solution with initial u0=0.2).
Case 2. Let ψk(t,u)=ut. The zero solution is uniformly Lipschitz stable with M=1.2 and δ=2 (see Figure 3) but it is not globally uniformly Lipschitz stable.
Theorem 1. [9] Let the following conditions be satisfied:
1. Conditions (H1), (H2) are fulfilled.
2. There exists a function V(t,x)∈Λ(R+,Rn) and
(i) the inequalities
b(||x||)≤V(t,x)≤a(||x||),x∈Sρ,t∈R+ |
hold, where b∈K([0,ρ]), a∈M([0,ρ]), ρ>0;
(ii) the inequality D+V(t,x)≤g(t,V(t,x)),t∈[0,s1]⋃∞k=1(tk,sk),x∈Sρ holds where D+V(t,x)=limsuph→0+1h{V(t,x)−V(t−h,x−hH(t,x))};
(iii) for any k=1,2,… the inequality
V(t,Ψk(t,y))≤ψk(t,V(sk−0,y)),t∈(sk,tk+1],y∈Sρ |
holds.
3. The zero solution of (3.2) is uniformly Lipschitz stable (globally uniformly Lipschitz stable).
Then the zero solution of (3.1) is uniformly Lipschitz stable (globally uniformly Lipschitz stable).
Corollary 1. Let (H1) be satisfied and the inequalities
n∑i=1xiHi(t,x)≤0,t∈[0,s1]∞⋃k=0[tk,sk],x∈Rn:x=(x1,x2,…,xn)T | (3.5) |
and
n∑i=1(Ψk,i(t,x))2≤n∑i=1x2i,x∈Rn,t∈(sk,tk+1],k=1,2,… | (3.6) |
hold.
Then the zero solution of (3.1) is globally uniformly Lipschitz stable.
Proof. For V(t,x)=0.5∑ni=1x2i, x=(x1,x2,…,xn)T from (3.5) it follows that condition 2(ⅱ) of Theorem 1 is satisfied with g(t,u)≡0 and from (3.6) it follows that 2(ⅲ) of Theorem 1 is fulfilled with ψk(t,u)=u. Then the comparison scalar equation is the following equation:
u′=0 for t∈∪∞k=0(tk,sk+1],u(t)=u(sk−0) for t∈(sk,tk],k=1,2,…,u(t0)=u0. | (3.7) |
The solution of the IVP (3.7) is u(t)=u0 for t≥t0 and therefore the zero solution of (3.7) is globally uniformly Lipschitz stable with M=1. According to Theorem 1 the zero solution of (3.1) is globally uniformly Lipschitz stable.
We will study stability properties of several different types of NNNI (2.1).
We will introduce the following assumptions:
Assumption A2. The neuron activation functions are locally Lipschitz, i.e. there exist positive numbers Li, i=1,2,…,n and λ>0 such that |fi(u)−fi(v)|≤Li|u−v|, i=1,2,…,n for u,v∈R: |u−x∗i|≤λ,|v−x∗i|≤λ where the point x∗∈Rn is the equilibrium from assumption A1.
Assumption A3. There exist positive numbers Mi,j, i,j=1,2,…,n such that |ai,j(t)|≤Mi,j for t∈[0,s1]⋃∪∞k=1(tk,sk+1].
Assumption A4. There exists a number B>0 such that ci(t)≥B for t∈[0,s1]⋃∪∞k=1(tk,sk+1], i=1,2,…,n and the inequality
B≥0.5maxi=1,2,…,nn∑j=1Mij Lj+0.5n∑i=1maxj=1,2…,nMijLj | (4.1) |
holds.
Assumption A5. For any k=1,2,… and x∈Rn: ||x−x∗||≤λ the inequalities
n∑i=1(Φk,i(t,xi)−Φk,i(t,x∗i))2≤n∑i=1(xi−x∗i)2, t∈[sk,tk] |
hold where x∗ is the equilibrium point from assumption A1 and λ is the number defined in A2.
Remark 8. Assumption A5 and the inequality (4.1) are connected with the strength of the self-regulating negative feedback of each neuron which dominates its own contribution to the entire network. Such phenomena are observed in Lotka-Volterra systems modeling interacting population systems.
Theorem 2. Let assumptions A1 - A5 be satisfied. Then the equilibrium x∗ of NNNI (2.1) is globally uniformly Lipschitz stable.
Proof. From the definition of the functions Fj, assumption A2 and the inequality 2ab≤a2+b2, a,b∈R we get the inequalities
n∑i=1xiHi(t,x)=−n∑i=1ci(t)x2i+n∑i=1xin∑j=1aij(t)Fj(xj)≤−n∑i=1ci(t)x2i+n∑i=1n∑j=1Mij|xi| |fj(xj+x∗j)−fj(x∗j)|≤−n∑i=1ci(t)x2i+n∑i=1n∑j=1Mij|xi| Lj|xj|≤−n∑i=1ci(t)x2i+0.5n∑i=1n∑j=1Mij Ljx2i+0.5n∑i=1n∑j=1Mij Ljx2j≤(−B+0.5maxi=1,2,…,nn∑j=1Mij Lj+0.5n∑i=1maxj=1,2…,nMijLj)n∑j=1x2j. | (4.2) |
From assumption A5 it follows that the inequality (3.5) holds.
From assumption A5 it follows that the inequality (3.6) holds for Ψk,i(t,yi)=Φk,i(t,yi+x∗i)−Φk,i(t,x∗i).
According to Corollary 1 the zero solution of (2.1) globally uniformly Lipschitz stable. This proves the claim.
EXAMPLE 3. Let n=3, t0=0, sk=(k−0.5k)π, tk=kπ,k=1,2,…. Consider the NNNI
x′i(t)=−|csc(t)|xi(t)+3∑j=1aij(t)cos(xj(t))+0.5π|csc(t)| for t∈∪∞k=0(tk,sk+1], i=1,2,3xi(t)=Φk,i(t,xi(sk−0)) for t∈(sk,tk], k=1,2,…xi(t0)=x0i | (4.3) |
with Φk,i(t,u)=sinkt u+(1−sinkt)0.5π and A=aij(t) is given by
A(t)=(−0.1sint0.40.3−t25t2+10.3t5t+1t10t+1−0.2cost−0.1sint). |
The point x∗=(0.5π,0.5π,0.5π) is the equilibrium point of NNNI (4.3).
Note assumptions A2, A4 are satisfied with L=1 and B=0.9. Next assumption A5 is fulfilled because
3∑i=1(Φk,i(t,xi)−Φk,i(t,x∗i))2=3∑i=1(sinkt xi+(1−sinkt)0.5π−0.5π)2=3∑i=1(sinkt(xi−0.5π))2≤3∑i=1(xi−0.5π)2, t∈[sk,tk] |
Also, |aij|≤Mij, i,j=1,2,3, t≥0 where M={Mij}, is given by
M=(0.10.40.30.20.30.20.10.20.1) |
and ∑3i=1maxjMij+maxi∑3j=1Mij=(0.4+0.3+0.2)+max(0.4,0.9,0.6)=0.9+0.9=1.8, i.e. inequality (4.1) holds. Therefore, according to Theorem 2 the equilibrium of NNNI(4.3) is globally uniformly Lipschitz stable.
Consider the system (4.3) without any impulses. The right side part is not defined for all t≥t0 because of the function csc(t). The appropriate perturbing of the model by non-instantaneous impulses can induce existence of the solution and stability properties of the equilibrium.
There are many types of activation functions which are not Lipschitz. In this case, we assume:
Assumption A6. There exists a function ξ∈C([0,s1]⋃∪∞k=1(tk,sk+1],R) such that for any point x: x−x∗∈S(λ) and t∈[0,s1]⋃∪∞k=1(tk,sk+1] the inequality
(xi−x∗i)(n∑j=1ai,j(t)(fj(xj)−fj(x∗j)))≤ξ(t)n∑j=1(xj−x∗j)2, i=1,2,…,n | (4.4) |
holds where x∗ is the equilibrium from condition A1.
Theorem 3. Let the following conditions be satisfied:
1. Assumptions A1 and A6 are satisfied.
2. The functions ci∈C(∪∞k=0(tk,sk+1],R+), i=1,,…,n.
3. The impulsive functions Φk,i∈C([sk,tk]×R,R) and there exists a positive number λ such that for k=1,2,… and x∈Rn:x−x∗∈S(λ) the inequality
n∑i=1(Φk,i(t,xi)−Φk,i(t,x∗i))2≤μk(t)n∑i=1(xi−x∗i)2,t∈[sk,tk] |
holds where x∗ is the equilibrium from condition 1.
4. There exists a function p∈C1([0,∞),R+):0<α≤p(t)≤β,t≥0 such that p(t)≤p(sk) for t∈[sk,tk] and the zero solution of the scalar non-instantaneous impulsive equation (3.2) with g(t,u)=2γ(t)u and ψk(t,u)=μk(t)u is uniformly Lipschitz stable (globally uniformly Lipschitz stable) where
γ(t)=−mini=1,2,…,nci(t)+nξ(t)+0.5p′(t)p(t),t∈∪∞k=0(tk,sk+1]. |
Then the equilibrium x∗ of NNNI (2.1) is uniformly Lipschitz stable (globally uniformly Lipschitz stable).
Proof. Consider the Lyapunov function V(t,x)=p(t)∑ni=1x2i, x=(x1,x2,…,xn). Then condition 2(ⅰ) of Theorem 1 is satisfied with b(u)=α√u and a(u)=K√u, K=max{1,β}. From Remark 5 it follows that b∈K([0,λ]) and a∈M([0,λ]). From assumption A6 it follows that the function F in NNNI (2.4) satisfies the inequality
ui(n∑j=1ai,j(t)Fj(uj))≤ξ(t)n∑i=1u2i(t), t∈∪∞k=0(tk,sk+1], u∈S(λ). | (4.5) |
From inequality (4.5) with ui=yi−x∗i∈S(λ) we obtain for the Dini derivative among the system (2.4):
D+V(t,y)=p′(t)n∑i=1y2i+2p(t)n∑i=1yi(y′i)=p′(t)n∑i=1y2i+2p(t)n∑i=1(−ci(t)y2i+yin∑j=1ai,j(t)Fj(yj))≤p′(t)n∑i=1y2i−2p(t)n∑i=1(mini=1,2,…,nci(t))y2i+2p(t)n∑i=1yin∑j=1ai,j(t)Fj(yj)≤p′(t)n∑i=1y2i−2p(t)n∑i=1(mini=1,2,…,nci(t))y2i+2p(t)n∑i=1ξ(t)n∑j=1y2j=2p(t)(−mini=1,2,…,nci(t)+nξ(t)−0.5p′(t)p(t))n∑i=1y2i(t)≤2γ(t)p(t)n∑i=1y2i(t)=2γ(t)V(t,y), t∈∪∞k=0(tk,sk+1]. | (4.6) |
Therefore, Condition 2(ⅱ) of Theorem 1 is satisfied with g(t,u)=2γ(t)u.
Then for t∈(sk,tk], y∈Sλ the Lyapunov function satisfies
V(t,ϕk(t,yi))=p(t)n∑i=1ϕ2k,i(t,y)≤μk(t)p(t)n∑i=1y2i≤μk(t)p(si)n∑i=1y2i=μk(t)V(sk−0,y). |
Therefore, Condition 2(ⅲ) of Theorem 1 is satisfied with ψk(t,u)=μk(t)u.
According to Theorem 1 the zero solution of NNNI (2.4) is uniformly Lipschitz stable (globally uniformly Lipschitz stable)
This proves the claim of Theorem 3.
Remark 9. The conditions in Theorem 2 and Theorem 3 differ. The inequality for the impulsive functions in condition 3 of Theorem 3 is weaker than Assumption A5. At the same time if Assumption A2 is satisfied then Assumption A6 follows. However the inequality in condition 4 of Theorem 3 is more restrictive than the inequality in Assumption A3. Both Theorems are applicable in different cases.
We will illustrate the difference between the applications of Theorem 2 and Theorem 3 on an example.
EXAMPLE 4. To simplify the calculations we consider the simplest neural network of two units, i.e. n=2. Let t0∈[0,1), sk=k, tk=k+0.5, k=1,2,… and consider the nonlinear NNNI
x′1(t)=−e−tx1(t)−0.5t+1f1(x1(t))+1t+1f2(x2(t))+0.5e−t−1t+1x′2(t)=−1t+1x2(t)−1t+1f2(x2(t))+1t+1 for t∈∪∞k=0(tk,sk+1],x1(t)=1tx1(sk−0)+0.5t−1tx2(t)=1tx2(sk−0) for t∈(sk,tk], k=1,2,…. | (4.7) |
where the activation functions are the Probit function f1(u)=√2 erf−1(2u−1) and the cosine function f2(u)=cos(u), c1(t)=e−t, c2(t)=1t+1, I1(t)=0.5e−t−1t+1, I2(t)=1t+1, Φ1,k(t,u)=1tu−0.5t−1t, Φ2,k(t,u)=ut and
A(t)=(−0.5t+11t+10−1t+1 ). |
Then the equilibrium is x∗1=0.5,x∗2=0.
Theorem 3 is not applicable because assumptions A4 and A2 for c1(t),c2(t) and f1(u) are not satisfied.
Now using (u−0.5)f1(u)≥(u−0.5)2, u∈R (see Figure 4) it follows that
−0.5t+1(u−0.5)f1(u)≤−0.5t+1(u−0.5)2≤0.5t+1(u−0.5)2, u∈R. | (4.8) |
According to inequalities (4.8) and |cos(u)−1|≤|u|, u∈R (see Figure 5) we obtain
(x1−0.5)(a1,1(t)(f1(x1)−f1(0.5))+a1,2(t)(f2(x2)−f2(0)))=(x1−0.5)(a1,1(t)f1(x1)+a1,2(t)(f2(x2)−1))=−0.5t+1(x1−0.5)f1(x1)+1t+1(x1−0.5)(f2(x2)−1))≤0.5t+1(x1−0.5)2+1t+1(x1−0.5)(cos(x2)−1)≤0.5t+1(x1−0.5)2+1t+1 |(x1−0.5)||x2|≤1t+1(x1−0.5)2+0.5t+1x22≤1t+1((x1−0.5)2+x22). | (4.9) |
From inequality u(cos(u)−1)≤u2 (see Figure 6) it follows that
x2(a2,1(t)(f1(x1)−f1(0.5))+a2,2(t)(f2(x2)−f2(0)))=−1t+1x2(cos(x2)−1)≤0.5t+1x22≤1t+1((x1−0.5)2+x22). | (4.10) |
Therefore, assumption A6 is satisfied with ξ(t)=1t+1.
Consider the function p(t)=1.1+cos(t), 0.1<p(t)≤1.1 for t≥0. The inequality p(t)≤p(sk) holds for t∈[sk,tk].
Also, applying the inequality e−t≤11+t (see Figure 7) we obtain the inequality
γ(t)=−mini=1,2ci(t)+2ξ(t)+0.5p′(t)p(t)=−e−t+211+t+0.5sin(t)1.1+cos(t)≤1+0.5sin(t)1.1+cos(t)≤1.05, t∈(tk,sk], k=1,2,…. | (4.11) |
Condition 3 of Theorem 3 is satisfied with μk(t)=1t2 because ∑2i=1(Φk,i(t,xi)−Φk,i(t,x∗i))2=1t2((x1−0.5)2+x22), t∈[sk,tk].
Then the scalar non-instantaneous impulsive equation (3.2) will be reduced to
u′=2.1u for t∈[t0,1]∪∞k=1(k+0.5,k+1],u(t)=1t2u(sk−0)) for t∈(k,k+0.5],k=1,2,…,u(t0)=u0. | (4.12) |
The zero solution of (4.12) is uniformly Lipschitz stable (see Figure 8) with M=11 and δ=1. According to Theorem 3 the equilibrium of (4.7) is uniformly Lipschitz stable.
Several sufficient conditions are given for stability of the equilibrium point of neural networks with non-instantaneous impulses. We study the general case when
− the self-regulating parameters of all units are variable in time;
− the functions of the connection between two neurons in the network are variable in time bounded functions;
− the activation functions are non-Lipschitz (usually in the literature only the case of Lipschitz activation functions is studied).
The investigation is based on applying the Lyapunov approach and various types of Lyapunov functions are used. These piecewise continuous Lyapunov functions are used as an apparatus to study stability properties of neural networks with non-instantaneous for different types of activation functions. Many examples with different types of activation functions such as cos function, Logit function, continuous Tan-Sigmoid function, are provided to illustrate the practical application of our results.
All authors declare no conflicts of interest in this paper.
[1] |
K. L. Prem, Y. Liu, T. W. Russell, A. J. Kucharski, R. M. Eggo, N. Davies, et.al., The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: A modelling study, Lancet Public Health, 5 (2020), e261–e270. https://doi.org/10.1016/S2468-2667(20)30073-6 doi: 10.1016/S2468-2667(20)30073-6
![]() |
[2] |
B. L. Dickens, J. R. Koo, J. T. Lim, M. Park, S. Quaye, H. Sun, et. al., Modelling lockdown and exit strategies for COVID-19 in Singapore, Lancet Regional Health – Western Pacific, 1 (2020) 100004. https://doi.org/10.1016/j.lanwpc.2020.100004 doi: 10.1016/j.lanwpc.2020.100004
![]() |
[3] |
E. Mahase, Covid-19: Novavax vaccine efficacy is 86% against UK variant and 60% against South African variant, BMJ, 372 (2021), n296. https://doi.org/10.1136/bmj.n296 doi: 10.1136/bmj.n296
![]() |
[4] |
G. Persad, M. E. Peek, E. J. Emanuel, Fairly prioritizing groups for access to COVID-19 vaccines, JAMA, 324 (2020), 1601–1602. https://doi.org/10.1001/jama.2020.18513 doi: 10.1001/jama.2020.18513
![]() |
[5] |
K. Liu, Y. Lou, Optimizing COVID-19 vaccination programs during vaccine shortages, Infect Dis. Model, 7 (2022), 286–298. https://doi.org/10.1016/j.idm.2022.02.002 doi: 10.1016/j.idm.2022.02.002
![]() |
[6] |
P. C. Jentsch, M. Anand, C. T. Bauch, Prioritising COVID-19 vaccination in changing social and epidemiological landscapes: A mathematical modelling study, Lancet Infect. Dis., 3099 (2021), 00057-8. https://doi.org/10.1016/S1473-3099(21)00057-8 doi: 10.1016/S1473-3099(21)00057-8
![]() |
[7] |
G. Persad, E. J. Emanuel, S. Sangenito, A. Glickman, S. Phillips, E. A. Largent, Public perspectives on COVID-19 vaccine prioritization, JAMA Netw. Open, 4 (2021), e217943. https://doi.org/10.1001/jamanetworkopen.2021.7943 doi: 10.1001/jamanetworkopen.2021.7943
![]() |
[8] | E. Rumpler, M. J. Feldman, M. T. Bassett, M. Lipsitch, Equitable COVID-19 vaccine prioritization: Front-line workers or 65–74 year olds?, preprint, medRxiv, no. 2022.02.03.22270414. https://doi.org/10.1101/2022.02.03.22270414 |
[9] |
L. A. C. Chapman, P. Shukla, I. Rodríguez-Barraquer, P. B. Shete, T. M. León, K. Bibbins-Domingo, et al., Risk factor targeting for vaccine prioritization during the COVID-19 pandemic, Sci. Rep., 12 (2022), 3055. https://doi.org/10.1038/s41598-022-06971-5 doi: 10.1038/s41598-022-06971-5
![]() |
[10] |
S. Epstein, K. Ayers, B. K. Swenor, COVID-19 vaccine prioritisation for people with disabilities, Lancet Public Health, 6 (2021), e361. https://doi.org/10.1016/S2468-2667(21)00093-1 doi: 10.1016/S2468-2667(21)00093-1
![]() |
[11] |
S. Han, J. Cai, J. Yang, J. Zhang, Q. Wu, W. Zheng, et al., Time-varying optimization of COVID-19 vaccine prioritization in the context of limited vaccination capacity, Nat. Commun., 12 (2021), 4673. https://doi.org/10.1038/s41467-021-24872-5 doi: 10.1038/s41467-021-24872-5
![]() |
[12] |
E. K. Lee, Z. L. Li, Y. K. Liu, J. LeDuc, Strategies for vaccine prioritization and mass dispensing, Vaccines (Basel), 9 (2021), 506. https://doi.org/10.3390/vaccines9050506 doi: 10.3390/vaccines9050506
![]() |
[13] |
H. Tatapudi, R. Das, T. K. Das, Impact of vaccine prioritization strategies on mitigating COVID-19: An agent-based simulation study using an urban region in the United States, BMC Med. Res. Methodol., 21 (2021), 272. https://doi.org/10.1186/s12874-021-01458-9 doi: 10.1186/s12874-021-01458-9
![]() |
[14] |
R. Strodel, L. Dayton, H. M. Garrison-Desany, G. Eber, C. Beyrer, J. Arscott, et al., COVID-19 vaccine prioritization of incarcerated people relative to other vulnerable groups: An analysis of state plans, PLoS One, 16 (2021), e0253208. https://doi.org/10.1371/journal.pone.0253208 doi: 10.1371/journal.pone.0253208
![]() |
[15] |
K. M. Bubar, K. Reinholt, S. M. Kissler, M. Lipsitch, S. Cobey, Y. H. Grad, et al., Model-informed COVID-19 vaccine prioritization strategies by age and serostatus, Science, 371 (2021), 916–921. https://doi.org/10.1126/science.abe6959 doi: 10.1126/science.abe6959
![]() |
[16] |
J. H. Buckner, G. Chowell, M. R. Springborn, Dynamic prioritization of COVID-19 vaccines when social distancing is limited for essential workers, Proc. Natl. Acad. Sci. U. S. A., 118 (2021), e2025786118. https://doi.org/10.1073/pnas.2025786118 doi: 10.1073/pnas.2025786118
![]() |
[17] |
H. Nishiura, K. Iwata, A simple mathematical approach to deciding the dosage of vaccine against pandemic H1N1 influenza, Euro. Surveill., 14 (2009), 1–5. https://doi.org/10.2807/ese.14.45.19396-en doi: 10.2807/ese.14.45.19396-en
![]() |
[18] |
M. A. Billah, M. M. Miah, M. N. Khan, Reproductive number of coronavirus: A systematic review and meta-analysis based on global level evidence, PLoS One, 15 (2020), e0242128. https://doi.org/10.1371/journal.pone.0242128 doi: 10.1371/journal.pone.0242128
![]() |
[19] |
A. T. Levin, W. P. Hanage, N. Owusu-Boaitey, K. B. Cochran, S. P. Walsh, G. Meyerowitz-Katz, Assessing the age specificity of infection fatality rates for COVID-19: Systematic review, meta-analysis, and public policy implications, Eur. J. Epidemiol., 35 (2020), 1123–1138. https://doi.org/10.1007/s10654-020-00698-1 doi: 10.1007/s10654-020-00698-1
![]() |
[20] |
L. Munasinghe, Y. Asai, H. Nishiura, Quantifying heterogeneous contact patterns in Japan: A social contact survey, Theor. Biol. Med. Model., 16 (2019), 6. https://doi.org/10.1186/s12976-019-0102-8 doi: 10.1186/s12976-019-0102-8
![]() |
[21] |
E. Mahase, Covid-19: Reports from Israel suggest one dose of Pfizer vaccine could be less effective than expected, BMJ, 372 (2021), n217. https://doi.org/10.1136/bmj.n217 doi: 10.1136/bmj.n217
![]() |
[22] | H. Nishiura, Tracking Public Health and Social Measures, in World Health Organization, 2021, work in progress. |
[23] |
N. Dagan, N. Barda, E. Kepten, O. Miron, S. Perchik, M. A. Katz, et al., BNT162b2 mRNA Covid-19 vaccine in a nationwide mass vaccination setting, N. Engl. J. Med., 384 (2021), 1412–1423. https://doi.org/10.1056/NEJMoa2101765 doi: 10.1056/NEJMoa2101765
![]() |
[24] |
S. J. B. Hanley, E. Yoshioka, Y. Ito, R. Kishi, HPV vaccination crisis in Japan, Lancet, 385 (2015), 2571. https://doi.org/10.1016/S0140-6736(15)61152-7 doi: 10.1016/S0140-6736(15)61152-7
![]() |
[25] |
S. Moore, E. M. Hill, M. J. Tildesley, L. Dyson, M. J. Keeling, Vaccination and non-pharmaceutical interventions for COVID-19: A mathematical modelling study, Lancet Infect. Dis., 3099 (2021), 793–802. https://doi.org/10.1016/S1473-3099(21)00143-2 doi: 10.1016/S1473-3099(21)00143-2
![]() |
[26] |
J. Wise, Covid-19: The E484K mutation and the risks it poses, BMJ, 372 (2021), n359. https://doi.org/10.1136/bmj.n359 doi: 10.1136/bmj.n359
![]() |
[27] |
T. Burki, Understanding variants of SARS-CoV-2, Lancet, 397 (2021), 462. https://doi.org/10.1016/S0140-6736(21)00298-1 doi: 10.1016/S0140-6736(21)00298-1
![]() |
[28] |
D. A. Collier, A. De Marco, I. A. T. M. Ferreira, B. Meng, R. Datir, A. C. Walls, et. al., Sensitivity of SARS-CoV-2 B.1.1.7 to mRNA vaccine-elicited antibodies, Nature, 593 (2021), 136–141. https://doi.org/10.1038/s41586-021-03412-7 doi: 10.1038/s41586-021-03412-7
![]() |
[29] |
D. Planas, T. Bruel, L. Grzelak, F. Guivel-Benhassine, I. Staropoli, F. Porrot, et al., Sensitivity of infectious SARS-CoV-2 B.1.1.7 and B.1.351 variants to neutralizing antibodies, Nat. Med., 27 (2021), 917–924. https://doi.org/10.1038/s41591-021-01318-5 doi: 10.1038/s41591-021-01318-5
![]() |
[30] |
A. Muik, A. K. Wallisch, B. Sänger, K. A. Swanson, J. Mühl, W. Chen, et al., Neutralization of SARS-CoV-2 lineage B.1.1.7 pseudovirus by BNT162b2 vaccine–elicited human sera, Science, 371 (2021), 1152–1153. https://doi.org/10.1126/science.abg6105 doi: 10.1126/science.abg6105
![]() |
[31] |
P. Wang, M. S. Nair, L. Liu, S. Iketani, Y. Luo, Y. Guo, et al., Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7, Nature, 593 (2021), 130–135. https://doi.org/10.1038/s41586-021-03398-2 doi: 10.1038/s41586-021-03398-2
![]() |
[32] |
D. Zhou, W. Dejnirattisai, P. Supasa, C. Liu, A. J. Mentzer, H. M. Ginn, et al., Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera, Cell, 184 (2021), 2348–2361. https://doi.org/10.1016/j.cell.2021.02.037 doi: 10.1016/j.cell.2021.02.037
![]() |
[33] |
K. Leung, M. Jit, G. M. Leung, J. T. Wu, The allocation of COVID-19 vaccines and antivirals against emerging SARS-CoV-2 variants of concern in East Asia and Pacific region: A modelling study, Lancet Regional Health – Western Pacific, 21 (2022), 100389. https://doi.org/10.1016/j.lanwpc.2022.100389 doi: 10.1016/j.lanwpc.2022.100389
![]() |
[34] |
C. C. John, V. Ponnusamy, S. K. Chandrasekaran, R. Nandakumar, A survey on mathematical, machine learning and deep learning models for COVID-19 transmission and diagnosis, IEEE Rev. Biomed. Eng., 15 (2022), 325–340. https://doi.org/10.1109/RBME.2021.3069213 doi: 10.1109/RBME.2021.3069213
![]() |
[35] |
S. M. Saadat, Z. R. Tehrani, J. Logue, M. Newman, M. B. Frieman, A. D. Harris, et al., Binding and neutralization antibody titers after a single vaccine dose in health care workers previously infected with SARS-CoV-2, JAMA, 325 (2021), 1467–1469. https://doi.org/10.1001/jama.2021.3341 doi: 10.1001/jama.2021.3341
![]() |
[36] | COVID-19 Dashboard, the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University, 2021. Available From: https://gisanddata.maps.arcgis.com/apps/opsdashboard/index.html#/bda7594740fd40299423467b48e9ecf6 (accessed Mar. 30, 2021). |
[37] | COVID-19 Advisory Board, Ministry of Health, Labor and Welfare (in Japanese), From https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/0000121431_00216.html (accessed Apr. 30, 2021). |
[38] |
I. Locatelli, B. Trächsel, V. Rousson, Estimating the basic reproduction number for COVID-19 in Western Europe, PLoS One, 16 (2021), 1–9. https://doi.org/10.1371/journal.pone.0248731 doi: 10.1371/journal.pone.0248731
![]() |
[39] |
Z. Zhuang, S. Zhao, Q. Lin, P. Cao, Y. Lou, L. Yang, et al., Preliminary estimates of the reproduction number of the coronavirus disease (COVID-19) outbreak in Republic of Korea and Italy by 5 March 2020, Int. J. Infect. Dis., 95 (2020), 308–310. https://doi.org/10.1016/j.ijid.2020.04.044 doi: 10.1016/j.ijid.2020.04.044
![]() |
[40] |
M. Al-Raeei, The basic reproduction number of the new coronavirus pandemic with mortality for India, the Syrian Arab Republic, the United States, Yemen, China, France, Nigeria and Russia with different rate of cases, Clin. Epidemiol. Glob. Heal., 9 (2021), 147–149. https://doi.org/10.1016/j.cegh.2020.08.005 doi: 10.1016/j.cegh.2020.08.005
![]() |
1. | Antonia Dobreva, 2021, 2321, 0094-243X, 030006, 10.1063/5.0040098 | |
2. | Ravi Agarwal, Snezhana Hristova, Donal O’Regan, Lyapunov Functions and Lipschitz Stability for Riemann–Liouville Non-Instantaneous Impulsive Fractional Differential Equations, 2021, 13, 2073-8994, 730, 10.3390/sym13040730 | |
3. | Snezhana Hristova, Stepan Tersian, Radoslava Terzieva, Lipschitz Stability in Time for Riemann–Liouville Fractional Differential Equations, 2021, 5, 2504-3110, 37, 10.3390/fractalfract5020037 | |
4. | Mohammed Said Souid, Krasimira Ivanova, 2022, 2459, 0094-243X, 030034, 10.1063/5.0083544 | |
5. | Hao Deng, Chuandong Li, Yinuo Wang, Asymptotic stability of non‐instantaneous impulsive systems and T‐S fuzzy non‐instantaneous impulsive control for nonlinear systems, 2023, 1751-8644, 10.1049/cth2.12448 |