# Exact Solutions of Travelling Wave Model via Dynamical System Method

Authors: Heng Wang; Longwei Chen; Hongjiang Liu
Journal: Abstract and Applied Analysis
 (2016)

---

## Abstract

By using the method of dynamical system, the exact travelling wave solutions of the coupled nonlinear Schrödinger-Boussinesq equations are studied. Based on this method, the bounded exact travelling wave solutions are obtained which contain solitary wave solutions and periodic travelling wave solutions. The solitary wave solutions and periodic travelling wave solutions are expressed by the hyperbolic functions and the Jacobian elliptic functions, respectively. The results show that the presented findings improve the related previous conclusions. Furthermore, the numerical simulations of the solitary wave solutions and the periodic travelling wave solutions are given to show the correctness of our results.


---

## Body

## 1. Introduction

In laser and plasma physics, the significant problems under interactions between a nonlinear real Boussinesq field and a nonlinear complex Schrödinger field have been raised [1]. In particular, the investigation of the coupled Schrödinger-Boussinesq equations has attracted much attention of physicists and mathematicians. In [2], the existence of the solutions for the equations was investigated. In [3], Guo and Du studied the local and global well-posedness of the periodic boundary value problem for the nonlinear Schrödinger-Boussinesq system. In [4], the approximate solutions and conservation law for the nonlinear Schrödinger-Boussinesq equations have been studied. The study of the coupled nonlinear Schrödinger-Boussinesq equations plays a crucial role in the study of nonlinear scientific fields.In this paper, we consider the following coupled nonlinear Schrödinger-Boussinesq equations [5]:(1)
i
E
t
+
E
x
x
+
γ
E
=
N
E
,
3
N
t
t
-
N
x
x
x
x
+
3
N
2
x
x
+
δ
N
x
x
=
E
2
x
x
, where γ, δ are real parameters. E is complex function and N is real function. Equations (1) are known to describe various physical processes in laser and plasma physics, such as Langmuir field amplitude, modulational instabilities, and intense electromagnetic waves [6]. The study of travelling wave solutions of the coupled Schrödinger-Boussinesq equations has also attracted much attention of physicists and mathematicians. In [7], Farah and Pastor used the (G
′
/
G)-expansion method to construct travelling wave solutions for the equations. In [8], Chen and Xu used the F-expansion method to obtain some periodic wave solutions for (1). In [9], Cai et al. studied coupled equations (1) by the modified expansion method and so on. However, we notice that the previous authors did not study the nonlinear dynamics of (1) and did not find all possible travelling wave solutions. Therefore, it is essential to study the nonlinear dynamics of (1) and find all possible travelling wave solutions of (1). Here, we use the approach of dynamical system to solve (1) and to give some bounded travelling wave solutions of (1) [10, 11].The rest of this paper is built up as follows. In Section2, we give the description of the dynamical system method. In Section 3, we apply this method to solve (1) and numerical simulations are conducted for the solitary wave solutions and the periodic travelling wave solutions to (1) with the aid of the Maple software. Finally, a conclusion is given in Section 4.


## 2. The Dynamical System Method

In this section, we describe the dynamical system method for finding travelling wave solutions of nonlinear wave equations. Suppose an(
n
+
1
)-dimensional nonlinear partial differential equation is given as follows:(2)
P
t
,
x
i
,
u
t
,
u
x
i
,
u
x
i
x
i
,
u
x
i
x
j
,
u
t
t
,
…
=
0
,
i
,
j
=
1,2
,
…
,
n
. The main steps of the dynamical system method are as follows.Step 1 (reduction of (2)).
Making a transformationu
(
t
,
x
1
,
x
2
,
…
,
x
n
)
=
ϕ
(
ξ
),  ξ
=
∑
i
=
1
n
‍
k
i
x
i
-
c
t, (2) can be reduced to a nonlinear ordinary differential equation(3)
D
ξ
,
ϕ
,
ϕ
ξ
,
ϕ
ξ
ξ
,
ϕ
ξ
ξ
ξ
,
…
=
0
, where k
i are nonzero constant and c is the wave speed. Integrating several times for (3), if it can be reduced to the second-order nonlinear ordinary differential equation(4)
E
ϕ
,
ϕ
ξ
,
ϕ
ξ
ξ
=
0
, then let ϕ
ξ
=
d
ϕ
/
d
ξ
=
y; (4) can be reduced to a two-dimensional dynamical system(5)
d
ϕ
d
ξ
=
y
,
d
y
d
ξ
=
f
ϕ
,
y
, where f
ϕ
,
y is an integral expression or a fraction. If f
(
ϕ
,
y
) is a fraction such as f
(
ϕ
,
y
)
=
F
(
ϕ
,
y
)
/
g
(
ϕ
) and g
(
ϕ
s
)
=
0, d
y
/
d
ξ does not exist when ϕ
=
ϕ
s. Then we will make a transformation d
ξ
=
g
(
ϕ
)
d
ζ; system (5) can be rewritten as(6)
d
ϕ
d
ζ
=
g
ϕ
y
,
d
y
d
ζ
=
F
ϕ
,
y
, where ζ is a parameter. If (2) can be reduced to system (5) or (6), then we can go on to the next step.Step 2 (discussion of bifurcations of phase portraits of system (5)).
If system (5) is an integral system, systems (5) and (6) can be reduced to the differential equation(7)
d
y
d
ϕ
=
f
ϕ
,
y
y
,
d
y
d
ϕ
=
F
ϕ
,
y
g
ϕ
y
=
f
ϕ
,
y
y
, and then systems (5) and (6) have the same first integral (i.e., Hamiltonian) as follows:(8)
H
ϕ
,
y
=
h
, where h is an integral constant. According to the first integral, we can get all kinds of phase portraits in the parametric space. Because the phase orbits that defined the vector fields of system (5) (or system (6)) determine all their travelling wave solutions of (2), we can investigate the bifurcations of phase portraits of system (5) (or system (6)) to seek the travelling wave solutions of (2). Usually, a periodic orbit always corresponds to a periodic wave solution; a homoclinic orbit always corresponds to a solitary wave solution; a heteroclinic orbit (or the so-called connecting orbit) always corresponds to kink (or antikink) wave solution. When we find all phase orbits, we can get the value of h or its range.Step 3 (calculation of the first equation of system (5)).
Afterh is determined, we can get the following relationship from (8):(9)
y
=
y
ϕ
,
h
;that is, d
ϕ
/
d
ξ
=
y
(
ϕ
,
h
). If expression (9) is an integral expression, then substituting it into the first term of (5) and integrating we obtain(10)
∫
ϕ
0
ϕ
d
φ
y
φ
,
h
=
∫
0
ξ
d
τ
,where ϕ
(
0
) and 0 are initial constants. Usually, the initial constants can be taken by a root of (9) or inflection points of the travelling waves. Taking proper initial constants and integrating (10), through the Jacobian elliptic functions [12], we can obtain the exact travelling wave solutions of (2).From the above description of the “three-step method,” we can see that solutions of (2) can be obtained by studying and solving the dynamical system simplified by (2). Therefore, this approach is called dynamical system method. The different nonlinear wave equations correspond to different dynamical systems. The different dynamical systems correspond to different travelling wave solutions. This is the whole process of the dynamical system method.


## 3. Travelling Wave Solutions of (1)

### 3.1. The Reduced Dynamical System of (1)

Following the procedure described in Section2, we solve (1) by using the dynamical system method.Use the transformation(11)
E
x
,
t
=
ϕ
ξ
e
i
η
,
N
x
,
t
=
ψ
ξ
,
ξ
=
k
x
-
c
t
,
η
=
p
x
+
l
t
, where k, c, p, and l are travelling wave parameters. Substituting (11) into the first equation of (1), canceling e
i
η, and separating the real and imaginary parts, we have(12)
2
k
p
-
c
ϕ
′
=
0
,
k
2
ϕ
′
′
-
p
2
+
l
-
γ
ϕ
-
ϕ
ψ
=
0
. Obviously, from (12), we know that if ϕ
′
=
0, (1) has the trivial solution. Otherwise, (12) must be satisfied:(13)
2
k
p
-
c
=
0
. Substituting (11) into the second equation of (1) and integrating twice (integral constant is zero), we have(14)
k
4
ψ
′
′
-
δ
k
2
+
3
c
2
ψ
-
3
k
2
ψ
2
+
k
2
ϕ
2
=
0
. Therefore, (1) is reduced to(15)
2
k
p
-
c
=
0
,
k
2
ϕ
′
′
-
p
2
+
l
-
γ
ϕ
-
ϕ
ψ
=
0
,
k
4
ψ
′
′
-
δ
k
2
+
3
c
2
ψ
-
3
k
2
ψ
2
+
k
2
ϕ
2
=
0
. It is very difficult to solve these equations by some ordinary methods because of the coupling of (1), so we consider the special transformation in subtle ways:(16)
ϕ
=
n
ψ
. Here, n is a constant to be determined later. Although in the ordinary course of events there is not linear relation between ϕ and ψ, the transformation is a key point to look for the travelling wave solutions of (1). Substituting (16) into (15), (15) is changed into(17)
2
k
p
-
c
=
0
,
k
2
ψ
′
′
-
p
2
+
l
ψ
-
ψ
2
=
0
,
β
k
3
ψ
′
′
-
c
ψ
+
α
k
2
ψ
2
-
m
2
k
ψ
2
=
0
. Compared with coefficients of the second equation and the third equation of (17), they follow(18)
n
=
±
2
,
p
2
+
l
-
γ
=
δ
k
2
+
3
c
2
k
2
,
2
k
p
-
c
=
0
. A result of the freedom of these parameters is consistency; under condition (18), (17) is simplified to the following equation:(19)
k
2
ψ
′
′
-
p
2
+
l
-
γ
ψ
-
ψ
2
=
0
. Suppose that k
≠
0 and write that A
=
p
2
+
l
-
γ
/
k
2, B
=
1
/
k
2. Thus, (19) has the following form:(20)
ψ
′
′
-
A
ψ
-
B
ψ
2
=
0
, which corresponds to the two-dimensional Hamiltonian system(21)
d
ψ
d
ξ
=
y
,
d
y
d
ξ
=
A
ψ
+
B
ψ
2
.


### 3.2.  The Bifurcations of Phase Portraits of System (21)

We divided by (21) and solve the Bernoulli one-order differential equation; we have the first integral(22)
H
ψ
,
y
=
1
2
y
2
-
1
2
A
ψ
2
-
1
3
B
ψ
3
.Let the integral constant beh; that is, H
(
ψ
,
y
)
=
h. Now we consider the phase portraits of (21). Let the right hand terms of system (21) be zeros; that is, y
=
0 and A
ψ
+
B
ψ
2
=
0; we obtain that system (21) has two equilibrium points S
(
-
A
/
B
,
0
) and O
(
0,0
). For H
(
ψ
,
y
)
=
1
/
2
y
2
-
1
/
2
A
ψ
2
-
1
/
3
B
ψ
3
=
h, we write h
0
=
H
(
0,0
)
=
0, h
1
=
H
(
-
A
/
B
,
0
)
=
-
A
3
/
6
B
2. Because only bounded travelling waves are meaningful to a physical model, we just pay our attention to the bounded solutions of (1) which are physically acceptable. In addition, because of B
=
1
/
k
2
>
0, here we just consider the phase portraits of (21) when B
>
0. With the change of the parameter group of A and B, the system has different phase portraits for (21) which are shown in Figure 1 (drawn by software Maple).Figure 1
The bifurcations of phase portraits of (21).From Figure1, we summarize crucial conclusions as follows:(1)
whenA
=
0, O is a cusp; when A
>
0 (<
0), O is a saddle point (center point) and S is a center (saddle point);
(2)
ifA
≠
0, system (21) has a unique homoclinic orbit Γ which is asymptotic to the saddle and enclosing the center;
(3)
ifA
≠
0, there is a family of periodic orbits which are enclosing the center and filling up the interior of the homoclinic orbit Γ.In the first image of Figure1, we take A
=
0, B
>
0. In the second image of Figure 1, we take A
<
0, B
>
0. In the third image of Figure 1, we take A
>
0, B
>
0.


### 3.3. The Exact Travelling Wave Solutions of (1)

We can get the following relationship from (22):(23)
y
=
A
ψ
2
+
2
3
B
ψ
3
+
2
h
;then substituting it into the first term of (21) and integrating we use the Maple software and the Jacobian elliptic functions; we have the following parametric representation.(
1
) When A
>
0, B
>
0, and h
=
h
0
=
0 (see the third image of Figure 1), there exists a smooth solitary solution which corresponds to a smooth homoclinic orbit Γ of (21) defined by H
(
ψ
,
y
)
=
h
0
=
0; we have the parametric representation:(24)
ψ
ξ
=
-
3
A
+
3
A
tanh
2
⁡
A
/
2
ξ
2
B
.(
2
) When A
<
0, B
>
0, and h
=
h
1
=
-
A
3
/
6
B
2 (see the second image of Figure 1), there exists a smooth solitary solution which corresponds to a smooth homoclinic orbit Γ of (21) defined by H
(
ψ
,
y
)
=
h
1; we have the parametric representation:(25)
ψ
ξ
=
A
B
-
3
A
2
B
sech
2
⁡
A
2
ξ
.(
3
) When A
>
0 (A
<
0), B
>
0, and h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)) (see the second image and third image of Figure 1), similarly, there exist periodic travelling wave solutions which correspond to the family of periodic orbits Γ
h of (21) defined by H
(
ψ
,
y
)
=
h, h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)); we have the following parametric representation:(26)
ψ
ξ
=
z
3
+
z
2
-
z
3
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
2
-
z
3
z
1
-
z
3
,where the parameters z
1, z
2, and z
3 and z
1
>
z
2
>
z
3 are defined by y
2
=
2
h
+
A
ψ
2
+
2
/
3
B
ψ
3
=
2
/
3
B
(
z
1
-
ψ
)
(
z
2
-
ψ
)
(
ψ
-
z
3
).By using the above results and considering condition (18), according to (11), we obtain the travelling wave solutions of (1) as follows. By using the numerical simulation method, the 3D graphics of bounded solutions of (1) are also shown in Figures 2–7.Figure 2
The 3D graphics of|
E
1
| (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.1).Figure 3
The 3D graphics ofN
1 (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.1).Figure 4
The 3D graphics of|
E
2
| (n
=
2, k
=
1, p
=
1
/
11, l
=
-
27, c
=
2
1
/
11, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.5).Figure 5
The 3D graphics ofN
2 (n
=
2, k
=
1, p
=
1
/
11, l
=
-
27, c
=
2
1
/
11, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.5).Figure 6
The 3D graphics of|
E
3
| (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, h
1
=
-
243
/
4, -
5
≤
x
≤
5, and 0
≤
t
≤
5).Figure 7
The 3D graphics ofN
3 (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, h
1
=
-
243
/
4, -
5
≤
x
≤
5, and 0
≤
t
≤
5).(
1
) When A
>
0, B
>
0, and h
=
h
0
=
0,(27)
E
1
x
,
t
=
-
3
A
n
2
B
+
3
A
n
2
B
tanh
2
⁡
A
2
ξ
e
i
η
,
N
1
x
,
t
=
-
3
A
n
2
B
+
3
A
n
2
B
tanh
2
⁡
A
2
ξ
.(
2
) When A
<
0, B
>
0, and h
=
h
1
=
-
A
3
/
6
B
2,(28)
E
2
x
,
t
=
n
A
B
-
3
n
A
2
B
sech
2
⁡
A
2
ξ
e
i
η
,
N
2
x
,
t
=
n
A
B
-
3
n
A
2
B
sech
2
⁡
A
2
ξ
.(
3
) When A
>
0 (A
<
0), B
>
0, and h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)),(29)
E
3
x
,
t
=
n
z
3
+
n
z
2
-
z
3
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
2
-
z
3
z
1
-
z
3
e
i
η
,
N
3
x
,
t
=
z
3
+
z
2
-
z
3
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
2
-
z
3
z
1
-
z
3
.From Figures2–7, it is easy to see that E
1, N
1, E
2, and N
2 are solitary wave solutions which are expressed by the hyperbolic functions. E
3 and N
3 are periodic travelling wave solution which is expressed by Jacobian elliptic functions. Note that our solutions are different from the given ones in [5, 8, 9].


## 3.1. The Reduced Dynamical System of (1)

Following the procedure described in Section2, we solve (1) by using the dynamical system method.Use the transformation(11)
E
x
,
t
=
ϕ
ξ
e
i
η
,
N
x
,
t
=
ψ
ξ
,
ξ
=
k
x
-
c
t
,
η
=
p
x
+
l
t
, where k, c, p, and l are travelling wave parameters. Substituting (11) into the first equation of (1), canceling e
i
η, and separating the real and imaginary parts, we have(12)
2
k
p
-
c
ϕ
′
=
0
,
k
2
ϕ
′
′
-
p
2
+
l
-
γ
ϕ
-
ϕ
ψ
=
0
. Obviously, from (12), we know that if ϕ
′
=
0, (1) has the trivial solution. Otherwise, (12) must be satisfied:(13)
2
k
p
-
c
=
0
. Substituting (11) into the second equation of (1) and integrating twice (integral constant is zero), we have(14)
k
4
ψ
′
′
-
δ
k
2
+
3
c
2
ψ
-
3
k
2
ψ
2
+
k
2
ϕ
2
=
0
. Therefore, (1) is reduced to(15)
2
k
p
-
c
=
0
,
k
2
ϕ
′
′
-
p
2
+
l
-
γ
ϕ
-
ϕ
ψ
=
0
,
k
4
ψ
′
′
-
δ
k
2
+
3
c
2
ψ
-
3
k
2
ψ
2
+
k
2
ϕ
2
=
0
. It is very difficult to solve these equations by some ordinary methods because of the coupling of (1), so we consider the special transformation in subtle ways:(16)
ϕ
=
n
ψ
. Here, n is a constant to be determined later. Although in the ordinary course of events there is not linear relation between ϕ and ψ, the transformation is a key point to look for the travelling wave solutions of (1). Substituting (16) into (15), (15) is changed into(17)
2
k
p
-
c
=
0
,
k
2
ψ
′
′
-
p
2
+
l
ψ
-
ψ
2
=
0
,
β
k
3
ψ
′
′
-
c
ψ
+
α
k
2
ψ
2
-
m
2
k
ψ
2
=
0
. Compared with coefficients of the second equation and the third equation of (17), they follow(18)
n
=
±
2
,
p
2
+
l
-
γ
=
δ
k
2
+
3
c
2
k
2
,
2
k
p
-
c
=
0
. A result of the freedom of these parameters is consistency; under condition (18), (17) is simplified to the following equation:(19)
k
2
ψ
′
′
-
p
2
+
l
-
γ
ψ
-
ψ
2
=
0
. Suppose that k
≠
0 and write that A
=
p
2
+
l
-
γ
/
k
2, B
=
1
/
k
2. Thus, (19) has the following form:(20)
ψ
′
′
-
A
ψ
-
B
ψ
2
=
0
, which corresponds to the two-dimensional Hamiltonian system(21)
d
ψ
d
ξ
=
y
,
d
y
d
ξ
=
A
ψ
+
B
ψ
2
.


## 3.2.  The Bifurcations of Phase Portraits of System (21)

We divided by (21) and solve the Bernoulli one-order differential equation; we have the first integral(22)
H
ψ
,
y
=
1
2
y
2
-
1
2
A
ψ
2
-
1
3
B
ψ
3
.Let the integral constant beh; that is, H
(
ψ
,
y
)
=
h. Now we consider the phase portraits of (21). Let the right hand terms of system (21) be zeros; that is, y
=
0 and A
ψ
+
B
ψ
2
=
0; we obtain that system (21) has two equilibrium points S
(
-
A
/
B
,
0
) and O
(
0,0
). For H
(
ψ
,
y
)
=
1
/
2
y
2
-
1
/
2
A
ψ
2
-
1
/
3
B
ψ
3
=
h, we write h
0
=
H
(
0,0
)
=
0, h
1
=
H
(
-
A
/
B
,
0
)
=
-
A
3
/
6
B
2. Because only bounded travelling waves are meaningful to a physical model, we just pay our attention to the bounded solutions of (1) which are physically acceptable. In addition, because of B
=
1
/
k
2
>
0, here we just consider the phase portraits of (21) when B
>
0. With the change of the parameter group of A and B, the system has different phase portraits for (21) which are shown in Figure 1 (drawn by software Maple).Figure 1
The bifurcations of phase portraits of (21).From Figure1, we summarize crucial conclusions as follows:(1)
whenA
=
0, O is a cusp; when A
>
0 (<
0), O is a saddle point (center point) and S is a center (saddle point);
(2)
ifA
≠
0, system (21) has a unique homoclinic orbit Γ which is asymptotic to the saddle and enclosing the center;
(3)
ifA
≠
0, there is a family of periodic orbits which are enclosing the center and filling up the interior of the homoclinic orbit Γ.In the first image of Figure1, we take A
=
0, B
>
0. In the second image of Figure 1, we take A
<
0, B
>
0. In the third image of Figure 1, we take A
>
0, B
>
0.


## 3.3. The Exact Travelling Wave Solutions of (1)

We can get the following relationship from (22):(23)
y
=
A
ψ
2
+
2
3
B
ψ
3
+
2
h
;then substituting it into the first term of (21) and integrating we use the Maple software and the Jacobian elliptic functions; we have the following parametric representation.(
1
) When A
>
0, B
>
0, and h
=
h
0
=
0 (see the third image of Figure 1), there exists a smooth solitary solution which corresponds to a smooth homoclinic orbit Γ of (21) defined by H
(
ψ
,
y
)
=
h
0
=
0; we have the parametric representation:(24)
ψ
ξ
=
-
3
A
+
3
A
tanh
2
⁡
A
/
2
ξ
2
B
.(
2
) When A
<
0, B
>
0, and h
=
h
1
=
-
A
3
/
6
B
2 (see the second image of Figure 1), there exists a smooth solitary solution which corresponds to a smooth homoclinic orbit Γ of (21) defined by H
(
ψ
,
y
)
=
h
1; we have the parametric representation:(25)
ψ
ξ
=
A
B
-
3
A
2
B
sech
2
⁡
A
2
ξ
.(
3
) When A
>
0 (A
<
0), B
>
0, and h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)) (see the second image and third image of Figure 1), similarly, there exist periodic travelling wave solutions which correspond to the family of periodic orbits Γ
h of (21) defined by H
(
ψ
,
y
)
=
h, h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)); we have the following parametric representation:(26)
ψ
ξ
=
z
3
+
z
2
-
z
3
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
2
-
z
3
z
1
-
z
3
,where the parameters z
1, z
2, and z
3 and z
1
>
z
2
>
z
3 are defined by y
2
=
2
h
+
A
ψ
2
+
2
/
3
B
ψ
3
=
2
/
3
B
(
z
1
-
ψ
)
(
z
2
-
ψ
)
(
ψ
-
z
3
).By using the above results and considering condition (18), according to (11), we obtain the travelling wave solutions of (1) as follows. By using the numerical simulation method, the 3D graphics of bounded solutions of (1) are also shown in Figures 2–7.Figure 2
The 3D graphics of|
E
1
| (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.1).Figure 3
The 3D graphics ofN
1 (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.1).Figure 4
The 3D graphics of|
E
2
| (n
=
2, k
=
1, p
=
1
/
11, l
=
-
27, c
=
2
1
/
11, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.5).Figure 5
The 3D graphics ofN
2 (n
=
2, k
=
1, p
=
1
/
11, l
=
-
27, c
=
2
1
/
11, γ
=
-
25, δ
=
-
3, -
5
≤
x
≤
5, and 0
≤
t
≤
0.5).Figure 6
The 3D graphics of|
E
3
| (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, h
1
=
-
243
/
4, -
5
≤
x
≤
5, and 0
≤
t
≤
5).Figure 7
The 3D graphics ofN
3 (n
=
2, k
=
p
=
1, c
=
2, l
=
-
17, γ
=
-
25, δ
=
-
3, h
1
=
-
243
/
4, -
5
≤
x
≤
5, and 0
≤
t
≤
5).(
1
) When A
>
0, B
>
0, and h
=
h
0
=
0,(27)
E
1
x
,
t
=
-
3
A
n
2
B
+
3
A
n
2
B
tanh
2
⁡
A
2
ξ
e
i
η
,
N
1
x
,
t
=
-
3
A
n
2
B
+
3
A
n
2
B
tanh
2
⁡
A
2
ξ
.(
2
) When A
<
0, B
>
0, and h
=
h
1
=
-
A
3
/
6
B
2,(28)
E
2
x
,
t
=
n
A
B
-
3
n
A
2
B
sech
2
⁡
A
2
ξ
e
i
η
,
N
2
x
,
t
=
n
A
B
-
3
n
A
2
B
sech
2
⁡
A
2
ξ
.(
3
) When A
>
0 (A
<
0), B
>
0, and h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)),(29)
E
3
x
,
t
=
n
z
3
+
n
z
2
-
z
3
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
2
-
z
3
z
1
-
z
3
e
i
η
,
N
3
x
,
t
=
z
3
+
z
2
-
z
3
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
2
-
z
3
z
1
-
z
3
.From Figures2–7, it is easy to see that E
1, N
1, E
2, and N
2 are solitary wave solutions which are expressed by the hyperbolic functions. E
3 and N
3 are periodic travelling wave solution which is expressed by Jacobian elliptic functions. Note that our solutions are different from the given ones in [5, 8, 9].


## 4. Conclusion

To summarize, by using the dynamical system method, the bounded exact travelling wave solutions (solitary wave solutions and periodic wave solutions) have been obtained for the coupled nonlinear Schrödinger-Boussinesq equations. The dynamical system method is a good method to obtain exact solutions, which can not only obtain exact solutions but also understand nonlinear dynamics of travelling wave equations. We show that the hyperbolic function solutions and the Jacobian elliptic function solutions we found in this paper are different from the solutions presented by other authors before. The results enrich the diversity of wave structures of the coupled nonlinear Schrödinger-Boussinesq equations.Furthermore, if there is not the condition thatB
>
0, system (21) has another case. When A
>
0 (A
<
0), B
<
0, there exist periodic travelling wave solutions corresponding to the family of periodic orbits Γ
h of (21) defined by H
(
ψ
,
y
)
=
h, h
∈
(
h
1
,
0
) (h
∈
(
0
,
h
1
)); we have following parametric representation:(30)
ψ
ξ
=
z
1
-
z
1
-
z
2
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
1
-
z
2
z
1
-
z
3
,where the parameters z
1, z
2, and z
3 and z
1
>
z
2
>
z
3 are defined by y
2
=
2
h
+
A
ψ
2
+
2
/
3
B
ψ
3
=
2
/
3
|
B
|
(
z
1
-
ψ
)
(
ψ
-
z
2
)
(
ψ
-
z
3
). Then we obtain additional travelling wave solutions of (1) as follows:(31)
E
x
,
t
=
n
z
1
-
n
z
1
-
z
2
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
1
-
z
2
z
1
-
z
3
e
i
η
,
N
x
,
t
=
z
1
-
z
1
-
z
2
sn
2
6
B
z
1
-
z
3
6
ξ
,
z
1
-
z
2
z
1
-
z
3
.



---

*Source: 9290734-2016-01-20.xml*