About this Journal Submit a Manuscript Table of Contents
Mathematical Problems in Engineering
VolumeΒ 2011Β (2011), Article IDΒ 810217, 23 pages
doi:10.1155/2011/810217
Research Article

Trigonometric Function Periodic Wave Solutions and Their Limit Forms for the KdV-Like and the PC-Like Equations

Department of Mathematics, South China University of Technology, Guangzhou 510640, China

Received 15 February 2011; Accepted 26 May 2011

Academic Editor: MingΒ Li

Copyright Β© 2011 Liu Zhengrong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We use the bifurcation method of dynamical systems to study the periodic wave solutions and their limit forms for the KdV-like equation 𝑒 𝑑 + π‘Ž ( 1 + 𝑏 𝑒 ) 𝑒 𝑒 π‘₯ + 𝑒 π‘₯ π‘₯ π‘₯ = 0 , and PC-like equation 𝑣 𝑑 𝑑 βˆ’ 𝑣 𝑑 𝑑 π‘₯ π‘₯ βˆ’ ( π‘Ž 1 𝑣 + π‘Ž 2 𝑣 2 + π‘Ž 3 𝑣 3 ) π‘₯ π‘₯ = 0 , respectively. Via some special phase orbits, we obtain some new explicit periodic wave solutions which are called trigonometric function periodic wave solutions because they are expressed in terms of trigonometric functions. We also show that the trigonometric function periodic wave solutions can be obtained from the limits of elliptic function periodic wave solutions. It is very interesting that the two equations have similar periodic wave solutions. Our work extend previous some results.

1. Introduction

Many authors have investigated the KdV-like equation 𝑒 𝑑 + π‘Ž ( 1 + 𝑏 𝑒 ) 𝑒 𝑒 π‘₯ + 𝑒 π‘₯ π‘₯ π‘₯ = 0 , ( 1 . 1 ) and the PC-like equation 𝑣 𝑑 𝑑 βˆ’ 𝑣 𝑑 𝑑 π‘₯ π‘₯ βˆ’ ξ€· π‘Ž 1 𝑣 + π‘Ž 2 𝑣 2 + π‘Ž 3 𝑣 3 ξ€Έ π‘₯ π‘₯ = 0 . ( 1 . 2 ) For example, Dey [1, 2] studied the exact Himiltonian density and the conservation laws, and gave two kink solutions for (1.1). Zhang et al. [3, 4] gave some solitary wave solutions and singular wave solutions for (1.1) by using two different methods. Yu [5] got an exact kink soliton for (1.1) by using homogeneous balance method. Grimshaw et al. [6] studied the large-amplitude solitons for (1.1). Fan [7, 8] gave some bell-shaped soliton solutions, kink-shaped soliton, and Jacobi periodic solutions for (1.1) by using algebraic method. Tang et al. [9] investigated solitary waves and their bifurcations for (1.1) by employing bifurcation method of dynamical systems. Peng [10] used the modified mapping method to get some solitary wave solutions composed of hyperbolic functions, periodic wave solutions composed of Jacobi elliptic functions, and singular wave solution composed of triangle functions for (1.1). Chow et al. [11] described the interaction between a soliton and a breather for (1.1) by using the Hirota bilinear method. Kaya and Inan [12] studied solitary wave solutions for (1.1) by using Adomian decomposition method. Yomba [13] used Fan's subequation method to construct exact traveling wave solutions composed of hyperbolic functions or Jacobi elliptic functions for (1.1).

Zhang and Ma [14] gave some explicit solitary wave solutions composed of hyperbolic functions by using solving algebraic equations for (1.2). Li and Zhang [15] used bifurcation method of dynamical system to study the bifurcation of traveling wave solutions and construct solitary wave solutions for (1.2). Kaya [16] discussed the exact and numerical solitary wave solutions by using a decomposition method for (1.2). Rafei et al. [17] gave numerical solutions by using He's method for (1.2).

Recently, many authors have presented some useful methods to deal with the problems in equations, for instance [1830].

In this paper, we use the bifurcation method mentioned above to study the periodic wave solutions for (1.1) and (1.2). Through some special phase orbits, we obtain new expressions of periodic wave solutions which are composed of trigonometric functions sin  πœ‰ or cos  πœ‰ . These solutions are called trigonometric function periodic wave solutions. We also check the correctness by using the software Mathematica.

In Section 2, we will state our results for (1.1). In Section 3, we will state our results for (1.2). In Sections 4, and 5, we will give derivations for our main results. Some discussions and the orders for testing the correctness of the solutions will be given in Section 6.

2. Trigonometric Function Periodic Wave Solutions for (1.1)

In this section, we state our main results for (1.1). In order to state these results conveniently, we give some preparations. For given constant 𝑐 β‰  0 , on π‘Ž βˆ’ 𝑏 plane we define some lines and regions as follows.( 1 )When 𝑐 < 0 , we define lines 𝑙 1 𝑙 : 𝑏 = 0 , 2 π‘Ž : 𝑏 = βˆ’ , 𝑙 6 𝑐 3 : 𝑏 = βˆ’ 3 π‘Ž , 𝑙 1 6 𝑐 4 : π‘Ž = 0 , ( 2 . 1 ) and regions 𝐴 𝑖 ( 𝑖 = 1 - 8 ) , as Figure 1(a).

fig1
Figure 1: The locations of the lines 𝑙 𝑖 , π‘˜ 𝑖 ( 𝑖 = 1 , 2 , 3 , 4 ) and the regions 𝐴 𝑗 , 𝐡 𝑗 ( 𝑗 = 1 , 2 , … , 8 ) for given constant 𝑐 β‰  0 .
( 2 )When 𝑐 > 0 , we define lines π‘˜ 1 π‘˜ : 𝑏 = 0 , 2 π‘˜ : π‘Ž = 0 , 3 : 𝑏 = βˆ’ 3 π‘Ž , π‘˜ 1 6 𝑐 4 π‘Ž : 𝑏 = βˆ’ , 6 𝑐 ( 2 . 2 ) and regions 𝐡 𝑖 ( 𝑖 = 1 - 8 ) , as Figure 1(b).

Using the lines and regions in Figure 1, we narrate our results as follows.

Proposition 2.1. For arbitrary given constant 𝑐 β‰  0 , let πœ‰ = π‘₯ βˆ’ 𝑐 𝑑 . ( 2 . 3 ) Then, (1.1) has the following periodic wave solutions. (1)When c < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 1 or 𝐴 5 , the expression of the periodic wave solution is 𝑒 1 ( πœ‰ ) = 6 𝑐 √ π‘Ž + π‘Ž ξ‚€ √ ( π‘Ž + 6 𝑏 𝑐 ) c o s  , βˆ’ 𝑐 πœ‰ ( 2 . 4 ) which has the following limit forms.( i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 1 and ( π‘Ž , 𝑏 ) tends to the line 𝑙 1 , 𝑒 1 ( πœ‰ ) tends to the periodic blow-up solution 𝑒 ∘ 1 ( πœ‰ ) = 6 𝑐 π‘Ž  ξ‚€ √ 1 + c o s βˆ’ 𝑐 πœ‰  ξ‚„ ( 2 . 5 ) (see Figure 2).( i i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 5 and ( π‘Ž , 𝑏 ) tends to the line 𝑙 1 , 𝑒 1 ( πœ‰ ) tends to the periodic blow-up solution 𝑒 βˆ— 1 ( πœ‰ ) = 6 𝑐 π‘Ž  ξ‚€ √ 1 βˆ’ c o s βˆ’ 𝑐 πœ‰  ξ‚„ ( 2 . 6 ) (see Figure 3).( i i i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 1 or 𝐴 5 , and ( π‘Ž , 𝑏 ) tends to 𝑙 2 , 𝑒 1 ( πœ‰ ) tends to the trivial solution 𝑒 ( πœ‰ ) = 6 𝑐 / π‘Ž .( 2 )When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 2 , or when 𝑐 > 0 and ( π‘Ž , 𝑏 ) ∈ 𝐡 5 , the expression of the periodic wave solution is 𝑒 2 𝛼 ( πœ‰ ) = 0 ξ€· 𝑀 c o s 0 πœ‰ ξ€Έ + 𝛽 0 𝑝 0 ξ€· 𝑀 c o s 0 πœ‰ ξ€Έ + π‘ž 0 , ( 2 . 7 ) where Ξ” = 3 π‘Ž ( 3 π‘Ž + 1 6 𝑏 𝑐 ) , ( 2 . 8 ) 𝛼 0 ξ‚€ √ = βˆ’ 3 π‘Ž + Ξ”  ξ‚™ π‘Ž ξ‚€ √ π‘Ž βˆ’ Ξ”  4 π‘Ž 2 𝑏 2 , 𝛽 ( 2 . 9 ) 0 √ = βˆ’ 3 π‘Ž + 2 4 𝑏 𝑐 + Ξ” 2 π‘Ž 𝑏 2 𝑝 , ( 2 . 1 0 ) 0 = ξ‚™ π‘Ž ξ‚€ √ π‘Ž βˆ’ Ξ”  π‘ž π‘Ž 𝑏 , ( 2 . 1 1 ) 0 1 = βˆ’ ξ‚€ √ π‘Ž 𝑏 π‘Ž + Ξ”  𝑀 , ( 2 . 1 2 ) 0 = ξƒŽ √ 3 π‘Ž + 1 6 𝑏 𝑐 + Ξ” . 8 𝑏 ( 2 . 1 3 ) The solution 𝑒 2 ( πœ‰ ) has the following limit forms. ( i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 2 and ( π‘Ž , 𝑏 ) tends to 𝑙 3 , the 𝑒 2 ( πœ‰ ) tends to the peak-shaped solitary wave solution 𝑒 ∘ 2 ξ€· ( πœ‰ ) = 4 𝑐 3 + 2 𝑐 πœ‰ 2 ξ€Έ π‘Ž ξ€· βˆ’ 9 + 2 𝑐 πœ‰ 2 ξ€Έ ( 2 . 1 4 ) (see Figure 4).( i i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 2 and ( π‘Ž , 𝑏 ) tends to 𝑙 2 , 𝑒 2 ( πœ‰ ) tends to the trivial solution 𝑒 ( πœ‰ ) = 0 .( i i i )When 𝑐 > 0 , ( π‘Ž , 𝑏 ) ∈ 𝐡 5 and ( π‘Ž , 𝑏 ) tends to π‘˜ 1 , the 𝑒 2 ( πœ‰ ) tends to the periodic blow-up solution 𝑒 βˆ— 2 𝑐 ξ‚€ ( πœ‰ ) = 2 s i n 2 ξ‚€ √   𝑐 πœ‰ / 2 βˆ’ 3 π‘Ž s i n 2 ξ‚€ √  𝑐 πœ‰ / 2 ( 2 . 1 5 ) (see Figure 5).( 3 )When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 6 , or when 𝑐 > 0 and ( π‘Ž , 𝑏 ) ∈ 𝐡 1 , the expressions of the solution is 𝑒 3 𝛼 ( πœ‰ ) = 1 ξ€· 𝑀 c o s 1 πœ‰ ξ€Έ + 𝛽 1 𝑝 1 ξ€· 𝑀 c o s 1 πœ‰ ξ€Έ + π‘ž 1 , ( 2 . 1 6 ) where 𝛼 1 = ξ‚€ √ βˆ’ 3 π‘Ž + Ξ”  ξ‚™ π‘Ž ξ‚€ √ π‘Ž + Ξ”  4 π‘Ž 2 𝑏 2 , 𝛽 1 = √ 3 π‘Ž + 2 4 𝑏 𝑐 βˆ’ Ξ” 2 π‘Ž 𝑏 2 , 𝑝 1 = ξ‚™ π‘Ž ξ‚€ √ π‘Ž + Ξ”  , π‘ž π‘Ž 𝑏 1 = √ π‘Ž βˆ’ Ξ” , 𝑀 π‘Ž 𝑏 1 = ξƒŽ √ 3 π‘Ž + 1 6 𝑏 𝑐 βˆ’ Ξ” . 8 𝑏 ( 2 . 1 7 ) The solution 𝑒 3 ( πœ‰ ) has the following limit forms.( i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 6 and ( π‘Ž , 𝑏 ) tends to 𝑙 3 , the 𝑒 3 ( πœ‰ ) tends to the canyon-shaped solitary wave (see Figure 6) solution 𝑒 ∘ 2 ( πœ‰ ) . ( i i )When 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 6 and ( π‘Ž , 𝑏 ) tends to 𝑙 2 , 𝑒 3 ( πœ‰ ) tends to the trivial solution 𝑒 ( πœ‰ ) = 0 . ( i i i )When 𝑐 > 0 , ( π‘Ž , 𝑏 ) ∈ 𝐡 1 and ( π‘Ž , 𝑏 ) tends to π‘˜ 1 , the 𝑒 3 ( πœ‰ ) tends to the periodic blow-up wave solution 𝑒 βˆ— 1 ( πœ‰ ) (see Figure 3).

fig2
Figure 2: The limiting precess of 𝑒 1 ( πœ‰ ) when 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 1 , and ( π‘Ž , 𝑏 ) tends to the line 𝑙 1 , where π‘Ž = 4 and 𝑐 = βˆ’ 1 .
fig3
Figure 3: The limiting precess of 𝑒 1 ( πœ‰ ) when 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 5 , and ( π‘Ž , 𝑏 ) tends to the line 𝑙 1 , where π‘Ž = βˆ’ 4 and 𝑐 = βˆ’ 1 .
fig4
Figure 4: The limiting precess of 𝑒 2 ( πœ‰ ) when 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 2 , and ( π‘Ž , 𝑏 ) tends to the line 𝑙 3 , where π‘Ž = 4 and 𝑐 = βˆ’ 1 .
fig5
Figure 5: The limiting precess of 𝑒 2 ( πœ‰ ) when 𝑐 > 0 , ( π‘Ž , 𝑏 ) ∈ 𝐡 5 , and ( π‘Ž , 𝑏 ) tends to the line π‘˜ 1 , where π‘Ž = βˆ’ 2 and 𝑐 = 1 .
fig6
Figure 6: The limiting precess of 𝑒 3 ( πœ‰ ) when 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ 𝐴 6 , and ( π‘Ž , 𝑏 ) tends to the line 𝑙 3 , where π‘Ž = βˆ’ 9 and 𝑐 = βˆ’ 1 .

Remark 2.2. Note that if 𝑒 = πœ‘ ( πœ‰ ) is a solution of (1.1), then 𝑒 = πœ‘ ( πœ‰ + π‘Ÿ ) also is solution of (1.1), where π‘Ÿ is a arbitrary constant. According to this fact and the results listed in Proposition 2.1, the following nine functions also are periodic wave solutions of (1.1).
(1) When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 1 or 𝐴 5 , the functions are 𝑒 1 1 ( πœ‰ ) = 6 𝑐 √ π‘Ž βˆ’ π‘Ž ξ‚€ √ ( π‘Ž + 6 𝑏 𝑐 ) c o s  , 𝑒 βˆ’ 𝑐 πœ‰ 1 2 ( πœ‰ ) = 6 𝑐 √ π‘Ž + ξ‚€ √ π‘Ž ( π‘Ž + 6 𝑏 𝑐 ) s i n  , 𝑒 βˆ’ 𝑐 πœ‰ 1 3 ( πœ‰ ) = 6 𝑐 √ π‘Ž βˆ’ ξ‚€ √ π‘Ž ( π‘Ž + 6 𝑏 𝑐 ) s i n  . βˆ’ 𝑐 πœ‰ ( 2 . 1 8 )
(2) When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 2 or when 𝑐 > 0 and ( π‘Ž , 𝑏 ) ∈ 𝐡 5 , the functions are 𝑒 2 1 ( πœ‰ ) = βˆ’ 𝛼 0 ξ€· 𝑀 c o s 0 πœ‰ ξ€Έ + 𝛽 0 βˆ’ 𝑝 0 ξ€· 𝑀 c o s 0 πœ‰ ξ€Έ + π‘ž 0 , 𝑒 2 2 𝛼 ( πœ‰ ) = 0 ξ€· 𝑀 s i n 0 πœ‰ ξ€Έ + 𝛽 0 𝑝 0 ξ€· 𝑀 s i n 0 πœ‰ ξ€Έ + π‘ž 0 , 𝑒 2 3 ( πœ‰ ) = βˆ’ 𝛼 0 ξ€· 𝑀 s i n 0 πœ‰ ξ€Έ + 𝛽 0 βˆ’ 𝑝 0 ξ€· 𝑀 s i n 0 πœ‰ ξ€Έ + π‘ž 0 . ( 2 . 1 9 )
(3) When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 6 , or when 𝑐 > 0 and ( π‘Ž , 𝑏 ) ∈ 𝐡 1 , the functions are 𝑒 3 1 ( πœ‰ ) = βˆ’ 𝛼 1 ξ€· 𝑀 c o s 1 πœ‰ ξ€Έ + 𝛽 1 βˆ’ 𝑝 1 ξ€· 𝑀 c o s 1 πœ‰ ξ€Έ + π‘ž 1 , 𝑒 3 2 𝛼 ( πœ‰ ) = 1 ξ€· 𝑀 s i n 1 πœ‰ ξ€Έ + 𝛽 1 𝑝 1 ξ€· 𝑀 s i n 1 πœ‰ ξ€Έ + π‘ž 1 , 𝑒 3 3 ( πœ‰ ) = βˆ’ 𝛼 1 ξ€· 𝑀 s i n 1 πœ‰ ξ€Έ + 𝛽 1 βˆ’ 𝑝 1 ξ€· 𝑀 s i n 1 πœ‰ ξ€Έ + π‘ž 1 . ( 2 . 2 0 )

Remark 2.3. In the given parametric regions, the solutions 𝑒 𝑖 ( πœ‰ ) , 𝑒 1 𝑖 ( πœ‰ ) , 𝑒 2 𝑖 ( πœ‰ ) , 𝑒 3 𝑖 ( πœ‰ ) ( 𝑖 = 1 , 2 , 3 ) , and 𝑒 ∘ 2 ( πœ‰ ) are nonsingular. The solutions 𝑒 ∘ 1 ( πœ‰ ) , 𝑒 βˆ— 1 ( πœ‰ ) , and 𝑒 βˆ— 2 ( πœ‰ ) are singular. The relationships of singular solutions and nonsingular solutions are displayed in the Proposition 2.1.

3. Trigonometric Function Periodic Wave Solutions for (1.2)

In this section, we state our main results for (1.2). For given π‘Ž 1 and 𝑐 ( π‘Ž 1 β‰  𝑐 2 ) , on π‘Ž 2 βˆ’ π‘Ž 3 plane we define some rays and regions as follows.( 1 )When 𝑐 2 < π‘Ž 1 , we define curves Ξ“ 1 : π‘Ž 2 > 0 , π‘Ž 3 Ξ“ = 0 , 2 : π‘Ž 2 > 0 , π‘Ž 3 = 2 π‘Ž 2 2 9 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , Ξ“ 3 : π‘Ž 2 > 0 , π‘Ž 3 = π‘Ž 2 2 4 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , Ξ“ 4 : π‘Ž 2 = 0 , π‘Ž 3 Ξ“ > 0 , 5 : π‘Ž 2 < 0 , π‘Ž 3 = π‘Ž 2 2 4 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , Ξ“ 6 : π‘Ž 2 < 0 , π‘Ž 3 = 2 π‘Ž 2 2 9 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , Ξ“ 7 : π‘Ž 2 < 0 , π‘Ž 3 Ξ“ = 0 , 8 : π‘Ž 2 = 0 , π‘Ž 3 < 0 , ( 3 . 1 ) and region π‘Š 𝑖 as the domain surrounded by Ξ“ 𝑖 and Ξ“ 𝑖 + 1 ( 𝑖 = 1 - 7 ) , π‘Š 8 as the domain surrounded by Ξ“ 8 and Ξ“ 1 (see Figure 7(a)).

fig7
Figure 7: The locations of the rays Ξ“ 𝑖 , 𝐿 𝑖 and the regions π‘Š 𝑖 , Ξ© 𝑖 ( 𝑖 = 1 , 2 , … , 8 ) for given π‘Ž 1 and 𝑐 .
( 2 )When 𝑐 2 > π‘Ž 1 , we define curves 𝐿 1 : π‘Ž 2 > 0 , π‘Ž 3 𝐿 = 0 , 2 : π‘Ž 2 = 0 , π‘Ž 3 𝐿 > 0 , 3 : π‘Ž 2 < 0 , π‘Ž 3 𝐿 = 0 , 4 : π‘Ž 2 < 0 , π‘Ž 3 = 2 π‘Ž 2 2 9 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , 𝐿 5 : π‘Ž 2 < 0 , π‘Ž 3 = π‘Ž 2 1 4 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , 𝐿 6 : π‘Ž 2 = 0 , π‘Ž 3 𝐿 < 0 , 7 : π‘Ž 2 > 0 , π‘Ž 3 = π‘Ž 2 2 4 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , 𝐿 8 : π‘Ž 2 > 0 , π‘Ž 3 = 2 π‘Ž 2 2 9 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ , ( 3 . 2 ) and region Ξ© 𝑖 as the domain surrounded by 𝐿 𝑖 and 𝐿 𝑖 + 1 ( 𝑖 = 1 - 7 ) , Ξ© 8 as the domain surrounded by 𝐿 8 and 𝐿 1 (see Figure 7(b)).

Using the rays and regions above, we state our results as follows.

Proposition 3.1. For given parameter π‘Ž 1 and constant 𝑐 satisfying 𝑐 2 β‰  π‘Ž 1 , let πœ‰ = π‘₯ βˆ’ 𝑐 𝑑 . Then, (1.2) has the following periodic wave solutions. ( 1 )When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 1 or π‘Š 6 , the expression of the periodic wave solution is 𝑣 1 𝑅 ( πœ‰ ) = 0 𝑅 1 + 𝑅 2 ξ€· 𝑅 c o s 3 πœ‰ ξ€Έ , ( 3 . 3 ) where 𝑅 0 ξ€· 𝑐 = 2 2 βˆ’ π‘Ž 1 ξ€Έ , 𝑅 1 = 2 π‘Ž 2 3 , 𝑅 2 = 1 3  1 8 π‘Ž 3 ξ€· 𝑐 2 βˆ’ π‘Ž 1 ξ€Έ + 4 π‘Ž 2 2 , 𝑅 3 = ξƒŽ π‘Ž 1 βˆ’ 𝑐 2 𝑐 2 . ( 3 . 4 ) For π‘Ž 2 β‰  0 , the periodic wave solution 𝑣 1 ( πœ‰ ) has the following limit forms.( i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 1 and ( π‘Ž 2 , π‘Ž 3 ) tends to the ray Ξ“ 1 , v 1 ( πœ‰ ) tends to the periodic blow-up solution 𝑣 ∘ 1 3 ξ€· 𝑐 ( πœ‰ ) = 2 βˆ’ π‘Ž 1 ξ€Έ π‘Ž 2 ξ‚€ √ 1 + c o s ξ‚€ ξ‚€ π‘Ž 1 βˆ’ 𝑐 2  πœ‰ . / | 𝑐 |   ( 3 . 5 ) The limiting process is similar to that in Figure 2.( i i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 6 and ( π‘Ž 2 , π‘Ž 3 ) tends to the ray Ξ“ 7 , 𝑣 1 ( πœ‰ ) tends to the periodic blow-up solution 𝑣 βˆ— 1 3 ξ€· 𝑐 ( πœ‰ ) = 2 βˆ’ π‘Ž 1 ξ€Έ π‘Ž 2 ξ‚€ √ 1 βˆ’ c o s ξ‚€ ξ‚€ π‘Ž 1 βˆ’ 𝑐 2  πœ‰ . / | 𝑐 |   ( 3 . 6 ) The limiting process is similar to that in Figure 3.( i i i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 1 and ( π‘Ž 2 , π‘Ž 3 ) tends to the curve Ξ“ 2 , or ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 6 and ( π‘Ž 2 , π‘Ž 3 ) tends to the curve Ξ“ 6 , 𝑣 1 ( πœ‰ ) tends to the trivial solution 𝑣 ( πœ‰ ) = 3 ( 𝑐 2 βˆ’ π‘Ž 1 ) / π‘Ž 2 .(2)When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 5 , or when 𝑐 2 > π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 1 , the expression of the periodic wave solution is 𝑣 2 ( πœ‰ ) = 𝑆 0 βˆ’ 2 𝑆 1 βˆ’ 𝑆 2 + 𝑆 3 ξ€· 𝑆 c o s 4 πœ‰ ξ€Έ , ( 3 . 7 ) where 𝑆 0 = βˆ’ π‘Ž 2 + √ πœ” 2 π‘Ž 3 , 𝑆 1 = βˆ’ π‘Ž 2 2 + 4 π‘Ž 3 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ + π‘Ž 2 √ πœ” π‘Ž 2 3 , 𝑆 2 = 2 3 π‘Ž 3 ξ‚€ βˆ’ π‘Ž 2 √ + 3 πœ”  , 𝑆 3 = 2 3 π‘Ž 3 ξ‚™ π‘Ž 2 ξ‚€ π‘Ž 2 √ + 3 πœ”  , 𝑆 4 = ξ‚™ βˆ’ 𝑆 1 π‘Ž 3 2 𝑐 2 , ( 3 . 8 ) π‘Ž πœ” = 2 2 βˆ’ 4 π‘Ž 3 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ . ( 3 . 9 ) The periodic wave solution 𝑣 2 ( πœ‰ ) has the following limit forms. ( i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 5 , and ( π‘Ž 2 , π‘Ž 3 ) tends to the curve Ξ“ 6 , 𝑣 2 ( πœ‰ ) tends to the trivial solution 𝑣 ( πœ‰ ) = 0 .( i i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 5 , and ( π‘Ž 2 , π‘Ž 3 ) tends to the curve Ξ“ 5 , the 𝑣 2 ( πœ‰ ) tends to the canyon-shaped solitary wave solution 𝑣 ∘ 2 2 ξ€· π‘Ž ( πœ‰ ) = 1 βˆ’ 𝑐 2 ξ€Έ ξ€Ί 1 2 𝑐 2 βˆ’ 9 𝑐 2 ξ€· π‘Ž βˆ’ 2 1 βˆ’ 𝑐 2 ξ€Έ πœ‰ 2 ξ€» π‘Ž 2 ξ€Ί 9 𝑐 2 ξ€· π‘Ž + 2 1 βˆ’ 𝑐 2 ξ€Έ πœ‰ 2 ξ€» . ( 3 . 1 0 ) The limiting process is similar to that in Figure 6.( i i i )When 𝑐 2 > π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 1 , and ( π‘Ž 2 , π‘Ž 3 ) tends to the ray 𝐿 1 , 𝑣 2 ( πœ‰ ) tends to the periodic blow-up wave solution 𝑣 βˆ— 2 π‘Ž ( πœ‰ ) = 1 βˆ’ 𝑐 2 2 π‘Ž 2 ⎑ ⎒ ⎒ ⎣ 1 + 3 t a n 2 βŽ› ⎜ ⎜ ⎝ ξƒŽ 𝑐 2 βˆ’ π‘Ž 1 4 𝑐 2 πœ‰ ⎞ ⎟ ⎟ ⎠ ⎀ βŽ₯ βŽ₯ ⎦ . ( 3 . 1 1 ) The limiting process is similar to that in Figure 2.(3)When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 2 , or when 𝑐 2 > π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 2 , the expression of the periodic wave solution is 𝑣 3 ( πœ‰ ) = 𝑇 0 + 2 𝑇 1 βˆ’ 𝑇 2 + 𝑇 3 ξ€· 𝑇 c o s 4 πœ‰ ξ€Έ , ( 3 . 1 2 ) where 𝑇 0 = βˆ’ π‘Ž 2 βˆ’ √ πœ” 2 π‘Ž 3 , 𝑇 1 = βˆ’ π‘Ž 2 2 + 4 π‘Ž 3 ξ€· π‘Ž 1 βˆ’ 𝑐 2 ξ€Έ βˆ’ π‘Ž 2 √ πœ” π‘Ž 2 3 , 𝑇 2 = 2 3 π‘Ž 3 ξ‚€ π‘Ž 2 √ + 3 πœ”  , 𝑇 3 = 2 3 π‘Ž 3 ξ‚™ π‘Ž 2 ξ‚€ π‘Ž 2 √ βˆ’ 3 πœ”  , 𝑇 4 = ξ‚™ βˆ’ 𝑇 1 π‘Ž 3 2 𝑐 2 . ( 3 . 1 3 ) The periodic wave solution 𝑣 3 ( πœ‰ ) has the following limit forms:( i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 2 , and ( π‘Ž 2 , π‘Ž 3 ) tends to the curve Ξ“ 2 , 𝑣 3 ( πœ‰ ) tends to the trivial solution 𝑣 ( πœ‰ ) = 0 .( i i )When 𝑐 2 < π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 2 , and ( π‘Ž 2 , π‘Ž 3 ) tends to the curve Ξ“ 3 , the 𝑣 3 ( πœ‰ ) tends to the peak-shaped solitary wave solution 𝑣 ∘ 2 ( πœ‰ ) . The limiting process is similar to that in Figure 4.( i i i )When 𝑐 2 > π‘Ž 1 , ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 2 , and ( π‘Ž 2 , π‘Ž 3 ) tends to the ray 𝐿 3 , the 𝑣 3 ( πœ‰ ) tends to the periodic blow-up wave solution 𝑣 βˆ— 2 ( πœ‰ ) . The limiting process is similar to that in Figure 5.

Remark 3.2. Similar to the reason in Remark 2.2, the following nine functions also are periodic wave solutions of (1.2).
( 1 )When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 1 or π‘Š 6 , the functions are 𝑣 1 1 𝑅 ( πœ‰ ) = 0 𝑅 1 βˆ’ 𝑅 2 ξ€· 𝑅 c o s 3 πœ‰ ξ€Έ , 𝑣 1 2 ( 𝑅 πœ‰ ) = 0 𝑅 1 + 𝑅 2 ξ€· 𝑅 s i n 3 πœ‰ ξ€Έ , 𝑣 1 3 𝑅 ( πœ‰ ) = 0 𝑅 1 βˆ’ 𝑅 2 ξ€· 𝑅 s i n 3 πœ‰ ξ€Έ . ( 3 . 1 4 ) ( 2 ) When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 5 or when 𝑐 2 > π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 1 , the functions are 𝑣 2 1 ( πœ‰ ) = S 0 + 2 𝑆 1 𝑆 2 + 𝑆 3 ξ€· 𝑆 c o s 4 πœ‰ ξ€Έ , 𝑣 2 2 ( πœ‰ ) = 𝑆 0 βˆ’ 2 𝑆 1 βˆ’ 𝑆 2 + 𝑆 3 ξ€· 𝑆 s i n 4 πœ‰ ξ€Έ , 𝑣 2 3 ( πœ‰ ) = 𝑆 0 + 2 𝑆 1 𝑆 2 + 𝑆 3 ξ€· 𝑆 s i n 4 πœ‰ ξ€Έ . ( 3 . 1 5 ) ( 3 ) When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 2 or when 𝑐 2 > π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 2 , the functions are 𝑣 3 1 ( πœ‰ ) = 𝑇 0 βˆ’ 2 𝑇 1 𝑇 2 + 𝑇 3 ξ€· 𝑇 c o s 4 πœ‰ ξ€Έ , 𝑣 3 2 ( πœ‰ ) = 𝑇 0 + 2 𝑇 1 βˆ’ 𝑇 2 + 𝑇 3 ξ€· 𝑇 s i n 4 πœ‰ ξ€Έ , 𝑣 3 3 ( πœ‰ ) = 𝑇 0 βˆ’ 2 𝑇 1 𝑇 2 + 𝑇 3 ξ€· 𝑇 s i n 4 πœ‰ ξ€Έ . ( 3 . 1 6 )

Remark 3.3. In the given regions, the solutions 𝑣 𝑖 ( πœ‰ ) , 𝑣 1 𝑖 ( πœ‰ ) , 𝑣 2 𝑖 ( πœ‰ ) , 𝑣 3 𝑖 ( πœ‰ ) ( 𝑖 = 1 , 2 , 3 ) , and 𝑣 ∘ 2 ( πœ‰ ) are nonsingular. The solutions 𝑣 ∘ 1 ( πœ‰ ) , 𝑣 βˆ— 1 ( πœ‰ ) , and 𝑣 βˆ— 2 ( πœ‰ ) are singular. The relationships of nonsingular solutions and singular solutions are displayed in Proposition 3.1.

4. The Derivation on Proposition 2.1

In order to derive the Proposition 2.1, letting 𝑐 be a constant and substituting 𝑒 = πœ‘ ( πœ‰ ) with πœ‰ = π‘₯ βˆ’ 𝑐 𝑑 into (1.1), we have βˆ’ 𝑐 πœ‘ ξ…ž + π‘Ž πœ‘ πœ‘ ξ…ž + π‘Ž 𝑏 πœ‘ 2 πœ‘ ξ…ž + πœ‘ ξ…ž ξ…ž ξ…ž = 0 . ( 4 . 1 )

Integrating (4.1) once and letting the integral constant be zero, it follows that π‘Ž βˆ’ 𝑐 πœ‘ + 2 πœ‘ 2 + π‘Ž 𝑏 3 πœ‘ 3 + πœ‘ ξ…ž ξ…ž = 0 . ( 4 . 2 )

Letting πœ‘ ξ…ž = 𝑦 , yields the following planar system: πœ‘ ξ…ž = 𝑦 , 𝑦 ξ…ž π‘Ž = 𝑐 πœ‘ βˆ’ 2 βˆ’ π‘Ž 𝑏 3 πœ‘ 3 . ( 4 . 3 )

Obviously, system (4.3) has the first integral 6 𝑦 2 βˆ’ 6 𝑐 πœ‘ 2 + 2 π‘Ž πœ‘ 3 + π‘Ž 𝑏 πœ‘ 4 = β„Ž . ( 4 . 4 )

Let πœ‘ 1 = √ βˆ’ 3 π‘Ž βˆ’ Ξ” , πœ‘ 4 π‘Ž 𝑏 2 = √ βˆ’ 3 π‘Ž + Ξ” , 4 π‘Ž 𝑏 ( 4 . 5 ) where Ξ” is defined in (2.8). Then, it is easy to see that system (4.3) has three singular points ( πœ‘ 1 , 0 ) , ( 0 , 0 ) and ( πœ‘ 2 , 0 ) when Ξ” > 0 , two singular points ( ( βˆ’ 3 / 4 𝑏 ) , 0 ) and ( 0 , 0 ) when Ξ” = 0 , unique singular point ( 0 , 0 ) when Ξ” < 0 .

Let 𝑒 𝑖 and 𝑓 𝑖 ( 𝑖 = 1 , 2 , 3 ) be, respectively, 𝑒 1 = √ βˆ’ π‘Ž βˆ’ π‘Ž 2 + 6 π‘Ž 𝑏 𝑐 , 𝑓 π‘Ž 𝑏 1 = √ βˆ’ π‘Ž + π‘Ž 2 + 6 π‘Ž 𝑏 𝑐 , 𝑒 π‘Ž 𝑏 2 = 1  √ 4 π‘Ž 𝑏 βˆ’ π‘Ž + ξ‚™ Ξ” βˆ’ 2 π‘Ž ξ‚€ √ π‘Ž βˆ’ Ξ”  ξƒͺ , 𝑓 2 = 1  √ 4 π‘Ž 𝑏 βˆ’ π‘Ž + ξ‚™ Ξ” + 2 π‘Ž ξ‚€ √ π‘Ž βˆ’ Ξ”  ξƒͺ , 𝑒 3 1 = βˆ’  √ 4 π‘Ž 𝑏 π‘Ž + ξ‚™ Ξ” + 2 π‘Ž ξ‚€ √ π‘Ž + Ξ”  ξƒͺ , 𝑓 3 1 = βˆ’  √ 4 π‘Ž 𝑏 π‘Ž + ξ‚™ Ξ” βˆ’ 2 π‘Ž ξ‚€ √ π‘Ž + Ξ”  ξƒͺ . ( 4 . 6 )

Using the qualitative analysis of dynamical systems, we obtain the bifurcation phase portraits of system (4.3) and the locations of 𝑒 𝑖 and 𝑓 𝑖 ( 𝑖 = 1 , 2 , 3 ) as Figures 8 and 9.

810217.fig.008
Figure 8: When 𝑐 < 0 , the bifurcation phase portraits of system (4.3) and the locations of 𝑒 𝑖 and 𝑓 𝑖 ( 𝑖 = 1 , 2 , 3 ) .
810217.fig.009
Figure 9: When 𝑐 > 0 , the bifurcation phase portraits of system (4.3) and the locations of 𝑒 𝑖 and 𝑓 𝑖 ( 𝑖 = 1 , 2 , 3 ) .

It is easy to test that the closed orbit passing ( 𝑒 𝑖 , 0 ) passes ( 𝑓 𝑖 , 0 ) ( 𝑖 = 1 , 2 , 3 ) . Thus, using the phase portraits in Figures 8 and 9, we derive 𝑒 𝑖 ( πœ‰ ) ( 𝑖 = 1 , 2 , 3 ) as follows.

( 1 )When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 1 or 𝐴 5 , the closed orbit passing the points ( 𝑒 1 , 0 ) and ( 𝑓 1 , 0 ) has expression ξ‚™ 𝑦 = Β± π‘Ž 𝑏 6 πœ‘  βˆ’ 𝑒 1 𝑓 1 + ξ€· 𝑒 1 + 𝑓 1 ξ€Έ πœ‘ βˆ’ πœ‘ 2 , w h e r e 𝑒 1 ≀ πœ‘ ≀ 𝑓 1 . ( 4 . 7 ) Substituting (4.7) into d πœ‘ / 𝑦 = d πœ‰ , we have d πœ‘  βˆ’ 𝑒 1 𝑓 1 + ξ€· 𝑒 1 + 𝑓 1 ξ€Έ πœ‘ βˆ’ πœ‘ 2 = ξ‚™ π‘Ž 𝑏 6 d πœ‰ . ( 4 . 8 ) Integrating (4.8) along the closed orbit and noting that 𝑒 = πœ‘ ( πœ‰ ) , we obtain the solution 𝑒 1 ( πœ‰ ) as (2.4).( 2 ) When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 2 or when 𝑐 > 0 and ( π‘Ž , 𝑏 ) ∈ 𝐡 5 , the closed orbit passing the points ( 𝑒 2 , 0 ) and ( 𝑓 2 , 0 ) has expression ξ‚™ 𝑦 = Β± π‘Ž 𝑏 6 ξ€· πœ‘ βˆ’ πœ‘ 1 ξ€Έ  βˆ’ 𝑒 2 𝑓 2 + ξ€· 𝑒 2 + 𝑓 2 ξ€Έ πœ‘ βˆ’ πœ‘ 2 , w h e r e 𝑒 2 ≀ πœ‘ ≀ 𝑓 2 . ( 4 . 9 ) Substituting (4.9) into d πœ‘ / 𝑦 = d πœ‰ , we get d πœ‘ ξ€· πœ‘ βˆ’ πœ‘ 1 ξ€Έ  βˆ’ 𝑒 2 𝑓 2 + ξ€· 𝑒 2 + 𝑓 2 ξ€Έ πœ‘ βˆ’ πœ‘ 2 = ξ‚™ π‘Ž 𝑏 6 d πœ‰ . ( 4 . 1 0 ) Along the closed orbit integrating (4.10) and noting that 𝑒 = πœ‘ ( πœ‰ ) , we get the solution 𝑒 2 ( πœ‰ ) as (2.7).( 3 ) When 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 6 or when 𝑐 > 0 and ( π‘Ž , 𝑏 ) ∈ 𝐡 1 , the closed orbit passing the points ( 𝑒 3 , 0 ) and ( 𝑓 3 , 0 ) has expression ξ‚™ 𝑦 = Β± π‘Ž 𝑏 6 ξ€· πœ‘ 2 ξ€Έ  βˆ’ πœ‘ βˆ’ 𝑒 3 𝑓 3 + ξ€· 𝑒 3 + 𝑓 3 ξ€Έ πœ‘ βˆ’ πœ‘ 2 , w h e r e 𝑒 3 ≀ πœ‘ ≀ 𝑓 3 . ( 4 . 1 1 ) Substituting (4.11) into d πœ‘ / 𝑦 = d πœ‰ , it follows that d πœ‘ ξ€· πœ‘ 2 ξ€Έ  βˆ’ πœ‘ βˆ’ 𝑒 3 𝑓 3 + ξ€· 𝑒 3 + 𝑓 3 ξ€Έ πœ‘ βˆ’ πœ‘ 2 = ξ‚™ π‘Ž 𝑏 6 d πœ‰ . ( 4 . 1 2 ) Similarly, along the closed orbit integrating (4.12), we obtain 𝑒 3 ( πœ‰ ) as (2.16). From the expressions of these solutions, we get their limit forms. This completes the derivation on Proposition 2.1.

5. The Derivation on Proposition 3.1

In this section, we give derivation on Proposition 3.1. Let 𝑣 = πœ“ ( πœ‰ ) with πœ‰ = π‘₯ βˆ’ 𝑐 𝑑 , where 𝑐 is a constant. Thus, (1.2) becomes 𝑐 2 πœ“ ξ…ž ξ…ž βˆ’ 𝑐 2 πœ“ ξ…ž ξ…ž ξ…ž ξ…ž βˆ’ ξ€· π‘Ž 1 πœ“ + π‘Ž 2 πœ“ 2 + π‘Ž 3 πœ“ 3 ξ€Έ ξ…ž ξ…ž = 0 . ( 5 . 1 )

Integrating (5.1) twice and letting integral constant be zero, we get 𝑐 2 ξ€· πœ“ βˆ’ πœ“ ξ…ž ξ…ž ξ€Έ = π‘Ž 1 πœ“ + π‘Ž 2 πœ“ 2 + π‘Ž 3 πœ“ 3 . ( 5 . 2 )

Letting πœ“ ξ…ž = 𝑦 , we have the planar system πœ“ ξ…ž = 𝑦 , 𝑐 2 𝑦 ξ…ž = ξ€· 𝑐 2 βˆ’ π‘Ž 1 ξ€Έ πœ“ βˆ’ π‘Ž 2 πœ“ 2 βˆ’ π‘Ž 3 πœ“ 3 . ( 5 . 3 )

It is easy to see that system (5.3) has the first integral 𝑐 2 𝑦 2 + πœ“ 2 ξ‚΅ π‘Ž 3 2 πœ“ 2 + 2 π‘Ž 2 3 πœ“ + π‘Ž 1 βˆ’ 𝑐 2 ξ‚Ά = β„Ž , ( 5 . 4 ) and three singular points ( 0 , 0 ) , ( πœ“ 1 , 0 ) , and ( πœ“ 2 , 0 ) , where πœ“ 1 = βˆ’ π‘Ž 2 βˆ’ √ πœ” 2 π‘Ž 3 , πœ“ 2 = βˆ’ π‘Ž 2 + √ πœ” 2 π‘Ž 3 ( 5 . 5 ) and πœ” is defined in (3.9).

Let π‘š 𝑖 and 𝑛 𝑖 ( 𝑖 = 1 , 2 , 3 ) be, respectively, π‘š 1 = βˆ’ 2 π‘Ž 2 βˆ’  2 ξ€· π‘Ž 2 2 βˆ’ 9 π‘Ž 1 π‘Ž 3 + 9 π‘Ž 3 𝑐 2 ξ€Έ 3 π‘Ž 3 , 𝑛 1 = βˆ’ 2 π‘Ž 2 +  2 ξ€· π‘Ž 2 2 βˆ’ 9 π‘Ž 1 π‘Ž 3 + 9 π‘Ž 3 𝑐 2 ξ€Έ 3 π‘Ž 3 , π‘š 2 = βˆ’ π‘Ž 2 √ βˆ’ 3 ξ‚™ πœ” βˆ’ 2 π‘Ž 2 ξ‚€ π‘Ž 2 √ + 3 πœ”  6 π‘Ž 3 , 𝑛 2 = βˆ’ π‘Ž 2 √ βˆ’ 3 ξ‚™ πœ” + 2 π‘Ž 2 ξ‚€ π‘Ž 2 √ + 3 πœ”  6 π‘Ž 3 , π‘š 3 = βˆ’ π‘Ž 2 √ + 3 ξ‚™ πœ” βˆ’ 2 π‘Ž 2 ξ‚€ π‘Ž 2 √ βˆ’ 3 πœ”  6 π‘Ž 3 , 𝑛 3 = βˆ’ π‘Ž 2 √ + 3 ξ‚™ πœ” + 2 π‘Ž 2 ξ‚€ π‘Ž 2 √ βˆ’ 3 πœ”  6 π‘Ž 3 . ( 5 . 6 )

Similarly, using the qualitative analysis of dynamical systems, we get the bifurcation phase portraits of system (5.3) and the locations of π‘š 𝑖 and 𝑛 𝑖 ( 𝑖 = 1 , 2 , 3 ) as Figures 10 and 11.

810217.fig.0010
Figure 10: When 𝑐 2 < π‘Ž 1 , the bifurcation phase portraits of system (5.3) and the locations of π‘š 𝑖 and 𝑛 𝑖 ( 𝑖 = 1 , 2 , 3 ) .
810217.fig.0011
Figure 11: When 𝑐 2 > π‘Ž 1 , the bifurcation phase portraits of system (5.3) and the locations of π‘š 𝑖 and 𝑛 𝑖 ( 𝑖 = 1 , 2 , 3 ) .

It is easy to test that the closed orbit passing ( π‘š 𝑖 , 0 ) passes ( 𝑛 𝑖 , 0 ) ( 𝑖 = 1 , 2 , 3 ) . Thus, using the phase portraits in Figures 10 and 11, we derive 𝑣 𝑖 ( πœ‰ ) ( 𝑖 = 1 , 2 , 3 ) as follows.( 1 )When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 1 or π‘Š 6 , the closed orbit passing the points ( π‘š 1 , 0 ) and ( 𝑛 1 , 0 ) has expression ξ‚™ 𝑦 = Β± π‘Ž 3 2 𝑐 2 πœ“  βˆ’ π‘š 1 𝑛 1 + ξ€· π‘š 1 + 𝑛 1 ξ€Έ πœ“ βˆ’ πœ“ 2 , w h e r e π‘š 1 ≀ πœ“ ≀ 𝑛 1 . ( 5 . 7 ) Substituting (5.7) into d πœ“ / 𝑦 = d πœ‰ , we have d πœ“ πœ“  βˆ’ π‘š 1 𝑛 1 + ξ€· π‘š 1 + 𝑛 1 ξ€Έ πœ“ βˆ’ πœ“ 2 = ξ‚™ π‘Ž 3 2 𝑐 2 d πœ‰ . ( 5 . 8 ) Integrating (5.8) along the closed orbit and noting that 𝑣 = πœ“ ( πœ‰ ) , we get the solution 𝑣 1 ( πœ‰ ) as (3.3). ( 2 ) When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 5 , or when 𝑐 2 > π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 1 , the closed orbit passing the points ( π‘š 2 , 0 ) and ( 𝑛 2 , 0 ) has expression ξ‚™ 𝑦 = Β± π‘Ž 3 2 𝑐 2 ξ€· πœ“ 2 ξ€Έ  βˆ’ πœ“ βˆ’ π‘š 2 𝑛 2 + ξ€· π‘š 2 + 𝑛 2 ξ€Έ πœ“ βˆ’ πœ“ 2 , w h e r e π‘š 2 ≀ πœ“ ≀ 𝑛 2 . ( 5 . 9 ) From d πœ“ / 𝑦 = d πœ‰ and (5.9), it follows that d πœ“ ξ€· πœ“ 2 ξ€Έ  βˆ’ πœ“ βˆ’ π‘š 2 𝑛 2 + ξ€· π‘š 2 + 𝑛 2 ξ€Έ πœ“ βˆ’ πœ“ 2 = ξ‚™ π‘Ž 3 2 𝑐 2 d πœ‰ . ( 5 . 1 0 ) Integrating (5.10) along the closed orbit, we get 𝑣 2 ( πœ‰ ) as (3.7).( 3 ) When 𝑐 2 < π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ π‘Š 2 , or when 𝑐 2 > π‘Ž 1 and ( π‘Ž 2 , π‘Ž 3 ) ∈ Ξ© 2 , the closed orbit passing the points ( π‘š 3 , 0 ) and ( 𝑛 3 , 0 ) has expression ξ‚™ 𝑦 = Β± π‘Ž 3 2 𝑐 2 ξ€· πœ“ βˆ’ πœ“ 1 ξ€Έ  βˆ’ π‘š 3 𝑛 3 + ξ€· π‘š 3 + 𝑛 3 ξ€Έ πœ“ βˆ’ πœ“ 2 , w h e r e π‘š 3 ≀ πœ“ ≀ 𝑛 3 . ( 5 . 1 1 ) Substituting (5.11) into d πœ“ / 𝑦 = d πœ‰ , we have d πœ“ ξ€· πœ“ βˆ’ πœ“ 1 ξ€Έ  βˆ’ π‘š 3 𝑛 3 + ξ€· π‘š 3 + 𝑛 3 ξ€Έ πœ“ βˆ’ πœ“ 2 = ξ‚™ π‘Ž 3 2 𝑐 2 d πœ‰ . ( 5 . 1 2 ) Integrating (5.12) along the closed orbit, we obtain 𝑣 3 ( πœ‰ ) as (3.12). From the expressions of these solutions, we get their limiting properties. This completes the derivation on Proposition 3.1.

6. Discussions and Testing Orders

In this paper, Using the special closed orbits, we have obtained trigonometric function periodic wave solutions for (1.1) and (1.2), respectively. Their limit forms have been given. From these expressions, an interesting phenomena has been seen, that is, (1.1) and (1.2) have similar periodic wave solutions. Our work has extended previous results on periodic wave solutions.

Now, we point out that the trigonometric function periodic wave solutions can be obtained from the limits of the elliplic function periodic wave solution. For given real number πœ‡ , let πœ‡ 1 = 1 βŽ› ⎜ ⎜ ⎝ 4 ξ‚€ √ 1 2 π‘Ž 𝑏 βˆ’ 4 π‘Ž ( 2 + 𝑏 πœ‡ ) + 1 + i 3  π‘Ž 𝐹 0 2 𝐹 ξ‚€ √ + 2 i i + 3  𝐹 ⎞ ⎟ ⎟ ⎠ , πœ‡ 2 = 1 βŽ› ⎜ ⎜ ⎝ 4 ξ‚€ √ 1 2 π‘Ž 𝑏 βˆ’ 4 π‘Ž ( 2 + 𝑏 πœ‡ ) + 1 βˆ’ i 3  π‘Ž 𝐹 0 2 𝐹 ξ‚€ √ βˆ’ 2 i i + 3  𝐹 ⎞ ⎟ ⎟ ⎠ , πœ‡ 3 = 1 ξ‚΅ 6 π‘Ž 𝑏 βˆ’ 2 π‘Ž ( 2 + 𝑏 πœ‡ ) βˆ’ 4 π‘Ž 𝐹 0 2 𝐹 ξ‚Ά , + 2 𝐹 ( 6 . 1 ) where 𝐹 0 1 = ξ€· 8 βˆ’ 6 𝑏 πœ‡ + 1 5 𝑏 2 πœ‡ 2 + 1 0 𝑏 3 πœ‡ 3 ξ€Έ , 𝐹 0 2 = ξ€· ξ€· βˆ’ 9 𝑏 𝑐 + π‘Ž βˆ’ 2 + 𝑏 πœ‡ + 𝑏 2 πœ‡ 2 , 𝐹 ξ€Έ ξ€Έ 0 3 = ξ‚™ π‘Ž 3 ξ‚€ 8 𝐹 3 0 2 ξ€· + π‘Ž βˆ’ 5 4 𝑏 𝑐 ( βˆ’ 1 + 𝑏 πœ‡ ) + π‘Ž 𝐹 0 1 ξ€Έ 2  , ξ€· 𝐹 = 5 4 π‘Ž 2 𝑏 𝑐 ( βˆ’ 1 + 𝑏 πœ‡ ) βˆ’ π‘Ž 3 𝐹 0 1 + 𝐹 0 3 ξ€Έ 1 / 3 . ( 6 . 2 )

Assume that 𝑐 < 0 , ( π‘Ž , 𝑏 ) ∈ ( 𝐴 1 ) , and πœ‘ 1 < πœ‡ < 𝑒 1 . It is easy to check that πœ‡ 𝑖 ( 𝑖 = 1 , 2 , 3 ) are real and satisfy πœ‡ < 𝑒 1 < πœ‘ 2 < 𝑓 1 < πœ‡ 1 < πœ‘ 3 < πœ‡ 2 < 0 < πœ‡ 3 < πœ‘ 4 . ( 6 . 3 )

There are two closed orbits 𝑙 1 πœ‡ and 𝑙 2 πœ‡ (see Figure 12). The closed orbit 𝑙 1 πœ‡ passes the points ( πœ‡ , 0 ) and ( πœ‡ 1 , 0 ) . The closed orbit 𝑙 2 πœ‡ passes the points ( πœ‡ 2 , 0 ) and ( πœ‡ 3 , 0 ) .

810217.fig.0012
Figure 12: The locations of 𝑙 1 πœ‡ and 𝑙 2 πœ‡ when 𝑐 < 0 and ( π‘Ž , 𝑏 ) ∈ 𝐴 1 .

On πœ‘ βˆ’ 𝑦 plane, the expression of 𝑙 1 πœ‡ is 𝑦 2 = π‘Ž 𝑏 6 ξ€· πœ‡ 3 πœ‡ βˆ’ πœ‘ ξ€Έ ξ€· 2 πœ‡ βˆ’ πœ‘ ξ€Έ ξ€· 1 ξ€Έ βˆ’ πœ‘ ( πœ‘ βˆ’ πœ‡ ) , w h e r e πœ‡ ≀ πœ‘ ≀ πœ‡ 1 . ( 6 . 4 )

Substituting (6.4) into d πœ‘ / 𝑦 = d πœ‰ and integrating it along 𝑙 1 πœ‡ , we have 𝑔 s n βˆ’ 1 ξ‚™ ( s i n 𝑧 , π‘˜ ) = π‘Ž 𝑏 6 | | πœ‰ | | , ( 6 . 5 ) where 2 𝑔 =  ξ€· πœ‡ 3 βˆ’ πœ‡ 1 πœ‡ ξ€Έ ξ€· 2 ξ€Έ , ξ„Ά ξ„΅ ξ„΅ ⎷ βˆ’ πœ‡ π‘˜ = ξ€· πœ‡ 3 βˆ’ πœ‡ 2 πœ‡ ξ€Έ ξ€· 1 ξ€Έ βˆ’ πœ‡ ξ€· πœ‡ 3 βˆ’ πœ‡ 1 πœ‡ ξ€Έ ξ€· 2 ξ€Έ , ξ„Ά ξ„΅ ξ„΅ ⎷ βˆ’ πœ‡ s i n 𝑧 = ξ€· πœ‡ 3 βˆ’ πœ‡ 1 ξ€Έ ( πœ‘ βˆ’ πœ‡ ) ξ€· πœ‡ 1 πœ‡ βˆ’ πœ‡ ξ€Έ ξ€· 3 ξ€Έ . βˆ’ πœ‘ ( 6 . 6 )

Solving (6.5) for πœ‘ and noting that 𝑒 = πœ‘ ( πœ‰ ) , we obtain an elliptic function periodic wave solution 𝑒 πœ‡ ξ€· πœ‡ ( πœ‰ ) = 3 βˆ’ πœ‡ 1 ξ€Έ + πœ‡ 3 ξ€· πœ‡ 1 ξ€Έ βˆ’ πœ‡ s n 2 ( πœ‚ πœ‰ , π‘˜ ) πœ‡ 3 βˆ’ πœ‡ 1 + ξ€· πœ‡ 1 ξ€Έ βˆ’ πœ‡ s n 2 , ( πœ‚ πœ‰ , π‘˜ ) ( 6 . 7 ) where ξƒŽ πœ‚ = ξ€· πœ‡ π‘Ž 𝑏 3 βˆ’ πœ‡ 1 πœ‡ ξ€Έ ξ€· 2 ξ€Έ βˆ’ πœ‡ . 2 4 ( 6 . 8 )

Letting πœ‡ β†’ 𝑒 1 βˆ’ 0 , it follows that πœ‡ 1 β†’ 𝑓 1 , πœ‡ 2 β†’ 0 , πœ‡ 3 β†’ 0 , π‘˜ β†’ 0 , √ πœ‚ β†’ ( π‘Ž 𝑏 𝑒 1 𝑓 1 ) / 2 4 and s n 2 ( πœ‚ πœ‰ , π‘˜ ) β†’ s n 2 ( √ ( π‘Ž 𝑏 𝑒 1 𝑓 1 / 2 4 ) πœ‰ , 0 ) = s i n 2 ( √ ( π‘Ž 𝑏 𝑒 1 𝑓 1 / 2 4 ) πœ‰ ) .

Therefore, in (6.7) letting πœ‡ β†’ 𝑒 1 βˆ’ 0 , we obtain the trigonometric function periodic wave solution 𝑒 𝑒 ( πœ‰ ) = 1 𝑓 1 𝑓 1 + ξ€· 𝑒 1 βˆ’ 𝑓 1 ξ€Έ s i n 2 ξ‚΅  ξ€· π‘Ž 𝑏 𝑒 1 𝑓 1 ξ€Έ πœ‰ ξ‚Ά = / 2 4 βˆ’ 6 𝑐 √ βˆ’ π‘Ž + √ π‘Ž ( π‘Ž + 6 𝑏 𝑐 ) βˆ’ 2 π‘Ž ( π‘Ž + 6 𝑏 𝑐 ) s i n 2 √ ξ‚€ ξ‚€  πœ‰  = | 𝑐 | / 2 6 𝑐 √ π‘Ž βˆ’ π‘Ž ξ‚€ √ ( π‘Ž + 6 𝑏 𝑐 ) c o s  | 𝑐 | πœ‰ = 𝑒 1 1 ( πœ‰ ) . ( 6 . 9 )

Via Remark 2.2 and 𝑒 1 1 ( πœ‰ ) , further we get 𝑒 1 2 ( πœ‰ ) , 𝑒 1 3 ( πœ‰ ) and 𝑒 1 ( πœ‰ ) . Similarly, we can derive others trigonometric function periodic wave solutions.

We also have tested the correctness of each solution by using the software Mathematica. Here, we list two testing orders. Others testing orders are similar.( 1 )The orders for testing 𝑒 1 ( πœ‰ ) u = 6 c √ a + a  √ ( a + 6 b c ) c o s ξ‚„ βˆ’ c ( x βˆ’ c t ) ( 6 . 1 0 ) Simplify [ D [ u , t ] + a ( 1 + b u ) D [ u , x ] u + D [ u , { x , 3 } ] ] . ( 2 ) The orders for testing 𝑣 1 ( πœ‰ ) R 0 ξ€· = 2 βˆ’ a 1 + c 2 ξ€Έ , R 1 = 2 a 2 3 , R 2 = ξƒŽ 2 a 3 ξ€· βˆ’ a 1 + c 2 ξ€Έ + 4 a 2 2 9 , R 3 = ξƒŽ a 1 βˆ’ c 2 c 2 , R v = 0 R 1 + R 2 ξ€Ί R c o s 3 ξ€» , D [ ] , [ ] ( x βˆ’ c t ) v t t = v , { t , 2 } v t t x x = D v t t , { x , 2 } ( 6 . 1 1 ) Simplify [ v t t βˆ’ v t t x x βˆ’ D [ a 1 v + a 2 v 2 + a 3 v 3 , { x , 2 } ] ] .

Acknowledgment

Research is supported by the National Natural Science Foundation of China (no. 10871073) and the Research Expences of Central Universities for students.

References

  1. B. Dey, β€œDomain wall solutions of KdV-like equations with higher order nonlinearity,” Journal of Physics A, vol. 19, no. 1, pp. L9–L12, 1986. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  2. B. Dey, β€œKdV like equations with domain wall solutions and their Hamiltonians,” in Solitons, pp. 188–194, Springer, Berlin, Germany, 1988.
  3. J. F. Zhang, β€œNew solitary wave solution of the combined KdV and mKdV equation,” International Journal of Theoretical Physics, vol. 37, no. 5, pp. 1541–1546, 1998. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  4. J. F. Zhang, F. M. Wu, and J. Q. Shi, β€œSimple solition solution method for the combined KdV and mKdV equation,” International Journal of Theoretical Physics, vol. 39, no. 6, pp. 1697–1702, 2000. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  5. J. Yu, β€œExact solitary wave solutions to a combined KdV and mKdV equation,” Mathematical Methods in the Applied Sciences, vol. 23, no. 18, pp. 1667–1670, 2000. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  6. R. Grimshaw, D. Pelinovsky, E. Pelinovsky, and A. Slunyaev, β€œGeneration of large-amplitude solitons in the extended Korteweg-de Vries equation,” Chaos, vol. 12, no. 4, pp. 1070–1076, 2002. View at Publisher Β· View at Google Scholar Β· View at PubMed Β· View at Zentralblatt MATH
  7. E. G. Fan, β€œMultiple travelling wave solutions of nonlinear evolution equations using a unified algebraic method,” Journal of Physics A, vol. 35, no. 32, pp. 6853–6872, 2002. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  8. E. G. Fan, β€œUniformly constructing a series of explicit exact solutions to nonlinear equations in mathematical physics,” Chaos, Solitons and Fractals, vol. 16, no. 5, pp. 819–839, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  9. M. Y. Tang, R. Q. Wang, and Z. J. Jing, β€œSolitary waves and their bifurcations of KdV like equation with higher order nonlinearity,” Science in China A, vol. 45, no. 10, pp. 1255–1267, 2002. View at Zentralblatt MATH
  10. Y.-Z. Peng, β€œNew exact solutions to the combined KdV and mKdV equation,” International Journal of Theoretical Physics, vol. 42, no. 4, pp. 863–868, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  11. K. W. Chow, R. H. J. Grimshaw, and E. Ding, β€œInteractions of breathers and solitons in the extended Korteweg-de Vries equation,” Wave Motion, vol. 43, no. 2, pp. 158–166, 2005. View at Publisher Β· View at Google Scholar Β· View at MathSciNet
  12. D. Kaya and I. E. Inan, β€œA numerical application of the decomposition method for the combined KdV—MKdV equation,” Applied Mathematics and Computation, vol. 168, no. 2, pp. 915–926, 2005. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  13. E. Yomba, β€œThe extended Fan's sub-equation method and its application to KdV-MKdV, BKK and variant Boussinesq equations,” Physics Letters A, vol. 336, no. 6, pp. 463–476, 2005. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  14. W. G. Zhang and W. X. Ma, β€œExplicit solitary-wave solutions to generalized Pochhammer-Chree equations,” Applied Mathematics and Mechanics, vol. 20, no. 6, pp. 625–632, 1999. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  15. J. B. Li and L. J. Zhang, β€œBifurcations of traveling wave solutions in generalized Pochhammer-Chree equation,” Chaos, Solitons and Fractals, vol. 14, no. 4, pp. 581–593, 2002. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  16. D. Kaya, β€œThe exact and numerical solitary-wave solutions for generalized modified Boussinesq equation,” Physics Letters A, vol. 348, no. 3–6, pp. 244–250, 2006. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  17. M. Rafei, D. D. Ganji, H. R. M. Daniali, and H. Pashaei, β€œApplication of homotopy perturbation method to the RLW and generalized modified Boussinesq equations,” Physics Letters A, vol. 364, no. 1, pp. 1–6, 2007. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  18. R. Liu, β€œSome new results on explicit traveling wave solutions of k(m,n) equation,” Discrete and Continuous Dynamical Systems B, vol. 13, no. 3, pp. 633–646, 2010. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  19. M. Li, S. C. Lim, and S. Y. Chen, β€œExact solution of impulse response to a class of fractional oscillators and its stability,” Mathematical Problems in Engineering, vol. 2011, Article ID 657839, 9 pages, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  20. J.-B. Li, β€œExact traveling wave solutions and dynamical behavior for the (n+1)-dimensional multiple sine-Gordon equation,” Science in China A, vol. 50, no. 2, pp. 153–164, 2007. View at Publisher Β· View at Google Scholar Β· View at MathSciNet
  21. C. Cattani and G. Pierro, β€œComplexity on acute myeloid leukemia mRNA transcript variant,” Mathematical Problems in Engineering, vol. 2011, Article ID 379873, 16 pages, 2011. View at Publisher Β· View at Google Scholar
  22. C. Cattani and A. Kudreyko, β€œApplication of periodized harmonic wavelets towards solution of eigenvalue problems for integral equations,” Mathematical Problems in Engineering, vol. 2010, Article ID 570136, 8 pages, 2010. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  23. Z. R. Liu and Z. Y. Ouyang, β€œA note on solitary waves for modified forms of Camassa-Holm and Degasperis-Procesi equations,” Physics Letters A, vol. 366, no. 4-5, pp. 377–381, 2007. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  24. E. G. Bakhoum and C. Toma, β€œSpecific mathematical aspects of dynamics generated by coherence functions,” Mathematical Problems in Engineering, vol. 2011, Article ID 436198, 10 pages, 2011. View at Publisher Β· View at Google Scholar
  25. M.-Y. Tang and W.-L. Zhang, β€œFour types of bounded wave solutions of CH-γ equation,” Science in China A, vol. 50, no. 1, pp. 132–152, 2007. View at Publisher Β· View at Google Scholar Β· View at MathSciNet
  26. Q.-D. Wang and M.-Y. Tang, β€œNew explicit periodic wave solutions and their limits for modified form of Camassa-Holm equation,” Acta Mathematicae Applicatae Sinica, vol. 26, no. 3, pp. 513–524, 2010. View at Publisher Β· View at Google Scholar
  27. E. G. Bakhoum and C. Toma, β€œMathematical transform of traveling-wave equations and phase aspects of quantum interaction,” Mathematical Problems in Engineering, vol. 2010, Article ID 695208, 15 pages, 2010. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  28. R. Liu, β€œSeveral new types of solitary wave solutions for the generalized Camassa-Holm-Degasperis-Procesi equation,” Communications on Pure and Applied Analysis, vol. 9, no. 1, pp. 77–90, 2010. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  29. M. Li, β€œFractal time series—a tutorial review,” Mathematical Problems in Engineering, vol. 2010, Article ID 157264, 26 pages, 2010. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  30. S. Y. Chen, H. Tong, Z. Wang, S. Liu, M. Li, and B. Zhang, β€œImproved generalized belief propagation for vision processing,” Mathematical Problems in Engineering, vol. 2011, Article ID 416963, 12 pages, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH