
Citation: Qiangqiang Zhang, Yiyan Han, Chuandong Li, Le You. Constraint impulsive consensus of nonlinear multi-agent systems with impulsive time windows[J]. AIMS Mathematics, 2020, 5(4): 3682-3701. doi: 10.3934/math.2020238
[1] | Itxaso Filgueira-Garro, Carolina González-Ferrero, Diego Mendiola, María R. Marín-Arroyo . Effect of cultivar and drying methods on phenolic compounds and antioxidant capacity in olive (Olea europaea L.) leaves. AIMS Agriculture and Food, 2022, 7(2): 250-264. doi: 10.3934/agrfood.2022016 |
[2] | Ramadan ElGamal, Omar A. Hamed, Ahmed M. Rayan, Chuanping Liu, Sameh Kishk, Salim Al-Rejaie, Gamal ElMasry . Effect of convective and vacuum drying on some physicochemical and phytochemical characteristics of peppermint leaves. AIMS Agriculture and Food, 2025, 10(1): 17-39. doi: 10.3934/agrfood.2025002 |
[3] | Hathairat Chokthaweepanich, Chuleeporn Chumnanka, Sribud Srichaijaroonpong, Rungnapa Boonpawa . Effect of harvesting age and drying condition on andrographolide content, antioxidant capacity, and antibacterial activity in Andrographis paniculata (Burm.f.) Nees. AIMS Agriculture and Food, 2023, 8(1): 137-150. doi: 10.3934/agrfood.2023007 |
[4] | Hesti Kurniasari, Wahyudi David, Laras Cempaka, Ardiansyah . Effects of drying techniques on bioactivity of ginger (Zingiber officinale): A meta-analysis investigation. AIMS Agriculture and Food, 2022, 7(2): 197-211. doi: 10.3934/agrfood.2022013 |
[5] | María Quelal, Elena Villacrés, Karla Vizuete, Alexis Debut . Physicochemical characterization of sangorache natural colorant extracts (Amaranthus quitensis L.) prepared via spray- and freeze-drying. AIMS Agriculture and Food, 2023, 8(2): 343-358. doi: 10.3934/agrfood.2023019 |
[6] | Ebrahim Falahi, Zohre Delshadian, Hassan Ahmadvand, Samira Shokri Jokar . Head space volatile constituents and antioxidant properties of five traditional Iranian wild edible plants grown in west of Iran. AIMS Agriculture and Food, 2019, 4(4): 1034-1053. doi: 10.3934/agrfood.2019.4.1034 |
[7] | Febiani Dwi Utari, Dessy Agustina Sari, Laeli Kurniasari, Andri Cahyo Kumoro, Mohamad Djaeni, Ching-Lik Hii . The enhancement of sappanwood extract drying with foaming agent under different temperature. AIMS Agriculture and Food, 2023, 8(1): 214-235. doi: 10.3934/agrfood.2023012 |
[8] | Nur Fajriani Suaib, Didah Nur Faridah, Dede Robiatul Adawiyah, Nuri Andarwulan . Semiquantification of volatile compounds and identification of potential volatile markers and dry aroma from robusta second-crack roasted coffee processed from several post-harvest processing. AIMS Agriculture and Food, 2025, 10(1): 74-96. doi: 10.3934/agrfood.2025005 |
[9] | Esteban Largo-Avila, Carlos Hernán Suarez-Rodríguez, Jorge Latorre Montero, Madison Strong, Osorio-Arias Juan . The influence of hot-air mechanical drying on the sensory quality of specialty Colombian coffee. AIMS Agriculture and Food, 2023, 8(3): 789-803. doi: 10.3934/agrfood.2023042 |
[10] | Saima Latif, Muhammad Sohaib, Sanaullah Iqbal, Muhammad Hassan Mushtaq, Muhammad Tauseef Sultan . Comparative evaluation of nutritional composition, phytochemicals and sensorial attributes of lyophilized vs conventionally dried Grewia asiatica fruit pulp powder. AIMS Agriculture and Food, 2025, 10(1): 247-265. doi: 10.3934/agrfood.2025013 |
In the paper, we consider fractional Klein-Gordon-Schrödinger (KGS) equations with the fractional Laplacian
i∂tu−ν(−Δ)α2u+uf1(|u|2,ϕ)=0, | (1.1) |
∂ttϕ−Δϕ+ϕ−f2(|u|2,ϕ)=0, | (1.2) |
where, f1(|u|2,ϕ)=∂f(|u|2,ϕ)∂|u|2, f2(|u|2,ϕ)=∂f(|u|2,ϕ)∂ϕ,f:R+×R⟼R, u(⋅,t), ϕ(⋅,t) are complex and real valued unknown functions, and the fractional Laplace operator [1,2,3,4] is defined by
(−Δ)α2u(x)=Cn,α∫Rnu(x)−u(y)|x−y|n+α/2dy. |
Let ϕt=v. Then, the fractional KGS system (1.1)–(1.2) can be expressed as the following form
iut−ν(−Δ)α2u=−uf1(|u|2,ϕ), | (1.3) |
ϕt=v, | (1.4) |
vt−Δϕ+ϕ=f2(|u|2,ϕ). | (1.5) |
The fractional KGS system (1.3)–(1.5) has the following mass and energy conserved laws
Mass: M(t)=∫+∞−∞|u(x,t)|2dx=M(0),Energy: E(t)=∫+∞−∞12(v2+|∇ϕ|2+ϕ2)+ν|(−Δ)α4u|2−f(|u|2,ϕ)dx=E(0). |
When α=2, the fractional KGS system (1.1) and (1.2) reduces to the integer KGS system [5,6,7,8,9,10,11,12]. When α≠2, the fractional KGS system [13,14,15,16,17] describes various important physical phenomena. Recently, some attentions have been paid to fractional KGS system with fractional Laplacian, which are fractional version of classical KGS equations and consider long-range interactions. In general, the analytic solution of the fractional KGS system (1.1), (1.2) can not be derived due to nonlocality and nonlinearity. Therefore, numerical method plays an important role in the study of the fractional KGS system (1.1), (1.2). Various numerical methods [13,14,15,16,17] such as the finite difference method, pseudo-spectral method, and symplectic method have been studied the fractional KGS system with f(|u|2,ϕ)=|u|2ϕ, and the stability and convergence of the numerical methods have been discussed.
It is not difficult to find that there are some deficiencies in the numerical schemes of the fractional KGS system. The first deficiencie is that the above finite difference methods are based on fractional center difference scheme, and approximate solutions converge to the exact solution at the rate O(τ2+h2). To the best of the authors'knowledge, there are few research on higher-order scheme for the fractional KGS system. The second deficiencie that the above numerical schemes are only based on one dimension fractional KGS system, and high dimension fractional KGS systems are rarely studied. In addition, there exist few reports on numerical scheme for general fractional KGS system in Eqs (1.1), (1.2) with f(|u|2,ϕ)≠|u|2ϕ.
The main goal of this paper is to construct structure-preserving scheme [18,19,20,21,22] for solving one dimension and two dimension fractional KGS equations with fractional Laplacian operator. Some numerical schemes were proposed to approximate the fractional Laplacian operator. The fractional center difference scheme based on fractional Laplace operator was first developed in [23]. Based on this work, some high-order schemes were constructed, and the results have be applied to some fractional difference equations. In addition, finite element scheme and Fourier spectral scheme of fractional Laplace operator have been studied with some special boundary conditions. In this paper, we use high-order difference scheme for fractional Laplace operator with zero boundary condition. By using some useful lemmas, we prove that the scheme can preserve the mass and energy conservation laws.
All structure-preserving schemes for the nonlinear fractional KGS system (f(|u|2,ϕ)≠|u|2ϕ) face how to efficiently solve a large nonlinear system at each time step. In this paper, we consider structure-preserving scheme for general fractional KGS equations (1.3)–(1.5). The high order central differences are used for the space direction, with the Crank-Nicolson scheme applied to the time direction. However, the numerical scheme is nonlinear scheme, and it takes too much time in the numerical simulation. Recently, scalar auxiliary variable scheme was introduce for gradient flows, and the numerical scheme only requiring solving decoupled linear systems at each time step. As far as we know, there exist few reports on scalar auxiliary variable scheme for fractional KGS system (1.3)–(1.5). In this paper, we show a new scalar auxiliary variable scheme to solve two dimension fractional KGS system (1.3)–(1.5).
In this paper, we first show numerical scheme to solve one dimension fractional KGS system (1.3)–(1.5), and discussed conservation, existence and uniqueness, stability and convergence of the numerical scheme. Second, we show numerical scheme for two dimension fractional KGS system by high order finite difference scheme in span and Crank-Nicolson scheme in time, and discussed conservation of the numerical scheme. Third, we show scalar auxiliary variable scheme to solve two dimension fractional KGS equations. Finally, some numerical examples are given, confirm theoretical results and demonstrate the efficiency of the numerical schemes.
The outline of the paper is as follows. In Section 2, the structure-preserving scheme is proposed for one dimension fractional KGS system, and convergence and stability of the numerical scheme is proved. In Section 3, the high conservative difference scheme is proposed for two dimension fractional KGS system. In Section 4, the numerical experiments are given, and the results verify the efficiency of the conservative difference scheme. Finally, a conclusion and some discussions are given in Section 5.
In the section, we show structure-preserving scheme for one dimension fractional KGS system with f(|u|2,ϕ)≠|u|2ϕ, and prove that the scheme can preserve mass and energy conservation laws. Moreover, we show that the arising scheme is uniquely solvable, and approximate solutions converge to the exact solution at the rate O(τ2+h4).
The fractional KGS system (1.1)–(1.2) contains a fractional Schrödinger equation and a classic Klein-Gordon equation, and we consider boundary condition
u(x,t)=0,x∈R/Ω;ϕ(x,t)=0,x∈∂Ω; Ω=(a,b). |
Let M,N. Then, choose time-step τ=T/N and mesh size h=(b−a)/M. Denote
xj=a+jh, tn=nτ, j=0, 1, 2,⋯, M, n=0, 1, 2, ⋯ |
Then
unj=u(xj,tn), Unj≈u(xj,tn), ϕnj=ϕ(xj,tn), Φnj≈ϕ(xj,tn). |
Define
Ωh={xj|1≤j≤M−1}, Ωτ={tn|1≤n≤N−1},ˉΩh={xj|0≤j≤M}, ˉΩτ={tn|0≤n≤N}. |
Suppose w={wnj; j=0, 1, 2⋯ M, n=0, 1, 2⋯ N} be a discrete function in ˉΩh×ˉΩτ, and
Z0h={w=wj|w0=wM=0, j=0, 1, 2, ⋯ M}. |
For convenience, we define the finite difference operators
(wnj)t=wn+1j−wnjτ, (wnj)ˆt=wn+1j−wn−1j2τ, wn+12j=wn+1j+wnj2, ˜wn+12j=wn+1j+wn−1j2. |
For any grid function u={uj}, v={vj}, define the discrete inner product, L2-norm and Lp-norm as
⟨u,v⟩=hM−1∑j=1uj¯vj, ‖u‖2=⟨u,u⟩,‖u‖plph=hM−1∑j=1|uj|p, 1≤p<+∞, ‖u‖l∞h=sup0<j<M−1|uj|. |
For 0≤δ≤1, we also define the fractional Sobolev norm ‖u‖Hδ and semi-norm |u|Hδ as
‖u‖2Hδ=∫π/h−π/h(1+|k|2δ)|ˆu(k)|2dk, |u|2Hδ=∫π/h−π/h|k|2δ|ˆu(k)|2dk, |
where
ˆu(k)=1√2πh∑j∈zuje−ikxj,uj=1√2π∫π/h−π/hˆu(k)eikxjdk. |
Lemma 1. [24] Let f(x)∈C6[xL,xR], 2<j<M−2. Then
f″(xj)=43f(xj+h)−2f(xj)+f(xj−h)h2−13f(xj+2h)−2f(xj)+f(xj−2h)h2+O(h4). |
When j=1,M−1, then
f″(x1)=76f(x2)−2f(x1)+f(x0)h2−16f(x3)−2f(x2)+f(x1)h2−112f″(x0)−1144f(4)(x0)+O(h4),f″(xM−1)=76f(xM)−2f(xM−1)+f(xM−2)h2−16f(xM−1)−2f(xM−2)+f(xM−3)h2−112f″(xM)−1144f(4)(xM)+O(h4). |
Lemma 2. [20] Suppose that u∈L1(R) and
u∈L4+α(R):={u|∫+∞−∞(1+|ξ|)4+α|ˆu(ξ)|dξ<∞}. |
Then for a fixed h, we can obtain high order scheme
43Δαhu(x)−13Δα2hu(x)=−(−Δ)α2u(x)+O(h4), |
where
Δαhu(x)=h−α∞∑k=−∞g(α)ku(x−kh), g(α)k=(−1)kΓ(α+1)Γ(α2−k+1)Γ(α2+k+1). |
Lemma 3. For any grid functions un∈Z0h, we can obtain
Im⟨43Δαhun+12−13Δα2hun+12,un+12⟩=0,Re⟨43Δαhun+12−13Δα2hun+12,unt⟩=12τ(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2−43‖ΔαhUn‖2+13‖Δα2hUn‖2), |
where Im(s),Re(s) mean taking the imaginary part and the real part of a complex number s, respectively.
Lemma 4. [25] (Diserete Sobolev inequality) For every 12≤δ≤1, there exist a constant C=C(δ)>0 independent of h>0 such that
‖u‖l∞h≤Cσ‖u‖Hδh, |
for all u∈l2h.
Lemma 5. [25] (Gagliardo-Nirenberg inequality) For any 14<δ0≤1, there exist a constant Cδ0=C(δ0)>0 independent of h>0 such that
‖u‖l4h≤Cδ0‖u‖δ0/δHδ‖u‖1−δ0/δ. |
Applying Crank-Nicolson scheme in time and higher order difference scheme in space to the fractional KGS system (1.3)–(1.5), we can obtain the following numerical scheme
i(Unj)t−νˆΔαhUn+12j=−12[f(|Un+1j|2,Φnj)−f(|Unj|2,Φnj)|Un+1j|2−|Unj|2+f(|Un+1j|2,Φn+1j)−f(|Unj|2,Φn+1j)|Un+1j|2−|Unj|2]Un+12j, | (2.1) |
(Φnj)t=Vn+12j, | (2.2) |
(Vnj)t−˜Δ2hΦn+12j+Φn+12j=12[f(|Unj|2,Φn+1j)−f(|Unj|2,Φnj)Φn+1j−Φnj+f(|Un+1j|2,Φn+1j)−f(|Un+1j|2,Φnj)Φn+1j−Φnj], | (2.3) |
U0j=u0(xj),Φ0j=ϕ0(xj),V0j=ϕ1(xj), | (2.4) |
Un0=UnM,Φn0=ΦnM, | (2.5) |
where
ˆΔαhUn+12j=43ΔαhUn+12j−13Δα2hUn+12j, |
˜Δ2hΦn+12j={76Δ2hΦn+121−112Δ2hΦn+122,j=1,43Δ2hΦn+12j−13Δ22hΦn+12j,2<j<M−2,76Δ2hΦn+12M−1−112Δ2hΦn+12M−2,j=M−1. |
Theorem 1. The scheme (2.1)–(2.5) is conservative in the sense
Mass: Mn=Mn−1⋯=M0,Energy: En=En−1⋯=E0, |
where
Mn=‖Un‖2,En=ν(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2)+12[43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+‖Φn+1‖2+‖Vn+1‖2]−hM−1∑j=1|f(|Un+1j|2,Φn+1j). |
Proof: Multiplying h(Un+1j+Unj)∗ to Eq (2.1) and summing it up for 1≤j≤M−1 can yield
i⟨Unt,2Un+12⟩−ν⟨ˆΔαhUn+12,2Un+12⟩=⟨−12[f(|Un+1|2,Φn)−f(|Un|2,Φn)|Un+1|2−|Un|2+f(|Un+1|2,Φn+1)−f(|Un|2,Φn+1)|Un+1|2−|Un|2]Un+12,2Un+12⟩. | (2.6) |
Taking the imaginary part of Eq (2.6) and noting that Lemma 3 yields ‖Un+1‖2=‖Un‖2. When n=0,1,2,⋯, we can obtain Mn=Mn−1⋯=M0.
Multiplying 2h(Un+1j−Unj)∗/τ to Eq (2.1) and summing it up for 1≤j≤M−1 can yield
i⟨Unt,2Unt⟩−ν⟨ˆΔαhUn+12,2Unt⟩=⟨−12[f(|Un+1|2,Φn)−f(|Un|2,Φn)|Un+1|2−|Un|2+f(|Un+1|2,Φn+1)−f(|Un|2,Φn+1)|Un+1|2−|Un|2]Un+12,2Unt⟩. | (2.7) |
It follows from Lemma 3 that
⟨ˆΔαhUn+12,2Unt⟩=1τ(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2−43‖ΔαhUn‖2+13‖Δα2hUn‖2). |
Taking the real part of Eq (2.7) and noting that Lemma 3 yields
ντ(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2−43‖ΔαhUn‖2+13‖Δα2hUn‖2)=h2τM−1∑j=1[f(|Un+1j|2,Φn+1j)−f(|Unj|2,Φn+1j)+f(|Un+1j|2,Φnj)−f(|Unj|2,Φnj)]. | (2.8) |
Multiplying 2h(Vn+1j−Vnj)/τ to Eq (2.2) and summing it up for 1≤j≤M−1 can yield
⟨Φnt,2Vnt⟩=⟨Vn+12,2Vnt⟩. | (2.9) |
Multiplying 2h(Φn+1j−Φnj)/τ to Eq (2.3) and summing it up for 1≤j≤M−1 can yield
⟨Vnt,2Φnt⟩−⟨˜Δ2hΦn+12,2Φnt⟩+⟨Φn+12,2Φnt⟩=⟨12[f(|Un|2,Φn+1)−f(|Un|2,Φn)Φn+1−Φn+f(|Un+1|2,Φn+1)−f(|Un+1|2,Φn)Φn+1−Φn],2Φnt⟩. | (2.10) |
It follows from literature [24] that
⟨˜Δ2hΦn+12,2Φnt⟩=−1τ[(43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2−43‖Δ2hΦn‖2+13‖Δ22hΦn‖2)+h6(|Δ2hΦn+112|2−|Δ2hΦn12|2)+(|Δ2hΦn+1M−12|2−|Δ2hΦnM−12|2)]. |
Noting that equation ⟨Φnt,Vnt⟩=⟨Vnt,Φnt⟩, we obtain
1τ[(43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2−43‖Δ2hΦn‖2+13‖Δ22hΦn‖2)+h6(|Δ2hΦn+112|2−|Δ2hΦn12|2)+(|Δ2hΦn+1M−12|2−|Δ2hΦnM−12|2)]+1τ(‖Φn+1‖2−‖Φn‖2)+1τ(‖Vn+1‖2−‖Vn‖2)=hτM−1∑j=1[f(|Un+1j|2,Φn+1j)−f(|Un+1j|2,Φnj)+f(|Unj|2,Φn+1j)−f(|Unj|2,Φnj)]. | (2.11) |
Combining Eqs (2.8) and (2.11), we obtain
ν(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2)+12[43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+‖Φn+1‖2+‖Vn+1‖2]−hM−1∑j=1|f(|Un+1j|2,Φn+1j)=ν(43‖ΔαhUn‖2−13‖Δα2hUn‖2)+12[43‖Δ2hΦn‖2−13‖Δ22hΦn‖2+h6(|Δ2hΦn12|2+|Δ2hΦnM−12|2)+‖Φn+1‖2+‖Vn‖2]−hM−1∑j=1f(|Unj|2,Φnj). |
When n=0,1,2,⋯, we can obtain En=En−1⋯=E0.
In [15], we show the theoretical analysis for f(|u|2,ϕ)=|u|2ϕ+|u|4ϕ. In similar analysis of literature [15], we can obtain the theoretical analysis for f(|u|2,ϕ)=|u|2l⋅ϕm when l,m satisfy certain conditions. For simplicity of notation, we only consider f(|u|2,ϕ)=|u|2ϕ in next theoretical analysis.
Theorem 2. The scheme (2.1)–(2.5) is bounded in the discrete l∞h.
Proof: Using Theorem 1 can obtain ‖Un‖=˜C and
En=ν(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2)+12[43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+‖Φn+1‖2+‖Vn+1‖2]−hM−1∑j=1|Un+1j|2Φn+1j. | (2.12) |
Noting that Young's inequality ab≤14a2+b2, we can obtain
hM−1∑j=1|Un+1j|2Φn+1j≤hM−1∑j=1|Un+1j|2|Φn+1j|≤hM−1∑j=1(|Un+1j|4+14|Φn+1j|2)=‖Un+1‖4l4h+14‖Φn+1‖2. | (2.13) |
Using Gagliardo-Nirenberg inequality and noting that 14<σ0<α4 can yield
‖Un+1‖4l4h≤ˆCδ0‖Un+1‖8δ0αHα2‖Un+1‖4−8δ0α≤ˆCδ0(ε‖Un+1‖2Hα2+C(ε)). | (2.14) |
It follows from [22] that there exists a constant 1≤Cα≤(2π)α such that
43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2=Cα|Un+1|2Hα2, | (2.15) |
43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2=C|Φn+1|2H1. | (2.16) |
Substituting Eqs (2.15), (2.16) into Eq (2.12) can yield
Cαν|Un+1|2Hα2+12[C|Φn+1|2H1+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+‖Φn+1‖2+‖Vn+1‖2]≤En+ˆCδ0(ε‖Un+1‖2Hα2+C(ε))+14‖Φn+1‖2. |
Noting that ‖Un‖2Hα2=|Un|2Hα2+‖Un‖2,En=E0, we have
(νCα−ˆCδ0ε)|Un+1|2Hα2+12[C|Φn+1|2H1+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+12‖Φn+1‖2+‖Vn+1‖2]≤ˆCδ0C(ε)+E0. |
When ε<νCαˆCδ0, there exist a constant C, such that
|Φn+1|2H1≤C,‖Φn‖2≤C,‖Un‖2Hα2≤C. |
It follows from Lemma 4 that ‖Un‖l∞h≤C,‖Φn‖l∞h≤C.
Theorem 3. The solution of the numerical scheme (2.1)–(2.5) exists and unique.
Proof: Let Pn=Un−ˆUn,Rn=Vn−ˆVn,Qn=Φn−ˆΦn. Then, it follows from Eqs (2.1)–(2.15) that
i(Pnj)t−νˆΔαhPn+12j=−(Un+12jΦn+12j−ˆUn+12jˆΦn+12j),1<j<M−1,1<n<N, | (2.17) |
(Qnj)t=Rn+12j,1<j<M−1,1<n<N, | (2.18) |
(Rnj)t−˜Δ2hQj+Qn+12j=12(|Unj|2+|Un+1j|2)−12(|ˆUnj|2+|ˆUn+1j|2),1<j<M−1,1<n<N. | (2.19) |
Multiplying h(Pn+1j+Pnj)∗ to Eq (2.17) and summing it up for 1<j<M−1 can yield
i⟨Pnt−νˆΔαhPn+12,Pn+12⟩=−⟨Un+12Φn+12−ˆUn+12ˆΦn+12,Pn+12⟩. | (2.20) |
Taking the imaginary part of Eq (2.20) can obtain
12τ(‖Pn+1‖2−‖Pn‖2)=Im⟨−(Un+12Φn+12−ˆUn+12ˆΦn+12),Pn+12⟩. | (2.21) |
Noting that Theorem 2, then we can obtain
12τ(‖Pn+1‖2−‖Pn‖2)≤C(‖Pn+1‖2+‖Pn‖2+‖Qn+1‖2+‖Qn‖2). | (2.22) |
Multiplying (Rn+1j−Rnj)/τ to Eq (2.18) and summing it up for 1<j<M−1 can yield
⟨Qnt,Rnt⟩=⟨Rn+12,Rnt⟩. | (2.23) |
Multiplying (Qn+1j−Qnj)/τ to Eq (2.19) and summing it up for 1<j<M−1 can yield
⟨Rnt,Qnt⟩+[(43‖Δ2hQn+1‖2−13‖Δ22hQn+1‖2−43‖Δ2hQn‖2+13‖Δ22hQn‖2)+h6(|Δ2hQn+112|2−|Δ2hQn12|2)+h6(|Δ2hQn+1M−12|2−|Δ2hQnM−12|2)]+⟨Qn+1,Qnt⟩=⟨12(|Un|2+|Un+1|2)−12(|ˆUn|2+|ˆUn+1|2),Qnt⟩. | (2.24) |
It follows from Theorem 2 that the Eq (2.24) can be expressed as the following form
⟨Rnt,Qnt⟩+1τ[(43‖Δ2hQn+1‖2−13‖Δ22hQn+1‖2−43‖Δ2hQn‖2+13‖Δ22hQn‖2)+h6(|Δ2hQn+112|2−|Δ2hQn12|2)+h6(|Δ2hQn+1M−12|2−|Δ2hQnM−12|2)]+1τ(‖Qn+1‖−‖Qn‖)≤C2(‖Pn‖2+‖Rn+1‖2+‖Rn‖). |
Combining Eq (2.23) can yield
1τ(‖Rn+1‖2−‖Rn‖)+1τ[(43‖Δ2hQn+1‖2−13‖Δ22hQn+1‖2−43‖Δ2hQn‖2+13‖Δ22hQn‖2)+h6(|Δ2hQn+112|2−|Δ2hQn12|2)+h6(|Δ2hQn+1M−12|2−|Δ2hQnM−12|2)]+1τ(‖Qn+1‖−‖Qn‖)≤C2(‖Pn‖2+‖Rn+1‖2+‖Rn‖). | (2.25) |
Adding Eqs (2.22) and (2.25) can yield
1τ(‖Rn+1‖2−‖Rn‖)+1τ(‖Pn+1‖2−‖Pn‖)+1τ(‖Qn+1‖−‖Qn‖)+1τ[(43‖Δ2hQn+1‖2−13‖Δ22hQn+1‖2−43‖Δ2hQn‖2+13‖Δ22hQn‖2)+h6(|Δ2hQn+112|2−|Δ2hQn12|2)+h6(|Δ2hQn+1M−12|2−|Δ2hQnM−12|2)]≤C2(‖Pn‖2+‖Pn+1‖2+‖Qn‖2+‖Qn+1‖2+‖Rn+1‖2+‖Rn‖). |
Let
Bn=‖Rn+1‖2+‖Pn+1‖2+‖Qn+1‖+43‖ΔαhQn+1‖2−13‖Δα2hQn+1‖2+h6|Δ2hQn+112|2−h6|Δ2hQn+1M−12|2. |
Then, we can obtain
Bn−Bn−1≤Cτ(Bn+Bn−1). |
It follows from Gronwall's inequality that
max1≤n≤N(‖Bn‖2)≤e4CT‖B0‖2=0. |
Noting that ‖Pn+1‖2+‖Qn+1‖2≤Bn, we can obtain ‖Pn+1‖2=0,‖Qn+1‖2=0.
Let unj=u(xj,tn), vnj=v(xj,tn), ϕnj=ϕ(xj,tn). Then, we define the local truncation error as
Rn1=i(unj)t−νˆΔαhun+12j+un+12jϕn+12j,1<j<M−1,1<n<N, | (2.26) |
Rn2=(ϕnj)t−vn+12j,1<j<M−1,1<n<N, | (2.27) |
Rn3=(vnj)t−˜Δ2hϕj+ϕn+12j−12(|unj|2+|un+1j|2),1<j<M−1,1<n<N. | (2.28) |
According to Taylor expansion, we obtain the following result.
Theorem 4. |Rnj|≤ˆC(τ2+h4) holds as τ,h→0.
Theorem 5. Suppose that the problem (1.1),(1.2) has a smooth solution, then the solution Un,Φn of difference scheme (2.1)–(2.5) converges to the true solution u,ϕ with order O(τ2+h4) by the‖⋅‖l∞h norm.
Proof: Let enj=Unj−unj,ηnj=Φnj−ϕnj, ξj=Vnj−vnj. Then, we obtain
Rn1=i(enj)t−νˆΔαhen+12j+Fn+12j,1<j<M−1,1<n<N, | (2.29) |
Rn2=(ηnj)t−ξn+12j,1<j<M−1,1<n<N, | (2.30) |
Rn3=(ξnj)t−˜Δ2hηj+ηn+12j−Gnj,1<j<M−1,1<n<N, | (2.31) |
where
Fn+12=Un+12Φn+12−un+12ϕn+12,Gn=12(|Un|2+|Un+1|2)−12(|un|2+|un+1|2)=12(Un¯Un+Un+1¯Un+1)−12(un¯un+un+1¯un+1)=12(en¯Un+un¯en+en+1¯Un+1+un+1¯en+1). |
Multiplying h(en+1j+enj)∗ to Eq (2.29) and summing it up for 1<j<M−1 can yield
⟨Rn1,2en+12⟩=i⟨ent,2en+12⟩−ν⟨ˆΔαhen+12,2en+12⟩+⟨Fn+12,2en+12⟩. | (2.32) |
Taking the imaginary part of Eq (2.32) can yield
1τ(‖en+1‖2−‖en‖2)≤‖Rn1‖2+C(‖en+1‖2+‖en‖2+‖ηn+1‖2+‖ηn‖2). | (2.33) |
Multiplying h(ξn+1j−ξnj)/τ to Eq (2.30) and summing it up for 1<j<M−1 can yield
⟨Rn2,ξnt⟩=⟨ηnt−ξn+12j,ξnt⟩. | (2.34) |
Multiplying h(ηn+1j−ηnj)/τ to Eq (2.31) and summing it up for 1<j<M−1 can yield
⟨Rn3,ηnt⟩=⟨ξnt,ηnt⟩+1τ[(43‖Δ2hηn+1‖2−13‖Δ22hηn+1‖2−43‖Δ2hηn‖2+13‖Δ22hηn‖2)+h6(|Δ2hηn+112|2−|Δ2hηn12|2)+(|Δ2hηn+1M−12|2−|Δ2hηnM−12|2)]+1τ(‖ηn+1‖−‖ηn‖)+⟨Gn,ηnt⟩. | (2.35) |
It follows from Cauchy-Schwarz inequality that
⟨G,ηnt⟩=⟨G,Rn2−ηn+12⟩≤C(‖ηn+1‖2+‖ηn‖2+‖en+1‖2+‖en‖2), | (2.36) |
⟨Rn3,ηnt⟩=⟨Rn3,Rn2−ηn+12⟩≤C(‖Rn3‖2+‖Rn2‖2+‖ηn+1‖2+‖ηn‖2). | (2.37) |
Substituting Eqs. (2.36) and (2.37) into Eq (2.35) and noting that Eqs (2.33) and (2.34) yields
1τ(‖en+1‖2−‖en‖2)+1τ(‖ξn+1‖2−‖ηn‖2)+1τ[(43‖Δαhηn+1‖2−13‖Δα2hηn+1‖2−43‖Δαhηn‖2+13‖Δα2hηn‖2)+h6(|Δαhηn+112|2−|Δαhηn12|2)+(|Δαhηn+1M−12|2−|ΔαhηnM−12|2)]+1τ(‖ηn+1‖2−‖ηn‖2)≤C(‖Rn1‖2+‖Rn2‖2+‖Rn3‖2+‖ηn+1‖2+‖ηn‖2+‖en+1‖2+‖en‖2). |
Let
Bn=‖en+1‖2+‖ξn+1‖2+43‖Δαhηn+1‖2−13‖Δα2hηn+1‖2+h6(|Δαhηn+112|2+|Δαhηn+1M−12|2)+‖ηn+1‖2. |
Then, we can obtain
Bn−Bn−1≤τ(|Rn1‖2+|Rn2‖2+|Rn3‖2)+Cτ(Bn+Bn−1). |
It follows from Gronwall's inequality that
max1≤n≤NBn≤(B0+τN∑l=1(|Rl1‖2+|Rl2‖2))e8CT≤(B0+ˆCT(τ2+h4)2)e8CT. |
Noting that B0=O(τ2+h4), we can obtain
‖en‖≤C(τ2+h4),‖ηn‖H1≤C(τ2+h4),‖ηn‖≤C(τ2+h4). | (2.38) |
Multiplying h(en+1j−enj)∗/τ to Eq (2.29) and summing it up for 1<j<M−1 can yield
⟨Rn1,en+1−en⟩=i⟨ent,en+1−en⟩−ν⟨ˆΔαhen+12,en+1−en⟩+⟨Fn+12,en+1−en⟩. | (2.39) |
It follows from Eq (2.32) that can obtain
Re⟨Fn+12,en+1−en⟩=τRe⟨Fn+12,−iνˆΔαhen+12+iFn+12−iRn1⟩=τIm⟨Fn+12,−νˆΔαhen+12⟩−τ⟨Fn+12,Rn1⟩, |
and
⟨ˆΔαhen+12,en+1−en⟩=ˆCατ(|en+1|2Hα/2−|en|2Hα/2). |
Noting that
|Im⟨Fn+12,−νˆΔαhen+12⟩|≤C(|Fn+12|2+|e|nHα2+|e|n+1Hα2), |
we can yield
|en+1|2Hα/2−|en|2Hα/2≤τ(|en+1|2Hα/2+|en|2Hα/2+|Fn+12|2Hα/2+Re⟨Rn1,en+1−en⟩). | (2.40) |
It follows from [15] that can obtain
|Fn+12|2Hα/2≤C(|en+1|2Hα/2+|en|2Hα/2+(τ2+h4)2), ‖Fn+12‖2≤C(τ2+h4)2. |
Thus, Eq (2.40) can be expressed as
|en+1|2Hα/2−|en|2Hα/2≤Cτ(|en+1|2Hα/2+|en|2Hα/2+(τ2+h2)2)+Re⟨Rn1,en+1−en⟩. | (2.41) |
Summing up the superscript n to N and then replacing N by n, we get
|en+1|2Hα/2≤Cn∑l=0|el|2Hα/2+CT(τ2+h4)2. |
It follows from Gronwall Inequality that
|en+1|Hα/2≤C(τ2+h4). |
Noting that Eq (2.38) can yield ‖e‖l∞h≤C(τ2+h4), ‖η‖l∞h≤C(τ2+h4).
Let U=[U1,U2,⋯,UM−1], V=[V1,V2,⋯,VM−1], Φ=[Φ1,Φ2,⋯,ΦM−1]. Then, we rewrite the numerical scheme (2.1)–(2.5) as the following vector form
iUnt−ν˜AUn+12=−12[f(|Un+1|2,Φn)−f(|Un|2,Φn)|Un+1|2−|Un|2+f(|Un+1|2,Φn+1)−f(|Un|2,Φn+1)|Un+1|2−|Un|2]Un+12,Φnt=Vn+12,Vnt−˜BΦn+12+Φn+12=12[f(|Un|2,Φn+1)−f(|Un|2,Φn)Φn+1−Φn+f(|Un+1j|2,Φn+1)−f(|Un+1|2,Φn)Φn+1−Φn],U0=u0,Φ0=ϕ0,V0=ϕ1, |
where, matrices ˜A, ˜B represent differential matrix of fractional Laplacian operator and classical Laplacian operator, respectively. In order to solve above nonlinear numerical scheme, we will use the following iterative algorithm
iUn+1(s+1)−Unτ−ν˜AUn+1(s+1)+Un2=−12[f(|Un+1(s)|2,Φn)−f(|Un|2,Φn)|Un+1(s)|2−|Un|2+f(|Un+1(s)|2,Φn+1(s))−f(|Un|2,Φn+1(s))|Un+1(s)|2−|Un|2]Un+1(s)+Un2,Φn+1(s+1)−Φnτ=Vn+1(s)+Vn2,Vn+1(s+1)−Vnτ−˜BΦn+1(s+1)+Φn2+Φn+1(s+1)+Φn2=12[f(|Un|2,Φn+1(s+1))−f(|Un|2,Φn)Φn+1(s+1)−Φn+f(|Un+1(s+1)|2,Φn+1(s+1))−f(|Un+1(s+1)|2,Φn)Φn+1(s+1)−Φn]. |
In this section, we first show Crank-Nicolson scheme in time and high central difference scheme in space, and the obtained scheme preserves mass and energy conservation laws. However, the obtained discrete system is nonlinear system, and it takes too much time in the numerical simulation for two dimension case. Then, we show a equivalent form of two dimension fractional KGS system by introducing some new auxiliary variables. The new system is discretized by the scalar auxiliary variable scheme, and a linear discrete system is obtained, which can preserve energy conservation law.
Now, we consider boundary condition
u(x,y,t)=0, (x,y)∈R2/Ω; ϕ(x,y,t)=0, (x,y)∈∂Ω, Ω=(xL,xR)×(yL,yR). |
Let
hx=(xR−xL)/M, hy=(yR−yL)/M, τ=T/N, |
where M, N be positive integers. Then,
xj=xL+jhx, yk=yL+khy, tn=nτ. |
Denote
Ωhx,hy={(xj,yk)|1≤j≤M−1, 1≤k≤M−1}, Ωτ={tn|1≤n≤N−1},ˉΩhx,hy={(xj,yk)|0≤j≤M, 0≤k≤M}, ˉΩτ={tn|0≤n≤N}. |
Then, the grid function can be defined by Unj,k≈u(xj,yk,tn), Φnj,k≈ϕ(xj,yk,tn), where unj,k=u(xj,yk,tn), ϕnj,k=ϕ(xj,yk,tn). For any grid u={uj,k}, ϕ={ϕj,k}, we can define
⟨u,v⟩=hxhyM−1∑j=1M−1∑k=1uj,k¯vj,k, ‖u‖2=⟨u,u⟩. |
Lemma 6. [26,27] Suppose that u∈L1(R2) and
u∈L4+α(R2):={u|∫+∞−∞(1+|ξ|)4+α|ˆu(ξ)|dξ<∞}. |
Then, for a fixed h=hx=hy, we can obtain high order scheme
43Δαhu(x,y)−13Δα2hu(x,y)=Δα2u(x,y)+O(h4), | (3.1) |
where
Δαhu(x,y)=h−2α∞∑k=−∞∞∑l=−∞g(α)k,lu(x−kh,y−lh), |
and g(α)k,l are Fourier expansion coefficients of generation function
ρ(x,y)=[4sin2(x2)+4sin2(y2)]α2, |
which can be calculated as
g(α)k,l=1(2π)2∫∫[−π,π]2ρ(x,y)e−i(kx+ly)dxdy. |
Lemma 7. Let
˜Δ2hx,hyΦ(xj,yk)=˜Δ2hxΦ(xj,yk)+˜Δ2hyΦ(xj,yk), |
where
˜Δ2hxΦ(xj,yk)={76Δ2hΦ(x1,yk)n+12−112Δ2hΦ(x2,yk)n+12,j=1,43Δ2hΦ(xj,yk)n+12−13Δ22hΦ(xj,yk)n+12,2<j<M−2,76Δ2hΦ(xM−1,yk)n+12−112Δ2hΦ(xM−2,yk)n+12,j=M−1, |
˜Δ2hyΦ(xj,yk)={76Δ2hΦ(xj,y1)n+12−112Δ2hΦ(xj,y2)n+12,k=1,43Δ2hΦ(xj,yk)n+12−13Δ22hΦ(xj,yk)n+12,2<k<M−2,76Δ2hΦ(xj,yM−1)n+12−112Δ2hΦ(xj,yM−2)n+12,k=M−1. |
Then, we obtain
ΔΦ(xj,yk)=˜Δ2hx,hyΦ(xj,yk)+O(h4x+h4y). |
Applying Crank-Nicolson scheme in time and higher order difference scheme in space to two dimension fractional KGS system (1.3)–(1.5), we can obtain the following numerical scheme
i(Unj,k)t−νˆΔαhUn+12j,k=−12[f(|Un+1j,k|2,Φnj,k)−f(|Unj,k|2,Φnj,k)|Un+1j,k|2−|Unj,k|2+f(|Un+1j,k|2,Φn+1j,k)−f(|Unj,k|2,Φn+1j,k)|Un+1j,k|2−|Unj,k|2]Un+12j,k, | (3.2) |
j,k=1,2,⋯,M−1,n=0,1,⋯,(Φnj,k)t=Vn+12j,k,j,k=1,2,⋯,M−1,n=0,1,⋯, | (3.3) |
(Vnj,k)t−˜Δ2hx,hyΦn+12j,k+Φn+12j,k=12[f(|Unj,k|2,Φn+1j,k)−f(|Unj,k|2,Φnj,k)Φn+1j,k−Φnj,k+f(|Un+1j,k|2,Φn+1j,k)−f(|Un+1j,k|2,Φnj,k)Φn+1j,k−Φnj,k], | (3.4) |
j,k=1,⋯,M−1,n=0,1,⋯,U0j,k=u0(xj,yk),Φ0j,k=ϕ0(xj,yk),V0j,k=ϕ1(xj,yk),j,k=1,2,⋯,M−1, | (3.5) |
Un0,k=UnM,k,Φn0,k=ΦnM,k, Unj,0=Unj,M,Φnj,M=Φnj,M, n=0,1,⋯, | (3.6) |
where
ˆΔαhUn+12j,k=43ΔαhUn+12j,k−13Δα2hUn+12j,k. |
Let
˜Δ2hx,hyΦj,k=(Dx⊗I+I⊗Dy)Φ, 43ΔαhUj,k−13Δα2hUj,k=AU. |
Then, we can obtain
iUt−νAUn+12=−12[f(|Un+1|2,Φn)−f(|Un|2,Φn)|Un+1|2−|Un|2+f(|Un+1|2,Φn+1)−f(|Un|2,Φn+1)|Un+1|2−|Un|2]Un+12,n=0,1,⋯, | (3.7) |
Φnt=Vn+12,n=0,1,⋯, | (3.8) |
Vnt−(Dx⊗I+I⊗Dy)Φn+12+Φn+12=12[f(|Un|2,Φn+1)−f(|Un|2,Φn)Φn+1−Φn+f(|Un+1|2,Φn+1)−f(|Un+1|2,Φn)Φn+1−Φn],n=0,1,⋯, | (3.9) |
where, ⊗ represents kronecker product of matrices, Dx, Dy are differential matrix of x, y direction, and I is identity matrix.
Lemma 8. [19] Let A∈Rn×n have eigenvalues {λj}nj=1, and let B∈Rm×m have eigenvalues{μj}mj=1. Then the m×n eigenvalues of A⊗B are
λ1μ1,⋯,λ1μm,λ2μ1,⋯,λ2μm,⋯,λnμ1,⋯,λnμm. |
Lemma 9. [19] For matrixs A and B, (A⊗B)T=AT⊗BT.
Let Λh=Dx⊗I+I⊗Dy. Then, it follows from Lemmas 8 and 9 that the matrix Λh and A are symmetric positive matrixs. Moreover, there exists fractional symmetric positive difference quotient operator denoted by Λ12h and A12 such that
−⟨Λhv,v⟩=⟨Λ12hv,Λ12hv⟩, −⟨Av,v⟩=⟨A12v,A12v⟩. |
It follows from Lemmas 8 and 9 that the numerical scheme (3.7)–(3.9) preserves mass and energy conservation laws.
For f(|u|2,ϕ)=|u|2ϕ, it follows from Eqs (3.2)–(3.6) that we can obtain the following numerical scheme
i(Unj,k)t−νˆΔαhUn+12j,k=−Un+12j,kΦn+12j,k,j,k=1,2,⋯,M−1,n=0,1,⋯, | (3.10) |
(Φnj,k)t=Vn+12j,k,j,k=1,2,⋯,M−1,n=0,1,⋯, | (3.11) |
(Vnj,k)t−˜Δ2hx,hyΦn+12j,k+Φn+12j,k=12(|Unj,k|2+|Un+1j,k|2),j,k=1,⋯,M−1,n=0,1,⋯, | (3.12) |
U0j,k=u0(xj,yk),Φ0j,k=ϕ0(xj,yk),V0j,k=ϕ1(xj,yk),j,k=1,2,⋯,M−1, | (3.13) |
Un0,k=UnM,k,Φn0,k=ΦnM,k, Unj,0=Unj,M,Φnj,M=Φnj,M, n=0,1,⋯. | (3.14) |
We can prove that the resulting scheme (3.10)–(3.14) can preserve the mass and energy conservation laws. However, the above numerical scheme is nonlinear scheme. In order to construct linear scheme, we consider also the following finite difference scheme for fractional KGS system (1.3)–(1.5) with f(|u|2,ϕ)=|u|2ϕ
i(Unj,k)t−νˆΔαhUn+12j,k=−Un+12j,kΦn+12j,k,j,k=1,2,⋯,M−1,n=0,1,⋯, | (3.15) |
(Φnj,k)ˆt=˜Vnj,k,j,k=1,2,⋯,M−1,n=0,1,⋯, | (3.16) |
(Vnj,k)ˆt−˜Δ2hx,hy˜Φnj,k+˜Φnj,k=|Unj,k|2,j,k=2,⋯,M−2,n=0,1,⋯, | (3.17) |
U0j,k=u0(xj,yk),Φ0j,k=ϕ0(xj,yk),V0j,k=ϕ1(xj,yk),j,k=1,2,⋯,M−1, | (3.18) |
Un0,k=UnM,k,Φn0,k=ΦnM,k, Unj,0=Unj,M,Φnj,M=Φnj,M, n=0,1,⋯. | (3.19) |
We can also prove that the resulting scheme (3.15)–(3.19) can preserve the mass and energy conservation laws. However, the scheme is only conservation for f(|u|2,ϕ)=|u|2ϕ. For f(|u|2,ϕ)≠|u|2ϕ, we use scalar auxiliary variable scheme to obtain linearly implicit scheme. Let q=√⟨f(|u|2,ϕ),1⟩+C0. Then,
iut−ν(−Δα2)u=−u∂f(|u|2,ϕ)∂|u|2q√⟨f(|u|2,ϕ),1⟩+C0, | (3.20) |
ϕt=v, | (3.21) |
vt−Δϕ+ϕ=∂f(|u|2,ϕ)∂ϕq√⟨f(|u|2,ϕ),1⟩+C0, | (3.22) |
qt=⟨∂f(|u|2,ϕ)∂|u|2,2Re(u⋅ut)⟩+⟨∂f(|u|2,ϕ)∂ϕ,∂ϕ∂t⟩2√⟨f(|u|2,ϕ),1⟩+C0. | (3.23) |
Lemma 10. The fractional KGS system (3.20)–(3.23) has the following energy conserved laws
E(t)=∫+∞−∞v2+|∇ϕ|2+ϕ2+2|(−Δ)α4u|2+q2dxdy=E(0). |
Applying scalar auxiliary variable scheme in time and higher order difference scheme in space to the fractional KGS system (1.3)–(1.5), we can obtain the following numerical scheme
i(Unj,k)t−νˆΔαhUn+12j,k=−~Uj,kn+12∂f(|~Uj,kn+12|2,~Φj,kn+12)∂|~Uj,kn+12|2Qn+12j,k√⟨f(|~Uj,kn+12|2,~Φj,kn+12),1⟩+C0, | (3.24) |
j,k=1,2,⋯,M−1,n=0,1,⋯,(Φnj,k)t=Vn+12j,k,j=1,2,⋯,M−1,n=0,1,⋯, | (3.25) |
(Vnj,k)t−˜Δ2hx,hyΦn+12j,k+Φn+12j,k=∂f(|~Uj,kn+12|2,~Φj,kn+12)∂~Φj,kn+12Qn+12j,k√⟨f(|~Uj,kn+12|2,~Φj,kn+12),1⟩+C0, | (3.26) |
j,k=1,2,⋯,M−1,n=0,1,⋯,(Qnj,k)t=⟨∂f(|~Uj,kn+12|2,~Φj,kn+12)∂|~Un+12j,kn+12|2,2Re(~Uj,kn+12⋅(Unj,k)t)⟩+⟨∂f(|~Uj,k|2,~Φijn+12)∂~Φj,kn+12,(Φnj,k)t⟩2√⟨f(|~Uj,kn+12|2,~Φj,kn+12),1⟩+C0,j,k=1,⋯,M−1,n=0,1,⋯, | (3.27) |
where ~Uj,kn+12=3Unj,k−Un−1j,k2,~Φj,kn+12=3Φnj,k−Φn−1j,k2.
Theorem 6. The scheme (3.24)–(3.27) is conservative in the sense
Energy: En=En−1⋯=E0, |
where
En=43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2+43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+‖Φn+1‖2+‖Vn+1‖2−‖Qn+1‖2. |
Proof: Multiplying 2h(Un+1j,k−Unj,k)∗/τ to Eq (3.24) and summing it up for 1≤j,k≤M−1 can yield
i⟨Unt−νˆΔαhUn+12,2Un+12t⟩=−⟨˜Un+12∂f(|˜Un+12|2,˜Φn+12)∂|˜Un+12|2Qn+12√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,2Un+12t⟩. | (3.28) |
Taking the real part of Eq (3.28) yields
1τ(43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2−43‖ΔαhUn‖2+13‖Δα2hUn‖2)=Re⟨˜Un+12∂f(|˜Un+12|2,˜Φn+12)∂|˜Un+12|2Qn+12√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,2Un+12t⟩. | (3.29) |
Multiplying 2h(Vn+1j−Vnj)/τ to Eq (3.25) and summing it up for 1≤j,k≤M−1 can yield
⟨Φnt,2Vnt⟩=⟨Vn+12,2Vnt⟩. | (3.30) |
Multiplying 2h(Φn+1j−Φnj)/τ to Eq (3.26) and summing it up for 1≤j,k≤M−1 can yield
⟨Vnt−˜Δ2hx,hyΦn+12+Φn+12,2Φnt⟩=⟨∂f(|˜Un+12|2,˜Φn+12)∂˜Φn+12Qn+12√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,2Φnt⟩. | (3.31) |
It is easy to check that
⟨Φn+12,Φnt⟩=1τ(‖Φn+1‖2−‖Φn‖2). |
Noting that equation ⟨Φnt,Vnt⟩=⟨Vnt,Φnt⟩, we obtain
1τ[(43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2−43‖Δ2hΦn‖2+13‖Δ22hΦn‖2)+h6(|Δ2hΦn+112|2−|Δ2hΦn12|2)+(|Δ2hΦn+1M−12|2−|Δ2hΦnM−12|2)]+1τ(‖Φn+1‖2−‖Φn‖2)+1τ(‖Vn+1‖2−‖Vn‖2)=⟨∂f(|˜Un+12|2,˜Φn+12)∂˜Φn+12Qn+12√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,Φnt⟩. | (3.32) |
Multiplying 2hQn+12j,k to Eq (3.27) and summing it up for 1≤j,k≤M−1 can yield
⟨Qnt,2Qn+12⟩=⟨⟨∂f(|˜Un+12|2,˜Φn+12)∂|~Un+12n+12|2,2Re(˜Un+12⋅Unt)⟩+⟨∂f(|˜U|2,˜Φn+12)∂˜Φn+12,Φnt⟩2√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,2Qn+12⟩. | (3.33) |
Combining Eqs (3.29), (3.32), (3.33), we obtain
43‖ΔαhUn+1‖2−13‖Δα2hUn+1‖2+43‖Δ2hΦn+1‖2−13‖Δ22hΦn+1‖2+h6(|Δ2hΦn+112|2+|Δ2hΦn+1M−12|2)+‖Φn+1‖2+‖Vn+1‖2−‖Qn+1‖2=43‖ΔαhUn‖2−13‖Δα2hUn‖2+43‖Δ2hΦn‖2−13‖Δ22hΦn‖2+h6(|Δ2hΦn12|2+|Δ2hΦnM−12|2)+‖Φn‖2+‖Vn‖2−‖Qn‖2. |
When n=0,1,2,⋯,n, we can obtain En=En−1⋯=E0.
In above subsection, we construct some structure-preserving schemes to solve two dimension fractional KGS equations. In similar method of one dimension case, we can obtain the iterative algorithms of numerical scheme (3.2)–(3.6) and numerical scheme (3.10)–(3.14). Now, we consider the iterative algorithm of numerical scheme (3.15)–(3.19) and numerical scheme (3.24)–(3.27). Let U=[U11,U12,⋯,UM−1M−1], V=[V11,V12,⋯,VM−1M−1], Φ=[Φ11,Φ12,⋯,ΦM−1M−1]. Then, we rewrite the numerical scheme (3.15)–(3.19) and numerical scheme (3.24)–(3.27) as the following vector form
iUn+1−Unτ−νAUn+1+Un2=−Un+1+Un2Φn+1+Φn2, | (3.34) |
Φn+1−Φn−12τ=Vn+1+Vn−12, | (3.35) |
Vn+1−Vn−12τ−BΦn+1+Φn−12=|Un|2, | (3.36) |
and
iUnt−νAUn+12=−˜Un+12∂f(|˜Un+12|2,˜Φn+12)∂|˜Un+12|2Qn+12√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0, | (3.37) |
Φnt=Vn+12, | (3.38) |
Vnt−BΦn+12+Φn+12=∂f(|˜Un+12|2,˜Φn+12)∂˜Φn+12Qn+12√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0, | (3.39) |
Qnt=⟨∂f(|˜Un+12|2,˜Φn+12)∂|~Un+12n+12|2,2Re(˜Un+12⋅(Un)t)⟩+⟨∂f(|˜U|2,˜Φn+12)∂˜Φn+12,(Φn)t⟩2√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0, | (3.40) |
where, matrice B represents differential matrix of two dimension Laplacian operator.
Consider numerical scheme (3.34)–(3.36), if (Un, Φn, Vn), n=0,1,2,⋯ are known, then Φn+1, Vn+1 of numerical scheme (3.34)–(3.36) is solve by the following linear equations
Φn+1−τVn+1=Φn−1+τVn−1,Vn+1−τBΦn+1=τBΦn−1+Vn−1+2τ|Un|2 |
Then, we can obtain Un+1 of numerical scheme (3.34)–(3.36) by solve linear equations
iUn+1−ντ2AUn+1+τ4Un+1(Φn+1+Φn)=iUn+ντ2AUn−τ4Un(Φn+1+Φn). |
Consider numerical scheme (3.37)–(3.40), the numerical solution Un+1, Φn+1, Vn+1 of numerical scheme (3.37)–(3.40) is solve by the following linear equations
iUn+1−ντ2AUn+1−τ2˜Un+12∂f(|˜Un+12|2,˜Φn+12)∂|˜Un+12|2Qn+1√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0=iUn+ντ2AUn−τ2˜Un+12∂f(|˜Un+12|2,˜Φn+12)∂|˜Un+12|2Qn√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,Φn+1−τVn+1=Φn−τVn,Vn+1−τ2BΦn+1+τ2Φn+1−τ2∂f(|˜Un+12|2,˜Φn+12)∂˜Φn+12Qn+1√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0=Vn+τ2BΦn−τ2Φn−τ2∂f(|˜Un+12|2,˜Φn+12)∂˜Φn+12Qn√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0,Qn+1−⟨∂f(|˜Un+12|2,˜Φn+12)∂|~Un+12n+12|2,2Re(˜Un+12⋅(Un)t)⟩+⟨∂f(|˜U|2,˜Φn+12)∂˜Φn+12,(Φn+1)⟩2√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0=Qn+⟨∂f(|˜Un+12|2,˜Φn+12)∂|~Un+12n+12|2,2Re(˜Un+12⋅(Un)t)⟩+⟨∂f(|˜U|2,˜Φn+12)∂˜Φn+12,(Φn)⟩2√⟨f(|˜Un+12|2,˜Φn+12),1⟩+C0. |
In above sections, we study some numerical schemes to solve following one dimension and two dimension fractional KGS systems:
iut−ν(−Δ)α2u=−u∂f(|u|2,ϕ)∂|u|2, | (4.1) |
ϕtt−Δϕ+ϕ=∂f(|u|2,ϕ)∂ϕ, | (4.2) |
and
iut−ν(−Δ)α2u=−uϕ, | (4.3) |
ϕtt−Δϕ+ϕ=|u|2. | (4.4) |
Recently, some structure-preserving schemes such as linearly implicit conservative scheme, symplectic scheme and multi-symplectic scheme have been designed and investigated for solving classical and fractional KGS system. However, in these works main system (4.3), (4.4) have been considered. As far as we know, there exist few studies on system (4.1), (4.2). In this paper, we consider not noly structure-preserving scheme of fractional KGS system (4.3), (4.4) but also fractional KGS system (4.1), (4.2). For two dimension case, we show linearly implicit conservative scheme (3.15)–(3.19) and fully implicit conservative scheme (3.10)–(3.14) to solve fractional KGS system (4.3), (4.4), and the fully implicit conservative scheme (3.10)–(3.14) is can obtain numerical result by above iterative algorithm. Moreover, we show also linearly implicit conservative scheme (3.24)–(3.27) and fully implicit conservative scheme (3.2)–(3.6) to solve fractional KGS system (4.1), (4.2). For one dimension case, we only give a fully implicit conservative scheme (2.1)–(2.5), in fact, linearly implicit(3.15)–(3.19) and (3.24)–(3.27) are still applicable to one dimensional case.
In this section, we give some numerical experiments to show the efficiency of the structure-preserving schemes. The first numerical example shows the numerical errors and convergence rates of the structure-preserving scheme, and check conservation property of the schemes. The second numerical example shows the numerical result for the general fractional KGS system. The third numerical example shows the numerical result for two dimension fractional KGS system.
When α=2, the KGS system (4.3), (4.4) has the following solitary wave solutions [5]
u(x,t,v)=3√24√1−v2sech212√1−v2(x−vt−x0)exp(i(vx+1−v2+v42(1−v2)t)),ϕ(x,t,v)=34(1−v2)sech212√1−v2(x−vt−x0), |
where v is the propagating velocity of the wave and x0 is the initial phase. We consider initial-value (v=0.8,x0=−10)
u0=3√24√1−v2sech212√1−v2(x−x0)exp(i(vx)),ϕ0=34(1−v2)sech212√1−v2(x−x0). |
In this example, we text errors, convergence orders and conservation of mass and energy of one dimension KGS system by numerical scheme (2.1)–(2.5). First, we show errors and convergence orders of numerical scheme (2.1)–(2.5) at time t=1. For α≠2, the numerical exact solutions are obtained by a very fine mesh and a small time step. Then, we fix the space mesh h=0.00001 and time step τ=0.00001 to test time convergence orders and space convergence orders by numerical scheme (2.1)–(2.5). The Figures 1–2 show time convergence orders and space convergence orders for different α, and it is found that the scheme is of order 2 in time, order 4 in space. From Figures 1–2, we can draw the observations: the approximate solution converge to the exact solution at the rate O(τ2+h4), and consistent with the theoretical estimates of Theorems 5–6. Second, we examine the conservation of mass and energy with x∈[−20,20],t∈[0,100], τ=0.01,h=0.1. Figure 3 shows relative residuals on the mass and energy errors for different values of α by numerical scheme (2.1)–(2.5). It is found that the numerical scheme preserves mass and energy conservation very well although energy varies with α.
Consider following one dimension fractional KGS system
i∂tu−12(−Δ)α2u=uϕ+γ|u|2u, | (4.5) |
∂ttϕ+ϕxx+ϕ=|u|2, | (4.6) |
u0=3√24√1−v2sech212√1−v2(x−x0)exp(i(vx)), | (4.7) |
ϕ0=34(1−v2)sech212√1−v2(x−x0). | (4.8) |
In this example, we simulate solitary wave and collisions of two solitary waves of fractional KGS system (4.5)–(4.8) by numerical scheme (2.1)–(2.5). First, we simulate solitary wave of numerical solutions for different orders α and difference parameter γ by the numerical scheme (2.1)–(2.5).
Figures 4–6 display solitary wave of the numerical solutions for difference value of γ=0.8, 1.5, 2 and same value of α=2. It is found that the parameter γ affects the propagation velocity of the solitary wave, and larger γ, the propagation of the soliton got slower.
Figures 7–9 display solitary wave of the numerical solutions for difference value of α=2, 1.8, 1.5 and same value of γ=1. It is found that the parameter α affects also the propagation velocity of the solitary wave, and smaller α, the propagation of the soliton got slower.
Second, we consider collisions of two solitary waves with x∈[−20,20], t∈[0,30], τ=0.01,h=0.1, and the initial data are chosen as (p1=10,p2=10,v1=0.8,v2=−0.8)
u0=u(x−p1,0,v1)+u(x−p2,0,v2),ϕ0=ϕ(x−p1,0,v1)+ϕ(x−p2,0,v2). |
Figures 10–12 display solitary wave of the numerical solutions for difference value of γ=0.8, 1.5, 2 and same value of α=2. Figures 13–15 display solitary wave of the numerical solutions for difference value of α=2, 1.8, 1.5 and same value of γ=1. It is found that the parameters α, γ affect also the propagation velocity of the solitary wave. When smaller α and larger γ, the propagation of the soliton got slower, the soliton changes faster and even a high oscillation appears.
In the two example, we show some numerical results of fractional KGS system (4.5)–(4.8) for solitary wave case and collisions of two solitary waves case. In the process of time evolution, the solitary wave moves towards the boundary gradually, and produces some small waves around the solitary wave. In this paper, we only consider spatial range [−20,20]. If we want to get a better numerical result, we can expand the space a bit, but it will take more calculation time.
Consider two dimension fractional KGS system
i∂tu−12(−Δ)α2u=−uϕ, | (4.9) |
∂ttϕ+Δϕ+ϕ=|u|2, | (4.10) |
u(x,y,0)=2ex2+2y2+e−(x2+2y2)⋅ei5/cosh(√4x2+y2), | (4.11) |
ϕ(x,y,0)=e−(x2+y2),ϕt(x,0)=e−(x2+y2)/2. | (4.12) |
In the example, we simulate solitary wave numerical solutions for different orders α by the scalar auxiliary variable scheme (3.24)–(3.27). Figures 16–18 show the surface plots of the nucleon density |u|2 and meson field ϕ for different time t=3,5,8 and same value of α=2, respectively.
Figures 19–21 show the surface plots of the nucleon density |u|2 and meson field ϕ for different time t=3,5,8 and same value of α=1.8, respectively. From Figures 16–21, we fine that the meson field change periodically, and the order α affects the shape of nucleon field. They also show that α affects the propagation velocity of the solitary wave.
In the paper, we study structure-preserving scheme to solve one dimension and two dimension space fractional KGS equations. First, we use the high central differences scheme in space and Crank-Nicolson scheme in time to discrete one dimension fractional KGS equations, which preserve mass and energy conservation laws of the fractional system. Then, we show that the arising scheme is uniquely solvable and approximate solutions converge to the exact solution at the rate O(τ2+h4). Second, we give the high central differences scheme in space, Crank-Nicolson scheme and scalar auxiliary variable scheme in time for two dimension fractional KGS equations, which preserve one or more analytical properties of the fractional system. Finally, the numerical experiments including some one dimensional and two dimensional fractional KGS systems are given to verify the correctness of theoretical results.
The authors declare there is no conflict of interest.
This work is supported by National Natural Science Foundation of China (No.12161070).
[1] |
H. G. Tanner, A. Jadbabaie, G. J. Pappas, Flocking in fixed and switching, IEEE Trans. Autom. Control, 52 (2007), 863-868. doi: 10.1109/TAC.2007.895948
![]() |
[2] |
Y. Cao, W. Yu, W. Ren, et al. An overview of recent progress in the study of distributed multi-agent coordination, IEEE Trans. Ind. Inform., 9 (2013), 427-438. doi: 10.1109/TII.2012.2219061
![]() |
[3] | X. Li, X. Yang, T. Huang, Persistence of delayed cooperative models: Impulsive control method, Appl. Math. Comput., 342 (2019), 130-146. |
[4] |
Z. Guan, Y. Wu, G. Feng, Consensus analysis based on impulsive systems in multiagent networks, IEEE Trans. Circuits Syst. I, 59 (2012), 170-178. doi: 10.1109/TCSI.2011.2158715
![]() |
[5] |
Y. Han, C. Li, W. Zhang, et al. Impulsive consensus of multiagent systems with limited bandwidth based on encoding-decoding, IEEE Trans. Cybern., 50 (2020), 1-12. doi: 10.1109/TCYB.2020.2984906
![]() |
[6] |
A. Jadbabaie, J. Lin, A. S. Morse, Coordination of groups of mobile autonomous agents using nearest neighbor rules, IEEE Trans. Autom. Control, 48 (2003), 988-1001. doi: 10.1109/TAC.2003.812781
![]() |
[7] |
H. Du, G. Wen, G. Chen, et al. A distributed finite-time consensus algorithm for higher-order leaderless and leader-following multiagent systems, IEEE Trans. Syst. Man Cybern., 47 (2017), 1625-1634. doi: 10.1109/TSMC.2017.2651899
![]() |
[8] |
B. Cui, C. Zhao, T. Ma, et al. Leaderless and leader-following consensus of multi-agent chaotic systems with unknown time delays and switching topologies, Nonlinear Anal. Hybrid Syst., 24 (2017), 115-131. doi: 10.1016/j.nahs.2016.11.007
![]() |
[9] |
D. Yang, X. Li, J. Qiu, Output tracking control of delayed switched systems via state-dependent switching and dynamic output feedback, Nonlinear Anal. Hybrid Syst., 32 (2019), 294-305. doi: 10.1016/j.nahs.2019.01.006
![]() |
[10] |
X. Tan, J. Cao, X. Li, et al. Leader-following mean square consensus of stochastic multi-agent systems with input delay via event-triggered control, IET Control Theory Appl., 12 (2018), 299-309. doi: 10.1049/iet-cta.2017.0462
![]() |
[11] |
Q. Zhang, S. Chen, C. Yu, Impulsive consensus problem of second-order multi-agent systems with switching topologies, Commun. Nonlinear Sci. Numer. Simulat., 17 (2012), 9-16. doi: 10.1016/j.cnsns.2011.04.007
![]() |
[12] |
Z. Xu, C. Li, Y. Han, Leader-following fixed-time quantized consensus of multi-agent systems via impulsive control, J. Franklin Inst., 356 (2019), 441-456. doi: 10.1016/j.jfranklin.2018.10.009
![]() |
[13] |
J. Hu, Y. Hong, Leader-following coordination of multi-agent systems with coupling time delays, Phys. A Stat. Mech. Appl., 374 (2007), 853-863. doi: 10.1016/j.physa.2006.08.015
![]() |
[14] |
F. Wang, Y. Yang, Leader-following exponential consensus of fractional order nonlinear multiagents system with hybrid time-varying delay: A heterogeneous impulsive method, Phys. A, 482 (2017), 158-172. doi: 10.1016/j.physa.2017.04.049
![]() |
[15] |
G. Wen, W. Yu, Y. Xia, et al. Distributed tracking of nonlinear multiagent systems under directed switching topology: An observer-based protocol, IEEE Trans. Syst. Man Cybern. Syst., 47 (2017), 869-881. doi: 10.1109/TSMC.2016.2564929
![]() |
[16] |
Z. Guan, Z. Liu, G. Feng, et al. Impulsive consensus algorithms for second-order multi-agent networks with sampled information, Automatica, 48 (2012), 1397-1404. doi: 10.1016/j.automatica.2012.05.005
![]() |
[17] | Y. Han, C. Li, Z. Zeng, et al. Exponential consensus of discrete-time non-linear multi-agent systems via relative state-dependent impulsive protocols, Neural Netw., 108 (2018) 192-201. |
[18] |
Y. Han, C. Li, Z. Zeng, Asynchronous event-based sampling data for impulsive protocol on consensus of non-linear multi-agent systems, Neural Netw., 115 (2019), 90-99. doi: 10.1016/j.neunet.2019.03.009
![]() |
[19] | X. Li, J. Shen, R. Rakkiyappan, Persistent impulsive effects on stability of functional differential equations with finite or infinite delay, Appl. Math. Comput., 329 (2018), 14-22. |
[20] |
X. Li, P. Li, Q. Wang, Input/output-to-state stability of impulsive switched systems, Syst. Control Lett., 116 (2018), 1-7. doi: 10.1016/j.sysconle.2018.04.001
![]() |
[21] |
X. Liu, C. Du, P. Lu, et al. Distributed event-triggered feedback consensus control with statedependent threshold for general linear multiagent systems, Internat. J. Robust Nonlinear Control, 27 (2017), 2589-2609. doi: 10.1002/rnc.3700
![]() |
[22] |
Z. Cao, C. Li, X. Wang, et al. Finite-time consensus of linear multi-agent system via distributed event-triggered strategy, J. Franklin Inst., 355 (2018), 1338-1350. doi: 10.1016/j.jfranklin.2017.12.026
![]() |
[23] |
T. Ma, Z. Zhang, C. Bing, Adaptive consensus of multi-agent systems via odd impulsive control, Neurocomputing, 321 (2018), 139-145. doi: 10.1016/j.neucom.2018.09.007
![]() |
[24] |
Y. Wang, J. Yi, Consensus in second-order multi-agent systems via impulsive control using position-only information with heterogeneous delays, IET Control Theory Appl., 9 (2015), 336-345. doi: 10.1049/iet-cta.2014.0425
![]() |
[25] |
F. Jiang, D. Xie, M. Cao, Dynamic consensus of double-integrator multi-agent systems with aperiodic impulsive protocol and time-varying delays, IET Control Theory Appl., 11 (2017), 2879-2885. doi: 10.1049/iet-cta.2016.1515
![]() |
[26] |
Q. Zhu, J. Cao, R. Rakkiyappan, Exponential input-to-state stability of stochastic Cohen-Grossberg neural networks with mixed delays, Nonlinear Dyn., 79 (2014), 1085-1098. doi: 10.1007/s11071-014-1725-2
![]() |
[27] |
Q. Zhu, pth Moment exponential stability of impulsive stochastic functional differential equations with Markovian switching, J. Franklin Inst., 351 (2014), 3965-3986. doi: 10.1016/j.jfranklin.2014.04.001
![]() |
[28] |
Q. Zhu, J. Cao, Stability analysis of Markovian jump stochastic BAM neural networks with impulse control and mixed time delays, IEEE Trans. Neural Netw. Learn. Syst., 23 (2012), 467-479. doi: 10.1109/TNNLS.2011.2182659
![]() |
[29] |
Y. Li, J. Lou, Z. Wang, et al. Synchronization of dynamical networks with nonlinearly coupling function under hybrid pinning impulsive controllers, J. Franklin Inst., 355 (2018), 6520-6530. doi: 10.1016/j.jfranklin.2018.06.021
![]() |
[30] |
C. Huang, J. Lu, D. W. C. Ho, et al. Stabilization of probabilistic Boolean networks via pinning control strategy, Inf. Sci., 510 (2020), 205-217. doi: 10.1016/j.ins.2019.09.029
![]() |
[31] |
X. Wang, C. Li, T. Huang, et al. Impulsive control and synchronization of nonlinear system with impulse time window, Nonlinear Dynam., 78 (2014), 2837-2845. doi: 10.1007/s11071-014-1629-1
![]() |
[32] |
Y. Feng, C. Li, Comparison system of impulsive control system with impulse time windows, J. Intell. Fuzzy Syst., 32 (2017), 4197-4204. doi: 10.3233/JIFS-16457
![]() |
[33] |
X. Wang, H. Wang, C. Li, et al. Synchronization of coupled delayed switched neural networks with impulsive time window, Nonlinear Dynam., 84 (2016), 1747-1757. doi: 10.1007/s11071-016-2602-y
![]() |
[34] |
X. Wang, J. Yu, C. Li, et al. Robust stability of stochastic fuzzy delayed neural networks with impulsive time window, Neural Netw., 67 (2015), 84-91. doi: 10.1016/j.neunet.2015.03.010
![]() |
[35] |
Y. Feng, C. Li, T. Huang, Periodically multiple state-jumps impulsive control systems with impulse time windows, Neurocomputing, 193 (2016), 7-13. doi: 10.1016/j.neucom.2016.01.059
![]() |
[36] |
T. Ma, Z. Zhang, C. Bing, Variable impulsive consensus of nonlinear multi-agent systems, Nonlinear Anal. Hybrid Syst., 31 (2019), 1-18. doi: 10.1016/j.nahs.2018.07.004
![]() |
[37] |
Li. L, C. Li, H. Li, An analysis and design for time-varying structures dynamical networks via state constraint impulsive control, Int. J. Control, 92 (2019), 2820-2828. doi: 10.1080/00207179.2018.1459861
![]() |
[38] |
Li. L, C. Li, H. Li, Fully state constraint impulsive control for non-autonomous delayed nonlinear dynamic systems, Nonlinear Anal. Hybrid Syst., 29 (2019), 383-394. doi: 10.1016/j.nahs.2018.03.008
![]() |
[39] | X. Liao, L. Wang, P, Yu, Stability of dynamical systems, Elsevier, 2007. |
1. | Sofia G. Florença, Ana C. Ferrão, Cristina A. Costa, Maria João Barroca, Aida Silva, Soraia I. Pedro, Ofélia Anjos, Raquel P. F. Guiné, Cultural Heritage on gastronomic usages of honey: Recipe analysis and cluster segmentation, 2025, 10, 2471-2086, 461, 10.3934/agrfood.2025023 |