#### Abstract

A general algebraic method based on the generalized Jacobi elliptic functions expansion method, the improved general mapping deformation method, and the extended auxiliary function method with computerized symbolic computation is proposed to construct more new exact solutions for coupled Schrödinger-Boussinesq equations. As a result, several families of new generalized Jacobi elliptic function wave solutions are obtained by using this method, some of them are degenerated to solitary wave solutions and trigonometric function solutions in the limited cases, which shows that the general method is more powerful than plenty of traditional methods and will be used in further works to establish more entirely new solutions for other kinds of nonlinear partial differential equations arising in mathematical physics.

#### 1. Introduction

Nonlinear partial differential equations (NLPDEs) are widely used to describe complex physical phenomena arising in the world around us and various fields of science. The investigation of exact solutions of NLPDEs plays an important role in the study of these phenomena such as the nonlinear dynamics and the mechanism behind the phenomena. With the development of soliton theory, many powerful methods for obtaining exact solutions of NLPDEs have been presented, such as homotopy perturbation method [1], nonperturbative method [2], homogeneous balance method [3], Bäcklund transformation [4], Darboux transformation [5], extended tanh-function method [6], extended -expansion method [7], method [8], exp-function method [9], sine-cosine method [10], Jacobi elliptic function method [11], extended Riccati equation rational expansion method [12], extended auxiliary function method [13], and other methods [14, 15].

In [16, 17], Hong proposed a generalized Jacobi elliptic functions expansion method to obtain generalized exact solutions of NLPDEs. In [18], Hong and Lu proposed an improved general mapping deformation method which is more general than many other algebraic expansion methods [19, 20]. The solution procedure of this method, by the help of Matlab or Mathematica, is of utmost simplicity, and this method can be easily extended to all kinds of NLPDEs. In this work, we will propose the general algebraic method which contained the two methods [16–18] to obtain several new families of exact solutions for the coupled Schrödinger-Boussinesq equations.

#### 2. Summary of the General Algebraic Method

Consider a given nonlinear evolution equation with one physical field in two variables and We seek the following formal solutions of the given system by a new intermediate transformation: where are constants to be determined later. are arbitrary functions with the variables and . The parameter can be determined by balancing the highest order derivative terms with the nonlinear terms in (1). And is a solution of the following ordinary differential equation: Substituting (3) and (2) into (1) and setting the coefficients of and to zero yield a set of algebraic equations for , , and . Using the Mathematica to solve the algebraic equations and substituting each of the solutions of the set, that is, each of the expressions of into (2), we can get the solutions of (1). In order to obtain some new general solutions of (3), we assume that (3) has the following solutions: where are functions of to be determined later and the four functions , and are expressed as follows: where , and are arbitrary constants which ensure the denominator unequal to zero, and is a solution of the following equations: where “” denotes , “” denotes , , and are arbitrary constants, and the four functions , and satisfy the following relations: And , and satisfy one of the following relations at the same time.

*Family 1. *When

*Family 2. *When

*Family 3. *When

*Family 4. *When

Substituting (4), (7) along with (8a)–(8d) into (3) separately yields four families of polynomial equations for , and . Setting the coefficients of , and to zero yields a set of overdetermined differential equations (ODEs) in , and , and , solving the ODEs by Mathematica and Wu elimination, we can obtain many exact solutions of (3) according to (4), (5), (6). If we let , , , , , , , , and , we have ; our method contains the improved general mapping deformation method [18].

*Remark 1. *Our method proposed here is more general than the method [8], the extended Riccati equation rational expansion method [12], the extended auxiliary function method [13], the generalized Jacobi elliptic functions expansion method [16, 17], and many other algebra expansion methods [6, 7, 11, 18–21].

*Remark 2. *Equations (2) and (3) can be extended to the following forms:
where can be an arbitrary positive integer. is usually a positive integer. If is a fraction or a negative integer, we make the following transformation: (a)when is a fraction, we let , then return to determine the balance constant again; (b)when is a negative integer, we suppose , then return to determine the balance constant again.

*Remark 3. *Notice that

We find a meanful conclusion that this general method implies a BT of (1) with the compatible conditions (4), (5), (6), (7), and (8a)–(8d).

In the following, we will use this method to solve the Schrödinger-Boussinesq equations.

#### 3. Exact Solutions to the Coupled Schrödinger-Boussinesq Equations

We consider the coupled Schrödinger-Boussinesq equations [21–28] These equations are known to describe various physical processes in laser and plasma, such as formation, Langmuir field amplitude, intense electromagnetic waves, and modulational instabilities [22–25]. The problem of the complete integrability of this system has been studied by Chowdhury et al. from the point of view of Painlevé analysis [25]. The solitary wave solutions for system (11) have been obtained in [26, 27]. The Jacobi doubly periodic wave solutions and a range of other solutions for this system have been investigated in [21, 28]. We are interested in searching new generalized Jacobi elliptic function solutions for (11) by using our method.

We consider the following transformations: where and are constants to be determined later and and are arbitrary constants.

Substituting (12) and (13) into (11) and integrating the second equation of system (11) twice, we have By balancing the highest derivative term with the nonlinear terms in (14), we obtain . Therefore, we assume that (14) have the following solutions: Substituting (3), (13), and (15) into (14) and setting the coefficients of and to zero yield a set of over-determined equations (ODEs) for , and . After solving the ODEs by Mathematica, we could determine the following solutions.

*Family 1. *

*Family 2. *

Substituting (4), (5), (6), (7) along with (8a)–(8d) and (16) into (14) separately yields an ODEs; after solving the ODEs by Mathematica and Wu elimination, we can obtain the following solutions of (14) according to (4), (5), (6), (16).

*Case 1. *

*Case 2. *

*Case 3. *

*Case 4. *

*Case 5. *

*Case 6. *

*Case 7. *

We obtain the following solutions of (11) according to (4), (5), (12), (13), (15), and (16) and Cases 1–7: where is an arbitrary constant, , and are defined as Cases 1–7, and and are defined as follows:

With the similar process,substituting (4), (5), (6), and (7) along with (8a)–(8d) and (17) into (14) separately yields ODEs; after solving the ODEs by Mathematica and Wu elimination, we can obtain the following solutions of (11) according to (4), (5), (12), (13), (15), and (17) and the solutions of (3) mentioned in [29]:

*Remark 4. *If we let , , , , and , then we have
If we let , , or , then solution is in full agreement with the solution , and , mentioned in [21].

*Remark 5. *The seven types of explicit solutions we obtained here to (11) are not shown in the previous literature to our knowledge. They are new exact solutions of (11). Notice that , , and when the modulus , and , , and when the modulus . Solutions are degenerated to solitary wave solutions when the modulus and to triangular functions solutions when the modulus .

We can give the numerical simulation of to show their physical properties (see Figures 1, 2, and 3).

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

#### 4. Conclusion

In this paper, we succeed to propose a general algebraic method approach for finding new exact solutions of the nonlinear evolution equations. By using this method and computerized symbolic computation, we have found abundant new exact solutions for the coupled Schrödinger-Boussinesq equations (11). More importantly, our method is much simpler and powerful to find new solutions to various kinds of nonlinear evolution equations, such as KdV equation, Boussinesq equation, and zakharov equation. We believe that this method should play an important role for finding exact solutions in the mathematical physics.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 61070231), the Outstanding Personal Program in Six Fields of Jiangsu (Grant no. 2009188), the Graduate Student Innovation Project of Jiangsu Province (Grant no. CXLX13_673), and the General Program of Innovation Foundation of Nanjing Institute of Technology (Grant no. CKJB201218).