Conditions | E0 | E1 | E∗ |
R10<1 and R20<1 | Yes (GAS) | No | No |
R10<1 and R20>1 | Yes (Unstable) | Yes (GAS) | No |
R10>1 and R20<1 | Yes (Unstable) | No | Yes (GAS) |
R10>1 and R20>1 | Yes (Unstable) | Yes (Unstable) | Yes (GAS) |
In this study, an accurate tool is provided for the evaluation of the effect of joint motion effect on gait stability. This quantitative gait evaluation method relies exclusively on the analysis of data acquired using acceleration sensors. First, the acceleration signal of lower limb motion is collected dynamically in real-time through the acceleration sensor. Second, an algorithm based on improved dynamic time warping (DTW) is proposed and used to calculate the gait stability index of the lower limbs. Finally, the effects of different joint braces on gait stability are analyzed. The experimental results show that the joint brace at the ankle and the knee reduces the range of motions of both ankle and knee joints, and a certain impact is exerted on the gait stability. In comparison to the ankle joint brace, the knee joint brace inflicts increased disturbance on the gait stability. Compared to the joint motion of the braced side, which showed a large deviation, the joint motion of the unbraced side was more similar to that of the normal walking process. In this paper, the quantitative evaluation algorithm based on DTW makes the results more intuitive and has potential application value in the evaluation of lower limb dysfunction, clinical training and rehabilitation.
Citation: Xuecheng Weng, Chang Mei, Farong Gao, Xudong Wu, Qizhong Zhang, Guangyu Liu. A gait stability evaluation method based on wearable acceleration sensors[J]. Mathematical Biosciences and Engineering, 2023, 20(11): 20002-20024. doi: 10.3934/mbe.2023886
[1] | Yan’e Wang , Zhiguo Wang, Chengxia Lei . Asymptotic profile of endemic equilibrium to a diffusive epidemic model with saturated incidence rate. Mathematical Biosciences and Engineering, 2019, 16(5): 3885-3913. doi: 10.3934/mbe.2019192 |
[2] | Yoichi Enatsu, Yukihiko Nakata . Stability and bifurcation analysis of epidemic models with saturated incidence rates: An application to a nonmonotone incidence rate. Mathematical Biosciences and Engineering, 2014, 11(4): 785-805. doi: 10.3934/mbe.2014.11.785 |
[3] | Pan Yang, Jianwen Feng, Xinchu Fu . Cluster collective behaviors via feedback pinning control induced by epidemic spread in a patchy population with dispersal. Mathematical Biosciences and Engineering, 2020, 17(5): 4718-4746. doi: 10.3934/mbe.2020259 |
[4] | Jummy F. David, Sarafa A. Iyaniwura, Michael J. Ward, Fred Brauer . A novel approach to modelling the spatial spread of airborne diseases: an epidemic model with indirect transmission. Mathematical Biosciences and Engineering, 2020, 17(4): 3294-3328. doi: 10.3934/mbe.2020188 |
[5] | Chang-Yuan Cheng, Shyan-Shiou Chen, Xingfu Zou . On an age structured population model with density-dependent dispersals between two patches. Mathematical Biosciences and Engineering, 2019, 16(5): 4976-4998. doi: 10.3934/mbe.2019251 |
[6] | Yicang Zhou, Zhien Ma . Global stability of a class of discrete age-structured SIS models with immigration. Mathematical Biosciences and Engineering, 2009, 6(2): 409-425. doi: 10.3934/mbe.2009.6.409 |
[7] | Zongwei Ma, Hongying Shu . Viral infection dynamics in a spatial heterogeneous environment with cell-free and cell-to-cell transmissions. Mathematical Biosciences and Engineering, 2020, 17(3): 2569-2591. doi: 10.3934/mbe.2020141 |
[8] | Yujie Sheng, Jing-An Cui, Songbai Guo . The modeling and analysis of the COVID-19 pandemic with vaccination and isolation: a case study of Italy. Mathematical Biosciences and Engineering, 2023, 20(3): 5966-5992. doi: 10.3934/mbe.2023258 |
[9] | Zhen Jin, Guiquan Sun, Huaiping Zhu . Epidemic models for complex networks with demographics. Mathematical Biosciences and Engineering, 2014, 11(6): 1295-1317. doi: 10.3934/mbe.2014.11.1295 |
[10] | Qian Yan, Xianning Liu . Dynamics of an epidemic model with general incidence rate dependent on a class of disease-related contact functions. Mathematical Biosciences and Engineering, 2023, 20(12): 20795-20808. doi: 10.3934/mbe.2023920 |
In this study, an accurate tool is provided for the evaluation of the effect of joint motion effect on gait stability. This quantitative gait evaluation method relies exclusively on the analysis of data acquired using acceleration sensors. First, the acceleration signal of lower limb motion is collected dynamically in real-time through the acceleration sensor. Second, an algorithm based on improved dynamic time warping (DTW) is proposed and used to calculate the gait stability index of the lower limbs. Finally, the effects of different joint braces on gait stability are analyzed. The experimental results show that the joint brace at the ankle and the knee reduces the range of motions of both ankle and knee joints, and a certain impact is exerted on the gait stability. In comparison to the ankle joint brace, the knee joint brace inflicts increased disturbance on the gait stability. Compared to the joint motion of the braced side, which showed a large deviation, the joint motion of the unbraced side was more similar to that of the normal walking process. In this paper, the quantitative evaluation algorithm based on DTW makes the results more intuitive and has potential application value in the evaluation of lower limb dysfunction, clinical training and rehabilitation.
Population migration is a common phenomenon. With the migration of population, infectious diseases can easily spread from one area to another, so it is meaningful to consider population migration when studying the spread of infectious diseases [1,2,3,4,5,6,7].
Wang and Mulone [2] established an SIS infectious disease model with standard incidence based on two patches. It is proved that the basic reproduction number is the threshold of the uniform persistence and disappearance of the disease. The dispersal rate of the population will make the infectious disease persist or disappear in all patches. There will be no the phenomena that infectious diseases persists in one patch but disappears in the other.
Sun et al. [3] put forward an SIS epidemic model with media effect in a two patches setting. Under the assumption that the migration matrix is irreducible, it is proved that if the basic reproduction number is greater than 1 then the system persists and solutions converge to an endemic equilibrium and that if the basic reproduction number is less than 1 then solutions tend to an equilibrium without disease.
Gao et al. [5,6] studied an SIS multi-patch model with variable transmission coefficients. Their results show that the basic reproduction number R0 is a threshold parameter of the disease dynamics.
All the patch models referenced above assume that the migration matrix is irreducible. The studies in which the migration matrix is reducible are few. Therefore, based on the case of two patches, we consider that the individuals can only migrate from one patch to the other. In this case, the migration matrix is reducible. It can characterize the phenomenon that individuals migrate in one direction between two regions, such as, from the rural patch to the urban one [8] and from a small community hospital to a large teaching hospital [4].
It is well known that the incidence rate plays an important role in the modeling of infectious disease. Considering the saturation phenomenon for numerous infected individuals, Capasso and Serio [9] first introduce a nonlinear bounded function g(I) to form the interaction term g(I)S in 1978. It can characterize the behavioral changes of individuals, such as wearing masks or reducing their social activities and direct contact with others with the increase of infectious individuals. After that, the saturation incidence rate has attracted much attention and various nonlinear types of incidence rate are employed. The most commonly used types are Holling type Ⅱ λSI1+αI [10,11,12] and βSI1+αS [13], Monod-Haldane type λSI1+αI2 [14], Beddington-DeAngelis type λSI1+αS+βI [15,16,17] and Crowley-Martin type λSI(1+αS)(1+βI) [18,19].
In this paper, we consider infectious disease transmission models with saturation incidence rate. The rest of this paper is organized as follows: In Section 2, we establish a two-patch SIS model with saturating contact rate and one-directing population dispersal. We discuss the existence of disease-free equilibrium, boundary equilibrium and endemic equilibrium and prove the global asymptotic stability of the equilibriums in Section 3. In Section 4, we perform simulations to illustrate the results and analyze the effect of the contact rate and population migration on epidemic transmission. Finally, we discuss in Section 5.
In the two patches, the population is divided into two states: susceptible and infective. Thus we can establish a two-patch SIS model with saturating contact rate and one-directing population dispersal
{dS1(t)dt=A1−d1S1−β1S1I11+α1I1−mS1+γ1I1,dS2(t)dt=A2−d2S2−β2S2I21+α2I2+mS1+γ2I2,dI1(t)dt=β1S1I11+α1I1−d1I1−mI1−γ1I1,dI2(t)dt=β2S2I21+α2I2−d2I2+mI1−γ2I2, | (2.1) |
where Si is the number of susceptible population in patch i (i=1,2), Ii is the number of infective population in patch i (i=1,2), Ai is the recruitment into patch i (i=1,2), di is the natural mortality rate, γi is the recovery rate of an infective individual in patch i (i=1,2), m is the migration rate form patch 1 to patch 2. Since the individuals can migrate from the first patch to the second, patch 1 is the source patch and patch 2 is the sink patch. The initial conditions is
Si(0)>0, Ii(0)≥0, i=1,2, I1(0)+I2(0)>0. | (2.2) |
Denote the population in patch i by Ni. Then Ni=Si+Ii. From system (2.1), the differential equations governing the evolution of N1 and N2 are
{dN1(t)dt=A1−(d1+m)N1,dN2(t)dt=A2−d2N2+mN1. | (2.3) |
Obviously, system (2.3) has a unique equilibrium (N∗1,N∗2)=(A1d1+m,A2d2+mA1d2(d1+m)) which is globally asymptotically stable for (2.3). So (2.1) is equivalent the following system
{dN1(t)dt=A1−(d1+−m)N1,dN2(t)dt=A2−d2N2+mN1,dI1(t)dt=β1(N1−I1)I11+α1I1−d1I1−mI1−γ1I1,dI2(t)dt=β2(N2−I2)I21+α2I2−d2I2+mI1−γ2I2. | (2.4) |
Because limt→∞Ni(t)→N∗i (i=1,2), system (2.4) leads to the following limit system
{dI1(t)dt=β1(N∗1−I1)I11+α1I1−d1I1−mI1−γ1I1,dI2(t)dt=β2(N∗2−I2)I21+α2I2−d2I2+mI1−γ2I2. | (2.5) |
Let Ω={(I1,I2)|0≤I1≤N∗1,0≤I2≤N∗2}. Then Ω is invariant region for system (2.5).
Define the basic reproduction number in the two patches respectively by R10=β1A1(d1+γ1+m)2=β1d1+γ1+mN∗1, R20=β2(d2+γ2)N∗2. The basic reproduction number R10 gives the expected secondary infections in the source patch produced by a primary infected individual in the source patch when the population is supposed to be in the disease-free equilibrium. The basic reproduction number R20 gives the expected secondary infections in the sink patch produced by a primary infected individual in the sink patch when the population is supposed to be in the disease-free equilibrium. Then we have the following theorem.
Theorem 3.1. For the system (2.5), we have
(i) The disease-free equilibrium E0:=(0,0) always exists;
(ii) The boundary equilibrium E1:=(0,β2N∗2−d2−γ2(d2+γ2)α2+β2) exists if R20>1;
(iii) There is a unique epidemic equilibrium E∗ if R10>1.
Proof. (ⅰ) can be easily proved.
Let
β1(N∗1−I1)I11+α1I1−d1I1−mI1−γ1I1=0, | (3.1) |
β2(N∗2−I2)I21+α2I2−d2I2+mI1−γ2I2=0. | (3.2) |
From Eq (3.1), we can have I1=0 always satisfies Eq (3.1). When I1=0, from Eq (3.2), we have
I2=β2N∗2−d2−γ2(d2+γ2)α2+β2. |
If R20>1, then I2=β2N∗2−d2−γ2(d2+γ2)α2+β2>0. So The boundary equilibrium E1:=(0,β2N∗2−d2−γ2(d2+γ2)α2+β2) exists if R20>1. The conclusion (ⅱ) is proved.
If R10>1, Eq (3.1) has a positive solution I∗1=β1N∗1−(d1+m+γ1)(d1+m+γ1)α1+β1. Solve Eq (3.2), we have
I2=(β2N∗2−d2−γ2+mα2I1)±√(β2N∗2−d2−γ2+mα2I1)2+4[(d2+γ2)α2+β2]mI12[(d2+γ2)α2+β2]. | (3.3) |
Substituting I∗1 into Eq (3.3), we have
I∗2=(β2N∗2−d2−γ2+mα2I∗1)±√(β2N∗2−d2−γ2+mα2I∗1)2+4[(d2+γ2)α2+β2]mI∗12[(d2+γ2)α2+β2]. |
Since I∗2≥0 is meaning only, we take
I∗2=(β2N∗2−d2−γ2+mα2I∗1)+√(β2N∗2−d2−γ2+mα2I∗1)2+4[(d2+γ2)α2+β2]mI∗12[(d2+γ2)α2+β2]. |
So if R10>1, there is a unique epidemic equilibrium E∗=(I∗1,I∗2), where I∗1=β1N∗1−(d1+m+γ1)(d1+m+γ1)α1+β1 and I∗2=(β2N∗2−d2−γ2+mα2I∗1)+√(β2N∗2−d2−γ2+mα2I∗1)2+4[(d2+γ2)α2+β2]mI∗12[(d2+γ2)α2+β2]. The conclusion (ⅲ) is proved.
This completes the proof of the theorem.
From the above analysis, we have the following theorem.
Theorem 3.2. For the system (2.5), we have
(i) If R10<1 and R20<1, there is the disease-free equilibrium E0 only;
(ii) If R10<1 and R20>1, there are the disease-free equilibrium E0 and the boundary equilibrium E1;
(iii) If R10>1 and R20<1, there are the disease-free equilibrium E0 and the epidemic equilibrium E∗;
(iv) If R10>1 and R20>1, there are the disease-free equilibrium E0, the boundary equilibrium E1 and the epidemic equilibrium E∗.
Remark 3.1. Define the basic reproduction number R0 of the system (2.5) by the spectral radius of the next generation matrix [20], we have
R0=ρ(β1d1+γ1+mN∗10−mβ2N∗2(d2+γ2)(d1+γ1+m)β2(d2+γ2)N∗2), |
where ρ(A) denotes the spectral radius of a matrix A. So from the above analysis, we know that R0=max{R10,R20}.
The next, we shall discuss the local stability of the disease-free equilibrium firstly. Then we discuss the global asymptotical stability.
Theorem 3.3. For the system (2.5), we have
(i) If R10<1 and R20<1, the disease-free equilibrium E0 is locally asymptotically stable;
(ii) If R10>1 or R20>1, the disease-free equilibrium E0 is unstable.
Proof. The linearized system of (2.5) at the equilibrium E0 is
{dI1(t)dt=(β1N∗1−d1−m−γ1)I1,dI2(t)dt=(β2N∗2−d2−γ2)I2+mI1. | (3.4) |
The associated characteristic equation of the linearized system of (3.4) at the equilibrium E0 is
F(λ)=|λ−(β1N∗1−(d1+m+γ1))0−mλ−(β2N∗2−d2−γ2)|=0 | (3.5) |
It is easy to see that the two eigenvalues of characteristic Eq (3.5) are
λ1=β1N∗1−(d1+m+γ1)=(R10−1)(d1+m+γ1) |
and
λ2=β2N∗2−d2−γ2=(R20−1)(d2+γ2). |
So, when R10<1 and R20<1, the disease-free equilibrium E0 is locally asymptotically stable; However, if R10>1 or R20>1, the disease-free equilibrium E0 is unstable.
Remark 3.2. From Theorem 3.3, we know that for the system (2.5), if R0<1 the disease-free equilibrium E0 is locally asymptotically stable; if R0>1, the disease-free equilibrium E0 is unstable.
Theorem 3.4. For the system (2.5), if R10<1 and R20<1, the disease-free equilibrium E0 is globally asymptotically stable.
Proof. Since Ii1+αiIi≤Ii for i=1,2, from system (2.5), we can obtain that
{dI1(t)dt≤(β1N∗1−d1−m−γ1)I1,dI2(t)dt≤(β2N∗2−d2−γ2)I2+mI1. | (3.6) |
Define an auxiliary linear system using the right hand side of (3.6) as follows
{dI1(t)dt=(β1N∗1−d1−m−γ1)I1,dI2(t)dt=(β2N∗2−d2−γ2)I2+mI1. |
It can be rewritten as
(I1I2)′=(β1N∗1−d1−m−γ10mβ2N∗2−d2−γ2)(I1I2). | (3.7) |
if R10<1 and R20<1, we can solve (3.7) and know that limt→∞I1(t)=0 and limt→∞I2(t)=0. By the comparison principle [21], we can conclude that when R10<1 and R20<1, all non-negative solutions of (2.5) satisfy limt→∞Ii(t)=0 for i=1,2. So the disease-free equilibrium E0 is globally asymptotically stable.
In this subsection, we will discuss the local stability of the boundary equilibrium firstly. Then discuss the global asymptotical stability.
Theorem 3.5. For the system (2.5), if R10<1 and R20>1, the boundary equilibrium E1 is globally asymptotically stable.
Proof. The Jacobian matrix at the boundary equilibrium E1 of system (2.5) is
J=(β1N∗1−d1−m−γ10m(d2+γ2)(1−R20)−(β2α2+α22)((d2+γ2)(R20−1)(d2+γ2)α2+β2)2(1+α2(d2+γ2)(R20−1)(d2+γ2)α2+β2)2). |
The two eigenvalues of the Jacobian matrix are
λ1=β1N∗1−(d1+m+γ1)=(R01−1)(d1+m+γ1) |
and
λ2=(d2+γ2)(1−R02)−(β2α2+α22)((d2+γ2)(R2−1)(d2+γ2)α2+β2)2(1+α2(d2+γ2)(R2−1)(d2+γ2)α2+β2)2. |
So, when R10<1 and R20>1, λ1<0 and λ2<0. That is the boundary equilibrium E1 is locally asymptotically stable.
For every (I1(0),I2(0))∈Ω, assume the solution of the system (2.5) with initial value (I1(0),I2(0)) is (I1(t),I2(t)). Since
dI1(t)dt=(β1+(d1+γ1)α1+mα1)(β1N∗1−d1−m−γ1β1+(d1+γ1)α1+mα1−I1)I11+α1I1, |
if R10<1, dI1(t)dt<0, then I1(t) is positive and decreasing and limt→∞I1(t)=0. So for sufficiently small positive number ϵ1, there exists a T, such that I1(T)=ϵ1 and I1(t)<ϵ1 when t>T.
The following, we prove that for any ϵ>0, there exists a T∗>T such that |I2(T∗)−β2N∗2−d2−γ2β2+(d2+γ2)α2|<ϵ. And because E1=(0,β2N∗2−d2−γ2β2+(d2+γ2)α2) is locally asymptotically stable, we have E1 is globally asymptotically stable.
Since
dI2(t)dt=(β2+(d2+γ2)α2)(β2N∗2−d2−γ2β2+(d2+γ2)α2−I2)I21+α2I2+mI1, |
if I2(T)<β2N∗2−d2−γ2β2+(d2+γ2)α2, then dI2(t)dt>0 for t>T. So I2(t) is increasing and there exists T∗1 such that |I2(T∗1)−β2N∗2−d2−γ2β2+(d2+γ2)α2|<ϵ;
if I2(T)>β2N∗2−d2−γ2β2+(d2+γ2)α2, there are two cases:
ⅰ) I2(t) is decreasing for t>T. In this case, there exists T∗2>T, such that |I2(T∗2)−β2N∗2−d2−γ2β2+(d2+γ2)α2|<ϵ;
ⅱ) There exists T1>T, such that dI2(T1)dt>0. That is
(β2+(d2+γ2)α2)(β2N∗2−d2−γ2β2+(d2+γ2)α2−I2(T1))I2(T1)1+α2I2(T1)+mI1(T1)>0. |
Since I1(t)<ϵ1 for t>T, we have
(β2+(d2+γ2)α2)(β2N∗2−d2−γ2β2+(d2+γ2)α2−I2(T1))I2(T1)1+α2I2(T1)+mϵ1>0. |
Since I2(T)>β2N∗2−d2−γ2β2+(d2+γ2)α2, we have
(β2+(d2+γ2)α2)(β2N∗2−d2−γ2β2+(d2+γ2)α2−I2(T1))α2+β2+(d2+γ2)α2β2N∗2−d2−γ2+mϵ1>0. |
So
I2(T1)−β2N∗2−d2−γ2β2+(d2+γ2)α2<α2+β2+(d2+γ2)α2β2N∗2−d2−γ2(d2+γ2)α2+β2mϵ1. |
If only ϵ1<(d2+γ2)α2+β2(α2+β2+(d2+γ2)α2β2N∗2−d2−γ2)mϵ, then |I2(T∗2)−β2N∗2−d2−γ2β2+(d2+γ2)α2|<ϵ. It is completed.
In this subsection, we will discuss the local stability of the epidemic equilibrium firstly and then discuss the global asymptotical stability.
Theorem 3.6. For the system (2.5), if R10>1, the epidemic equilibrium E∗ is locally asymptotically stable.
Proof. The Jacobian matrix at the epidemic equilibrium E∗ of system (2.5) is
J=(−β1I1∗1+α1I∗1+β1(N∗1−I∗1)(1+α1I∗1)2−d1−m−γ10mβ2N∗2−2β2I∗2−2β2α2I∗22−(d2+γ2)(1+α2I∗2)2(1+α2I∗2)2)=((d1+m+γ1)1−R10−(β1α1d1+m+γ1+α21)I∗21(1+α1I∗1)20mβ2N∗2−2β2I∗2−2β2α2I∗22−(d2+γ2)(1+α2I∗2)2(1+α2I∗2)2). |
The two eigenvalues of the Jacobian matrix are
λ1=(d1+m+γ1)1−R10−(β1α1d1+m+γ1+α21)I∗21(1+α1I∗1)2 |
and
λ2=β2N∗2−2β2I∗2−2β2α2I∗22−(d2+γ2)(1+α2I∗2)2(1+α2I∗2)2. |
It is easy to see that if R10>1, λ1<0. The next, we need only prove the second eigenvalue λ2<0 if R10>1. Since (1+α2I∗2)2>0, we need only prove β2N∗2−2β2I∗2−2β2α2I∗22−(d2+γ2)(1+α2I∗2)2<0. Let
G(I∗2)=β2N∗2−2β2I∗2−2β2α2I∗22−(d2+γ2)(1+α2I∗2)2. |
Since I∗2≥(β2N∗2−d2−γ2+mα1I∗1)(d2+γ2)α2+β2, so
G(I∗2)=β2N∗2−2β2I∗2−2β2α2I∗22−(d2+γ2)(1+α2I∗2)2≤β2N∗2−(β2+(d2+γ2)α2)(β2N∗2−d2−γ2+mα1I∗1)(d2+γ2)α2+β2−d2−γ2=−mα1I∗1<0. |
This completes the proof.
Theorem 3.7. For the systeml (2.5), if R10>1, the epidemic equilibrium E∗ is globally asymptotically stable.
Proof. Since E∗ is stable when R10>1, we need only prove E∗ is globally attractive.
Consider the equation
dI1(t)dt=I1(β1(N∗1−I1)1+α1I1−d1−m−γ1). |
Let
f1(I1)=β1(N∗1−I1)1+α1I1−d1−m−γ1. |
Then f′1(I1)=β1−1−N∗1α1(1+α1I1)2<0. So f1(I1) is a monotonic decreasing function for all I1>0. Furthermore, f1(0)>0, f1(N∗1)<0 and f1(I∗1)=0 when R10>1. That means if I1∈(0,I∗1), f1(I1)>0 and dI1(t)dt>0; if I1∈(I∗1,N∗1), f1(I1)<0 and dI1(t)dt<0. Hence limt→∞I1(t)=I∗1. By Eq (3.3), limt→∞I2(t)=I∗2. Thus E∗ is globally asymptotically stable.
The results about the existence and stability of equilibria are summarized in Table 1.
Conditions | E0 | E1 | E∗ |
R10<1 and R20<1 | Yes (GAS) | No | No |
R10<1 and R20>1 | Yes (Unstable) | Yes (GAS) | No |
R10>1 and R20<1 | Yes (Unstable) | No | Yes (GAS) |
R10>1 and R20>1 | Yes (Unstable) | Yes (Unstable) | Yes (GAS) |
Remark 3.3. From Theorems 3.5 and 3.7, we know that for the system (2.5), if R0>1, the infectious disease is uniformly persistent. However, the infectious disease is not always uniformly persistent in every patch. If R0>1, but R10<1, the disease is uniformly persistent in the sink patch, but is extinct in the source patch. If R10>1, the disease is always uniformly persistent in every patch. This is a different conclusion resulted by the reducible migration matrix.
Remark 3.4. From Theorem 3.5, in the case that R10<1 and R20>1, the infection does not persist in the source patch but is able to persist in the sink patch. So, in the early stage of the spread of infectious disease, the sink patch should assess the reproduction number R20 reasonably and take control measures timely to prevent the epidemic.
In this section, we carry on numerical simulations to verify the theoretical conclusions, reveal the influence of the migration rate form patch 1 to patch 2 on the basic reproduction number, the transmission scale and transmission speed, and discuss the influence of the parameters α1 and α2 that measure the inhibitory effect on the basic reproduction number, the transmission scale and transmission speed.
To numerically illustrate the theoretical results, we need to choose some parameter values (see Table 2).
Parameter | Description | Value |
A1 | the recruitment rate of the population in patch 1 | 0.018 (Figures 1 and 2) |
0.03 (Figure 3) | ||
A2 | the recruitment rate of the population in patch 2 | 0.0005 (Figure 1) |
0.004 (Figures 2 and 3) | ||
β1 | the transmission rate in patch 1 | 0.00001 |
β2 | the transmission rate in patch 2 | 0.00005 |
α1 | the parameter that measure the inhibitory effect in patch 1 | 0.02 |
α2 | the parameter that measure the inhibitory effect in patch 1 | 0.02 |
d1 | the death rate in patch 1 | 0.0003 |
d2 | the death rate in patch 2 | 0.0003 |
γ1 | the death rate in patch 1 | 0.0001 |
γ2 | the death rate in patch 2 | 0.0001 |
m | the migration rate form the patch 1 to the patch 2 | 0.00005 |
We verify the theoretical conclusions firstly. Denote the density of the infective individuals in patch 1 by i1(t)=I1(t)N∗1. Denote the density of the infective individuals in patch 2 by i2(t)=I2(t)N∗2. Figure 1 shows the evolution of the density of infective individuals in the two patches when R10=0.8889 and R20=0.7500. As predicted by the analytic calculation, the infectious disease in the two patches will disappear eventually. Figure 2 shows the evolution of the density of infective individuals in the two patches when R10=0.8889 and R20=1.8750. We can see the infectious disease will be endemic in patch 2 and the infectious disease in patch 1 will disappear eventually. Figure 3 shows the evolution of the density of infective individuals in the two patches when R10=1.4815 and R20=2.2917. We can see the infectious disease will be endemic in the two patches. And we can see the infectious disease will be endemic in the two patches if R10>1 from the subfigures (c) and (d) of Figure 4.
Second, we reveal the influence of the migration rate m on the transmission in Figure 5. With the increasing of m, the density of infective individuals in patch 1 i1 is decreasing, however the density of infective individuals in patch 2 i2 is increasing.
Third, we reveal the parameters α1 and α2 on the transmission scale and transmission speed. We can see that when α1 is increasing, the density of infective individuals in patch 1 i1 is decreasing from Figure 6 and When α2 is increasing, the density of infective individuals in patch 2 i1 is decreasing from Figure 7.
Many scholars have studied infectious disease transmission with population migration [1,2,3,4,5,6,7], assuming that the migration matrix is irreducible, and found that the propagation dynamics of infectious diseases is determined by the basic reproduction number of the system. When the basic reproduction number is less than 1, the infectious disease eventually becomes extinct; when the basic reproduction number is larger than 1, the infectious disease is epidemic eventually. Since the migration matrix is irreducible, all patches are a connected whole. In all patches, infectious diseases are either extinct or epidemic. That is there is not the phenomenon that infectious diseases are extinct in some patches but epidemic in the others.
Because the studies about the spread of infectious diseases with reducible migration matrix are rare, in this paper, we proposed a two-patch SIS model with saturating contact rate and one-directing population dispersal, discussed the global asymptotic stability of the disease-free equilibrium, the boundary equilibrium and the endemic equilibrium respectively, and revealed the influence of saturating contact rate and migration rate on basic reproduction number and the transmission scale. We have the following main conclusions:
1) If R10>1 then the system tends to a global endemic equilibrium in which infected individuals are present in both patches provided initially there were infected individuals in the source patch; If R10<1 and R20>1 then the system converges to an equilibrium with infected individuals only in the sink patch; If R10<1 and R20<1 then the system converges to the disease-free equilibrium.
2) When migration rate is increasing, the density of infective individuals in the source patch is decreasing; but the density of infective individuals in the sink patch is increasing;
3) With the increasing of the parameter αi (i=1,2) in saturating contact rate, the density of infective individuals in patch i (i=1,2) is decreasing.
The similar conclusions can be obtained for the two patch SI model
{dS1(t)dt=A1−d1S1−β1S1I11+α1I1−mS1,dS2(t)dt=A2−d2S2−β2S2I21+α2I2+mS1,dI1(t)dt=β1S1I11+α1I1−d1I1−mI1,dI2(t)dt=β2S2I21+α2I2−d2I2+mI1. |
We can generalize the current model in many aspects to increase realism. For instance, the infection rate can be given by β(I)SI. We can give the properties on function β(I) such that β(I) is decreasing and tends to 0 when I tends to infinity. The mortality rates of the susceptible and infected individuals are the same in the current model. In fact, the disease-induced death rate can not be neglected sometimes. So the disease-induced death rate can be considered. It is also significant to consider heterogeneous number of contacts for each individual on complex network. There are many paper on this topic [22,23]. One can investigate the multi-patch epidemic model with reducible migration matrix.
This work is supported by the National Natural Sciences Foundation of China (Nos.12001501, 12071445, 11571324, 61603351), Shanxi Province Science Foundation for Youths (201901D211216), the Fund for Shanxi '1331KIRT'.
All authors declare no conflicts of interest in this paper.
[1] |
G. M. Scalera, M. Ferrarin, M. Rabuffetti, Gait regularity assessed by wearable sensors: Comparison between accelerometer and gyroscope data for different sensor locations and walking speeds in healthy subjects, J. Biomech., 113 (2020), 110115. https://doi.org/10.1016/j.jbiomech.2020.110115 doi: 10.1016/j.jbiomech.2020.110115
![]() |
[2] |
Y. P. Demir, S. A. Yildirim, Different walk aids on gait parameters and kinematic analysis of the pelvis in patients with adult neuromuscular disease, Neurosciences, 24 (2019), 36–44. https://doi.org/10.17712/nsj.2019.1.20180316 doi: 10.17712/nsj.2019.1.20180316
![]() |
[3] |
Y. Wang, Z. Liu, J. Xu, W. Yan, Heterogeneous network representation learning approach for ethereum identity identification, IEEE Trans. Comput. Soc. Syst., 10 (2023), 890–899. https://doi.org/10.1109/TCSS.2022.3164719 doi: 10.1109/TCSS.2022.3164719
![]() |
[4] |
K. G. M. Quispe, W. S. Lima, D. M. Batista, E. Souto, MBOSS: A symbolic representation of human activity recognition using mobile sensors, Sensors, 18 (2018). https://doi.org/10.3390/s18124354 doi: 10.3390/s18124354
![]() |
[5] |
S. Bahadori, J. M. Williams, T. W. Wainwright, Lower limb kinematic, kinetic and spatial-temporal gait data for healthy adults using a self-paced treadmill, Data Brief, 34 (2021), 106613. https://doi.org/10.1016/j.dib.2020.106613 doi: 10.1016/j.dib.2020.106613
![]() |
[6] |
M. F. Antwi-Afari, H. Li, Fall risk assessment of construction workers based on biomechanical gait stability parameters using wearable insole pressure system, Adv. Eng. Inf., 38 (2018), 683–694. https://doi.org/10.1016/j.aei.2018.10.002 doi: 10.1016/j.aei.2018.10.002
![]() |
[7] |
H. Ohtsu, S. Yoshida, T. Minamisawa, N. Katagiri, T. Yamaguchi, T. Takahashi, et al., Does the balance strategy during walking in elderly persons show an association with fall risk assessment?, J. Biomech., 103 (2020), 109657. https://doi.org/10.1016/j.jbiomech.2020.109657 doi: 10.1016/j.jbiomech.2020.109657
![]() |
[8] |
J. R. Brickner, J. L. Garzon, K. A. Cimprich, Walking a tightrope: The complex balancing act of R-loops in genome stability, Mol. Cell, 82 (2022), 2267–2297. https://doi.org/10.1016/j.molcel.2022.04.014 doi: 10.1016/j.molcel.2022.04.014
![]() |
[9] |
J. C. Schrijvers, J. C. van den Noort, M. van der Esch, J. Dekker, J. Harlaar, Objective parameters to measure (in)stability of the knee joint during gait: A review of literature, Gait Posture, 70 (2019), 235–253. https://doi.org/10.1016/j.gaitpost.2019.03.016 doi: 10.1016/j.gaitpost.2019.03.016
![]() |
[10] |
J. N. Katz, K. R. Arant, R. F. Loeser, Diagnosis and treatment of hip and knee osteoarthritis: A Review, J. Am. Med. Assoc., 325 (2021), 568–578. https://doi.org/10.1001/jama.2020.22171 doi: 10.1001/jama.2020.22171
![]() |
[11] |
M. Wang, X. Wang, C. Peng, S. Zhang, Z. Fan, Z. Liu, Research on EMG segmentation algorithm and walking analysis based on signal envelope and integral electrical signal, Photonic Netwrk Commun., 37 (2019), 195–203. https://doi.org/10.1007/s11107-018-0809-1 doi: 10.1007/s11107-018-0809-1
![]() |
[12] |
S. Borel, P. Schneider, C. J. Newman, Video analysis software increases the interrater reliability of video gait assessments in children with cerebral palsy, Gait Posture, 33 (2011), 727–729. https://doi.org/10.1016/j.gaitpost.2011.02.012 doi: 10.1016/j.gaitpost.2011.02.012
![]() |
[13] |
S. Chakraborty, A. Nandy, T. Yamaguchi, V. Bonnet, G. Venture, Accuracy of image data stream of a markerless motion capture system in determining the local dynamic stability and joint kinematics of human gait, J. Biomech., 104 (2020), 109718. https://doi.org/10.1016/j.gaitpost.2011.02.012 doi: 10.1016/j.gaitpost.2011.02.012
![]() |
[14] |
C. N. Armitano, H. J. Bennett, J. A. Haegele, S. Morrison, Assessment of the gait-related acceleration patterns in adults with autism spectrum disorder, Gait Posture, 75 (2020), 155–162. https://doi.org/10.1016/j.gaitpost.2019.09.002 doi: 10.1016/j.gaitpost.2019.09.002
![]() |
[15] |
Y. Shi, H. Li, X. Fu, R. Luan, Y. Wang, N. Wang, et al., Self-powered difunctional sensors based on sliding contact-electrification and tribovoltaic effects for pneumatic monitoring and controlling, Nano Energy, 110 (2023), 108339. https://doi.org/10.1016/j.nanoen.2023.108339 doi: 10.1016/j.nanoen.2023.108339
![]() |
[16] |
C. Mei, F. Gao, Y. Li, A determination method for gait event based on acceleration sensors, Sensors, 19 (2019), 5499. https://doi.org/10.3390/s19245499 doi: 10.3390/s19245499
![]() |
[17] |
J. Taborri, E. Palermo, S. Rossi, P. Cappa, Gait partitioning methods: A systematic review, Sensors, 16 (2016), 66. https://doi.org/10.3390/s16010066 doi: 10.3390/s16010066
![]() |
[18] | H. Chen, F. Gao, C. Chen, T. Tian, Estimation of ankle angle based on multi-feature fusion with random forest, in 2018 37th Chinese Control Conference (CCC), IEEE, (2018), 5549–5553. |
[19] |
M. D. Gor-García-Fogeda, R. Cano de la Cuerda, M. Carratalá Tejada, I. M. Alguacil-Diego, F. Molina-Rueda, Observational gait assessments in people with neurological disorders: A systematic review, Arch. Phys. Med. Rehabil., 97 (2016), 131–140. https://doi.org/10.1016/j.apmr.2015.07.018 doi: 10.1016/j.apmr.2015.07.018
![]() |
[20] |
C. R. Brown, S. J. Hillman, A. M. Richardson, J. L. Herman, J. E. Robb, Reliability and validity of the visual gait assessment scale for children with hemiplegic cerebral palsy when used by experienced and inexperienced observers, Gait Posture, 27 (2008), 648–652. https://doi.org/10.1016/j.gaitpost.2007.08.008 doi: 10.1016/j.gaitpost.2007.08.008
![]() |
[21] |
C. H. Lee, S. H. Chen, B. C. Jiang, T. L. Sun, Estimating postural stability using improved permutation entropy via tug accelerometer data for community-dwelling elderly people, Entropy, 22 (2020), 354–365. https://doi.org/10.3390/e22101097 doi: 10.3390/e22101097
![]() |
[22] |
S. Majumder, M. J. Deen, Wearable IMU-based system for real-time monitoring of lower-limb joints, IEEE Sensors J., 21 (2020), 8267–8275. https://doi.org/10.1109/JSEN.2020.3044800 doi: 10.1109/JSEN.2020.3044800
![]() |
[23] |
J. Taborri, J. Keogh, A. Kos, A. Santuz, A. Umek, C. Urbanczyk, E. van der Kruk, S. Rossi, Sport biomechanics applications using inertial, force, and EMG sensors: A literature overview, Appl. Bionic. Biomech., 27 (2020), 65–78. https://doi.org/10.1155/2020/2041549 doi: 10.1155/2020/2041549
![]() |
[24] |
A. Rajkumar, F. Vulpi, S. R. Bethi, H. K. Wazir, P. Raghavan, V. Kapila, Wearable Inertial Sensors for Range of Motion Assessment, IEEE Sensors J., 20 (2020), 3777–3787. https://doi.org/10.1109/JSEN.2019.2960320 doi: 10.1109/JSEN.2019.2960320
![]() |
[25] |
J. Liu, T. Lockhart, S. Kim, Prediction of the spatio-temporal gait parameters using inertial sensor, J. Mech. Med. Biol., 18 (2018), 121–135. https://doi.org/10.1142/S021951941840002X doi: 10.1142/S021951941840002X
![]() |
[26] |
S. Bahadori, J. M. Williams, T. W. Wainwright, Lower limb kinematic, kinetic and spatial-temporal gait data for healthy adults using a self-paced treadmill, Data Brief, 34 (2021), 106613. https://doi.org/10.1016/j.dib.2020.106613 doi: 10.1016/j.dib.2020.106613
![]() |
[27] |
J. Soulard, J. Vaillant, R. Balaguier, N. Vuillerme, Spatio-temporal gait parameters obtained from foot-worn inertial sensors are reliable in healthy adults in single- and dual-task conditions, Sci. Rep., 11 (2021), 10229. https://doi.org/10.1038/s41598-021-88794-4 doi: 10.1038/s41598-021-88794-4
![]() |
[28] |
S. M. Moghadam, T. Yeung, J. Choisne, A comparison of machine learning models' accuracy in predicting lower-limb joints' kinematics, kinetics, and muscle forces from wearable sensors, Sci. Rep., 13 (2023), 5046. https://doi.org/10.1038/s41598-023-31906-z doi: 10.1038/s41598-023-31906-z
![]() |
[29] |
G. Salatino, E. Bergamini, T. Marro, P. Gentili, M. Iosa, D. Morelli, et al., Gait stability assessment in Down and Prader-Willi syndrome children using inertial sensors, Gait Posture, 49 (2016), S16. https://doi.org/10.1016/j.gaitpost.2016.07.046 doi: 10.1016/j.gaitpost.2016.07.046
![]() |
[30] |
J. Johansson, A. Nordström, P. Nordström, Greater fall risk in elderly women than in men is associated with increased gait variability during multitasking, J. Am. Med., 17 (2016), 535–540. https://doi.org/10.1016/j.jamda.2016.02.009 doi: 10.1016/j.jamda.2016.02.009
![]() |
[31] |
N. J. Nelms, C. E. Birch, D. H. Halsey, M. Blankstein, B. D. Beynnon, Assessment of early gait recovery after anterior approach compared to posterior approach total hip arthroplasty: A smartphone accelerometer-based study, J. Arthroplasty, 35 (2019), 125–138. https://doi.org/10.1016/j.arth.2019.09.030 doi: 10.1016/j.arth.2019.09.030
![]() |
[32] |
P. Tamburini, F. Storm, C. Buckley, M. C. Bisi, R. Stagni, C. Mazzà, Moving from laboratory to real life conditions: Influence on the assessment of variability and stability of gait, Gait Posture, 59 (2018), 248–252. https://doi.org/10.1016/j.gaitpost.2017.10.024 doi: 10.1016/j.gaitpost.2017.10.024
![]() |
[33] |
N. Muthukrishnan, J. J. Abbas, N. Krishnamurthi, A wearable sensor system to measure step-based gait parameters for parkinson's disease rehabilitation, Sensors, 20 (2020), 6417. https://doi.org/10.3390/s20226417 doi: 10.3390/s20226417
![]() |
[34] |
J. Y. Wang, D. W. Gong, H. C. Luo, W. B. Zhang, L. Zhang, H. Zhang, et al., Measurement of step angle for quantifying the gait impairment of parkinson's disease by wearable sensors: Controlled study, JMIR mHealth uHealth, 8 (2020), 10–25. https://doi.org/10.2196/16650 doi: 10.2196/16650
![]() |
[35] |
A. Nguyen, N. Roth, N. H. Ghassemi, J. Hannink, T. Seel, J. Klucken, et al., Development and clinical validation of inertial sensor-based gait-clustering methods in Parkinson's disease, J. Neuroeng. Rehabil., 16 (2019), 147–159. https://doi.org/10.1186/s12984-019-0548-2 doi: 10.1186/s12984-019-0548-2
![]() |
[36] |
P. Caliandro, C. Conte, C. Iacovelli, A. Tatarelli, S. F. Castiglia, G. Reale, et al., Exploring risk of falls and dynamic unbalance in cerebellar ataxia by inertial sensor assessment, Sensors, 19 (2019), 5571. https://doi.org/10.3390/s19245571 doi: 10.3390/s19245571
![]() |
[37] |
P. Tamburini, F. Storm, C. Buckley, M. C. Bisi, R. Stagni, C. Mazza, Moving from laboratory to real life conditions: Influence on the assessment of variability and stability of gait, Gait Posture, 59 (2018), 248–252. https://doi.org/10.1016/j.gaitpost.2017.10.024 doi: 10.1016/j.gaitpost.2017.10.024
![]() |
[38] |
M. Diopa, A. Rahmani, A. Belli, V. Gautheron, A. Geyssant, J. Cottalorda, Influence of speed variation and age on the asymmetry of ground reaction forces and stride parameters of normal gait in children, J. Pediatric Orthopaedics-Part B, 13 (2004), https://doi.org/10.1097/01202412-200409000-00005 doi: 10.1097/01202412-200409000-00005
![]() |
[39] |
J. Zhao, Y. Lv, Output-feedback robust tracking control of uncertain systems via adaptive learning, Int. J. Control Autom. Syst., 21 (2023), 1108–1118. https://doi.org/10.1007/s12555-021-0882-6 doi: 10.1007/s12555-021-0882-6
![]() |
[40] |
Y. Shi, L. Li, J. Yang, Y. Wang, S. Hao, Center-based transfer feature learning with classifier adaptation for surface defect recognition, Mech. Syst. Signal Process., 188 (2023), 110001. https://doi.org/10.1016/j.ymssp.2022.110001 doi: 10.1016/j.ymssp.2022.110001
![]() |
[41] |
Z. Liu, D. Yang, Y. Wang, M. Lu, R. Li, EGNN: Graph structure learning based on evolutionary computation helps more in graph neural networks, Appl. Soft Comput., 135 (2023), 110040. https://doi.org/10.1016/j.asoc.2023.110040 doi: 10.1016/j.asoc.2023.110040
![]() |
[42] |
C. Tian, Z. Xu, L. Wang, Y. Liu, Arc fault detection using artificial intelligence: Challenges and benefits, Math. Biosci. Eng., 20 (2023), 12404–12432. https://doi.org/10.3934/mbe.2023552 doi: 10.3934/mbe.2023552
![]() |
[43] |
A. S. Alharthi, S. U. Yunas, K. B. Ozanyan, Deep learning for monitoring of human gait: A review, IEEE Sensors J., 19 (2019), 9575–9591. https://doi.org/10.1109/JSEN.2019.2928777 doi: 10.1109/JSEN.2019.2928777
![]() |
[44] |
T. Yao, F. Gao, Q. Zhang, Y. Ma, Multi-feature gait recognition with DNN based on sEMG signals, Math. Biosci. Eng., 18 (2021), 3521–3542. 18(4). https://doi.org/10.3934/mbe.2021177 doi: 10.3934/mbe.2021177
![]() |
[45] |
E. Sansano, R. Montoliu, Ó. B. Fernández, A study of deep neural networks for human activity recognition, Comput. Intell., 36 (2020), 1113–1139. https://doi.org/10.1111/coin.12318 doi: 10.1111/coin.12318
![]() |
[46] |
J. N. Mogan, C. P. Lee, K. M. Lim, K. S. Muthu, Gait-ViT: Gait recognition with vision transformer, Sensors, 22 (2022), 362. https://doi.org/10.3390/s22197362 doi: 10.3390/s22197362
![]() |
[47] |
R. Moe-Nilssen, A new method for evaluating motor control in gait under real-life environmental conditions. Part 1: The instrument, Clin. Biomech., 13 (1998), 320–327. https://doi.org/10.1016/S0268-0033(98)00089-8 doi: 10.1016/S0268-0033(98)00089-8
![]() |
[48] |
M. I. Esfahani, M. A. Nussbaum, Using smart garments to differentiate among normal and simulated abnormal gaits, J. Biomech., 93 (2019), 70–76. https://doi.org/10.1016/j.jbiomech.2019.06.009 doi: 10.1016/j.jbiomech.2019.06.009
![]() |
[49] |
M. S. Jia, T. H. Li, J. Wang, Audio fingerprint extraction based on locally linear embedding for audio retrieval system, Electronics, 9 (2020), 238–253. https://doi.org/10.1016/j.imu.2018.10.002 doi: 10.1016/j.imu.2018.10.002
![]() |
[50] |
P. A. Semblantes, V. H. Andaluz, J. Lagla, F. A. Chicaiza, A. Acurio, Visual feedback framework for rehabilitation of stroke patients, Inf. Med. Unlocked, 13 (2018), 41–50. https://doi.org/10.1016/j.imu.2018.10.002 doi: 10.1016/j.imu.2018.10.002
![]() |
[51] |
C. J. Su, C. Y. Chiang, J. Y. Huang, Kinect-enabled home-based rehabilitation system using dynamic time warping and fuzzy logic, Appl. Soft Comput., 22 (2014), 652–666. https://doi.org/10.1016/j.asoc.2014.04.020 doi: 10.1016/j.asoc.2014.04.020
![]() |
[52] | J. Minsu, D. Kim, Y. Kim, K. Jaehong, Automated dance motion evaluation using dynamic time warping and Laban movement analysis, in IEEE International Conference on Consumer Electronics, (2017), 141–142. |
[53] |
R. Haghighi Osgouei, D. Soulsby, F. Bello, Rehabilitation exergames: Use of motion sensing and machine learning to quantify exercise performance in healthy volunteers, JMIR Rehabil. Assistive Technol., 7 (2020), e17289. https://doi.org/10.2196/17289 doi: 10.2196/17289
![]() |
[54] |
I. Hagoort, N. Vuillerme, T. Hortobágyi, C. J. Lamoth, Outcome-dependent effects of walking speed and age on quantitative and qualitative gait measures, Gait Posture, 93 (2022), 39–46. https://doi.org/10.1016/j.gaitpost.2022.01.001 doi: 10.1016/j.gaitpost.2022.01.001
![]() |
[55] |
Z. Yan, X. Xu, Y. Wang, T. Li, B. Ma, L. Yang, et al., Application of ultrasonic doppler technology based on wavelet threshold denoising algorithm in fetal heart rate and central nervous system malformation detection, World Neurosurg., 30 (2020), 168–179. https://doi.org/10.1016/j.wneu.2020.10.030 doi: 10.1016/j.wneu.2020.10.030
![]() |
[56] |
R. Takeda, S. Tadano, M. Todoh, M. Morikawa, M. Nakayasu, S. Yoshinari, Gait analysis using gravitational acceleration measured by wearable sensors, J. Biomech., 42 (2009), 223–233. https://doi.org/10.1016/j.jbiomech.2008.10.027 doi: 10.1016/j.jbiomech.2008.10.027
![]() |
[57] |
Y. Liu, G. Yin, The Delaunay triangulation learner and its ensembles, Comput. Stat. Data Anal., 152 (2020), 1121–1135. https://doi.org/10.1016/j.csda.2020.107030 doi: 10.1016/j.csda.2020.107030
![]() |
[58] |
X. Yu, S. Xiong, A dynamic time warping based algorithm to evaluate kinect-enabled home-based physical rehabilitation exercises for older people, Sensors, 19 (2019), 2882. https://doi.org/10.3390/s19132882 doi: 10.3390/s19132882
![]() |
[59] |
T. Seel, J. Raisch, T. Schauer, IMU-based joint angle measurement for gait analysis, Sensors, 14 (2014), 6891–6909. https://doi.org/10.3390/s140406891 doi: 10.3390/s140406891
![]() |
[60] |
V. B. Semwal, A. Gupta, P. Lalwani, An optimized hybrid deep learning model using ensemble learning approach for human walking activities recognition, J. Supercomput., 77 (2021), 12256–12279. https://doi.org/10.1007/s11227-021-03768-7 doi: 10.1007/s11227-021-03768-7
![]() |
[61] |
X. Yu, J. Jang, S. Xiong, A large-scale open motion dataset (KFall) and Benchmark algorithms for detecting pre-impact fall of the elderly using wearable inertial sensors, Front. Aging Neurosci., 13 (2021), 692865. https://doi.org/10.3389/fnagi.2021.692865 doi: 10.3389/fnagi.2021.692865
![]() |
[62] |
R. J. Kate, Using dynamic time warping distances as features for improved time series classification, Data Mining Knowl. Discovery, 30 (2015), 283–312. https://doi.org/10.1007/s10618-015-0418-x doi: 10.1007/s10618-015-0418-x
![]() |
[63] |
S. J. Dixon, R. S. Hinman, M. W. Creaby, G. Kemp, K. M. Crossley, Knee joint stiffness during walking in knee osteoarthritis, Arth. Care Res., 62 (2010), 38–44. https://doi.org/10.1002/acr.20012 doi: 10.1002/acr.20012
![]() |
[64] |
D. H. Ro, T. Kang, D. Han, D. Y. Lee, H. S. Han, M. C. Lee, Quantitative evaluation of gait features after total knee arthroplasty: Comparison with age and sex-matched controls, Gait Posture, 75 (2020), 78–84. https://doi.org/10.1016/j.gaitpost.2019.09.026 doi: 10.1016/j.gaitpost.2019.09.026
![]() |
Conditions | E0 | E1 | E∗ |
R10<1 and R20<1 | Yes (GAS) | No | No |
R10<1 and R20>1 | Yes (Unstable) | Yes (GAS) | No |
R10>1 and R20<1 | Yes (Unstable) | No | Yes (GAS) |
R10>1 and R20>1 | Yes (Unstable) | Yes (Unstable) | Yes (GAS) |
Parameter | Description | Value |
A1 | the recruitment rate of the population in patch 1 | 0.018 (Figures 1 and 2) |
0.03 (Figure 3) | ||
A2 | the recruitment rate of the population in patch 2 | 0.0005 (Figure 1) |
0.004 (Figures 2 and 3) | ||
β1 | the transmission rate in patch 1 | 0.00001 |
β2 | the transmission rate in patch 2 | 0.00005 |
α1 | the parameter that measure the inhibitory effect in patch 1 | 0.02 |
α2 | the parameter that measure the inhibitory effect in patch 1 | 0.02 |
d1 | the death rate in patch 1 | 0.0003 |
d2 | the death rate in patch 2 | 0.0003 |
γ1 | the death rate in patch 1 | 0.0001 |
γ2 | the death rate in patch 2 | 0.0001 |
m | the migration rate form the patch 1 to the patch 2 | 0.00005 |
Conditions | E0 | E1 | E∗ |
R10<1 and R20<1 | Yes (GAS) | No | No |
R10<1 and R20>1 | Yes (Unstable) | Yes (GAS) | No |
R10>1 and R20<1 | Yes (Unstable) | No | Yes (GAS) |
R10>1 and R20>1 | Yes (Unstable) | Yes (Unstable) | Yes (GAS) |
Parameter | Description | Value |
A1 | the recruitment rate of the population in patch 1 | 0.018 (Figures 1 and 2) |
0.03 (Figure 3) | ||
A2 | the recruitment rate of the population in patch 2 | 0.0005 (Figure 1) |
0.004 (Figures 2 and 3) | ||
β1 | the transmission rate in patch 1 | 0.00001 |
β2 | the transmission rate in patch 2 | 0.00005 |
α1 | the parameter that measure the inhibitory effect in patch 1 | 0.02 |
α2 | the parameter that measure the inhibitory effect in patch 1 | 0.02 |
d1 | the death rate in patch 1 | 0.0003 |
d2 | the death rate in patch 2 | 0.0003 |
γ1 | the death rate in patch 1 | 0.0001 |
γ2 | the death rate in patch 2 | 0.0001 |
m | the migration rate form the patch 1 to the patch 2 | 0.00005 |