
Citation: Weijie Liu, Yan Hao, Jihong Jiang, Cong Liu, Aihua Zhu, Jingrong Zhu, Zhen Dong. Biopolymeric flocculant extracted from potato residues using alkaline extraction method and its application in removing coal fly ash from ash-flushing wastewater generated from coal fired power plant[J]. AIMS Environmental Science, 2017, 4(1): 27-41. doi: 10.3934/environsci.2017.1.27
[1] | Kishor D. Kucche, Sagar T. Sutar, Kottakkaran Sooppy Nisar . Analysis of nonlinear implicit fractional differential equations with the Atangana-Baleanu derivative via measure of non-compactness. AIMS Mathematics, 2024, 9(10): 27058-27079. doi: 10.3934/math.20241316 |
[2] | M. Mallika Arjunan, Nabil Mlaiki, V. Kavitha, Thabet Abdeljawad . On fractional state-dependent delay integro-differential systems under the Mittag-Leffler kernel in Banach space. AIMS Mathematics, 2023, 8(1): 1384-1409. doi: 10.3934/math.2023070 |
[3] | Veliappan Vijayaraj, Chokkalingam Ravichandran, Thongchai Botmart, Kottakkaran Sooppy Nisar, Kasthurisamy Jothimani . Existence and data dependence results for neutral fractional order integro-differential equations. AIMS Mathematics, 2023, 8(1): 1055-1071. doi: 10.3934/math.2023052 |
[4] | Muneerah AL Nuwairan, Ahmed Gamal Ibrahim . The weighted generalized Atangana-Baleanu fractional derivative in banach spaces- definition and applications. AIMS Mathematics, 2024, 9(12): 36293-36335. doi: 10.3934/math.20241722 |
[5] | Miguel Vivas-Cortez, Muhammad Uzair Awan, Sehrish Rafique, Muhammad Zakria Javed, Artion Kashuri . Some novel inequalities involving Atangana-Baleanu fractional integral operators and applications. AIMS Mathematics, 2022, 7(7): 12203-12226. doi: 10.3934/math.2022678 |
[6] | Muhammad Farman, Ali Akgül, Kottakkaran Sooppy Nisar, Dilshad Ahmad, Aqeel Ahmad, Sarfaraz Kamangar, C Ahamed Saleel . Epidemiological analysis of fractional order COVID-19 model with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(1): 756-783. doi: 10.3934/math.2022046 |
[7] | Abdelatif Boutiara, Mohammed M. Matar, Jehad Alzabut, Mohammad Esmael Samei, Hasib Khan . On ABC coupled Langevin fractional differential equations constrained by Perov's fixed point in generalized Banach spaces. AIMS Mathematics, 2023, 8(5): 12109-12132. doi: 10.3934/math.2023610 |
[8] | A. M. Sayed Ahmed, Hamdy M. Ahmed, Nesreen Sirelkhtam Elmki Abdalla, Assmaa Abd-Elmonem, E. M. Mohamed . Approximate controllability of Sobolev-type Atangana-Baleanu fractional differential inclusions with noise effect and Poisson jumps. AIMS Mathematics, 2023, 8(10): 25288-25310. doi: 10.3934/math.20231290 |
[9] | Mdi Begum Jeelani, Abeer S. Alnahdi, Mohammed A. Almalahi, Mohammed S. Abdo, Hanan A. Wahash, M. A. Abdelkawy . Study of the Atangana-Baleanu-Caputo type fractional system with a generalized Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(2): 2001-2018. doi: 10.3934/math.2022115 |
[10] | Lahcene Rabhi, Mohammed Al Horani, Roshdi Khalil . Existence results of mild solutions for nonlocal fractional delay integro-differential evolution equations via Caputo conformable fractional derivative. AIMS Mathematics, 2022, 7(7): 11614-11634. doi: 10.3934/math.2022647 |
With the development of neuroscience, the controlling mechanism and mode of biological motion have been paid much attention by biologists [1,2,3,4], and rhythmic movement is a common mode of motion in biology. Rhythmic movement refers to periodic movement with symmetry of time and space, such as walking, running, jumping, flying, swimming and so on. Biologists have shown that rhythmic movement is not related to the consciousness of the brain, but to the self-excitation of the lower nerve centers. It is a spatiotemporal motion mode controlled by a central pattern generator located in the spinal cord of vertebrates or in the thoracic and abdominal ganglia of invertebrates [5]. They have the ability to automatically generate complex high dimensional control signals for the coordination of the muscles during rhythmic movements [6,7,8,9].
In engineering, CPG can be regarded as a distributed system consisting of a group of coupled nonlinear oscillators. The generation of rhythmic signals can be realized by phase coupling. Changing the coupling relationship of oscillators can produce spatiotemporal sequence signals with different phase relations, and realize different movement modes. CPG of animals lays a foundation for the research of bionic robots. For example, in [10,11] the gait control of quadruped robots based on CPG is studied. Mathematically, there are several common types of CPG oscillators systems, such as Hopf oscillators systems [12,13], Kimura oscillators systems, Rayleigh oscillators systems, Matsuoa oscillators systems and VDP oscillator systems [14,15], etc.
Quadrupedal gait is a kind of gait that people are very concerned. The gait of quadruped is an important type described by a symmetrical system [16,17,18]. For example, in [17,18], base on the symmetry property, the primary and secondary gait modes of quadruped are described, respectively. In animal gait movement, the legs are coupled with each other, and the coupling strength affects the complexity of animal gait. In this paper, the delay of leg signal is considered according to CPG model, the basic gait CPG model of a class of quadruped is constructed by using VDP oscillators, and the ranges of coupling strength between legs under four basic gaits are given. This paper is organized as follows. Firstly, a kind of delay CPG network system is constructed by using VDP oscillator. Secondly, the conditions of Hopf bifurcation in VDP-CPG network corresponding to the four basic gaits are given, and the coupling ranges between legs in four basic gaits are given. Finally, the theoretical results are supported by numerical simulations.
CPG, as the control center of rhythmic motion, is a kind of neural network that can generate the output of rhythmic mode without sensor feedback. It sends out motion instructions from the high-level center to control the initial state of rhythmic motion, and integrates the feedback information and perception information of CPG to regulate the motion organically. The CPG network in this paper adopts the following network structure [14].
In Figure 1, LF, RF, LH and RH represent the animal's left foreleg, right foreleg, left hind leg and right hind leg, respectively. The black arrows represent the leg raising sequence, and the numbers in the circles are the phase difference between other legs and LF leg. In order to generate the rhythmic signal of each leg, the VDP oscillator used in this paper can refer to [14], the equation is as follows.
{˙x=y,˙y=α(p2−x2)˙x−w2x, |
where x is the output signal from oscillator, α, p and w are variable parameters which can influence the character of oscillators. Commonly, the shape of the wave is affected by parameter α, and the amplitude of an output counts on the parameter p mostly. The output frequency is mainly relying on the parameter w when the amplitude parameter p is fixed. But the alteration of parameter p can lightly change the frequency of the signal, and α also can effect the output frequency.
Four-legged muscle groups are regarded as VDP oscillators for feedback motion signals, respectively. The animal's left foreleg, right foreleg, right hind leg and left hind leg are recorded as oscillator x1,x2,x3 and x4, respectively.
Then the oscillator of the ith leg is as follows
{˙xi=yi,˙yi=αi(p2i−x2ki)yi−w2ixki,i=1,2,3,4, |
where xki=xi+4∑j=1,j≠iKijxj denotes the coupling variable. Here Kij is the coupling coefficient, which represents strength of coupling from j oscillator to i oscillator.
Because the motion state of each leg depends on the motion state of the other three legs in the past short time, the time delay is introduced as follows
xki=xi(t)+4∑j=1,j≠iKijxj(t−τ). |
Assuming that the biological mechanism of each leg is similar and the degree of excitation or inhibition is the same between legs, and the excitation is positive coupling, then the inhibition is negative coupling. Therefore,
α1=α2=α3=α4=α, |
p1=p2=p3=p4=p, |
w1=w2=w3=w4=w, |
Kij={K,whenthejlegexcitestheileg,−K,whenthejlegrestrainstheileg.K>0. |
Thus, we study the following VDP-CPG system
{˙xi=yi,˙yi=α(p2−(xi(t)+4∑j=1,j≠iKijxj(t−τ))2)yi−w2(xi(t)+4∑j=1,j≠iKijxj(t−τ)), | (1) |
where i=1,2,3,4. It is clear that the origin (0, 0, 0, 0, 0, 0, 0, 0) is an equilibrium of Eq (1).
In this section, we construct a VDP-CPG network which is used for generation four basic gaits patterns (walk, trot, pace and bound). Then we analyze the conditions for four gait systems to produce Hopf bifurcation.
In order to analyses the four basic gaits, we make the following assumptions.
(H1) h<0,
(H2) 2s−h2>0,19m<K2<19,
(H3) K2<m,
where h=αp2,s=w2,m=4h2s−h44s2.
In walking gait, one leg is inhibited by the other three legs, then there are
Kij=−K,i,j=1,2,3,4,i≠j. |
So the VDP-CPG network in walking gait is as follows
{˙xi=yi,˙yi=α(p2−(xi(t)+4∑j=1,j≠i(−K)xj(t−τ))2)yi−w2(xi(t)+4∑j=1,j≠i(−K)xj(t−τ)). | (2) |
This is a symmetric system. We first explore the symmetry of system (2), then study the existence of Hopf bifurcation of system (2).
Let Yi=(xiyi)∈R2,i=1,2,3,4, system (2) can be written in block form as follows
˙Yi=MYi(t)+NYi+1(t−τ)+NYi+2(t−τ)+NYi+3(t−τ)+g(Yi(t)),i=1,2,3,4(mod4), | (3) |
where
M=(01−w2αp2),N=(00Kw20), |
g(xiyi)=(0−α(xi−Kxi+1(t−τ)−Kxi+2(t−τ)−Kxi+3(t−τ))2yi). |
Let Γ be a compact Lie group. It follows from [19], system ˙u(t)=G(ut) is said to be Γ− equivariant if G(γut)=γG(ut) for all γ∈Γ. Let Γ=D4 be the dihedral group of order 8, which is generated by the cyclic group Z4 of order 4 together with the flip of order 2. Denote by ρ the generator of the cyclic subgroup Z4 and k the flip. Define the action of D4 on R8 by
(ρU)i=Ui+1,(kU)i=U6−i,Ui∈R2,i=1,2,3,4(mod4). |
Then it is easy to get the following lemma.
Lemma 3.1. System (3) is D4− equivariant.
The linearization of Eq (3) at the origin is
˙Yi=MYi(t)+NYi+1(t−τ)+NYi+2(t−τ)+NYi+3(t−τ),i=1,2,3,4(mod4). | (4) |
The characteristic matrix of Eq (4) is given by
A(τ,λ)=(λI2−M−Ne−λτ−Ne−λτ−Ne−λτ−Ne−λτλI2−M−Ne−λτ−Ne−λτ−Ne−λτ−Ne−λτλI2−M−Ne−λτ−Ne−λτ−Ne−λτ−Ne−λτλI2−M), |
where I2 is a 2×2 identity matrix. This is a block circulant matrix, from [20], we have
det(A(τ,λ))=3∏j=0det(λI2−M−χjNe−λτ−(χj)2Ne−λτ−(χj)3Ne−λτ), |
where χj=eπj2i,i is the imaginary unit. The characteristic equation of Eq (4) at the zero solution is
Δ(τ,λ)=det(A(τ,λ))=Δ1(Δ2)3, | (5) |
with
Δ1=λ(λ−h)+s(1−3Ke−λτ)),Δ2=λ(λ−h)+s(1+Ke−λτ),h=αp2,s=w2. |
Lemma 3.2. If (H1) and (H2) hold, for the equation Δ1=0, we have the following results.
(1) when τ=0, all roots of equation Δ1=0 have negative real parts,
(2) when τ>0, there exist τj, such that when τ=τj(j=0,1,2,…), Δ1(±iβ)=0 holds,
(3) the transversality condition:
Re(dλdτ)|λ=iβ+,τ=τjwalk+>0,Re(dλdτ)|λ=iβ−,τ=τjwalk−<0, |
where
β=β±=√2s−h2±√(h2−2s)2−4s2(1−9K2)2, |
τj=τjwalk±=1β±(−arccoss−β2±3Ks+2jπ+2π),j=0,1,2,…. |
Proof. (1) When τ=0, equation Δ1=0 becomes λ(λ−h)+s(1−3K)=0, and the solution is obtained as follows
λ=h±√h2−4s(1−3K)2. |
By (H1) and (H2), the roots of equation Δ1=0 have negative real parts.
(2) When τ>0, let λ=iβ(β>0) be a root of Δ1=0. Substituting iβ into Δ1=0, then we have
−β2−iβh+s(1−3Ke−iβτ)=0. |
Separating the real and imaginary parts, we get the following form
{s−β2=3Kscos(βτ),βh=3Kssin(βτ). | (6) |
If (H2) holds, by solving the above equation, we have
β±=β=√2s−h2±√(h2−2s)2−4s2(1−9K2)2, | (7) |
τjwalk±=τj=1β(−arccoss−β23Ks+2jπ+2π),j=0,1,2,…. |
(3) Let λ(τ)=α(τ)+iβ(τ) be the root of equation Δ1=0, satisfying α(τj)=0 and β(τj)=β. Taking the derivative of the equation Δ1=0 with respect to τ, we can get
dλdτ=−3Ksλeλτ2λ−h+3Ksτeλτ. |
Then
Re(dλdτ)|λ=iβ,τ=τj=3Ksβhsin(βτj)−6Ksβ2cos(βτj)(−h+3Ksτcos(βτj))2+(2β−3Ksτjsin(βτj))2, |
by (6) and (7), we have
Re(dλdτ)|λ=iβ+,τ=τjwalk+>0,Re(dλdτ)|λ=iβ−,τ=τjwalk−<0, |
which means that the transversality condition holds at τjwalk±,j=0,1,2,….
The lemma 3.2 holds.
Lemma 3.3. For Δ2=0, we have the following results.
(1) if (H1) holds, when τ=0 all roots of equation Δ2=0 have negative real parts,
(2) if (H3) holds, when τ>0 equation Δ2=0 has no pure imaginary root.
Proof. (1) When τ=0, equation Δ2=0 becomes λ(λ−h)+s(1+K)=0, and the solution is obtained as follows
λ=h±√h2−4s(1+K)2. |
By (H1), the roots of equation Δ2=0 have negative real parts.
(2) When τ>0, let λ=iβ(β>0) be a root of Δ2=0. Substituting iβ into Δ2=0 then we have
−β2−iβh+s(1+Ke−iβτ)=0. |
The real and imaginary parts of the above equation are separated, then we obtain
{s−β2=−Kscos(βτ),βh=−Kssin(βτ). |
By solving the above equation, we have
β=√2s−h2±√(h2−2s)2−4s2(1−K2)2. |
By (H3), we obtain (h2−2s)2−4s2(1−K2)<0, then the formula above is not valid. So the lemma 3.3 holds.
From lemma 3.2 and 3.3, we have following theorem.
Theorem 3.1. If (H1), (H2) and (H3) hold, then we have the following results.
(1) all roots of Eq (5) have negative real parts for 0≤τ<τ0walk, and at least a pair of roots with positive real parts for τ∈(τ0walk,τ0walk+ε), for some ε>0,
(2) zero equilibrium of system (2) is asymptotically stable for 0≤τ<τ0walk, and unstable for τ∈(τ0walk,τ0walk+ε), for some ε>0,
(3) when τ=τ0walk, system (2) undergoes a Hopf bifurcation at zero equilibrium, where τ0walk=min{τ0walk+,τ0walk−}.
Remark 3.1. Near the critical value τ=τ0walk, the periodic solution of system (2) at the origin accords with walking gait.
In a trot, a leg on the same diagonal as the current leg stimulates the current leg, and two legs on the other diagonal suppress the current leg, thus
K12=−K,K13=K,K14=−K,K21=−K,K23=−K,K24=K,K31=K,K32=−K,K34=−K,K41=−K,K42=K,K43=−K. |
The VDP-CPG network for trotting is as follows.
{˙xi=yi,˙yi=αp2yi−w2(xi(t)+(−K)xi+1(t−τ)+Kxi+2(t−τ)+(−K)xi+3(t−τ))−α(xi(t)+(−K)xi+1(t−τ)+Kxi+2(t−τ)+(−K)xi+3(t−τ))2yi. | (8) |
This is also a symmetric system. Similarly, by lemma 3.1, we have
Lemma 3.4. System (8) is D4− equivariant.
The characteristic matrix of linearization of Eq (8) is given by
A1(τ,λ)=(λI2−M−Ne−λτNe−λτ−Ne−λτ−Ne−λτλI2−M−Ne−λτNe−λτNe−λτ−Ne−λτλI2−M−Ne−λτ−Ne−λτNe−λτ−Ne−λτλI2−M). |
This is a block circulant matrix, and we have
det(A1(τ,λ))=3∏j=0det(λI2−M−χjNe−λτ+(χj)2Ne−λτ−(χj)3Ne−λτ), |
with χj=eπj2i.
The characteristic equation of linearization of Eq (8) at zero solution is
Δ(τ,λ)=det(A1(τ,λ))=Δ3(Δ4)3, | (9) |
where
Δ3=λ(λ−h)+s(1+3Ke−λτ), |
Δ4=λ(λ−h)+s(1−Ke−λτ). |
Similarly, by lemma 3.2 and 3.3, we have following lemmas.
Lemma 3.5. For the equation Δ3=0, we have the following results.
(1) if (H1) holds, when τ=0, all roots of equation Δ3=0 have negative real parts,
(2) if (H2) holds, when τ>0, there exist τj, such that when τ=τj(j=0,1,2,…), Δ3(±iβ)=0 holds,
(3) the transversality condition:
Re(dλdτ)|λ=iβ+,τ=τjtrot+>0,Re(dλdτ)|λ=iβ−,τ=τjtrot−<0, |
where
β=β±=√2s−h2±√(h2−2s)2−4s2(1−9K2)2, |
τj=τjtrot±=1β±(arccoss−β2±−3Ks+2jπ),j=0,1,2,… |
Lemma 3.6. For Δ4=0, we have the following results.
(1) if (H1)and K<1 hold, when τ=0, all roots of equation Δ4=0 have negative real parts,
(2) if (H3) holds, when τ>0, equation Δ4=0 has no pure imaginary root.
From lemma 3.5 and 3.6, we have following theorem.
Theorem 3.2. If (H1), (H2) and (H3) hold, we have the following results.
(1) all roots of Eq (9) have negative real parts for 0≤τ<τ0trot, and at least a pair of roots with positive real parts for τ∈(τ0trot,τ0trot+ε), for some ε>0,
(2) zero equilibrium of Eq (8) is asymptotically stable for 0≤τ<τ0trot, and unstable for τ∈(τ0trot,τ0trot+ε), for some ε>0,
(3) when τ=τ0trot, system (8) undergoes a Hopf bifurcation at zero equilibrium, where τ0trot=min{τ0trot+,τ0trot−}
Remark 3.2. Near the critical value τ=τ0trot, the periodic solution of system (8) at the origin accords with trotting gait.
In a pace, the leg on the same side (left or right) of the current leg stimulates the current leg, and the other two legs inhibit the current leg, thus
K12=−K,K13=−K,K14=K,K21=−K,K23=K,K24=−K,K31=−K,K32=K,K34=−K,K41=K,K42=−K,K43=−K. |
Thus Eq (1) becomes the following VDP-CPG pacing system.
{˙xi=yi,˙yi=α(p2−(xi(t)+(−K)xi+1(t−τ)+(−K)xi+2(t−τ)+Kxi+3(t−τ))2)yi−w2(xi(t)+(−K)xi+1(t−τ)+(−K)xi+2(t−τ)+Kxi+3(t−τ)),i=1,3(mod4) | (10) |
{˙xi=yi,˙yi=α(p2−(xi(t)+(−K)xi−1(t−τ)+Kxi+1(t−τ)+(−K)xi+2(t−τ))2)yi−w2(xi(t)+(−K)xi−1(t−τ)+Kxi+1(t−τ)+(−K)xi+2(t−τ)),i=2,4(mod4) |
and the linearization of Eq (10) at the origin is
{˙xi=yi,˙yi=αp2yi−w2(xi(t)+(−K)xi+1(t−τ)+(−K)xi+2(t−τ)+Kxi+3(t−τ))i=1,3(mod4) | (11) |
{˙xi=yi,˙yi=αp2yi−w2(xi(t)+(−K)xi−1(t−τ)+Kxi+1(t−τ)+(−K)xi+2(t−τ)),i=2,4(mod4) |
the characteristic equation of system (11) is
|Rm−m−m+m−Rm+m−m−m+Rm−m+m−m−R|=Δ5(Δ6)3=0, | (12) |
where
Δ5=λ(λ−h)+s(1+3Ke−λτ),Δ6=λ(λ−h)+s(1−Ke−λτ). |
R=(λ−1w2λ−αp2),m+=(00Kw2e−λτ0),m−=(00−Kw2e−λτ0), |
Similarly, by theorem 3.1, we have following theorem.
Theorem 3.3. If (H1), (H2) and (H3) hold, we have the following results.
(1) all roots of Eq (12) have negative real parts for 0≤τ<τ0pace, and at least a pair of roots with positive real parts for τ∈(τ0pace,τ0pace+ε), for some ε>0,
(2) zero equilibrium of system (10) is asymptotically stable for 0≤τ<τ0pace, and unstable for τ∈(τ0pace,τ0pace+ε), for some ε>0,
(3) when τ=τ0pace, system (10) undergoes a Hopf bifurcation at zero equilibrium,
where
τ0pace=min{τ0pace+,τ0pace−}, |
τjpace±=1β±(arccoss−β2±−3Ks+2jπ),j=0,1,2,…, |
β±=√2s−h2±√(h2−2s)2−4s2(1−9K2)2. |
Remark 3.3. Near the critical value τ=τ0pace, the periodic solution of system (10) at the origin accords with pacing gait.
In a bound, legs on the same side (front or hind) as the current leg stimulate the current leg, and the other two legs inhibit the current leg, thus
K12=K,K13=−K,K14=−K,K21=K,K23=−K,K24=−K, |
K31=−K,K32=−K,K34=K,K41=−K,K42=−K,K43=K. |
Eq (1) becomes the following bounding VDP-CPG system.
{˙xi=yi,˙yi=α(p2−(xi(t)+Kxi+1(t−τ)+(−K)xi+2(t−τ)+(−K)xi+3(t−τ))2)yi−w2(xi(t)+Kxi+1(t−τ)+(−K)xi+2(t−τ)+(−K)xi+3(t−τ)),i=1,3(mod4) | (13) |
{˙xi=yi,˙yi=α(p2−(xi(t)+Kxi−1(t−τ)+(−K)xi+1(t−τ)+(−K)xi+2(t−τ))2)yi−w2(xi(t)+Kxi−1(t−τ)+(−K)xi+1(t−τ)+(−K)xi+2(t−τ)),i=2,4(mod4) |
and the linearization of Eq (13) at the origin is
{˙xi=yi,˙yi=αp2yi−w2(xi(t)+Kxi+1(t−τ)+(−K)xi+2(t−τ)+(−K)xi+3(t−τ)),i=1,3(mod4) | (14) |
{˙xi=yi,˙yi=αp2yi−w2(xi(t)+Kxi−1(t−τ)+(−K)xi+1(t−τ)+(−K)xi+2(t−τ)),i=2,4(mod4) |
the characteristic equation of system (14) is
|Rm+m−m−m+Rm−m−m−m−Rm+m−m−m+R|=Δ7(Δ8)3=0, | (15) |
where
Δ7=λ(λ−h)+s(1+3Ke−λτ),Δ8=λ(λ−h)+s(1−Ke−λτ). |
Similarly, by theorem 3.1, we have following theorem.
Theorem 3.4. If (H1), (H2) and (H3) hold, we have the following results.
(1) all roots of Eq (15) have negative real parts for 0≤τ<τ0bound, and at least a pair of roots with positive real parts for τ∈(τ0bound,τ0bound+ε), for some ε>0,
(2) zero equilibrium of system (13) is asymptotically stable for 0≤τ<τ0bound, and unstable for τ∈(τ0bound,τ0bound+ε), for some ε>0,
(3) when τ=τ0bound, system (13) undergoes a Hopf bifurcation at zero equilibrium,
where
τ0bound=min{τ0bound+,τ0bound−}, |
τjbound±=1β±(arccoss−β2±−3Ks+2jπ),j=0,1,2,…, |
β±=√2s−h2±√(h2−2s)2−4s2(1−9K2)2. |
Remark 3.4. Near the critical value τ=τ0bound , the periodic solution of system (13) at the origin accords with bounding gait.
In this section, the numerical simulation of model is carried out to verify the results obtained in the previous sections. Let α=−1.5, p=1,w=4, K=0.3, according to the calculation, we obtain the h=−1.5,s=16,m=0.1357,K2=0.09,19m=0.0151. Thus 2s−h2=29.7500>0,19m<K2<min{m,19} and the critical value τ0walk=0.7039, τ0trot=τ0pace=τ0bound=0.1103 are obtained. Basing on Theorem 3.2, we know the zero equilibrium is asymptotically stable when τ<τ0trot (shown in Figure 2a), when τ>τ0trot, the zero equilibrium of system (8) is unstable, and the periodic solution corresponding to the trot gait occurs (see Figure 2b). From theorem 3.3, we know the zero equilibrium is asymptotically stable when τ<τ0pace (shown in Figure 3a), when τ>τ0pace, the zero equilibrium of system (10) is unstable, and the periodic solution corresponding to the pace gait occurs (see Figure 3b). From theorem 3.4, we know the zero equilibrium is asymptotically stable when τ<τ0bound (shown in Figure 4a), when τ>τ0bound, the zero equilibrium of system (13) is unstable, and the periodic solution corresponding to the bound gait occurs (see Figure 4b).
In this paper, a kind of CPG network system is constructed by using VDP oscillators, and a VDP-CPG network system with four basic gaits (walk, trot, pace and bound) is presented. By studying the corresponding characteristic equations of four gaits systems, it is found that the conditions for the periodic solutions of four gaits systems are h<0,2s−h2>0 and 19m<K2<min{m,19} and the critical values τjwalk,τjtrot,τjpace andτjbound,j=0,1,2⋯. Thus, the range of coupling strength between legs in four gaits is 19m<K2<min{m,19}. Finally, the numerical simulations show that the gait systems (trot, pace and bound) produce corresponding gaits near the corresponding critical value.
This research is supported by the Fundamental Research Funds for the Central Universities (No.2572019BC12). The authors wish to express their gratitude to the editors and the reviewers for the helpful comments.
The authors declare there is no conflict of interest
[1] |
Lee JB, Park KK, Eum HM, et al. (2006) Desalination of a thermal power plant wastewater by membrane capacitive deionization. Desalination 196: 125-134. doi: 10.1016/j.desal.2006.01.011
![]() |
[2] |
Ahmaruzzaman M (2010) A review on the utilization of fly ash. Prog Energ Combust Sci 36: 327-363. doi: 10.1016/j.pecs.2009.11.003
![]() |
[3] |
Gitari WM, Petrik LF, Etchebers O, et al. (2008) Utilization of fly ash for treatment of coal mines wastewater: Solubility controls on major inorganic contaminants. Fuel 87: 2450-2462. doi: 10.1016/j.fuel.2008.03.018
![]() |
[4] |
Liu WJ, Hao Y, Jiang JH, et al. (2016) Production of a bioflocculant from Pseudomonas veronii L918 using the hydrolyzate of peanut hull and its application in the treatment of ash-flushing wastewater generated from coal fired power plant. Bioresour Technol 218: 318-325. doi: 10.1016/j.biortech.2016.06.108
![]() |
[5] |
Rubio B, Izquierdo MT, Mayoral MC, et al. (2007) Unburnt carbon from coal fly ashes as a precursor of activated carbon for nitric oxide removal. J Hazard Mater 143: 561-566. doi: 10.1016/j.jhazmat.2006.09.074
![]() |
[6] |
Dutta BK, Khanra S, Mallick D (2009) Leaching of elements from coal fly ash: Assessment of its potential for use in filling abandoned coal mines. Fuel 88: 1314-1323. doi: 10.1016/j.fuel.2009.01.005
![]() |
[7] |
Karwowska E, Wojtkowska M, Andrzejewska D (2015) The influence of metal speciation in combustion waste on the efficiency of Cu, Pb, Zn, Cd, Ni and Cr bioleaching in a mixed culture of sulfur-oxidizing and biosurfactant-producing bacteria. J Hazard Mater 299: 35-41. doi: 10.1016/j.jhazmat.2015.06.006
![]() |
[8] |
Hu Y, Zhang P, Li J, et al. (2015) Stabilization and separation of heavy metals in incineration fly ash during the hydrothermal treatment process. J Hazard Mater 299: 149-157. doi: 10.1016/j.jhazmat.2015.06.002
![]() |
[9] |
Yao ZT, Xia MS, Sarker PK, et al. (2014) A review of the alumina recovery from coal fly ash, with a focus in China. Fuel 120: 74-85. doi: 10.1016/j.fuel.2013.12.003
![]() |
[10] |
Guarienti M, Gianoncelli A, Bontempi E, et al. (2014) Biosafe inertization of municipal solid waste incinerator residues by COSMOS technology. J Hazard Mater 279: 311-321. doi: 10.1016/j.jhazmat.2014.07.017
![]() |
[11] |
Wang FH, Zhang F, Chen YJ, et al. (2015) A comparative study on the heavy metal solidification/stabilization performance of four chemical solidifying agents in municipal solid waste incineration fly ash. J Hazard Mater 300: 451-458. doi: 10.1016/j.jhazmat.2015.07.037
![]() |
[12] | Campbell A (2002) The potential role of aluminium in Alzheimer's disease. Nephrol Dial Transplant 17 (suppl 2): 17-20. |
[13] | Dearfield KL, Abernathy CO, Ottley MS, et al. (1988) Acrylamide: its metabolism, developmental and reproductive effects, genotoxicity, and carcinogenicity. Mutat Res 195(1): 45-77. |
[14] |
Salehizadeh H, Shojaosadati S (2001) Extracellular biopolymeric flocculants: recent trends and biotechnological importance. Biotechnol Adv 19: 371-385. doi: 10.1016/S0734-9750(01)00071-4
![]() |
[15] |
Liu WJ, Yuan HL, Yang JS, et al. (2009) Characterization of bioflocculants from biologically aerated filter backwashed sludge and its application in dying wastewater treatment. Bioresour Technol 100: 2629-2632. doi: 10.1016/j.biortech.2008.12.017
![]() |
[16] |
Yin YJ, Tian ZM, Tang W, et al. (2014) Production and characterization of high efficiency bioflocculant isolated from Klebsiella sp. ZZ-3. Bioresour Technol 171: 336-342. doi: 10.1016/j.biortech.2014.08.094
![]() |
[17] |
Liu WJ, He RN, Liu C (2016) An alkali-tolerant strain Microbacterium esteraromaticum C26 produces a high yield of cation-independent bioflocculant. AIMS Environ Sci 3: 408-419. doi: 10.3934/environsci.2016.3.408
![]() |
[18] | Sun P.F, Lin H, Wang G, et al. (2014) Preparation of a new-style composite containing a key bioflocculant produced by Pseudomonas aeruginosa ZJU1 and its flocculating effect on harmful algal blooms. J Hazard Mater 284: 215-221. |
[19] |
Belbahloul M, Zouhri A, Anouar A (2015) Bioflocculants extraction from Cactaceae and their application in treatment of water and wastewater. J Water Process Eng 7: 306-313. doi: 10.1016/j.jwpe.2015.07.002
![]() |
[20] |
Salehizadeh H, Yan N (2014) Recent advances in extracellular biopolymer flocculants. Biotechnol Adv 32: 1506-1522. doi: 10.1016/j.biotechadv.2014.10.004
![]() |
[21] |
Guo J, Lau AK, Zhang Y, et al. (2015) Characterization and flocculation mechanism of a bioflocculant from potato starch wastewater. Appl Microbiol Biotechnol 99: 5855-5861. doi: 10.1007/s00253-015-6589-y
![]() |
[22] |
Aljuboori AHR, Uemura Y, Osman NB, et al. (2014) Production of a bioflocculant from Aspergillus niger using palm oil mill effluent as carbon source. Bioresour Technol 171: 66-70. doi: 10.1016/j.biortech.2014.08.038
![]() |
[23] | Pu SY, Qin LL, Che JP, et al. (2014) Preparation and application of a novel bioflocculant by two strains of Rhizopus sp. using potato starch wastewater as nutrilite. Bioresour Technol 162: 184-191. |
[24] |
Zhong C, Xu A, Chen L, et al. (2014) Production of a bioflocculant from chromotropic acid waste water and its application in steroid estrogen removal. Colloids Surf B Biointerfaces 122: 729-737. doi: 10.1016/j.colsurfb.2014.08.006
![]() |
[25] |
Chen H, Zhong C, Berkhouse H, et al. (2016) Removal of cadmium by bioflocculant produced by Stenotrophomonas maltophilia using phenol-containing wastewater. Chemosphere 155: 163-169. doi: 10.1016/j.chemosphere.2016.04.044
![]() |
[26] |
Gang C, Zhang Y, Li C, et al. (2015) Production of a bioflocculant from methanol wastewater and its application in arsenite removal. Chemosphere 141: 274-281. doi: 10.1016/j.chemosphere.2015.08.009
![]() |
[27] |
Guo J, Yang C, Zeng G (2013) Treatment of swine wastewater using chemically modified zeolite and bioflocculant from activated sludge. Bioresour Technol 143: 289-297. doi: 10.1016/j.biortech.2013.06.003
![]() |
[28] |
Guo J, Jing M (2015) Bioflocculant from pre-treated sludge and its applications in sludge dewatering and swine wastewater pretreatment. Bioresour Technol 196: 736-740. doi: 10.1016/j.biortech.2015.07.113
![]() |
[29] |
Guo J, Yu J, Xin X, et al. (2015) Characterization and flocculation mechanism of a bioflocculant from hydrolyzate of rice stover. Bioresour Technol 177: 393-397. doi: 10.1016/j.biortech.2014.11.066
![]() |
[30] |
Liu WJ, Zhao CC, Jiang JH, et al. (2015) Bioflocculant production from untreated corn stover using Cellulosimicrobium cellulans L804 isolate and its application to harvesting microalgae. Biotechnol Biofuels 8: 1. doi: 10.1186/s13068-014-0179-6
![]() |
[31] |
Lu H, Gui Y, Zheng L, et al. (2013) Morphological, crystalline, thermal and physicochemical properties of cellulose nanocrystals obtained from sweet potato residue. Food Res Int 50: 121-128. doi: 10.1016/j.foodres.2012.10.013
![]() |
[32] |
Wu D (2016) Recycle Technology for Waste Residue in Potato Starch Processing: A Review. Procedia Environ Sci 31: 108-112. doi: 10.1016/j.proenv.2016.02.015
![]() |
[33] |
Liu C, Wang K, Jiang JH, et al. (2015) A novel bioflocculant produced by a salt-tolerant, alkaliphilic and biofilm-forming strain Bacillus agaradhaerens C9 and its application in harvesting Chlorella minutissima UTEX2341. Biochem Eng J 93: 166-172. doi: 10.1016/j.bej.2014.10.006
![]() |
[34] |
Krupa AND, Abigail MEA, Santhosh C, et al. (2016) Optimization of process parameters for the microbial synthesis of silver nanoparticles using 3-level Box–Behnken Design. Ecol Eng 87: 168-174. doi: 10.1016/j.ecoleng.2015.11.030
![]() |
[35] |
Masmoudi M, Besbes S, Chaabouni M, et al. (2008) Optimization of pectin extraction from lemon by-product with acidified date juice using response surface methodology. Carbohydr Polym 74: 185-192. doi: 10.1016/j.carbpol.2008.02.003
![]() |
[36] |
Zhang C, Mu T (2011) Optimisation of pectin extraction from sweet potato (Ipomoea batatas, Convolvulaceae) residues with disodium phosphate solution by response surface method. Int J Food Sci Technol 46: 2274-2280. doi: 10.1111/j.1365-2621.2011.02746.x
![]() |
[37] |
Wang S, Fang C, Wu J, et al. (2007) Optimization of pectin extraction assisted by microwave from apple pomace using response surface methodology. J Food Eng 78: 693-700. doi: 10.1016/j.jfoodeng.2005.11.008
![]() |
[38] |
Pinheiro E, Silva I, Gonzaga L, et al. (2008) Optimization of extraction of high-ester pectin from passion fruit peel (Passiflora edulis flavicarpa) with citric acid by using response surface methodology. Bioresour Technol 99: 5561-5566. doi: 10.1016/j.biortech.2007.10.058
![]() |
[39] |
Lei X, Yao C, Shao Z, et al. (2015) Effective harvesting of the microalgae Chlorella vulgaris via flocculation-flotation with bioflocculant. Bioresour Technol 198: 922-925. doi: 10.1016/j.biortech.2015.08.095
![]() |
[40] |
Liu WJ, Liu C, Yuan HL, et al. (2015) The mechanism of kaolin clay flocculation by a cation-independent bioflocculant produced by Chryseobacterium daeguense W6. AIMS Environ Sci 2: 169-179. doi: 10.3934/environsci.2015.2.169
![]() |
[41] |
Li L, Pan G (2013) A universal method for flocculating harmful algal blooms in marine and fresh waters using modified sand. Environ Sci Technol 47: 4555-4562. doi: 10.1021/es305234d
![]() |
[42] | Yuan Y, Zhang H, Pan G (2016) Flocculation of cyanobacterial cells using coal fly ash modified chitosan. Wat Res 97: 11-18. |
1. | Liqin Liu, Xiaoxiao Liu, Chunrui Zhang, REALIZATION OF NEURAL NETWORK FOR GAIT CHARACTERIZATION OF QUADRUPED LOCOMOTION, 2022, 12, 2156-907X, 455, 10.11948/20210005 | |
2. | Mingfang Chen, Kangkang Hu, Yongxia Zhang, Fengping Qi, Motion coordination control of planar 5R parallel quadruped robot based on SCPL-CPG, 2022, 14, 1687-8140, 168781402110709, 10.1177/16878140211070910 | |
3. | Zigen Song, Jiayi Zhu, Jian Xu, Gaits generation of quadruped locomotion for the CPG controller by the delay-coupled VDP oscillators, 2023, 111, 0924-090X, 18461, 10.1007/s11071-023-08783-2 | |
4. | Yangyang Han, Guoping Liu, Zhenyu Lu, Huaizhi Zong, Junhui Zhang, Feifei Zhong, Liyu Gao, A stability locomotion-control strategy for quadruped robots with center-of-mass dynamic planning, 2023, 24, 1673-565X, 516, 10.1631/jzus.A2200310 | |
5. | Zigen Song, Fengchao Ji, Jian Xu, Is there a user-friendly building unit to replicate rhythmic patterns of CPG systems? Synchrony transition and application of the delayed bursting-HCO model, 2024, 182, 09600779, 114820, 10.1016/j.chaos.2024.114820 | |
6. | Liqin Liu, Chunrui Zhang, A neural network model for goat gait, 2024, 21, 1551-0018, 6898, 10.3934/mbe.2024302 |