Discrete Dynamics in Nature and Society

Volume 2011, Article ID 293607, 24 pages

http://dx.doi.org/10.1155/2011/293607

## Extended Mixed Function Method and Its Application for Solving Two Classic Toda Lattice Equations

Center for Nonlinear Science Research, College of Mathematics, Honghe University, Yunnan, Mengzi 661100, China

Received 14 April 2011; Revised 26 May 2011; Accepted 15 June 2011

Academic Editor: Zhengqiu Zhang

Copyright © 2011 Weiguo Rui and Xinsong Yang. 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

The mixed function method is extended from the -dimensional space to the -dimensional one, even those forms of exact solution do not exist in -dimensional NDDEs. By using this extended method, the Toda lattice and -dimensional Toda lattice equations are studied. Some new solutions such as discrete solitary wave solutions, discrete kink and antikink wave solutions, and discrete breather soliton solutions are obtained, and their dynamic properties are discussed.

#### 1. Introduction

In [1], motivated by the structures of the exact solutions obtained from Darboux transformation, the authors constructed the mixed function method in (1 + 1)-dimensional space. Further, they studied the generalized Hybrid lattice equation and the two-component Volterra lattice equation by using this method, then they obtained some new exact solutions such as discrete solitary wave solutions and kink wave solutions. By this token, the mixed function method is as powerful as the exp method [2–7] and tanh method [8–10]. In this paper, we will extend this method from the (1 + 1)-dimensional space to the (2 + 1)-dimensional one and present some new forms on establishing different kinds of exact solutions, then this extended method will be applied to many higher-order nonlinear differential-difference equations (NDDEs). Moreover, we provide a simplified way to solve large numbers of high-power algebraic equations using a computer. Of course, these high-power algebraic equations mentioned above are all derived from high-order or high-power NDDEs. Therefore, the procedure of computation can be greatly simplified, and time is saved. We notice that this extended method is as powerful as the other symbolic computation methods such as tanh-function method, sine-cosine function method, exp-function method, Jacobian elliptic function method, ()-expansion method, and Adomian decomposition method, and these methods (please refer to [2–18] and references cited therein) are popular tools in the field of the nonlinear differential-difference equations.

We consider the following (2 + 1)-dimensional Toda lattice equation [19] and the Toda lattice equation [20] It is well known that there is another form for the (2 + 1)-dimensional Toda lattice equation and there are another two forms for the Toda lattice equation Under the transformation , (1.3) and (1.4) can be reduced to (1.1) and (1.2), respectively. For (1.1) and (1.2), there are important physical models and biological models which have been studied by many authors, see [7, 10, 21–39] and references cited therein. For example, the Toda lattice equation was used as DNA models [36–38], this equation can also be used as a model to describe the pressure pulse wave in aorta.

Recently, by using the exp-function method, Zhu [7] studied (2 + 1)-dimensional Toda lattice equation and obtained some exact solutions of exp-function type. By using tanh method, S. Zhang and H.-Q. Zhang [10] studied the same equation, they obtained some exact solutions of hyperbolic function type. By using the symbolic computation of hyperbolic tangent solutions, Baldwin et al. [11] studied both (1.1) and (1.2), and they obtained some exact solutions of tanh-function type. By using the modified hyperbolic function method, Zhi et al. [34] studied (1.2), they obtained some discrete soliton solutions of hyperbolic function type. In this paper, using the extended mixed function method, we discuss (1.1) and (1.2), and new solutions which are different from the results in [7, 10, 11, 34] are obtained.

The rest of this paper is organized as follows. In Section 2, we introduce the extended mixed function method. In Section 3, by using the extended method, we obtain new exact solutions of (1.1), and discuss their dynamic properties. In Section 4, we will obtain new exact solutions of (1.2), and also discuss their dynamic properties.

#### 2. The Extended Mixed Function Method

In this section, we begin formally introducing the extended mixed function method motivated by [1].

First, we extend the mixed function method from (1 + 1)-dimensional space to (2 + 1)-dimensional one, and later, we present some new forms of solutions on establishing different kinds of exact solutions of NDDEs. Without loss of generality, we suppose that a polynomial (2 + 1)-dimensional NDDEs and (1 + 1)-dimensional NDDEs have the following form: respectively, where in (2.1) and in (2.2), is a discrete variable, and and are two continuous variable. The general framework of the extended method for NDDEs is shown in the following.

*Step 1. *Corresponding to (2.1), as in [1], we suppose that its exact solutions has the following three expansions:
where is a positive integer which can be given by the homogeneous balance principle or the character of the idiographic NDDEs, the parameters are constants which need to be determined later, and .

Besides the above forms of solutions, we will add three new expansions. If (2.1) has no term , then we suppose their expansions as follows:
where is an arbitrary function which can be chosen as hyperbolic functions and trigonometric functions. In [1], we supposed that (2.2)'s exact solution is one of the following two kinds of expansion forms:
We add two new forms of solutions for (2.2) here
where is a positive integer which can be given by the homogeneous balance principle or the character of the idiographic NDDEs, the parameters are constants which need to be determined later, and .

*Step 2. *By using the homogeneous balance principle or according to character of the idiographic NDDEs, we determine the value of then substitute it in the expressions (2.3)–(2.12) of Step 1. Sometimes, we can directly assume that or . In [8], by using the homogeneous balance method, the balance number of the generalized hybrid equation has been obtained successfully.

*Step 3. *Substituting the presupposed solutions determined by Step 2 in the original equation (2.1) or (2.2), then setting the coefficients of all independent terms in or to zero, and we get a series of algebraic equations from which the corresponding undetermined constants are explicitly solved by the use of mathematical software *Maple* or *Mathematica*.

*Step 4. *Substituting the values of these constants given by Step 3 in the solutions presupposed by Step 2, thus, the exact solutions of the original equation (2.1) or (2.2) are obtained finally.

#### 3. New Exact Solutions of (2 + 1)-Dimensional Toda Lattice Equation and Their Dynamic Properties

In this section, by using the extended method shown in the Section 2, we discuss exact solutions of (2 + 1)-dimensional Toda lattice equation (1.1) and their dynamic properties.

##### 3.1. Exact Solutions of the Form (2.3)

Suppose that (1.1) has exact solutions of the form (2.3). By using the balance procedure, we easily obtain . Thus, we let where , , ,, ,,, are constants to be determined later, and . After substituting (3.1) into (1.1), multiplying both sides by the common denominator and dividing both sides by the common factor , it follows: where

In (3.2), setting the coefficients of all independent terms ,, to zero, we get a series of algebraic equations as follows: Equation (3.4) is a group of high-power algebraic equations, so the computational load may be heavy when we use computer to solve it. Notice that , we let As a result, the equation becomes , that is, , which is a low-power equation. Similarly, we obtain Notice that , thus, (3.4) can be reduced to a group of low-power algebraic equations as follows: Solving (3.7) yields and are arbitrary nonzero constants. The process of computation can be greatly simplified by the above settings, and time is saved.

From (3.1), (3.8), and (3.9), we obtain two exact solutions of (1.1) as follows: where are given by (3.9).

In fact, (3.10) and (3.11) are solutions of local discretization because the is a discrete variable while the are continuous variables. Here we still call them discrete exact solutions. In order to describe the dynamic properties of these two discrete soliton solutions intuitively, we plot their profile figures for some fixed parametric values, as shown in Figure 1. Setting , when , the solution (3.10) shows a shape of discrete kink soliton, see Figure 1(a); when , the solution (3.10) shows another shape of discrete kink soliton, see Figure 1(b). The shape shown by Figure 1(b) has jumping phenomenon, but this is not the case shown by Figure 1(a). Setting , when , the solution (3.11) shows a shape of discrete antikink soliton, see Figure 1(c); when , the solution (3.11) shows another shape of discrete antikink soliton, see Figure 1(d). From Figures 1(c) to 1(d), the waveforms have both discrete and continuous character, and the soliton is discrete along the -axes and is continuous along the -axes.

##### 3.2. Exact Solutions of the Form (2.5)

Suppose that (1.1) has exact solution of the form (2.5). By using the balance procedure, we obtain . For simplicity, here we only consider the case of . Suppose that where , , , , , , are constants to be determined later, and , the is an arbitrary constant. As in Section 3.1, after substituting (3.12) into (1.1), multiplying both sides of the result by the common denominator , it follows: where In (3.13), setting the coefficients of all independent exp-function terms to zero, it follows Let . Then Thus, (3.15) can be reduced to the following equations: Solving (3.17), we obtain and are arbitrary nonzero constants with .

Based on (3.12), (3.18), (3.19), and (3.20), we obtain four exact solutions of (1.1) as follows: where are given by (3.19), (3.20), and .

The dynamic properties of solutions (3.21), (3.22), (3.23), and (3.24) are similar to those of (3.10) and (3.11). As an example, we plot two profile figures of the solutions (3.21) and (3.22), see Figure 2.

##### 3.3. Exact Solutions of the Form (2.6)

Suppose that (1.1) has exact solution of the form (2.6). By using the balance procedure, we obtain . So, we suppose that (1.1)'s solution has the following form: where , , , , , , are constants to be determined later, is an arbitrary constant, and is an arbitrary function. As in Section 3.1, after substituting (3.25) into (1.1), multiplying both sides of the result by the common denominator , it follows: where

In (3.26), setting the coefficients of all independent exp-function terms to zero, we obtain Directly solving (3.28), we obtain and are arbitrary constants. Thus, we obtain a class of exact solutions of (1.1) as follows: where is an arbitrary function, for which let us set to be , , , , , ; respectively, we obtain a series of soliton solutions and breather solutions as follows:

In order to describe the dynamic properties of the above solutions intuitively, as examples, we draw the profile figures of the solutions (3.31), (3.34), and (3.35) in Figure 3. Setting , the solution (3.31) shows a shape of discrete soliton, see Figure 3(a); the solutions (3.34) and (3.35) show two shapes of discrete breather oscillations, see Figures 3(c) and 3(e). Setting , the solution (3.31) shows a shape of continuous double soliton, see Figure 3(b); the solutions (3.34) and (3.35) show two shapes of continuous breather oscillations, see Figures 3(d) and 3(f). From Figure 3, the solutions (3.30)–(3.35) have diplex dynamic characters, that is, the discrete character and the continuous character coexist. If the parameter is fixed, the profiles show discrete dynamic character, see Figures 3(a), 3(c), and 3(e); while the parameter is fixed, the profiles show continuous dynamic character, see Figures 3(b), 3(d), and 3(f).

##### 3.4. Exact Solutions of the Form (2.7)

Suppose that (1.1) has exact solution of the form (2.7). By using the balance procedure, we obtain . So, we suppose that (1.1) has exact solution as the following form: where , , , , , are constants to be determined later, is an arbitrary constant, and is an arbitrary function. As in Section 3.3, by using (3.36) and (1.1), we can obtain the exact solution of (1.1); here, we omit those processes and directly give the result as follows: Similarly, setting to be , respectively, we obtain a series of soliton solutions and breather solutions as follows:

The dynamic properties of solutions (3.38)–(3.43) are similar to those of solutions (3.30)–(3.35). As examples, we plot six profile figures of the solutions (3.39), (3.42), and (3.43), see Figure 4. Setting , the solution (3.39) shows a shape of discrete soliton, see Figure 4(a); the solutions (3.42) and (3.43) show two shapes of discrete breather oscillations, see Figures 4(c) and 4(e). Setting , (3.39) shows a shape of continuous double soliton, see Figure 4(b); the solutions (3.42) and (3.43) show two shapes of continuous breather oscillations, see Figures 4(d) and 4(f).

##### 3.5. Exact Solutions of the Form (2.8)

Suppose that (1.1) has exact solution of the form (2.8). By using the balance procedure, we obtain . So, we suppose that (1.1) has exact solution as the following form: where , ,,, , are constants to be determined later, is an arbitrary constant, and is an arbitrary function. As in Section 3.3, by using (3.44) and (1.1) and the expatiatory computation, we obtain . So we obtain a family of exact solutions with arbitrary function as follows: Similarly, assume the same functions as in Section 3.3, we obtain a series of soliton solutions and breather solutions as follows:

The dynamic properties of solutions (3.46)–(3.51) are similar to those of solutions (3.30)–(3.35). As examples, we plot six profile figures of the solutions (3.46), (3.47), and (3.48), see Figure 5. Setting , the solution (3.46) shows a shape of discrete kink soliton, see Figure 5(a); the solution (3.47) shows a shape of discrete soliton, see Figure 5(c); the solution (3.48) shows a shape of discrete breather oscillations, see Figure 5(e). Setting , the solution (3.46) shows a shape of continuous double kink soliton, see Figure 5(b); the solution (3.47) shows a shape of continuous double soliton, see Figure 5(c); the solution (3.48) shows a shape of continuous breather oscillations, see Figure 5(e).

#### 4. New Exact Solutions of Toda Lattice Equation and Its Dynamic Properties

In this section, using the extended method offered in Section 2, we discuss the exact solutions of the Toda lattice equation (1.2).

##### 4.1. Exact Solutions of the Form (2.9)

Suppose that (1.2) has exact solution of the form (2.9) as follows: After substituting (4.1) in (1.2), multiplying both sides by the common denominator , it follows: where In (4.2), setting the coefficients of all independent exp-function terms to zero, we obtain Let , then Thus, (4.4) can be reduced to Solving (4.6) yields and are arbitrary constants. From (4.1), (4.7), we obtain discrete kink or antikink soliton solution of (1.2) as follows: Similarly, by using the same method, we can obtain exact solutions as the form of (2.10); here, we omit them.

##### 4.2. Exact Solutions of the Form (2.11)

Suppose that (1.2) has exact solution of the form (2.11) as follows: where , ,, , , are constants to be determined later and . After substituting (4.9) in (1.2), multiplying both sides by the common denominator , it follows: where

,

,

,

,

,