We introduce the variational iteration method for solving the generalized Degasperis-Procesi equation. Firstly, according to the variational iteration, the Lagrange multiplier is found after making the correction functional. Furthermore, several approximations of which is converged to are obtained, and the exact solutions of Degasperis-Procesi equation will be obtained by using the traditional variational iteration method with a suitable initial approximation . Finally, after giving the perturbation item, the approximate solution for original equation will be expressed specifically.

1. Introduction

The theory of soliton has extensive applications in physics, mechanics, and combustion science. In recent years, many researchers studied the soliton theory in the fields of shock wave [1, 2], light scattering, quantum mechanics, atmospheric physics, neural networks, explosion, and combustion [3]. There are many new methods for searching the soliton solution of nonlinear evolution equations such as hyperbolic tangent function method [4], the homogeneous balance methods [5], Jacobi elliptic function expansion method [3], and pseudo-spectral method [6].

The variational iteration method (VIM) was developed, in 1999, by He [713]. The VIM gives rapidly convergent successive approximations of the exact solution if such a solution exists; otherwise, a few approximations can be used for numerical purposes. The Adomian decomposition method suffers from the complicated computational work needed for the derivation of Adomian polynomials for nonlinear terms. The VIM has no specific requirements, such as linearization, small parameters for nonlinear operators. Therefore, the VIM can overcome the foregoing restrictions and limitations of perturbation techniques, so that it provides us with a possibility to analyze strongly nonlinear problems. On the other hand, the VIM is capable of greatly reducing the size of calculation while still maintaining high accuracy of the numerical solution [14]. Moreover, the power of the method gives it a wider applicability in handling a huge number of analytical and numerical applications. The VIM was successfully applied to study a variety of differential equations. It is based on Lagrange multiplier, and it has the merits of simplicity and easy execution. As a result, it has been proved by many authors to be a powerful mathematical tool for addressing various kinds of linear and nonlinear problem. For example, this method was used for solving nonlinear wave equations and the Laplace equation by Wazwaz [14]. The VIM for solving linear systems of ODEs with constant coefficients was studied by Khojasteh Salkuyeh [15]. Helmholtz equation was researched by Momani and Abuasad [16]. Geng [17] introduced the piecewise VIM for solving Riccati differential equation. Fractional vibration equation was researched by Das [18]. Furthermore, higher order boundary value problems were researched by Xu [19],Noor, and Mohyud-Din [20]. Noor et al. [21] applied a modified He’s variational iteration method for solving singular fourth-order parabolic partial differential equations. The proposed modification is made by introducing He’s polynomials in the correction functional. Ghorbani and Saberi-Nadjafi [22] modified the VIM by constructing an initial trial function without unknown parameters. Sevimlican [23] constructed approximate Green’s function for a vector equation for the electric field by using VIM.

In this paper, we are concerned with the variational iterations method for solving the generalized Degasperis-Procesi equation. As a review, we will recall the VIM briefly in Section 2.

2. Variational Iteration Method

In this section, the basic concepts of variational iteration method (VIM) are introduced. Here, a description of method [715] is given to handle the general nonlinear problem. Consider the differential equation of the form where is a linear operator, is a nonlinear operator, and is the inhomogeneous term. According to He’s variational iteration method, we can construct a correction functional for (1) as follows: where is a general Lagrange multiplier, which can be identified optimally via variational theory [12, 24]. Here is considered as a restricted variation [14, 25] which means ; the subscript denotes the th approximations. The successive approximations , of the solution , can be obtained after using the obtained Lagrange multiplier and the zeroth approximation , which are selected from any function that satisfies the initial conditions. With determined, several approximations , follow. Consequently, the exaction solution may be obtained as In fact, the VIM depends on the suitable selection of the initial approximation . Moreover, we use a well-known, powerful tool to prove the convergence of the sequence obtained via the VIM and its rate. It is the Banach’s fixed point theorem that follows.

Theorem 1 (Banach’s fixed point theorem). Assume that is a Banach space and is a nonlinear mapping, and suppose that for some constant . Then has a unique fixed point. Furthermore, the sequence with an arbitrary choice of converges to the fixed point of .

According to Theorem 1, for the nonlinear mapping a sufficient condition for the convergence of the variational iteration method is strict contraction of . Furthermore, the sequence (2) converges to the fixed point of which is also the solution of problem (1). Some modifications to prove the convergence speed and to lengthen the interval of convergence for VIM series solution are suggested in [17, 2630].

3. The Variational Iteration of Generalized Degasperis-Procesi Equation

Degasperis and Procesi consider the following family of third-order dispersive conservation laws [31], where , , , , , and are real constants. In this family, only three equations satisfy asymptotic integrability conditions [31]. That is, if , , , , , and , (8) is the KdV equation If , , , , , and , (8) is the Camassa-Holm equation If , , , , , and , (8) is the Degasperis-Procesi equation It should be mentioned that both C-H and D-P equations are derived as members of a one-parameter family of asymptotic shallow water approximations to the Euler equations. It shows that the two equations are physically relevant; otherwise, the D-P equation would be of purely theoretical interest.

Variational iteration method for KdV-Burgers and Lax’s seventh-order KdV equations has been studied by Soliman [32]. In this paper, we consider the generalized Degasperis-Procesi equation that was proposed in [33]. is the generalized perturbation item. We suppose is a sufficiently smooth function of the variable.

Step 1. Make the independent variable transformation: Here, is an arbitrary complex number. is wave number; is wave velocity. Substituting (13) into (12), we have Here, is the derivative of with respect to ; that is, . .

Step 2. From [34], we find the special solution, when is identical to :

Remark 2. Notice that , , , , and , where . These solutions contain the other four types of forms named , , , and .

Step 3. Make the correction functional
Here, , , , and are considered as a restricted variation [35]. That is,

Step 4. Under the above condition, make the correct functional stationary with respect to ; noticing that , we have
For arbitrary , from the above relation, we obtain the Euler-Language equation:
Solve (19), we derive
Substituted (20) into (16), we have the integration form:

From the above solution procedure, we can see that the approximate solutions converge to its exact solution. That is, , is the approximate solution with arbitrary degree of accurate solitary wave of Degasperis-Procesi equation.

Step 5. Calculation of the approximate solution.
According to the integration form (21),we can calculate the approximate solution. Firstly, let (15) be the zero-order approximate solution: Substitute (15) into (21). We obtain the one-order approximate solution in which .
Then, substitute (23) into (21). We can obtain the second-order approximate solution :
Using the same method, we can get the higher order approximate solution.

4. The Optical Soliton Perturbation Solution and the Numerical Example

Specially, we set Then (12) change to We obtain the zero-order and the one-order approximate solution of (20)

We set we will obtain the following numerical example. See Tables 1 and 2.

5. Conclusion

By the analysis of structure on the left side of (8) and the properties of about the variable and the analytical variational iteration formula, we can prove that the sequence of functions of decided by (21) is uniform convergence. So the limit function of is the solution of the equation. Moreover, the zero-order approximate solution is the soliton of (8) in which ; it should be specially pointed out that the more accurate the identification of the multiplier is, the faster the approximations converge to their exact solution.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The work is supported by the Youth Foundation of National Natural Science Foundation of China (Grant no. 71101072).