Research article

Multiple periodic solutions of differential delay systems with 2k lags

  • Received: 15 December 2020 Accepted: 14 April 2021 Published: 21 April 2021
  • MSC : 34A34, 34K13

  • The quantity of 2(2k+1)-periodic solutions to a specific differential delay system with 2k lags is studied and resolved by variational methods. Several results are revealed and two examples are given to illustrate the application of the main results.

    Citation: Li Zhang, Huihui Pang, Weigao Ge. Multiple periodic solutions of differential delay systems with 2k lags[J]. AIMS Mathematics, 2021, 6(7): 6815-6832. doi: 10.3934/math.2021399

    Related Papers:

    [1] Yameng Duan, Wieslaw Krawcewicz, Huafeng Xiao . Periodic solutions in reversible systems in second order systems with distributed delays. AIMS Mathematics, 2024, 9(4): 8461-8475. doi: 10.3934/math.2024411
    [2] Zhifeng Lu, Fei Wang, Yujuan Tian, Yaping Li . Lag synchronization of complex-valued interval neural networks via distributed delayed impulsive control. AIMS Mathematics, 2023, 8(3): 5502-5521. doi: 10.3934/math.2023277
    [3] Kasbawati, Mariani, Nur Erawaty, Naimah Aris . A mathematical study of effects of delays arising from the interaction of anti-drug antibody and therapeutic protein in the immune response system. AIMS Mathematics, 2020, 5(6): 7191-7213. doi: 10.3934/math.2020460
    [4] Mouataz Billah Mesmouli, Amir Abdel Menaem, Taher S. Hassan . Effectiveness of matrix measure in finding periodic solutions for nonlinear systems of differential and integro-differential equations with delays. AIMS Mathematics, 2024, 9(6): 14274-14287. doi: 10.3934/math.2024693
    [5] Yanshou Dong, Junfang Zhao, Xu Miao, Ming Kang . Piecewise pseudo almost periodic solutions of interval general BAM neural networks with mixed time-varying delays and impulsive perturbations. AIMS Mathematics, 2023, 8(9): 21828-21855. doi: 10.3934/math.20231113
    [6] Huahai Qiu, Li Wan, Zhigang Zhou, Qunjiao Zhang, Qinghua Zhou . Global exponential periodicity of nonlinear neural networks with multiple time-varying delays. AIMS Mathematics, 2023, 8(5): 12472-12485. doi: 10.3934/math.2023626
    [7] Yan Yan . Multiplicity of positive periodic solutions for a discrete impulsive blood cell production model. AIMS Mathematics, 2023, 8(11): 26515-26531. doi: 10.3934/math.20231354
    [8] Ruizhi Yang, Dan Jin, Wenlong Wang . A diffusive predator-prey model with generalist predator and time delay. AIMS Mathematics, 2022, 7(3): 4574-4591. doi: 10.3934/math.2022255
    [9] Aziz Belmiloudi . Time-varying delays in electrophysiological wave propagation along cardiac tissue and minimax control problems associated with uncertain bidomain type models. AIMS Mathematics, 2019, 4(3): 928-983. doi: 10.3934/math.2019.3.928
    [10] Yun Xin, Hao Wang . Positive periodic solution for third-order singular neutral differential equation with time-dependent delay. AIMS Mathematics, 2020, 5(6): 7234-7251. doi: 10.3934/math.2020462
  • The quantity of 2(2k+1)-periodic solutions to a specific differential delay system with 2k lags is studied and resolved by variational methods. Several results are revealed and two examples are given to illustrate the application of the main results.



    Delay differential equations naturally appear in many fields of science and engineering. Such equations have been proposed as models for a variety of physiological processes and conditions including production of blood cells, respiration and cardiac arrhythemias. In recent years, the existence of positive periodic solutions for delay differential equations has been received considerable attention. J. Kaplan and J. Yorke [10] gave conditions for the existence of 4 or 6 periodic solutions of equation

    x(t)=ni=1f(x(ti)),xR (1.1)

    while n=1 and n=2, respectively. After then, lots of results are achieved for this equation in general cases [1,2,3,4,5,6,7,8,9,10,11,12]. W. Ge [4,5,6] proved several existence theorems of periodic solutions for (1.1) in the general case by use of the fixed-point theorem in cone. After transforming (1.1) into a Hamiltonian system and applying a theorem given by Mawhin and Willem [13], J.Li and X. He [11,12] proved a theorem for the existence of several periodic solutions under the condition xf(x)>0(xf(x)<0) for x0. G. Fei [1,2] studied the multiplicity of periodic solutions of (1.1) and proved theorems for the cases n=2k1 and n=2k by using the index theory. In his work, the way to construct the functional for the second case is quite different from that for the first case.

    Applying S1 index theory Z. Guo and J. Yu [8] gave a theorem for the multiple periodic solutions of the delay-differential system with one lag in the form

    x(t)=F(x(tl)),xRN,FC0(RN,RN). (1.2)

    And by applying the same index theory, in [9], they obtained results on the multiplicity of 2(n+1)-periodic solutions to the delay differential system

    x(t)=nt=1F(x(ti)),xRN,FC0(RN,RN). (1.3)

    when n=2k1. After then B. Zheng and Z. Guo studied (1.3) while n=2k and gave two criteria for the multiplicity of 2(2k+1)-periodic solutions. The theorems given in both[9] and [14] contain a condition to calculate the following function

    Ψ(A,B)=j=1Ψj(A,B),n=2k1,
    Ψ(A,B)=j=1,jΓΨj(A,B),ΓN+,n=2k,

    for symmetric matrices A and B, respectively. This is a hard task to be fulfilled since it is not definite for the problem how large the j must be to ensure Ψj(A,B)=0.

    In [7], we researched the same problem for system (1.3) when n=2k1 and by constructing a new functional we gave criteria for the multiplicity of periodic solutions only based on the eigenvalues of symmetric matrices A and B.

    In this paper we are to study the multiplicity of 2(2k+1)-periodic solutions for system (1.3) when n=2k, i.e.,

    x(t)=2ki=1F(x(ti)),xRN,FC0(RN,RN). (1.4)

    with the conditions FC1(RN,R),F(x)=F(x),F(0)=0,

    x(tk1)=x(t). (1.5)

    At the same time, assume there are real symmetric matrices A and A0 such that

    F(x)=Ax+o(|x|)as|x|,F(x)=A0x+o(|x|)as|x|0.

    In section 2, some notations and variation structure for system (1.4) are introduced. Meanwhile, a lemma about the relation between the periodic solutions of system(1.4) and the critical points of the functional Φ are also proved in section 2. Section 3 presents and proves the main results of this paper. Before the proof of the main results some lemmas about the calculation of the differential of functional Φ are proven. Two examples are given in section 4 to illustrate the applications of the main results.

    Firstly, we consider a linear space which consists of all the 2(2k+1)-periodic functions satisfying (1.5) in RN, i.e.

    H12={x(t)=i=1(aicos(2i1)π2k+1t+bisin(2i1)π2k+1t):ai,biRN,i=1(2i1)(|ai|2+|bi|2)<}

    Since our goal is to find the 2(2k+1)-periodic solutions for system (1.4) with condition (1.5), we discuss now the special requirement for a solution of system (1.4) in H12. After then, we are able to choose a suitable Hilbert space for our problem. Assume ni=mai=0 if n<m.

    If x(t) is a 2(2k+1)-periodic solution of (1.4) satisfying (1.5), then, x(tl)=l1i=0F(x(ti))2ki=l+1F(x(ti)),l=0,1,,2k. Hence,

    2kl=0(1)lx(tl)=0. (2.1)

    From the equation above, one has 2kl=0(1)lx(tl)=c. However, c=2kl=0(1)lx(tl2k1)=2kl=0(1)lx(tl)=c, which means c=0. Then,

    2kl=0(1)lx(tl)=0. (2.2)

    By Fourier's expansion theory and (1.5), it holds that

    x(t)=i=1(aicos(2i1)π2k+1t+bisin(2i1)π2k+1t)=l=02k+1i=1(al(2k+1)+icos(2l(2k+1)+2i1)π2k+1t+bl(2k+1)+isin(2l(2k+1)+2i1)π2k+1t).

    Substituting the above expression into (2.2) we have

    al(2k+1)+k+1=bl(2k+1)+k+1=0. (2.3)

    Then suppose

    X={xH12:x(t2k1)=x(t),2ki=0(1)ix(ti)=0}.

    For xX, define P:XX:

    xl=0(ki=1+2k+1i=k+2)(2l(2k+1)+2i1)(al(2k+1)+icos(2l(2k+1)+2i1)π2k+1t+bl(2k+1)+isin(2l(2k+1)+2i1)π2k+1t),

    and denote

    x,yX=2(2k+1)0(Px(t),y(t))dt,x,y=2(2k+1)0(x(t),y(t))dt,x,yX.

    Define

    xX=x,xX,x2=x,x=P1x,xX.

    It is obvious that (X,X) is a Hilbert space and P:XX is an invertible and self-adjoint operator. Define

    Z={z(t)}={(z1(t),z2(t),,z2k(t))}={(x(t),x(t1),,x(t2k+1)):xX}X2k.

    Therefore z1(t)=x(t)X. Furthermore, let

    X(i)={acos(2i1)π2k+1t+bsin(2i1)π2k+1t:a,bRN},
    Z(i)={(x(t),x(t1),,x(t2k+1)):x(t)X(i)}.

    Then,

    Z=i=1il(2k+1)+k+1Z(i)=l=0ki=1Z(l(2k+1)+i).

    Denote H(z)=2k1i=0F(x(ti))+F(2k1i=0(1)i+1x(ti)), then,

    H(z)=(H(z)z1,H(z)z2,,H(z)z2k)=(F(x(t))F(x(t2k)),F(x(t1))+F(x(t2k)),,F(x(t1))(1)iF(x(t2k)),,F(x(t2k+1))+F(x(t2k)))

    where x(t2k)=2k1i=0(1)i+1x(ti), zl=x(tl+1),l=1,2,,2k, and from (2.2),

    z2k+1=2k1i=0(1)i+1x(ti)=2ki=1(1)izi.

    Define Φ:ZR as a functional of z:

    Φ(z)=12Lz,z+G(z) (2.4)

    where

    Lz=P1(2k1i=1(1)i+1x(ti),,l1i=0(1)i+lx(ti)+2k1i=l+1(1)i+l+1x(ti),,2k2i=0(1)i+2k1x(ti)),G(z)=2(2k+1)0(2k1i=0F(x(ti))+F(2k1i=0(1)i+1x(ti)))dt,z,y=2(2k+1)0(ˆPz,y)dt,z,yZ,

    ˆP=(P,P,,P):ZZ. Then

    z=z,z=2kx,xX=2kx2X, (2.5)
    Φ(z)=122(2k+1)02k1l=0(l1i=0(1)i+lx(ti)+2k1i=l+1(1)i+l+1x(ti),x(tl))dt2(2k+1)0[2k1i=0F(x(ti))+F(2k1i=0(1)i+1x(ti))]dt=120i<l2k12(2k+1)0(1)i+l+1[(x(ti),x(tl))(x(ti),x(tl))]dt2(2k+1)0[2k1i=0F(x(ti))+F(2k1i=0(1)i+1x(ti))]dt.

    By Mawhin-Willem's theorem [13, Theorem 1.4], Φ is continuously differentiable and

    Φ(z),v=Lu,vZ2(2k+1)0(2k1l=0(F(x(tl))(1)lF(x(t2k)),y(tl)))dt=2k1l=02(2k+1)0(l1i=0(1)i+lx(ti)+2k1i=l+1(1)i+l+1x(ti)F(x(tl))(1)l+1F(x(t2k)),y(tl)))dt

    where z(t)=(x(t),x(t1),,x(t2k+1)),v(t)=(y(t),y(t1),,y(t2k+1)).

    Denote Φ(z)=(Φ1(z),Φ2(z),,Φ2k(z)), where

    Φl+1(z)=P1(l1i=0(1)i+lx(ti)+2k1i=l(1)i+l+1x(ti)F(x(tl))+(1)lF(x(t2k))),l=0,1,,2k1. (2.6)

    If a point zZ satisfies Φ(z)=0, then z is called a critical point of Φ. The following lemma reveals the relation between the 2(2k+1)-periodic solutions of system (1.4) and the critical points of Φ. Based on this lemma we can discuss the multiplicity of critical points of Φ instead of discussing that of 2(2k+1)-periodic solutions of system (1.4). In the following, denote x(t) as the first N components of z(t) in Z and set 1i=0x(ti)=0.

    Lemma 2.1 If xX, z=(x(t),x(t1),,x(t2k+1))Z, then the following propositions are equivalent to each other.

    1) x(t) is one of the solutions of system (1.4),

    2) z(t) satisfies Φ1(z)=0, i.e.

    2k1i=1(1)i+1x(ti)F(x(t))+F(x(t2k))=0, (2.7)

    3) z(t) is one of the critical points of Φ on Z, i.e., z(t) satisfies

    l1i=0(1)i+lx(ti)+2k1i=l+1(1)i+l+1x(ti)F(x(tl))+(1)lF(x(t2k))=0,l=0,1,,2k1(2.8l)

    Proof. 1) 2).

    Suppose x(t) is a solution of system (1.4).

    By the periodicity of x(t) and (1.5), one has

    x(tl)=l1i=0F(x(ti))2k1i=l+1F(x(ti))F(x(t2k))(2.9l)

    l=1,2,,2k1, x(t2k)=2k1i=0(1)i+1x(ti). Multiplying (1)l+1 to the equalities (2.9l) and summing them up, we have

    2k1l=1(1)l+1x(tl)=F(x(t))F(x(t2k)),

    which means (2.7) holds.

    2) 3).

    System (2.7) is the same as (2.80).

    From (2.7) one has

    2k1i=1(1)i+1x(ti1)F(x(t1))+F(x(t2k1))=0,

    i.e.,

    0=2k1i=2(1)ix(ti)+x(t2k)F(x(t1))F(x(t))=2k1i=2(1)ix(ti)+2k1i=0(1)i+1x(ti)F(x(t1))F(x(t)). (2.10)

    Then (2.7) and (2.10) yield

    x(t)+2k1i=2(1)ix(ti)F(x(t1))F(x(t2k))=0,

    which means (2.81).

    Suppose (2.8l) holds. Then,

    0=li=1(1)i+l+1x(ti)+2k1i=l+2(1)i+l+2x(ti)+(1)lx(t2k)F(x(tl1))+(1)l+1F(x(t))=li=1(1)i+l+1x(ti)+2k1i=l+2(1)i+l+2x(ti)+(1)l+12k1i=0(1)ix(ti)F(x(tl1))+(1)l+1F(x(t)). (2.11)

    Multiplying (2.7)with (1)l+1 and add it to (2.11), one can obtain (2.8(l+1)).

    3) 1).

    Summing the equalities from (2.81) to (2.8(2k1)), one has

    x(t)=2ki=1F(x(ti)),

    which implies x(t) is a solution of system (1.4).

    Denote ^Pi:ZZ(i), Z(i)={(x(t),x(t1),,x(t2k+1)):x(t)=aicos(2i1)π2k+1t+bisin(2i1)π2k+1t,ai,biRN},i(2l+1)(2k+1); Z((2l+1)(2k+1))=0,lN+{0}. Furthermore, let Pi=ij=1ˆPj:Zij=1Z(i),i=1,2,. Then, one has LPi=PiL.

    Let G(z)=(P1(F(x(t))F(x(t2k))),P1(F(x(t1))+F(x(t2k))),,P1(F(x(tl))+(1)l+1F(x(t2k))),,P1(F(x(t2k+1))F(x(t2k)))). Then Φ(z)=Lz+G(z). L is a bounded self-adjoint operator and G(z):ZX2k is compact since (Gz)(t) is differentiable with respect to t.

    Denote B,B0R2kN×2kN as

    B=(2k+12kA,,2k+12kA),B0=(2k+12kA0,,2k+12kA0)

    and for z(t)=(x(t),x(t1),,x(t2k+1)),

    Bz=(2k+12kAx(t),2k+12kAx(t1),,2k+12kAx(t2k+1)),B0z=(2k+12kA0x(t),2k+12kA0x(t1),,2k+12kA0x(t2k+1)). (2.12)

    Define δz(t)=T2kNz(t):ZZ, where

    T2kN=[IIIIIOOOOIOOOOIO] (2.13)

    is a 2kN×2kN matrix, I is a N×N unit matrix and O is a N×N zero matrix.

    From (1.5) and (2.2), one has δz(t)=z(t+1)andδ2k+1z(t)=z(t). Then, {1,δ,δ2,,δ4k+1} is a Lie group with δ4k+2=1, where 1 stands for the identity transform in Z.

    Now we have

    Φ(δz)=Φ(z),
    Φ(δz)=δΦ(z),

    i.e., with respect to the Lie group {1,δ,δ2,,δ4k+1}, Φ is invariant and Φ is δequivariant.

    In order to prove our results we shall apply the following lemma.

    Lemma 2.2 [2, Lemma 2.4]

    Assume there are two S1-invariant linear subspaces, Z+ and ZZ, and r>0 such that

    (a) Z+Z is a closed and finite codimension in Z,

    (b) ˆL(Z)Z with ˆL=LP1B or ˆL=LP1B0,

    (c) there exists cR such that

    Φ(z)c<Φ(0),zZSr={zZ:z=r},

    (d) there exists c0R such that

    infzZ+Φ(z)c0,

    (e) Φ satisfies (P.S)c-condition for c0cc.

    Then Φ possesses at least m=12[dim(ZZ+)codimZ(ZZ+)] different critical orbits in Φ1([c0,c]) provided m>0.

    Remark 2.1 (P.S)c-condition in (e) can be replaced by (P.S)-condition since (P.S)-condition implies (P.S)c-condition for each cR.

    Remark 2.2 If (P.S)c-condition in (e) is replaced by (P.S)-condition, (c) and (d) can be changed into

    (c) there exists cR such that

    Φ(z)c>Φ(0),zZSr={zZ:z=r},

    (d) there exists c0R such that

    infzZ+Φ(z)c0

    if Φ(z) is replaced by Φ(z).

    Remark 2.3 When Z=i=1Z(i), Z+(i)=Z+Z(i), Z(i)=ZZ(i), we have dimZ(i)=2N and

    m=12i=1[dim(Z+(i)Z(i))codimZ(i)(Z+(i)+Z(i))]=12i=1[dim(Z+(i)Z(i))(dimZ(i)dim(Z+(i)+Z(i)))]=12i=1[dimZ+(i)+dimZ(i)2N]. (2.14)

    Suppose {αj:j=1,2,,N} are the eigenvalues of A and {βj:j=1,2,,N} are the eigenvalues of A0, dj,ejRN are the eigenvectors of A and A0 with respect to αj and βj, respectively. Assume {d1,d2,,dN} and {e1,e2,,eN} are two orthogonal bases of RN.

    Lemma 3.1 If zZ(i),il(k+1), then there is γ=π2k+1 such that

    Lz,z=2k+12k(2i1)γtan(2i1)γ2z,z. (3.1)

    Proof. Denote || as the norm in RN, then

    z,z=2k2(2k+1)0|x(t)|2dt=2k(2k+1)(|ai|2+|bi|2). (3.2)

    It follows that

    Lz,z=0m<l2k1(1)m+l+12(2k+1)0[(x(tm),x(tl))(x(tm),x(tl))]dt=20m<l2k1(1)m+l+12(2k+1)0(x(tm),x(tl))dtj=lm=22k1j=12k1l=j(1)2lj+12(2k+1)0(x(tl+j),x(tl))dt=22k1j=1(1)j+1(2kj)2(2k+1)0(x(t+j),x(t))dt=22k1j=1(2kj)(1)j+1(2i1)γ(2k+1)(|a|2+|b|2)sinj(2i1)γ=2(2i1)γ(2k+1)2k1j=1(2kj)(1)j+1sinj(2i1)γ(|a|2+|b|2).

    Since

    2k1j=1(2kj)(1)j+1sinj(2i1)γ=12cos(2i1)γ22k1j=1(1)j+1(2kj)[sin(2j1)(2i1)γ2+sin(2j+1)(2i1)γ2]=12cos(2i1)γ2[2ksin(2i1)γ2+2kj=1(1)jsin(2j1)(2i1)γ2]=ktan(2i1)γ2+12cos(2i1)γ22kj=1(1)j[sinj(2i1)γcos(2i1)γ2cosj(2i1)γsin(2i1)γ2]=(k+12)tan(2i1)γ2,

    one has

    Lz,z=(2k+1)2(2i1)γtan(2i1)γ2(|ai|2+|bi|2)=2k+12k(2i1)γtan(2i1)γ2z,z.

    Lemma 3.1 is proved.

    Denote Ej={λej:λR}, Dj={λdj:λR}, j=1,2,,N, l0 and

    indαj={ki=1{l:0<[2l(2k+1)+(2i1)]γtan(2i1)γ2<αj},αj>0,0,αj=0,ki=1{l:0<[(2l+1)(2k+1)+2(k+1i)]γtan(2i1)γ2<αj},αj<0,
    indβj={ki=1{l:0<[2l(2k+1)+(2i1)]γtan(2i1)γ2<βj},βj>0,0,βj=0,ki=1{l:0<[(2l+1)(2k+1)+2(k+1i)]γtan(2i1)γ2<βj},βj<0,

    then RN=Nj=1Dj=Nj=1Ej. Suppose

    indA=Nj=1indαj,indA0=Nj=1indβj, (3.3)

    and

    XD,j(i)={x(t)=aicos(2i1)γ2t+bisin(2i1)γ2t,ai,biDj},
    XE,j(i)={x(t)=aicos(2i1)γ2t+bisin(2i1)γ2t,ai,biEj},

    therefore,

    X(i)=Nj=1XD,j(i)=Nj=1XE,j(i). (3.4)

    Denote

    x(t)=l=0ki=1[al(2k+1)+icos(2l(2k+1)+2i1)γt+bl(2k+1)+isin(2l(2k+1)+2i1)γt+al(2k+1)+k+1+icos((2l+1)(2k+1)+2i)γt+bl(2k+1)+k+1+isin((2l+1)(2k+1)+2i)γt].

    It follows from (3.1) and (3.2) that for z1(t)=x(t), one has

    Lz,z=γ(2k+1)2l=0ki=1[(2l(2k+1)+2i1)tan(2l(2k+1)+2i1)γ2(|al(2k+1)+i|2+|bl(2k+1)+i|2)+(2l(2k+1)+2(k+1+i)1)tan(2l(2k+1)+2(k+1+i)1)γ2(|al(2k+1)+k+1+i|2+|bl(2k+1)+k+1+i|2)]=γ(2k+1)2l=0ki=1[(2l(2k+1)+2i1)tan(2i1)γ2(|al(2k+1)+i|2+|bl(2k+1)+i|2)((2l+1)(2k+1)+2(k+1i))tan(2i1)γ2(|al(2k+1)+2k+2i|2+|bl(2k+1)+2k+2i|2)]. (3.5)

    and

    P1Ax,x=2(k+1)0(Ax,x)dt=(2k+1)αj(|ai|2+|bi|2),xXD,j(i),
    P1A0x,x=2(k+1)0(A0x,x)dt=(2k+1)βj(|ai|2+|bi|2),xXE,j(i),

    and

    ˆP1Bz,z=αj(2k+1)2(|ai|2+|bi|2),xXD,j,ˆP1B0z,z=βj(2k+1)2(|ai|2+|bi|2),xXE,j. (3.6)

    Denote

    Z={zZ:LzˆP1Bz,z<0},Z0={zZ:LzˆP1B0z,z<0},
    Z0={zZ:LzˆP1Bz,z=0},Z00={zZ:LzˆP1B0z,z=0},
    Z+={zZ:LzˆP1Bz,z>0},Z+0={zZ:LzˆP1B0z,z>0}.

    Then

    Z=ZZ0Z+=Z0Z00Z+0. (3.7)

    From (3.1) and Lemma 3.1 [7], we have the following lemmas.

    Lemma 3.2 ˆL(Z)Z, ˆL(Z0)Z0, ˆL(Z+)Z+, ˆL(Z0)Z0, ˆL(Z00)Z00, ˆL(Z+0)Z+0, (ˆL=LˆP1B or ˆL=LˆP1B0).

    Lemma 3.3 All the subspaces of Z,

    Z0+Z+,Z+0+Z,Z+0+Z+Z00,Z0+Z++Z00,Z+0+Z+Z0,
    Z0+Z++Z0,Z0+Z++Z00+Z0,Z+0+Z+Z00+Z0,

    are of finite codimensions in Z.

    Denote Γ+={αj>0:there arel0,i{1,2,,k}such thatαj=γ(2l(2k+1)+2i1)tan(2i1)γ2},Γ={αj<0:there arel0,i{1,2,,k}such thatαj=γ((2l+1)(2k+1)+2(k+1i))tan(2i1)γ2},Γ+0={βj>0:there arel0,i{1,2,,k}such thatβj=γ(2l(2k+1)+2i1)tan(2i1)γ2},Γ0={βj<0:there arel0,i{1,2,,k}such thatβj=γ((2l+1)(2k+1)+2(k+1i))tan(2i1)γ2},Γ=Γ+Γ,Γ0=Γ+0Γ0.

    For i{1,2,,k} and l0, denote

    η+={(l,i):there isαjΓ+such thatγ(2l(2k+1)+2i1)tan(2i1)γ2=αj},

    η={(l,i):there isαjΓsuch thatγ((2l+1)(2k+1)+2(k+1i))tan(2i1)γ2=αj},

    η+0={(l,i):there isβjΓ+0such thatγ(2l(2k+1)+2i1)tan(2i1)γ2=βj},

    η0={(l,i):there isβjΓ0such thatγ((2l+1)(2k+1)+2(k+1i))tan(2i1)γ2=βj}.

    Let D={Dj:αjΓ} and Π:RND be an orthogonal projection with ΠRN=D. Assume

    (A1) FC1(RN,R) satisfies (1.4) and F(x)=F(x),F(0)=0,

    (A2) there exist M and rC0(R+,R+) with r()=+, r(s)s0 as s, such that |F(x)12(Ax,x)|M+r(Πx) whenever x{Dj:αjΓ},

    (A±3) ±[F(x)12(Ax,x)]>0,as|x|,

    (A±4) ±[F(x)12(A0x,x)]>0,0<|x|1.

    By a standard argument as the proof of Lemma 2.1 [13] and Lemma 3.3 [8], we have

    Lemma 3.4 Assume (A1) and (A2) hold. Then, Φ(x) defined by (2.4) satisfies (P.S)-condition.

    Lemma 3.5 Suppose (A1) and (A2) hold. Then there is I>0 such that

    m=1i<il(2k+1)+k+1[dim(Z+(i)Z(i))codimZ(i)(Z+(i)+Z(i))]=1iIil(2k+1)+k+1[dim(Z+(i)Z(i))codimZ(i)(Z+(i)+Z(i))],

    if

    (Z+,Z){(Z+,Z0),(Z+0,Z),(Z,Z+0),(Z0,Z+),(Z++Z0,Z0),(Z+,Z0+Z00),(Z0,Z++Z0),(Z0+Z00,Z+),(Z+0,Z+Z0),(Z+0+Z00,Z),(Z+Z0,Z+0),(Z,Z+0+Z00),(Z++Z0,Z0+Z00),(Z+0+Z00,Z+Z0),(Z+Z0,Z+0+Z00),(Z0+Z00,Z++Z0)}.

    Now we give the main results of this paper.

    Theorem 3.1 Assume (A1) and (A2) hold. Then system (1.4) possesses at least

    m=max{ind(A)ind(A0)ηη+0,ind(A0)ind(A)η0η+}

    different 2(2k+1)-periodic orbits satisfying x(t(2k+1))=x(t) provided m>0.

    Corollary 3.1 Suppose (A1) and (A2) hold.

    ⅰ) If Γ+Γ=, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)η+0,ind(A0)ind(A)η0}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅱ) If Γ+0Γ0=, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)η,ind(A0)ind(A)η+}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅲ) If Γ+0Γ0=Γ+Γ=Γ+0Γ0=, then system (1.4) possesses at least

    m=max{ind(A)ind(A0),ind(A0)ind(A)}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    Corollary 3.1 can be directly obtained from Theorem 3.1.

    Theorem 3.2 Suppose (A1) and (A2) hold.

    ⅰ) If (A+3) holds, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)+η+η+0,ind(A0)ind(A)η0η+}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅱ) If (A3) holds, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)ηη+0,ind(A0)ind(A)η0+η}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅲ) If (A+4) holds, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)ηη+0,ind(A0)ind(A)+η+0η+}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅳ) If (A4) holds, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)η+η0,ind(A0)ind(A)η0η+}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅴ) If (A+3), (A4) hold, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)+η++η0,ind(A0)ind(A)η0η+}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅵ) If (A3), (A+4) hold, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)ηη+0,ind(A0)ind(A)+η+η+0}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅶ) If (A+3), (A+4) hold, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)+η+η0,ind(A0)ind(A)+η+0η+}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    ⅷ) If (A3), (A4) hold, then system (1.4) possesses at least

    m=max{ind(A)ind(A0)η+η0,ind(A0)ind(A)+η+η0}

    different 2(2k+1)-periodic orbits satisfying x(t2k1)=x(t) provided that m>0.

    Proof of Theorem 3.1.

    Suppose m=ind(A0)ind(A)η0η+>0.

    Let Z+=Z and Z=Z+0,zZ+. There is σ>0 such that

    12LzˆP1Bz,zσx22.

    Then,

    |F(x)12P1Ax,xX|σ2(2k+1)x22+M.

    Therefore,

    Φ(z)=12Lz,z+G(z)=12Lz,z2(2k+1)02ki=0F(x(ti))dt=12Lz,z(2k+1)2(2k+1)0F(x(t))dt=12Lz12ˆP1Bz,z(2k+1)2(2k+1)0[F(x(t))12P1Ax,xX]dt12σx222(2k+1)M+,x22,

    which implies that there is c0R such that

    Φ(z)c0,zZ+.

    At the same time, there are cR and γ>0 such that

    Φ(z)c.

    Suppose c0<c. Then conditions (c) and (d) in Lemma 2.2 are satisfied.

    Since Lemmas 2.1 and 3.2–3.4 imply (P.S)-condition as well as conditions (a) and (b) in Lemma 2.2 hold under the requirement of Theorem 3.1, we only need to compute m.

    From Lemma 3.5 we have

    m=121iIil(2k+1)+k+1[dimZ+(i)Z(i)codimZ(i)(Z+(i)+Z(i))]=121iIil(2k+1)+k+1[dimZ+(i)+dimZ(i)dimZ(i)]=121iIil(2k+1)+k+1[dimZ+(i)+dimZ(i)2N]. (3.8)

    Then,

    m=121iIil(2k+1)+k+1[dim(Z+(i)+dimZ0(i))2N]=12l(2k+1)+iIi{1,2,,k},l0[dimZ+(l(2k+1)+i)+dimZ0(l(2k+1)+i)2N]+12l(2k+1)+k+1+iIi{1,2,,k},l0[dimZ+(l(2k+1)+k+1+i)+dimZ0(l(2k+1)+k+1+i)2N]=l(2k+1)+iIi{1,2,,k},l0[dimZ+(l(2k+1)+i)+dimZ0(l(2k+1)+i)2N]+12l(2k+1)+2(k+1)iIi{1,2,,k},l0[dimZ+(l(2k+1)+2(k+1)i)+dimZ0(l(2k+1)+2(k+1)i)2N].

    Denote I1={l(2k+1)+iI:l0,i{1,2,,k}},I2={l(2k+1)+2(k+1)iI:l0,i{1,2,,k}}, and indA+=αj>0indαj,indA=αj<0indαj,indA+0=βj>0indβj,indA0=βj<0indβj, then,

    indA=indA++indA,indA0=indA+0+indA0.

    Obviously,

    l(2k+1)+iIi{1,2,,k},l0dimZ0(l(2k+1)+i)=2indA+0,

    l(2k+1)+iIi{1,2,,k},l0dimZ+(l(2k+1)+i)=l(2k+1)+iIi{1,2,,k},l0[2NdimZ(l(2k+1)+i)dimZ0(l(2k+1)+i)]=2NI12indA+2η+,

    l(2k+1)+2(k+1)iIi{1,2,,k},l0dimZ0(l(2k+1)+2(k+1)i)=l(2k+1)+2(k+1)iIi{1,2,,k},l0[2NdimZ00dimZ+0]=2NI2+2indA02η0.

    Therefore,

    m=indA+0+indA0indA+indAη+η0=indA0indAη+η0.

    Theorem 3.1 is proved.

    Proof of Theorem 3.2.

    Since the proof of Theorem 3.2 for each case is similar, we prove the theorem only for case (ⅱ).

    Without loss of generality, suppose

    indA0indAη0+η>max{0,indAindA0ηη+0},

    and denote Z+=Z++Z0,Z=Z0.

    From (A3), there is M>0 such that

    [F(x)12(Ax,x)]>M,xRN.

    Then, for zZ++Z0,

    Φ(z)=12Lz,z+G(z)=12Lz12ˆP1Bz,z+G(z)+12ˆP1Bz,zG(z)+12ˆP1Bz,z=(2k+1)2(2k+1)0[F(x(t))12(Ax,x)]dt>(2k+1)2(2k+1)0(M)dt=2M(2k+1)2,

    which implies that condition (d) in Lemma 2.2 holds.

    From (3.8), one has

    m=121iIil+k+1[dimZ+(i)+dimZ(i)2N]=121iIil(2k+1)+k+1[dimZ+(i)+dimZ0(i)+dimZ0(i)2N]=121iIil(2k+1)+k+1[dimZ+(i)+dimZ0(i)2N]+121iIil(2k+1)+k+1dimZ0(i)=indA0indAη+η0+12dimZ0=indA0indAη+η0+η++η=indA0indAη0+η.

    Theorem 3.2 is proved.

    Example 4.1 Consider the number of 6-periodic solutions of the following system

    x(t)=F(x(t1))F(x(t2)) (4.1)

    where x=(x1,x2)R2,

    F(x)=π123(9x21+22x1x2+9x22)(x431+2x432),|x|,
    F(x)=π123(29x21+223x1x2+7x22)+(x21+x22)32,|x|0.

    Obviously, FC1(R2,R),F(0)=0,F(x)=F(x),

    F(x)=[3π2311π6311π633π23][x1x2][43x13183x132],|x|
    F(x)=[29π63113π63113π637π63][x1x2]+[3(x21+x22)2x13(x21+x22)12x2],|x|0

    and

    A=[3π2311π6311π633π23],A0=[29π6311π611π67π63].

    The eigenvalues of A and A0 are

    α1=π33,α2=10π33andβ1=π33,β2=10π33,

    respectively. Since k=1, one has

    indα1={l0:0<13(6l+1)π13=π33(6l+1)<π33}=0,indα2={l0:0<π33(6l+5)<10π33}=1,indβ1={l0:0<π33(6l+5)<π33}=0,indβ2={l0:0<π33(6l+1)<10π33}=2.

    Therefore,

    indA=1,indA0=2.

    At the same time,

    η+=1,η=0,η+0=0,η0=0,

    and all the conditions of Theorem 3.2 (vi) are satisfied. Then system (4.1) has at least

    m=indA0indA+η+η+0=3

    different 6-periodic orbits satisfying x(t3)=x(t).

    Example 4.2 Let N=2. We discuss the multiplicity of 6-periodic solutions of system (4.1), where x=(x1,x2)R2,

    F(x)=(25π63x21+5π3x1x25π23x22x41x62)(1φ(|x|))+(5π63x215π3x1x2+5π63x22+x431+x852)φ(|x|).

    Then

    F(x)=A0x[4x316x52],|x|0,
    F(x)=Ax+[43x13185x352],|x|,

    where

    A0=[25π335π35π35π3],A=[5π335π35π35π33].

    A and A0 have their eigenvalues α1=10π33,α2=20π33, and β1=10π3,β2=10π33, respectively.

    So

    indα1={l0:0<π33(6l+5)<10π33}=1,indα2={l0:0<π33(6l+1)<20π33}=4,indβ1={l0:0<π33(6l+5)<10π3}=5,indβ2={l0:0<π33(6l+5)<10π33}=1,

    and then

    indA=3,indA0=6.

    On the other hand, we have

    η+=η=η+0=η0=0.

    Therefore,

    m=indAindA0=9,

    and system (4.1) has at least 9 different 6-periodic orbits satisfying x(t3)=x(t) by Theorem 3.1.

    Supported by R & D Program of Beijing Municipal Education Commission (No. KM202011417010), Academic Research Projects of Beijing Union University (No. ZB10202004), Scientific Research Projects of Beijing Union University (No. WZ10201902).

    The authors declare that there is no conflicts of interest regarding the publication of this paper.



    [1] G. Fei, Multiple periodic solutins of differential delay equations vai Hamiltonian systems(Ⅰ), Nonlinear Anal. Theor., 65 (2006), 25–39. doi: 10.1016/j.na.2005.06.011
    [2] G. Fei, Multiple periodic solutins of differential delay equations vai Hamiltonian systems(Ⅱ), Nonlinear Anal. Theor., 65 (2006), 40–58. doi: 10.1016/j.na.2005.06.012
    [3] W. Ge, Periodic solutions of the differential delay equation x(t)=f(x(t1)), Acta Math. Sin., 12 (1996), 113–121.
    [4] W. Ge, On the existence of periodic solutions of the differential delay equations with multiple lags, Acta Math. Appl. Sin., 17 (1994), 173–181.
    [5] W. Ge, Two existence theorems of periodic solutions for differential delay equations, Chinese Ann. Math., 15 (1994), 217–224.
    [6] W. Ge, Oscillatory periodic solutions of differential delay equations with multiple lags, Chinese Sci. Bull., 42 (1997), 444–447. doi: 10.1007/BF02882587
    [7] W. Ge, L. Zhang, Multiple periodic solutions of delay differential systems with 2k-1 lags via variational approach, Discrete Cont. Dyn. A, 36 (2016), 4925–4943. doi: 10.3934/dcds.2016013
    [8] Z. Guo, J. Yu, Multiplicity results for periodic solutions to delay differential equations via critical point theory, J. Differ. Equations, 218 (2005), 15–35. doi: 10.1016/j.jde.2005.08.007
    [9] Z. Guo, J. Yu, Multiplicity results on period solutions to higher dimensional differential equations with multiple delays, J. Dyn. Differ. Equ., 23 (2011), 1029–1052. doi: 10.1007/s10884-011-9228-z
    [10] J. Kaplan, J. Yorke, Ordinary differential equations which yield periodic solutions of differential delay equations, J. Math. Anal. Appl., 48 (1974), 317–324. doi: 10.1016/0022-247X(74)90162-0
    [11] J. Li, X. He, Multiple periodic solutions of differential delay equations created by asymptotically linear Hailtonian systems, Nonlinear Anal. Theor., 31 (1998), 45–54. doi: 10.1016/S0362-546X(96)00058-2
    [12] J. Li, X. He, Proof and gengeralization of Kaplan-Yorke' conjecture under the condition f(0)>0 on periodic solution of differential delay equations, Sci. China Ser. A, 42 (1999), 957–964. doi: 10.1007/BF02880387
    [13] J. Mawhin, M. Willem, Critical point theory and Hamiltonian systems, New York, Springer-Verlag, 1989.
    [14] B. Zheng, Z. Guo, Multiplicity results on periodic solutions to higher-dimensional differential equations with multiple delays, Rocky Mountain J. Math., 44 (2014), 1715–1744.
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2498) PDF downloads(94) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog