# Traveling Wave Solutions and Infinite-Dimensional Linear Spaces of Multiwave Solutions to Jimbo-Miwa Equation

Authors: Lijun Zhang; C. M. Khalique
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

The traveling wave solutions and multiwave solutions to (3 + 1)-dimensional Jimbo-Miwa equation are investigated in this paper. As a result, besides the
exact bounded solitary wave solutions, we obtain the existence of two families
of bounded periodic traveling wave solutions and their implicit formulas by
analysis of phase portrait of the corresponding traveling wave system. We
derive the exact 2-wave solutions and two families of arbitrary finiteN-wave
solutions by studying the linear space of its Hirota bilinear equation, which
confirms that the (3 + 1)-dimensional Jimbo-Miwa equation admits multiwave solutions of any order and is completely integrable.


---

## Body

## 1. Introduction

Various nonlinear partial differential equations (NLPDEs) have been proposed to model different kinds of phenomena in natural and applied sciences such as fluid dynamics, plasma physics, solid-state physics, optical fibers, acoustics, mechanics, biology, and mathematical finance. Obviously, it is of significant importance to study the solutions of such NLPDEs from both theoretical and practical points of view. However, the solution spaces of nonlinear equations are infinite-dimensional and contain diverse solution structures, so it is usually a difficult job to determine the solutions to nonlinear NLPDEs.A great idea to generate exact solutions of NLPDEs is to reduce the NLPDEs into some algebraic equations by assuming the solutions to have some special forms or satisfy some solvable simpler equations. This can be seen in, for example, the exp-function method [1], the tanh function method [2], the homogeneous balance method [3, 4], the auxiliary function method [5, 6], the sech-function method [7], the sine-cosine method [8, 9], the tanh-coth method [10], the Jacobi elliptic function method [11], the F-expansion method, and the extended F-expansion method [12]. Normally, it is not an easy task to solve these nonlinear algebraic equations reduced from NLPDEs because they involve many parameters. However, the rapid development of symbolic computation makes it relatively easy to solve these algebraic equations [13]. At the same time one should be very careful when applying these methods as different methods might give the same solutions. For example, solutions obtained by using the sech-function method, the tanh-coth method, and the exp-function method are actually the same [14, 15], because f
(
sech
2
x
)
=
f
(
1
-
tanh
2
x
)
=
f
(
4
/
(
e
x
+
e
-
x
)
2
), for any function f.Recently, the planar dynamical system theorem has been employed to study the traveling wave solutions of NLPDEs [16, 17]. The best advantage of this approach is that the boundedness, periodicity, the shapes of the solutions, and even the singular traveling wave solutions of NLPDEs can be recognized easily from the corresponding orbits of their phase portraits under various different parametric conditions. Also the exact solutions can be derived at the same time. It is worth pointing out that Hirota’s bilinear method [18] is also an amazing method to find exact solutions to some NLPDEs. By some independent variable transformation, various nonlinear equations of mathematical physics are transformed into Hirota’s bilinear equations [18, 19], which possess some specific properties and thus might be applied to study the solution sets of nonlinear differential equations. Recently, even some programs have been designed and some algorithms have been proposed on searching for integrable bilinear equations [20–26]. Based on Hirota’s bilinear form, soliton solutions were obtained by the Hirota perturbation technique [18], the multiple exp-function algorithm [27], and other methods [28–31]. Even some generalized bilinear form has been proposed recently by Ma [32].The(
3
+
1
)-dimensional Jimbo-Miwa equation
(1)
u
x
x
x
y
+
3
(
u
x
u
y
)
x
+
2
u
y
t
-
3
u
x
z
=
0
was first proposed and studied by Jimbo and Miwa [33]. Recently, this equation has attracted a great deal of attention which mainly focuses on its solutions, integrability properties and symmetries. Ma and Lee [34] proposed a direct approach to solve (1) by using rational function transformations. Li and Dai [35] applied the generalized Riccati equation method to look for its exact solutions. However, in [36], Kudryashov and Sinelshchikov pointed out some mistakes and proved that some solutions in [35] could be written in a uniform form and thus they are not new at all.The(
3
+
1
)
-dimensional Jimbo-Miwa equation can be transformed into the Hirota bilinear equation
(2)
(
D
x
3
D
y
+
2
D
t
D
y
-
3
D
x
D
z
)
f
·
f
=
0
through the dependent variable transformation u
=
2
(
ln
⁡
f
)
x.In this paper, we firstly study the bounded traveling wave solutions of the Jimbo-Miwa equation (1) by investigating the bifurcation and phase portraits of a planar cubic polynomial ordinary differential equation by using the planar dynamical system theory [16, 37]. We then explore the 1-wave and 2-wave solutions and the infinite-dimensional linear spaces of multiwave solutions to Jimbo-Miwa equation by employing Hirota’s bilinear method, thus confirming that Jimbo-Miwa equation possesses multiwave solutions of arbitrary order.


## 2. Traveling Wave Solutions to the(
3
+
1
)
-Dimensional Jimbo-Miwa Equation

To investigate the traveling wave solutions to the(
3
+
1
)
-dimensional Jimbo-Miwa equation, we make the traveling wave transformation u
(
x
,
y
,
z
,
t
)
=
u
(
ξ
), ξ
=
k
x
+
l
y
+
m
z
+
n
t, under which (1) is reduced to the nonlinear ordinary differential equation
(3)
k
3
l
d
4
u
d
ξ
4
+
6
k
2
l
d
u
d
ξ
d
2
u
d
ξ
2
+
(
2
l
n
-
3
k
m
)
d
2
u
d
ξ
2
=
0
.
Integrating (1) once with respect to ξ gives
(4)
k
3
l
d
3
u
d
ξ
3
+
3
k
2
l
(
d
u
d
ξ
)
2
+
(
2
l
n
-
3
k
m
)
d
u
d
ξ
=
g
,
where g is an arbitrary constant. Let v
=
d
u
/
d
ξ; then (1) becomes
(5)
k
3
l
d
2
v
d
ξ
2
+
3
k
2
l
v
2
+
(
2
l
n
-
3
k
m
)
v
=
g
,
which is a second-order nonlinear ordinary differential equation. We will study the solutions of (5) by planar dynamical system method and thus derive the traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation.


### 2.1. Bounded Solutions of (5)

First, we rewrite (5) in a simpler and more general form, namely,
(6)
d
2
v
d
ξ
2
=
a
v
2
+
b
v
+
c
,
where a
=
-
3
/
k, b
=
(
3
k
m
-
2
l
n
)
/
(
k
3
l
), and c
=
g
/
(
k
3
l
). We now study the bifurcation and exact solutions of (6). Let v
′
=
w; then (6) is equivalent to the dynamical system
(7)
v
′
=
w
,
w
′
=
a
v
2
+
b
v
+
c
,
which has the Hamiltonian
(8)
H
(
v
,
w
)
=
w
2
2
-
a
3
v
3
-
b
2
v
2
-
c
v
.
Clearly, the phase orbits defined by the vector fields of system (7) determine all solutions of (6). The bounded solutions of (6) correspond to the bounded phase orbits of system (7), which we now investigate. Along the orbit corresponding to H
(
v
,
w
)
=
h,
(9)
(
d
v
d
ξ
)
2
=
2
a
3
v
3
+
b
v
2
+
2
c
v
+
2
h
.
Consequently, the general formula of the solutions of (6) can be expressed as
(10)
∫
0
v
d
v
(
2
a
v
3
/
3
)
+
b
v
2
+
2
c
v
+
2
h
=
±
∫
ξ
0
ξ
d
ξ
.However, it is not easy to know the properties and the shapes of (10) which actually is determined by the parameters a, b, c, and h. Clearly, the abscissas of equilibrium points of system (7) are the zeros of a
v
2
+
b
v
+
c
=
0. Obviously, the system has no bounded orbits when b
2
-
4
a
c
<
0. We suppose that b
2
-
4
a
c
>
0 in order to study the bounded orbits of system (7). Denoting v
±
=
(
-
b
±
b
2
-
4
a
c
)
/
2
a, then E
+
(
v
+
,
0
) and E
-
(
v
-
,
0
) are two equilibrium points of system (7). By the theory of planar dynamical system, we know that E
+ is a saddle point and E
- is a center. Denote h
±
=
H
(
v
±
,
0
), and, by careful computation, we get
(11)
h
±
=
(
b
2
-
4
a
c
)
(
-
b
±
b
2
-
4
a
c
)
+
2
a
b
c
12
a
2
.
Obviously, h
-
<
h
<
h
+. H
(
v
,
w
)
=
h
+ corresponds to homoclinic orbits, and H
(
v
,
w
)
=
h
- corresponds to the center E
- and H
(
v
,
w
)
=
h, where h
+
<
h
<
h
- corresponds to a family of closed orbits surrounding the center E
-, which are surrounded by a homoclinic orbit. That is to say, (10) defines bounded solutions if and only if h
+
≤
h
<
h
-. To be exact, (10) defines a family of periodic solutions when h
+
<
h
<
h
-.Whenh
=
h
+, (10) defines a bounded solution which approaches v
+ as ξ goes to infinity. Actually,
(12)
2
a
3
v
3
+
b
v
2
+
2
c
v
+
2
h
+
=
2
a
3
(
v
-
v
+
)
2
(
v
-
v
0
)
,
where v
0
=
-
(
b
+
2
b
2
-
4
a
c
)
/
2
a, so (10) can be reduced to
(13)
∫
v
0
v
d
v
(
v
-
v
+
)
a
(
v
-
v
0
)
=
2
3
(
ξ
-
ξ
0
)
,
from which we can get the exact solution as
(14)
v
=
v
+
-
(
v
+
-
v
0
)
sech
2
(
a
(
v
+
-
v
0
)
6
(
ξ
-
ξ
0
)
)
.By further simplification, (14) becomes
(15)
v
=
v
+
-
3
b
2
-
4
a
c
2
a
sech
2
[
1
2
(
b
2
-
4
a
c
)
1
/
4
(
ξ
-
ξ
0
)
]
,
which is an exact bounded solution of (5).Thus, we have the following lemma.Lemma 1.
The general second-order ODE (6) has bounded solutions if and only if b
2
-
4
a
c
>
0. The bounded solutions can be expressed as (10) in an implicit form. In fact, provided h
-
<
h
<
h
+, (10) defines a family of bounded periodic solutions and h
=
h
+ defines a bounded solution which approaches v
+ as ξ goes to infinity and can be expressed explicitly as (15), where v
+
=
(
-
b
+
b
2
-
4
a
c
)
/
2
a and h
± is defined by (11).


### 2.2.  Bounded Traveling Wave Solutions to the(
3
+
1
)
-Dimensional Jimbo-Miwa Equation

According to the analysis and results in Section2.1, we know that (5) has only two kinds of bounded solutions, among which one is a family of periodic solutions and another is a family of solutions approaching a fixed number as ξ goes to infinity. Note that what we aim to study is the bounded traveling wave solutions to the (
3
+
1
)-dimensional Jimbo-Miwa equation which are determined by v
=
d
u
/
d
ξ and v satisfies (5). So we have to investigate how we can get the bounded solutions to (3) from the bounded solution of (5).Clearly,u
(
ξ
)
=
∫
ξ
0
ξ
‍
v
(
ξ
)
d
ξ, and v
(
ξ
) can be expressed implicitly as (10). By the geometry meaning of the integral and the properties of the solutions of (5), we get the traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation. To get the bounded solution, we choose the integral constant g to be 0; that is, c
=
0 in (15), and so
(16)
u
(
ξ
)
=
g
-
3
|
b
|
a
tanh
[
|
b
|
2
(
ξ
-
ξ
0
)
]
;
that is,
(17)
u
(
ξ
)
=
g
-
k
|
2
l
n
-
3
k
m
k
3
l
|
×
tanh
[
1
2
|
2
l
n
-
3
k
m
k
3
l
|
(
ξ
-
ξ
0
)
]
,
which is a family of exact bounded kink traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation, where g is an arbitrary constant.However, we may not get bounded solutions from the family of periodic solutions of (5). It is easy to see that if v
(
ξ
) is a periodic solution of (5), then u
(
ξ
)
=
∫
ξ
0
ξ
‍
v
(
ξ
)
d
ξ is bounded if and only if ∫
0
T
v
(
ξ
)
d
ξ
=
0, where T is the period of the function v
(
ξ
). Recall that the period of the function v
(
ξ
) which is defined by (10) with h
-
<
h
<
h
+ continuously depends on the parameters a, b, c, and h. So ∫
0
T
‍
v
(
ξ
)
d
ξ continuously depends on the parameters a, b, c, and h too. Let V
(
a
,
b
,
c
,
h
)
=
∫
0
T
‍
v
(
ξ
)
d
ξ; then V
(
a
,
b
,
c
,
h
) is a continuous function of a, b, c, and h. We now prove the existence of the root of V
(
a
,
b
,
c
,
h
)
=
0 to get the existence of the bounded periodic traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation.By the theory of planar dynamical system, whenh
-
<
h
<
h
+,
(18)
2
a
3
v
3
+
b
v
2
+
2
c
v
+
2
h
+
=
2
a
3
(
v
-
v
1
)
(
v
-
v
2
)
(
v
-
v
3
)
,
where v
1
>
v
+
>
v
2
>
v
-
>
v
3
>
v
0 in the case of a
>
0 and v
0
>
v
1
>
v
-
>
v
2
>
v
+
>
v
3 in the case of a
<
0. The periodic solution satisfies v
+
>
v
2
>
v
(
ξ
)
>
v
3
>
v
0 when a
>
0 and v
0
>
v
1
>
v
(
ξ
)
>
v
2
>
v
+ when a
<
0. Clearly, v
0
>
0 when a
>
0 and 0
>
b
>
-
4
a
c
/
3. So v
(
ξ
)
>
0 and thus V
(
a
,
b
,
c
,
h
)
>
0 when a
>
0 and 0
>
b
>
-
4
a
c
/
3. However, v
(
ξ
)
<
v
+
<
0 and thus V
(
a
,
b
,
c
,
h
)
<
0 when a
>
0
,
c
>
0, and b
>
0. There must exist at least one zero of V
(
a
,
b
,
c
,
h
)
=
0 in the region a
>
0 of the parameter space since V
(
a
,
b
,
c
,
h
) is a continuous function of a, b, c, and h. The same happens in the region a
<
0 of the parameter space. So, we know that there exist at least two families of bounded periodic traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation. Thus we have the following theorem.Theorem 2.
The(
3
+
1
)
-dimensional Jimbo-Miwa equation has two types of bounded traveling wave solutions as given below.(1)
The(
3
+
1
)
-dimensional Jimbo-Miwa equation has a family of exact bounded kink traveling wave solutions
(19)
u
(
x
,
y
,
z
,
t
)
=
g
-
k
|
2
l
n
-
3
k
m
k
3
l
|
×
tanh
[
1
2
|
2
l
n
-
3
k
m
k
3
l
|
(
k
x
+
l
y
+
m
z
+
n
t
-
ξ
0
)
]
,
whereξ
0 and g are two arbitrary constants.
(2)
The(
3
+
1
)-dimensional Jimbo-Miwa equation has at least two families of bounded periodic traveling wave solutions which are determined implicitly by (10) and u
(
ξ
)
=
∫
ξ
0
ξ
‍
v
(
ξ
)
d
ξ, where ξ
=
k
x
+
l
y
+
m
z
+
n
t and ξ
0 is an arbitrary constant.


## 2.1. Bounded Solutions of (5)

First, we rewrite (5) in a simpler and more general form, namely,
(6)
d
2
v
d
ξ
2
=
a
v
2
+
b
v
+
c
,
where a
=
-
3
/
k, b
=
(
3
k
m
-
2
l
n
)
/
(
k
3
l
), and c
=
g
/
(
k
3
l
). We now study the bifurcation and exact solutions of (6). Let v
′
=
w; then (6) is equivalent to the dynamical system
(7)
v
′
=
w
,
w
′
=
a
v
2
+
b
v
+
c
,
which has the Hamiltonian
(8)
H
(
v
,
w
)
=
w
2
2
-
a
3
v
3
-
b
2
v
2
-
c
v
.
Clearly, the phase orbits defined by the vector fields of system (7) determine all solutions of (6). The bounded solutions of (6) correspond to the bounded phase orbits of system (7), which we now investigate. Along the orbit corresponding to H
(
v
,
w
)
=
h,
(9)
(
d
v
d
ξ
)
2
=
2
a
3
v
3
+
b
v
2
+
2
c
v
+
2
h
.
Consequently, the general formula of the solutions of (6) can be expressed as
(10)
∫
0
v
d
v
(
2
a
v
3
/
3
)
+
b
v
2
+
2
c
v
+
2
h
=
±
∫
ξ
0
ξ
d
ξ
.However, it is not easy to know the properties and the shapes of (10) which actually is determined by the parameters a, b, c, and h. Clearly, the abscissas of equilibrium points of system (7) are the zeros of a
v
2
+
b
v
+
c
=
0. Obviously, the system has no bounded orbits when b
2
-
4
a
c
<
0. We suppose that b
2
-
4
a
c
>
0 in order to study the bounded orbits of system (7). Denoting v
±
=
(
-
b
±
b
2
-
4
a
c
)
/
2
a, then E
+
(
v
+
,
0
) and E
-
(
v
-
,
0
) are two equilibrium points of system (7). By the theory of planar dynamical system, we know that E
+ is a saddle point and E
- is a center. Denote h
±
=
H
(
v
±
,
0
), and, by careful computation, we get
(11)
h
±
=
(
b
2
-
4
a
c
)
(
-
b
±
b
2
-
4
a
c
)
+
2
a
b
c
12
a
2
.
Obviously, h
-
<
h
<
h
+. H
(
v
,
w
)
=
h
+ corresponds to homoclinic orbits, and H
(
v
,
w
)
=
h
- corresponds to the center E
- and H
(
v
,
w
)
=
h, where h
+
<
h
<
h
- corresponds to a family of closed orbits surrounding the center E
-, which are surrounded by a homoclinic orbit. That is to say, (10) defines bounded solutions if and only if h
+
≤
h
<
h
-. To be exact, (10) defines a family of periodic solutions when h
+
<
h
<
h
-.Whenh
=
h
+, (10) defines a bounded solution which approaches v
+ as ξ goes to infinity. Actually,
(12)
2
a
3
v
3
+
b
v
2
+
2
c
v
+
2
h
+
=
2
a
3
(
v
-
v
+
)
2
(
v
-
v
0
)
,
where v
0
=
-
(
b
+
2
b
2
-
4
a
c
)
/
2
a, so (10) can be reduced to
(13)
∫
v
0
v
d
v
(
v
-
v
+
)
a
(
v
-
v
0
)
=
2
3
(
ξ
-
ξ
0
)
,
from which we can get the exact solution as
(14)
v
=
v
+
-
(
v
+
-
v
0
)
sech
2
(
a
(
v
+
-
v
0
)
6
(
ξ
-
ξ
0
)
)
.By further simplification, (14) becomes
(15)
v
=
v
+
-
3
b
2
-
4
a
c
2
a
sech
2
[
1
2
(
b
2
-
4
a
c
)
1
/
4
(
ξ
-
ξ
0
)
]
,
which is an exact bounded solution of (5).Thus, we have the following lemma.Lemma 1.
The general second-order ODE (6) has bounded solutions if and only if b
2
-
4
a
c
>
0. The bounded solutions can be expressed as (10) in an implicit form. In fact, provided h
-
<
h
<
h
+, (10) defines a family of bounded periodic solutions and h
=
h
+ defines a bounded solution which approaches v
+ as ξ goes to infinity and can be expressed explicitly as (15), where v
+
=
(
-
b
+
b
2
-
4
a
c
)
/
2
a and h
± is defined by (11).


## 2.2.  Bounded Traveling Wave Solutions to the(
3
+
1
)
-Dimensional Jimbo-Miwa Equation

According to the analysis and results in Section2.1, we know that (5) has only two kinds of bounded solutions, among which one is a family of periodic solutions and another is a family of solutions approaching a fixed number as ξ goes to infinity. Note that what we aim to study is the bounded traveling wave solutions to the (
3
+
1
)-dimensional Jimbo-Miwa equation which are determined by v
=
d
u
/
d
ξ and v satisfies (5). So we have to investigate how we can get the bounded solutions to (3) from the bounded solution of (5).Clearly,u
(
ξ
)
=
∫
ξ
0
ξ
‍
v
(
ξ
)
d
ξ, and v
(
ξ
) can be expressed implicitly as (10). By the geometry meaning of the integral and the properties of the solutions of (5), we get the traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation. To get the bounded solution, we choose the integral constant g to be 0; that is, c
=
0 in (15), and so
(16)
u
(
ξ
)
=
g
-
3
|
b
|
a
tanh
[
|
b
|
2
(
ξ
-
ξ
0
)
]
;
that is,
(17)
u
(
ξ
)
=
g
-
k
|
2
l
n
-
3
k
m
k
3
l
|
×
tanh
[
1
2
|
2
l
n
-
3
k
m
k
3
l
|
(
ξ
-
ξ
0
)
]
,
which is a family of exact bounded kink traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation, where g is an arbitrary constant.However, we may not get bounded solutions from the family of periodic solutions of (5). It is easy to see that if v
(
ξ
) is a periodic solution of (5), then u
(
ξ
)
=
∫
ξ
0
ξ
‍
v
(
ξ
)
d
ξ is bounded if and only if ∫
0
T
v
(
ξ
)
d
ξ
=
0, where T is the period of the function v
(
ξ
). Recall that the period of the function v
(
ξ
) which is defined by (10) with h
-
<
h
<
h
+ continuously depends on the parameters a, b, c, and h. So ∫
0
T
‍
v
(
ξ
)
d
ξ continuously depends on the parameters a, b, c, and h too. Let V
(
a
,
b
,
c
,
h
)
=
∫
0
T
‍
v
(
ξ
)
d
ξ; then V
(
a
,
b
,
c
,
h
) is a continuous function of a, b, c, and h. We now prove the existence of the root of V
(
a
,
b
,
c
,
h
)
=
0 to get the existence of the bounded periodic traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation.By the theory of planar dynamical system, whenh
-
<
h
<
h
+,
(18)
2
a
3
v
3
+
b
v
2
+
2
c
v
+
2
h
+
=
2
a
3
(
v
-
v
1
)
(
v
-
v
2
)
(
v
-
v
3
)
,
where v
1
>
v
+
>
v
2
>
v
-
>
v
3
>
v
0 in the case of a
>
0 and v
0
>
v
1
>
v
-
>
v
2
>
v
+
>
v
3 in the case of a
<
0. The periodic solution satisfies v
+
>
v
2
>
v
(
ξ
)
>
v
3
>
v
0 when a
>
0 and v
0
>
v
1
>
v
(
ξ
)
>
v
2
>
v
+ when a
<
0. Clearly, v
0
>
0 when a
>
0 and 0
>
b
>
-
4
a
c
/
3. So v
(
ξ
)
>
0 and thus V
(
a
,
b
,
c
,
h
)
>
0 when a
>
0 and 0
>
b
>
-
4
a
c
/
3. However, v
(
ξ
)
<
v
+
<
0 and thus V
(
a
,
b
,
c
,
h
)
<
0 when a
>
0
,
c
>
0, and b
>
0. There must exist at least one zero of V
(
a
,
b
,
c
,
h
)
=
0 in the region a
>
0 of the parameter space since V
(
a
,
b
,
c
,
h
) is a continuous function of a, b, c, and h. The same happens in the region a
<
0 of the parameter space. So, we know that there exist at least two families of bounded periodic traveling wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation. Thus we have the following theorem.Theorem 2.
The(
3
+
1
)
-dimensional Jimbo-Miwa equation has two types of bounded traveling wave solutions as given below.(1)
The(
3
+
1
)
-dimensional Jimbo-Miwa equation has a family of exact bounded kink traveling wave solutions
(19)
u
(
x
,
y
,
z
,
t
)
=
g
-
k
|
2
l
n
-
3
k
m
k
3
l
|
×
tanh
[
1
2
|
2
l
n
-
3
k
m
k
3
l
|
(
k
x
+
l
y
+
m
z
+
n
t
-
ξ
0
)
]
,
whereξ
0 and g are two arbitrary constants.
(2)
The(
3
+
1
)-dimensional Jimbo-Miwa equation has at least two families of bounded periodic traveling wave solutions which are determined implicitly by (10) and u
(
ξ
)
=
∫
ξ
0
ξ
‍
v
(
ξ
)
d
ξ, where ξ
=
k
x
+
l
y
+
m
z
+
n
t and ξ
0 is an arbitrary constant.


## 3.N
-Wave Solutions Linear Subspace of the (
3
+
1
)
-Dimensional Jimbo-Miwa Equation

In this section we study theN
-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation by the linear superposition principle to Hirota bilinear equations which was proposed firstly by Ma and Fan [30].Let wave variablesη
i
=
k
i
·
x, 1
≤
i
≤
N, where k
i
=
{
k
1
,
i
,
k
2
,
i
,
k
3
,
i
,
k
4
,
i
}
T, x
=
{
x
,
y
,
z
,
t
}
T. According to the linear superposition principle to Hirota bilinear equations [30, 31], we can get the following sufficient and necessary conditions for span
⁡
{
e
η
1
,
…
,
e
η
N
} being a subspace of the Hirota bilinear equation (2).Theorem 3.
Letk
i
=
{
k
1
,
i
,
k
2
,
i
,
k
3
,
i
,
k
4
,
i
}
T, x
=
{
x
,
y
,
z
,
t
}
T, and f
N
=
∑
i
=
1
N
‍
ϵ
i
e
k
i
·
x. Then, for any constants {
ϵ
i
,
1
≤
i
≤
N
}, f
N solves the Hirota bilinear equation (2) if and only if
(20)
(
k
1
,
i
-
k
1
,
j
)
3
(
k
2
,
i
-
k
2
,
j
)
+
2
(
k
4
,
i
-
k
4
,
j
)
(
k
2
,
i
-
k
2
,
j
)
-
3
(
k
1
,
i
-
k
1
,
j
)
(
k
3
,
i
-
k
3
,
j
)
=
0
for any 1
≤
j
<
i
≤
N.It follows from Theorem3 that f
N is an N-wave solution to the Hirota bilinear equation (2) if (20) holds for any 1
≤
j
<
i
≤
N. However, the corresponding u
(
x
,
y
,
z
,
t
)
=
2
(
f
N
)
x
/
f
N is not necessarily an N-wave solution to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) even if {
η
1
,
η
2
,
…
,
η
N
} is linear independent. In fact,
(21)
u
(
x
,
y
,
z
,
t
)
=
2
(
f
N
)
x
f
N
=
2
∑
i
=
1
N
ϵ
i
k
1
,
i
e
k
i
·
x
∑
i
=
1
N
ϵ
i
e
k
i
·
x
=
2
k
1,1
+
∑
i
=
2
N
k
1
,
i
ϵ
i
′
e
(
k
i
-
k
1
)
·
x
1
+
∑
i
=
2
N
ϵ
i
′
e
(
k
i
-
k
1
)
·
x
=
2
k
1,1
+
2
∑
i
=
2
N
k
1
,
i
′
ϵ
i
′
e
k
i
′
·
x
1
+
∑
i
=
2
N
ϵ
i
′
e
k
i
′
·
x
,
where ϵ
i
′
=
ϵ
i
/
ϵ
1, k
1
,
i
′
=
k
1
,
i
-
k
1,1 and k
i
′
=
k
i
-
k
1. Obviously, the function above is an N
-
1-wave solution to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) if {
η
i
′,   i
=
2,3
,
…
,
N
} is independent; that is, {
k
i
′,   i
=
2
,
…
,
N
} is linear independent, solves (20), and satisfies the dispersion relations. Here, k
i
′
=
{
k
1
,
i
′
,
k
2
,
i
′
,
k
3
,
i
′
,
k
4
,
i
′
}
T satisfies the dispersion relations meaning that k
1
,
i
′
3
k
2
,
i
′
+
2
k
4
,
i
′
k
2
,
i
′
-
3
k
1
,
i
′
k
3
,
i
′
=
0 holds.Recall that the goal of this paper is to investigate the multiwave solutions to the(
3
+
1
)
-dimensional Jimbo-Miwa equation. Combining the conclusion of Theorem 3 with the analysis above, we can now state the following theorem on the N-wave solution to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1).Theorem 4.
Letk
i
=
{
k
1
,
i
,
k
2
,
i
,
k
3
,
i
,
k
4
,
i
}
T, x
=
{
x
,
y
,
z
,
t
}
T, 1
≤
i
≤
N. For any constants g and {
ϵ
i
,
1
≤
i
≤
N
}, suppose that
(22)
u
(
x
,
y
,
z
,
t
)
=
2
∑
i
=
1
N
k
1
,
i
ϵ
i
e
k
i
·
x
1
+
∑
i
=
1
N
ϵ
i
e
k
i
·
x
+
g
.
Then u
(
x
,
y
,
z
,
t
) is an N-wave solution to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1), if {
k
1
,
k
2
,
…
,
k
N
} is a linear independent solution set of equation (20) satisfying the dispersion relations.To get otherN-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation, we study the linear independent solution sets to (20) which are required to satisfy the dispersion relations. Actually, (20) is a system possessing N
(
N
+
1
)
/
2 coupled equations (plus the dispersion conditions). It is usually not so easy to get the solutions of (20). Fortunately, the number of the equations is 1 when N
=
1 and it is 3 when N
=
2, which might make it easy to get the solutions in these two cases. Clearly, 1-wave solution is the traveling wave solution. Let us check what kinds of traveling wave solutions we can obtain by the Hirota bilinear method first.For the caseN
=
1 in (22), we only need to find the independent solutions k
1 of the dispersion relation
(23)
(
k
1,1
)
3
k
2,1
+
2
k
4,1
k
2,1
-
3
k
1,1
k
3,1
=
0
.That is to say, we get the traveling wave solutions to the(
3
+
1
)
-dimensional Jimbo-Miwa equation from (20). Suppose that l
i, i
=
1,2
,
3,4 are any arbitrary constants, l
1
≠
0 and l
3
=
l
2
(
l
1
3
+
2
l
4
)
/
(
3
l
1
); then k
1
=
{
l
1
,
l
2
,
l
3
,
l
4
}
T satisfies (23), and thus
(24)
u
(
x
,
y
,
z
,
t
)
=
2
l
1
ϵ
e
k
i
·
x
1
+
ϵ
e
k
i
·
x
+
g
,
where ϵ and g are any arbitrary constants, is a 1-wave solution, that is, traveling wave solution to the (
3
+
1
)
-dimensional Jimbo-Miwa equation.By further computation, (24) could be rewritten as
(25)
u
(
x
,
y
,
z
,
t
)
=
g
-
2
l
1
1
+
ϵ
e
k
i
·
x
.
For ϵ
>
0,
(26)
2
l
1
1
+
ϵ
e
k
i
·
x
=
l
1
[
1
-
tanh
1
2
(
k
i
·
x
+
ξ
0
)
]
and, for ϵ
<
0,
(27)
2
l
1
1
+
ϵ
e
k
i
·
x
=
l
1
[
1
-
coth
⁡
1
2
(
k
i
·
x
+
ξ
0
)
]
,
where ξ
0
=
ln
⁡
|
ϵ
| is an arbitrary constant. Consequently, (24) can be rewritten as
(28)
u
(
x
,
y
,
z
,
t
)
=
g
-
l
1
tanh
1
2
(
k
i
·
x
+
ξ
0
)
,
(29)
u
(
x
,
y
,
z
,
t
)
=
g
-
l
1
coth
⁡
1
2
(
k
i
·
x
+
ξ
0
)
.
By parametric transformation, we can find that the traveling wave solution (28) is exactly the same as solution (19) in Section 2. However, (29) is a family of unbounded solutions.


### 3.1. 2-Wave Solutions to the(
3
+
1
)
-Dimensional Jimbo-Miwa Equation

To obtain the 2-wave solutions to the(
3
+
1
)
-dimensional Jimbo-Miwa equation, we now study the solutions to the following coupled algebraic equations:
(30)
k
1,1
3
k
2,1
+
2
k
4,1
k
2,1
-
3
k
1,1
k
3,1
=
0
,
k
1,2
3
k
2,2
+
2
k
4,2
k
2,2
-
3
k
1,2
k
3,2
=
0
,
(
k
1,2
-
k
1,1
)
3
(
k
2,2
-
k
2,1
)
+
2
(
k
4,2
-
k
4,1
)
(
k
2,2
-
k
2,1
)
-
3
(
k
1,2
-
k
1,1
)
(
k
3,2
-
k
3,1
)
=
0
.From (30), we get
(31)
k
3,1
=
k
1,1
3
+
2
k
4,1
3
k
1,1
k
2,1
,
(32)
k
3,2
=
k
1,2
3
+
2
k
4,2
3
k
1,2
k
2,2
,
(33)
k
3,2
-
k
3,1
=
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
3
(
k
1,2
-
k
1,1
)
×
(
k
2,2
-
k
2,1
)
.
Clearly, from (31)–(33), we have
(34)
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
3
(
k
1,2
-
k
1,1
)
(
k
2,2
-
k
2,1
)
=
k
1,2
3
+
2
k
4,2
3
k
1,2
k
2,2
-
k
1,1
3
+
2
k
4,1
3
k
1,1
k
2,1
,
from which we get
(35)
k
2,2
=
(
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
(
(
k
1,1
3
+
2
k
4,1
)
(
k
1,2
-
k
1,1
)
-
k
1,1
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
×
(
(
k
1,2
3
+
2
k
4,2
)
(
k
1,2
-
k
1,1
)
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
)
×
k
1,1
k
1,2
k
2,1
.
Now, substituting (35) into (32), we obtain
(36)
k
3,2
=
(
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
(
(
k
1,1
3
+
2
k
4,1
)
(
k
1,2
-
k
1,1
)
-
k
1,1
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
×
(
(
k
1,2
3
+
2
k
4,2
)
(
k
1,2
-
k
1,1
)
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
)
×
(
k
1,2
3
+
2
k
4,2
)
3
k
1,2
2
k
1,1
k
2,1
.
Consequently, we get the 2-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation and we can now state the following theorem.Theorem 5.
Letk
i
=
{
k
1
,
i
,
k
2
,
i
,
k
3
,
i
,
k
4
,
i
}
T, x
=
{
x
,
y
,
z
,
t
}
T, i
=
1,2. For any constants g and {
ϵ
i
,
i
=
1,2
},
(37)
u
(
x
,
y
,
z
,
t
)
=
2
k
1,1
ϵ
1
e
k
1
·
x
+
k
1,2
ϵ
2
e
k
2
·
x
1
+
ϵ
1
e
k
1
·
x
+
ϵ
2
e
k
2
·
x
+
g
is a 2-wave solution to the (
3
+
1
)-dimensional Jimbo-Miwa equation (1), where k
1,1, k
1,2, k
2,1, k
4,1, and k
4,2 are arbitrary constants and the constants k
3,1, k
2,2, and k
3,2 are determined by (31), (35), and (36), respectively.


### 3.2.
N-Wave Solutions to the (
3
+
1
)
-Dimensional Jimbo-Miwa Equation

To get theN-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation, we need to investigate the independent solution set {
k
1
,
k
2
,
…
,
k
N
} of (20), which satisfies the dispersion relations. Generally, it is very difficult to get the solution if N is greater than 3 in which case the number of (20) is N
(
N
+
1
)
/
2. However, it might be possible to solve these equations by assuming some special relations between these parameters [30, 31].Following the idea in [30, 31], suppose that k
i
=
{
a
1
k
i
n
1
,
a
2
k
i
n
2
,
a
3
k
i
n
3
,
a
4
k
i
n
4
}
T,   i
=
1,2
,
…
,
N, where a
j
′
s and n
j
′
s
,
j
=
1,2
,
3,4, are parameters to be determined later. To get the possible solutions to (20), we choose (
n
1
,
n
2
,
n
3
,
n
4
) to be (
1,2
,
4,3
) or (
1
,
-
1,1
,
3
).For the case when(
n
1
,
n
2
,
n
3
,
n
4
)
=
(
1,2
,
4,3
), substituting k
i
=
{
a
1
k
i
n
1
,
a
2
k
i
n
2
,
a
3
k
i
n
3
,
a
4
k
i
n
4
}
T, i
=
1,2
,
…
,
N, into (20) gives
(38)
a
1
3
a
2
(
k
i
-
k
j
)
3
(
k
i
2
-
k
j
2
)
+
2
a
4
a
2
(
k
i
3
-
k
j
3
)
(
k
i
2
-
k
j
2
)
-
3
a
1
a
4
(
k
i
-
k
j
)
(
k
i
4
-
k
j
4
)
=
0
.
Obviously, (38) holds for arbitrary values of k
i and k
j if and only if a
1
,
a
2
,
a
3, and a
4 satisfy the two equations
(39)
2
a
2
(
a
1
3
+
2
a
4
)
-
3
a
1
a
3
=
0
,
a
2
(
a
1
3
-
4
a
4
)
=
0
.
Solving (39), we get a
4
=
(
1
/
4
)
a
1
3, a
3
=
a
1
2
a
2, where a
1 and a
2 are free parameters. It is easy to see that k
i
=
{
a
1
k
i
,
a
2
k
i
2
,
a
3
k
i
4
,
a
4
k
i
3
}
T, i
=
1,2
,
…
,
N, satisfy the dispersion condition. Consequently, we know that the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) possesses N-wave solution for any arbitrary positive integer N. Thus, we have the following theorem.Theorem 6.
For any arbitrary positive integerN and for any different k
i, i
=
1,2
,
…
,
N,
(40)
u
(
x
,
y
,
z
,
t
)
=
2
a
1
∑
i
=
1
N
k
i
ϵ
i
e
a
1
k
i
x
+
a
2
k
i
2
y
+
a
1
2
a
2
k
i
4
z
+
(
1
/
4
)
a
1
3
k
i
3
t
1
+
∑
i
=
1
N
ϵ
i
e
a
1
k
i
x
+
a
2
k
i
2
y
+
a
1
2
a
2
k
i
4
z
+
(
1
/
4
)
a
1
3
k
i
3
t
+
g
is a family of N-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1), where ϵ
i, i
=
1,2
,
…
,
N, and a
1 and a
2 are arbitrary constants.For the case when(
n
1
,
n
2
,
n
3
,
n
4
)
=
(
1
,
-
1,1
,
3
), substituting k
i
=
{
a
1
k
i
,
a
2
k
i
-
1
,
a
3
k
i,a
4
k
i
3
}
T, i
=
1,2
,
…
,
N, into (20) gives
(41)
a
1
3
a
2
(
k
i
-
k
j
)
3
(
k
i
-
1
-
k
j
-
1
)
+
2
a
4
a
2
(
k
i
3
-
k
j
3
)
(
k
i
-
1
-
k
j
-
1
)
-
3
a
1
a
4
(
k
i
-
k
j
)
(
k
i
3
-
k
j
3
)
=
0
.
Clearly, (41) holds for arbitrary values of k
i and k
j if and only if a
1
,
a
2
,
a
3, and a
4 satisfy
(42)
2
a
2
(
a
1
3
-
a
4
)
-
3
a
1
a
3
=
0
,
a
2
(
a
1
3
+
2
a
4
)
=
0
.
Solving (42), we obtain a
4
=
-
a
1
3
/
2
,
a
3
=
a
1
2
a
2, where a
1 and a
2 are free parameters. Unfortunately, unlike the case when (
n
1
,
n
2
,
n
3
,
n
4
)
=
(
1,2
,
4,3
) above, it is easy to see that k
i
=
{
a
1
k
i
,
a
2
k
i
-
1
,
a
3
k
i
,
a
4
k
i
3
}
T, i
=
1,2
,
…
,
N, do not satisfy the dispersion condition any more. However, if we let k
0
=
{
a
1
,
a
2
,
a
3
,
a
4
}
T, then, for any i
=
1,2
,
…
,
N,
(43)
k
i
′
=
k
i
-
k
0
=
{
a
1
(
k
i
-
1
)
,
a
2
(
k
i
-
1
-
1
)
,
a
3
(
k
i
-
1
)
,
a
4
(
k
i
3
-
1
)
}
T
satisfies the dispersion relations because {
k
0
,
k
1
,
k
2
,
…
,
k
N
} is a solution set of (20).Consequently, besides the family ofN-wave solution (40), the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) has another group of N-wave solution for any arbitrary positive integer N and so we have the following theorem.Theorem 7.
For any arbitrary positive integerN and for any different k
i, i
=
1,2
,
…
,
N,
(44)
u
(
x
,
y
,
z
,
t
)

=

2
a
1

∑
i
=
1
N

(
k
i
-
1
)
ϵ
i
e
a
1
(
k
i
-
1
)
x
+
a
2
(
k
i
-
1
-
1
)
y
+
a
1
2
a
2
(
k
i
-
1
)
z
-
(
1
/
2
)
a
1
3
(
k
i
3
-
1
)
t
1
+

∑
i
=
1
N

ϵ
i
e
a
1
(
k
i
-
1
)
x
+
a
2
(
k
i
-
1
-
1
)
y
+
a
1
2
a
2
(
k
i
-
1
)
z
-
(
1
/
2
)
a
1
3
(
k
i
3
-
1
)
t
+
g
is a family of N-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1), where ϵ
i, i
=
1,2
,
…
,
N, and a
1 and a
2 are arbitrary constants.


## 3.1. 2-Wave Solutions to the(
3
+
1
)
-Dimensional Jimbo-Miwa Equation

To obtain the 2-wave solutions to the(
3
+
1
)
-dimensional Jimbo-Miwa equation, we now study the solutions to the following coupled algebraic equations:
(30)
k
1,1
3
k
2,1
+
2
k
4,1
k
2,1
-
3
k
1,1
k
3,1
=
0
,
k
1,2
3
k
2,2
+
2
k
4,2
k
2,2
-
3
k
1,2
k
3,2
=
0
,
(
k
1,2
-
k
1,1
)
3
(
k
2,2
-
k
2,1
)
+
2
(
k
4,2
-
k
4,1
)
(
k
2,2
-
k
2,1
)
-
3
(
k
1,2
-
k
1,1
)
(
k
3,2
-
k
3,1
)
=
0
.From (30), we get
(31)
k
3,1
=
k
1,1
3
+
2
k
4,1
3
k
1,1
k
2,1
,
(32)
k
3,2
=
k
1,2
3
+
2
k
4,2
3
k
1,2
k
2,2
,
(33)
k
3,2
-
k
3,1
=
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
3
(
k
1,2
-
k
1,1
)
×
(
k
2,2
-
k
2,1
)
.
Clearly, from (31)–(33), we have
(34)
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
3
(
k
1,2
-
k
1,1
)
(
k
2,2
-
k
2,1
)
=
k
1,2
3
+
2
k
4,2
3
k
1,2
k
2,2
-
k
1,1
3
+
2
k
4,1
3
k
1,1
k
2,1
,
from which we get
(35)
k
2,2
=
(
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
(
(
k
1,1
3
+
2
k
4,1
)
(
k
1,2
-
k
1,1
)
-
k
1,1
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
×
(
(
k
1,2
3
+
2
k
4,2
)
(
k
1,2
-
k
1,1
)
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
)
×
k
1,1
k
1,2
k
2,1
.
Now, substituting (35) into (32), we obtain
(36)
k
3,2
=
(
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
(
(
k
1,1
3
+
2
k
4,1
)
(
k
1,2
-
k
1,1
)
-
k
1,1
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
×
(
(
k
1,2
3
+
2
k
4,2
)
(
k
1,2
-
k
1,1
)
-
k
1,2
[
(
k
1,2
-
k
1,1
)
3
+
2
(
k
4,2
-
k
4,1
)
]
)
-
1
)
×
(
k
1,2
3
+
2
k
4,2
)
3
k
1,2
2
k
1,1
k
2,1
.
Consequently, we get the 2-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation and we can now state the following theorem.Theorem 5.
Letk
i
=
{
k
1
,
i
,
k
2
,
i
,
k
3
,
i
,
k
4
,
i
}
T, x
=
{
x
,
y
,
z
,
t
}
T, i
=
1,2. For any constants g and {
ϵ
i
,
i
=
1,2
},
(37)
u
(
x
,
y
,
z
,
t
)
=
2
k
1,1
ϵ
1
e
k
1
·
x
+
k
1,2
ϵ
2
e
k
2
·
x
1
+
ϵ
1
e
k
1
·
x
+
ϵ
2
e
k
2
·
x
+
g
is a 2-wave solution to the (
3
+
1
)-dimensional Jimbo-Miwa equation (1), where k
1,1, k
1,2, k
2,1, k
4,1, and k
4,2 are arbitrary constants and the constants k
3,1, k
2,2, and k
3,2 are determined by (31), (35), and (36), respectively.


## 3.2.
N-Wave Solutions to the (
3
+
1
)
-Dimensional Jimbo-Miwa Equation

To get theN-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation, we need to investigate the independent solution set {
k
1
,
k
2
,
…
,
k
N
} of (20), which satisfies the dispersion relations. Generally, it is very difficult to get the solution if N is greater than 3 in which case the number of (20) is N
(
N
+
1
)
/
2. However, it might be possible to solve these equations by assuming some special relations between these parameters [30, 31].Following the idea in [30, 31], suppose that k
i
=
{
a
1
k
i
n
1
,
a
2
k
i
n
2
,
a
3
k
i
n
3
,
a
4
k
i
n
4
}
T,   i
=
1,2
,
…
,
N, where a
j
′
s and n
j
′
s
,
j
=
1,2
,
3,4, are parameters to be determined later. To get the possible solutions to (20), we choose (
n
1
,
n
2
,
n
3
,
n
4
) to be (
1,2
,
4,3
) or (
1
,
-
1,1
,
3
).For the case when(
n
1
,
n
2
,
n
3
,
n
4
)
=
(
1,2
,
4,3
), substituting k
i
=
{
a
1
k
i
n
1
,
a
2
k
i
n
2
,
a
3
k
i
n
3
,
a
4
k
i
n
4
}
T, i
=
1,2
,
…
,
N, into (20) gives
(38)
a
1
3
a
2
(
k
i
-
k
j
)
3
(
k
i
2
-
k
j
2
)
+
2
a
4
a
2
(
k
i
3
-
k
j
3
)
(
k
i
2
-
k
j
2
)
-
3
a
1
a
4
(
k
i
-
k
j
)
(
k
i
4
-
k
j
4
)
=
0
.
Obviously, (38) holds for arbitrary values of k
i and k
j if and only if a
1
,
a
2
,
a
3, and a
4 satisfy the two equations
(39)
2
a
2
(
a
1
3
+
2
a
4
)
-
3
a
1
a
3
=
0
,
a
2
(
a
1
3
-
4
a
4
)
=
0
.
Solving (39), we get a
4
=
(
1
/
4
)
a
1
3, a
3
=
a
1
2
a
2, where a
1 and a
2 are free parameters. It is easy to see that k
i
=
{
a
1
k
i
,
a
2
k
i
2
,
a
3
k
i
4
,
a
4
k
i
3
}
T, i
=
1,2
,
…
,
N, satisfy the dispersion condition. Consequently, we know that the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) possesses N-wave solution for any arbitrary positive integer N. Thus, we have the following theorem.Theorem 6.
For any arbitrary positive integerN and for any different k
i, i
=
1,2
,
…
,
N,
(40)
u
(
x
,
y
,
z
,
t
)
=
2
a
1
∑
i
=
1
N
k
i
ϵ
i
e
a
1
k
i
x
+
a
2
k
i
2
y
+
a
1
2
a
2
k
i
4
z
+
(
1
/
4
)
a
1
3
k
i
3
t
1
+
∑
i
=
1
N
ϵ
i
e
a
1
k
i
x
+
a
2
k
i
2
y
+
a
1
2
a
2
k
i
4
z
+
(
1
/
4
)
a
1
3
k
i
3
t
+
g
is a family of N-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1), where ϵ
i, i
=
1,2
,
…
,
N, and a
1 and a
2 are arbitrary constants.For the case when(
n
1
,
n
2
,
n
3
,
n
4
)
=
(
1
,
-
1,1
,
3
), substituting k
i
=
{
a
1
k
i
,
a
2
k
i
-
1
,
a
3
k
i,a
4
k
i
3
}
T, i
=
1,2
,
…
,
N, into (20) gives
(41)
a
1
3
a
2
(
k
i
-
k
j
)
3
(
k
i
-
1
-
k
j
-
1
)
+
2
a
4
a
2
(
k
i
3
-
k
j
3
)
(
k
i
-
1
-
k
j
-
1
)
-
3
a
1
a
4
(
k
i
-
k
j
)
(
k
i
3
-
k
j
3
)
=
0
.
Clearly, (41) holds for arbitrary values of k
i and k
j if and only if a
1
,
a
2
,
a
3, and a
4 satisfy
(42)
2
a
2
(
a
1
3
-
a
4
)
-
3
a
1
a
3
=
0
,
a
2
(
a
1
3
+
2
a
4
)
=
0
.
Solving (42), we obtain a
4
=
-
a
1
3
/
2
,
a
3
=
a
1
2
a
2, where a
1 and a
2 are free parameters. Unfortunately, unlike the case when (
n
1
,
n
2
,
n
3
,
n
4
)
=
(
1,2
,
4,3
) above, it is easy to see that k
i
=
{
a
1
k
i
,
a
2
k
i
-
1
,
a
3
k
i
,
a
4
k
i
3
}
T, i
=
1,2
,
…
,
N, do not satisfy the dispersion condition any more. However, if we let k
0
=
{
a
1
,
a
2
,
a
3
,
a
4
}
T, then, for any i
=
1,2
,
…
,
N,
(43)
k
i
′
=
k
i
-
k
0
=
{
a
1
(
k
i
-
1
)
,
a
2
(
k
i
-
1
-
1
)
,
a
3
(
k
i
-
1
)
,
a
4
(
k
i
3
-
1
)
}
T
satisfies the dispersion relations because {
k
0
,
k
1
,
k
2
,
…
,
k
N
} is a solution set of (20).Consequently, besides the family ofN-wave solution (40), the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) has another group of N-wave solution for any arbitrary positive integer N and so we have the following theorem.Theorem 7.
For any arbitrary positive integerN and for any different k
i, i
=
1,2
,
…
,
N,
(44)
u
(
x
,
y
,
z
,
t
)

=

2
a
1

∑
i
=
1
N

(
k
i
-
1
)
ϵ
i
e
a
1
(
k
i
-
1
)
x
+
a
2
(
k
i
-
1
-
1
)
y
+
a
1
2
a
2
(
k
i
-
1
)
z
-
(
1
/
2
)
a
1
3
(
k
i
3
-
1
)
t
1
+

∑
i
=
1
N

ϵ
i
e
a
1
(
k
i
-
1
)
x
+
a
2
(
k
i
-
1
-
1
)
y
+
a
1
2
a
2
(
k
i
-
1
)
z
-
(
1
/
2
)
a
1
3
(
k
i
3
-
1
)
t
+
g
is a family of N-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1), where ϵ
i, i
=
1,2
,
…
,
N, and a
1 and a
2 are arbitrary constants.


## 4. Concluding Remarks

The dynamical system theory was employed to study the traveling wave solutions of the(
3
+
1
)
-dimensional Jimbo-Miwa equation (1). The multiwave solutions were investigated by studying the linear space of the corresponding Hirota bilinear equation. The exact formulas of two families of multiwave solutions of any order were obtained as well. Wazwaz [38] employed Hirota’s bilinear method and derived 1-wave and 2-wave solutions to this equation and stated that it is completely integrable and it admits multiple-soliton solutions of any order. In this paper we explicitly obtained two different families of N-wave solutions to the (
3
+
1
)
-dimensional Jimbo-Miwa equation (1) given by (40) and (44) and this confirmed the statement given in [38] that the Jimbo-Miwa equation (1) admits multiple-wave solutions of any order N.



---

*Source: 963852-2014-04-09.xml*