1. Introduction
As we all know, it is difficult for planar polynomial differential systems to characterize their centers and isochronous centers. Hence, up to now the sufficient and necessary condition determining an isochronous center can only be found by making some appropriate analytic changes of coordinates which let the original system be reduced to a linear system. This kind of appropriate analytic change is very difficult to obtain, so only a handful of isochronous systems have been investigated. Several classes of known isochronous systems studied are as follows: quadratic isochronous centers (see [1]); isochronous centers of a linear center perturbed by third, fourth and fifth degree homogeneous polynomials (see [2–4]); complex polynomial systems (see [5]); reversible systems (see [6, 7]); and isochronous centers of cubic systems with degenerate infinity (see [8]).
In [9], a class of ninth degree system with four isochronous centers was investigated. In this paper, we consider the following ninth degree system:
(1)dxdt=-λy(x2+y2)4+y(p40x4+p31x3y+p22x2y2+p04y4)+y(p60x6+p42x4y2+p24x2y4+p06y6),dydt=λx(x2+y2)4+x(q40x4+q22x2y2+q13xy3+q04y4)+x(q60x6+q42x4y2+q24x2y4+q06y6),
where
(2)p40=-q04=12(-10λ+4A32+8A50-5β),p31=-q13=-8B50,p22=-q22=10λ-6A32-4A50+5β,p04=-q40=12(-2λ-β),p60=-q06=12(12λ-4A32-8A50+7β),p42=-q24=12(-4λ+8A32+16A50-9β),p24=-q42=-32(4λ-4A32-8A50+5β),p06=-q60=12(4λ+β), λβ≠0.
The paper is organized as follows. In Section 2, the preliminary methods to calculate the focal values (or Lyapunov constants) and period constants which are necessary in Sections 3 and 5 are presented. In Section 3, two appropriate transformations, changing system (1) to a class of Z4-equivariant fifth degree system in which the focal values are calculated, are provided. Based on it, the condition and proof that the infinity and four elementary singular points of (1) can be centers are gained. The bifurcations of 21 limit circles are the subject of Section 4. At last, in Section 5, the period constants with simpler expressions are given; moreover, the sufficient and necessary condition and proof are given to illustrate that the infinity and four elementary singular points of (1) become isochronous centers.
2. Some Preliminary Results
In [10, 11], the complex system of complex center is
(3)dzdT=z+∑k=2∞Zk(z,w)=Z(z,w),dwdT=-w-∑k=2∞Wk(z,w)=-W(z,w),
where
(4)Zk(z,w)=∑α+β=kaαβzαwβ, Wk(z,w)=∑α+β=kbαβwαzβ,
and the authors gave two recursive algorithms to determine necessary conditions for a center. We now restate the definitions and algorithms.
Theorem 1 (see [11]).
For system (3), we can derive successively the terms of the following formal series:
(5)M(z,w)=∑α+β=0∞cαβzαwβ,
such that
(6)∂(MZ)∂z-∂(MW)∂w=∑m=1∞(m+1)μm(zw)m,
where c00=1, for all ckk∈R, k=1,2,…, and for any integer m, μm is determined by the following recursive formulae:
(7)c00=1;when (α=β>0) or α<0, or β<0, cαβ=0;elsecαβ=1β-α∑k+j=3α+β+2[(α+1)ak,j-1-(β+1)bj,k-1]cα-k+1,β-j+1,μm=∑k+j=32m+2(ak,j-1-bj,k-1)cm-k+1,m-j+1.
Theorem 2 (see [12]).
For system (3), we can derive uniquely the following formal series:
(8)f(z,w)=z+∑k+j=2∞ckj′zkwj,g(z,w)=w+∑k+j=2∞dkj′wkzj,
where ck+1,k′=dk+1,k′=0, k=1,2,…, such that
(9)dfdT=f(z,w)+∑j=1∞pj′zj+1wj,dgdT=-g(z,w)-∑j=1∞qj′wj+1zj,
and when k-j-1≠0, ckj′ and dkj′ are determined by the following recursive formulae:
(10)ckj′=1j+1-k∑α+β=3k+j+1[(k-α+1)aα,β-1 -(j-β+1)bβ,α-1]ck-α+1,j-β+1′,dkj′=1j+1-k∑α+β=3k+j+1[(k-α+1)bα,β-1 -(j-β+1)aβ,α-1]dk-α+1,j-β+1′,
and for any positive integer j, pj′ and qj′ are determined by the following recursive formulae:
(11)pj′=∑α+β=32j+2[(j-α+2)aα,β-1 -(j-β+1)bβ,α-1]cj-α+2,j-β+1′,qj′=∑α+β=32j+2[(j-α+2)bα,β-1 -(j-β+1)aβ,α-1]dj-α+2,j-β+1′.
In the above expression, we have let c10′=d10′=1, c01′=d01′=0, and if α<0 or β<0, let aαβ=bαβ=cαβ′=dαβ′=0.
3. Singular Point Quantities and Center Conditions
It can be easily checked that system (1) have four elementary critical points (±1,0), (0,±1) and the infinity. After introducing the method to calculate the focal values and period constants of system (1), we try to make some appropriate transformations so as to carry out our investigation.
By means of the Bendixson homeomorphous transformation
(12)u=xx2+y2, v=yx2+y2,
and time transformation
(13)dt=(x2+y2)4dτ,
system (1) can be transformed into the following real system:
(14)dudτ=-λv+(-2λ+2A32+4A50-52β)u2v+(2λ+β2)v3+(3λ-2A32-4A50+32β)u4v+8B50u3v2+(2λ-2A32+4A50+β)u2v3-(λ+β2)v5,dvdτ=λu-(2λ+β2)u3+(2λ-2A32-4A50+52β)uv2+(λ+β2)u5+(-2λ+2A32-4A50-β)u3v2-8B50u2v3+(-3λ+2A32+4A50-32β)uv4.
Then, the infinity and four elementary focuses (±1,0), (0,±1) of (1) become the origin and (±1,0), (0,±1) of system (14), respectively. Thus, we have the following theorems.
Theorem 3.
System (14) is a class of Z4-equivariant differential polynomial system of fifth degree about point (-1,0).
Proof.
The system is invariant by the transformation of rotation
(15)u=xcosπ2-ysinπ2, v=xcosπ2+ysinπ2.
Then system (14) is a Z4-equivariant polynomial system about point (-1,0).
Now, we will consider the (0,0) of system (14). Let
(16)z=x+iy, w=x-iy, T=it.
System (14) can be transformed into the following complex system:
(17)dzdT=λz-(λ+a32+2a50+14β)wz2+(a32-λ+2a50-3β4)w3+a50z5+(λ-a32-a50+12β)w4z+a32w2z3,dwdT=-λw+(λ+b32+2b50+14β)zw2-(b32-λ+2b50-3β4)z3-b50w5-(λ-b32-b50+12β)z4w-b32z2w3,
where
(18)a50=A50+B50i2, b50=A50-B50i2,a32=A32-2B50i2, b32=A32+2B50i2.
Theorem 4.
The first two singular point quantities at the origin of system (17) are as follows:
(19)u1=0, u2=2iB50λ.
Theorem 5.
For system (17), the first two singular point quantities are zero if and only if the following condition holds:
(20)B50=0.
Theorem 6.
For system (14), all the singular point quantities at the origin are zero if and only if the first two singular point quantities are zero; that is, the condition in Theorem 5 holds. Correspondingly, the condition in Theorem 5 is the center condition of the origin.
Proof.
If condition (20) holds, system (14) becomes
(21)dudτ=-λv+(-2λ+2A32+4A50-52β)u2v+(2λ+β2)v3+(3λ-2A32-4A50+32β)u4v+(2λ-2A32+4A50+β)u2v3-(λ+β2)v5,dvdτ=λu-(2λ+β2)u3+(2λ-2A32-4A50+52β)uv2+(λ+β2)u5+(-2λ+2A32-4A50-β)u3v2+(-3λ+2A32+4A50-32β)uv4.
Remark (21) is symmetric with v axis.
System (14) is a class of Z4-equivariant differential polynomial system, so that the four symmetric elementary focus points of system (14) have the same topological structure; then they have the same center conditions and isochronous center conditions and so on. Without loss of generality, we only need to consider the case of (1,0).
Let
(22)x=u-1, y=v.
Then, system (14) can be changed into
(23)dxdt=-βy+(8λ-4A32-8A50+β)xy+(16λ-10A32-20A50+132β)x2y+(12λ-8A32-16A50+6β)x3y+(3λ-2A32-4A50+32β)x4y+8B50y2+8B50x3y2+(4λ-2A32+4A50+32β)y3+(4λ-4A32+8A50+2β)xy3+24B50xy2+24B50x2y2+(2λ-2A32+4A50+β)x2y3-(λ+12β)y5,dydt=βx+(4λ+72β)x2+(8λ+92β)x3+(5λ+52β)x4+(λ+12β)x5+(32β-8A50)y2+(4A32-4λ-16A50-12β)xy2+(6A32-6λ-12A50-3β)x2y2+(-2λ+2A32-4A50-β)x3y2+8B50y3+16B50xy3+8B50x2y3+(-3λ+2A32+4A50-32β)y4+(-3λ+2A32+4A50-32β)xy4.
By noting that
(24)z=x+iy, w=x-iy, T=it
system (23) can be transformed into its complex system:
(25)dzdT=βz+(3λ-2a32+34β)w2+(3λ-3a32-2a50+54β)w3+(λ-a32-a50+12β)w4+(2λ-8a50+52β)wz+(6λ-3a32-6a50+3β)w2z+(4λ-4a32-4a50+2β)w3z+(λ-a32-a50+12β)w4z+(2a32-λ+8a50+14β)z2+(5a32-λ-2a50+14β)wz2+3a32w2z2+(a32+10a50)z3+2a32wz3+a32w2z3+5a50z4+a50z5,dwdT=-βw-(3λ-2b32+34β)z2-(3λ-3b32-2b50+54β)z3-(λ-b32-b50+12β)z4-(2λ-8b50+52β)wz-(6λ-3b32-6b50+3β)z2w-(4λ-4b32-4b50+2β)z3w-(λ-b32-b50+12β)z4w-(2b32-λ+8b50+14β)w2-(5b32-λ-2b50+14β)zw2-3b32w2z2-(b32+10b50)w3-2b32zw3-b32z2w3-5b50w4-b50w5.
Through computing, we have the following theorem.
Theorem 7.
The first five singular point quantities at the origin of system (25) are
(26)u1=-4iB50(4λ-2A32-12A50+5β)β2,u2=8iB50(β-4A50)(2A32-28A50+5β)β3,u3=8iB50(β-4A50)2(-4480A502+560B502 +1992A50β-215β2)(β5)-1u4=-64iB50(β-4A50)2f175β7,u5=-16384iA503B50(β-4A50)2g1827307548828125β9,
where
(27)f=67737600A504-63598080A503β+22017616A502β2-3335040A50β3+186875β4,g=16146722095087449600A503-9886197819204477680A502β+2031858784197168536A50β2-142685023001803215β3.
Theorem 8.
For system (25), the first five singular point quantities are zero if and only if one of the following conditions holds:
(28)λ=12(A32-4A50), β=4A50;(29)B50=0.
Then we can get the following theorem.
Theorem 9.
For system (25), all the singular point quantities at the origin are zero if and only if the first five singular point quantities are zero; that is, one of the conditions in Theorem 8 holds. Correspondingly, the conditions in Theorem 8 are the center conditions of the origin.
Proof.
Under condition (28), the system (25) becomes
(30)dxdt=-4A50y-20A50xy-(2A32+26A50)x2y-(2A32+16A50)x3y-(12A32+4A50)x4y+8B50y2+24B50xy2+24B50x2y2+8B50x3y2+2A50y3+(8A50-2A32)xy3+(4A50-A32)x2y3-12A32y5,dydt=4A50x+2(A32+3A50)x2+2(2A32+A50)x3+52A32x4+12A32x5-2A50y2+2(A32-5A50)xy2+3(A32-4A50)x2y2+(A32-4A50)x3y2+8B50y3+16B50xy3+8B50x2y3+(12A32+4A50)y4+(12A32+4A50)xy4.
System (30) has an integrating factor
(31)u(x,y)=1(1+2x+x2+y2)3.
Under condition (29), system (25) becomes
(32)dxdt=-βy+(8λ-4A32-8A50+β)xy+(16λ-10A32-20A50+132β)x2y+(12λ-8A32-16A50+6β)x3y+(3λ-2A32-4A50+32β)x4y+(4λ-2A32+4A50+32β)y3+(4λ-4A32+8A50+2β)xy3+(2λ-2A32+4A50+β)x2y3-(λ+12β)y5,dydt=βx+(4λ+72β)x2+(8λ+92)βx3+(5λ+52β)x4+(λ+12β)x5+(-8A50+32β)y2+(-4λ+4A32-16A50-12β)xy2+(-6λ+6A32-12A50-3β)x2y2+(-2λ+2A32-4A50-β)x3y2+(-3λ+2A32+4A50-32β)y4+(-3λ+2A32+4A50-32β)xy4,
where (32) is symmetric with x axis.
So we get the following theorem.
Theorem 10.
The system (14) has five centers (0,0), (±1,0), (0,±1) when B50=0. Namely, the system (1) has five centers the infinity and (±1,0), (0,±1) when B50=0.
The system (14) has four centers (±1,0), (0,±1) when λ=(1/2)(A32-4A50), β=4A50, B50≠0. Namely, the system (1) has four centers (±1,0), (0,±1) when λ=(1/2)(A32-4A50), β=4A50, and B50≠0.
4. Bifurcation of Limit Circles
Theorem 11.
For system (14), the (1,0) is a fifth order weak focus if and only if
(33)λ=12(20A50-5β),A32=12(28A50-5β),B502=4480A502-1992A50β+215β2560,f=0, g≠0.
Proof.
When u1=u2=u3=u4=0, we could find suitable β and A50 to satisfy
(34)B502>0, u4=0,(35)Resultant[f,g,A50]=16309978866916163079012825599795663295988359261138540000000000000000000000000000000000β12.
So the (1,0) is a fifth order weak focus.
Theorem 12.
Assume that the conditions in Theorem 11 hold. Then the 5 limit cycles can bifurcate from the (1,0) of System (14), which lie in the neighborhood of the origin.
Proof.
If the conditions in Theorem 11 hold, we can obtain that the Jacobian of the u1,u2,u3,u4 with respect to λ,A32,B50,A50(36)J=D(u1,u2,u3,u4)D(λ,A32,B50,A50) =-2097152B505(-4A50+β)4175β17 ×(203212800A504-192864000A503β +67884512A50β2-10506964A50β3 +603755β4);Resultant[f,J,A50]=-661345270173473293422286219719475200000000β16.
So J≠0. Applying the theory in [13], the result is completed.
Theorem 13.
Assume that the conditions in Theorem 11 hold. Then 1 limit cycle can bifurcate from the (0,0) of System (14), which lie in the neighborhood of the origin.
Since System (14) is Z4-equivariant, taking into account Theorems 12 and 13, we conclude the following.
Theorem 14.
Assume that the conditions in Theorem 11 hold. Then 5 limit cycles can bifurcate from each of the four fine focus points (±1,0), (0,±1) of System (14), and one limit circle could be created from the origin at the same time; namely, 21 limit cycles can bifurcate from the System (14) in all, five of which are located in the neighborhood of each of the four fine focus points (±1,0), (0,±1) and one of which is located in the neighborhood of origin.
In fact, there are 21 limit cycles for system (1); 5 limit cycles can bifurcate from each of the four fine focus points (±1,0), (0,±1) of System (1), and one limit circle could be created from the infinity at the same time; namely, 21 limit cycles can bifurcate from the system (1) in all. It is a very interesting result.
5. Period Constants and Isochronous Center Conditions
In this section, we focus on the calculation of period constants and isochronous center conditions of the origin and (1,0) for system (17) and (25). According to Theorems 2, we begin by computing period constants through recursive formulae. Under the center conditions given in Section 3, we discuss three different relevant cases: (a) about the origin, case condition (20) holds.
Theorem 15.
Under condition (20), the origin of system (17) is a complex isochronous center if and only if the following conditions holds
(37)A32=B50=0, A50=-2λ, β=-4λ.
Proof.
If condition (20) holds,substituting B50=0 into the recursive formulae in Theorem 2, we can easily get the first four period constants
(38)τ1=-4λ+2A32+4A50-β2λ,τ2=-64λ2+16λA50+20λβ+3β28λ2,τ3=-(4λ+β)(80λ2+46λβ+9β2)32λ3,τ4=3β(4λ+β)(104λ2+40λβ+9β2)512λ4.
Clearly, τ1=τ2=τ3=τ4=0 implies that condition (37) holds.
On the other hand, case condition (37) holds, and system (17) is changed into
(39)dzdT=-z(z+1)(z-1)(z2+1),dwdT=(-1+w)w(1+w)(1+w2).
By means of transformation
(40)ξ=z(z4-1)-1/4, η=w(w4-1)-1/4,
system (39) is reduced to a linear system.
Thus, the origin of system (17) is a complex isochronous center.
(b) about (1,0) of system (25).
When condition (28) holds, substituting λ=(1/2)(A32-4A50), β=4A50, into the recursive formulae in Theorem 2, we easily obtain the first two period constants
(41)τ1=-A32-2A50A32-4A50, τ2=-2.
Thus (1,0) could not be an isochronous center when condition (28) holds.
Theorem 16.
(
1,0
)
of system (25) is a complex isochronous center if and only if condition (37) holds.
Proof.
When condition (29) holds, substituting B50 into the recursive formulae in Theorem 2, we easily obtain the first four period constants:
(42)τ1=-4λ+10A32-4A50+β2,τ2=(80λ-13392λ2+80A50+17856λA50 -5952A502-20β-15624λβ +10416A50β-4557β2)(200)-1,τ3=18000(-237312λ2+2939328λ3 +208896λA50-5878656λ2A50 -46592A502-870912A503 -299104λβ+5143824λ2β +162816A50β-6858432A50λβ +2286144A502β-103592β2 +3000564λβ2-2000376A50β2 +583443β3+3919104λA502),τ4=140000(-51200λ2+28097664λ3+17356032λ4 +57600λA50-29498368λ2A50 -46282752λ3A50+6545472λA502 +46282752λ2A502+421504A503 -20570112λA503+3428352A504 -70400λβ+40046432λ2β +40497408λ3β+39600A50β -11200A502-80994816λ2A50β +996112A502β+53996544λA502β -11999232A503β-24200β2 +17793272λβ2+35435232λ2β2 -3856888A50β2-47246976λA50β2 -24098656λA50β+15748992A502β2 +2329754β3+13780368λβ3 -9186912A50β3+2009637β4).
Denoting by
(43)g1=R(τ2,τ3,λ), g2=R(τ2,τ4,λ),
where R(τi,τj,A) is the resultant of τi,τj with respect to A, we have R(g1,g2,β)=0. From polynomioalremainder[g2,g1,A50]=0, we could get that β=2A50. When β=2A50, τ2=τ3=τ4=0 deduce A32=B50=0, A50=-2λ, β=-4λ.
Besides, when condition (37) holds, system (25) is changed into
(44)dzdT=14z(z+1)(z+2)(z2+2z+2),dwdT=-14(1+w)w(2+w)(2+2w+w2).
Letting that
(45)ξ=z(2+z)(2+2z+z2)(z+1)4,η=w(2+w)(2+2w+w2)(w+1)4,
system (44) is reduced to a linear system.
From Theorems 15 and 16, we can conclude the following theorem.
Theorem 17.
The infinity and four elementary singular points (±1,0), (0,±1) are five isochronous centers of (1) if and only if condition (37) holds; namely A32=B50=0, A50=-2λ, β=-4λ.