Nonlinear phenomena occur widely in various scientific fields, such as fluid mechanics, solid state physics, etc. As an important mathematical model to describe nonlinear phenomena, nonlinear partial differential equation has been widely concerned by many scholars in recent years, and the exact solutions of the PDE, including the soliton solution and wave solution, has always been a subject of interest to mathematical physicists. Recently, the methods of solving nonlinear partial differential equations mainly include Backlund-transformation method, non-locally symmetric method, Lie group method, and invariant subspace method, etc.
The invariant subspace method was first proposed by Galaktionov et al., and then extended by many scholars and widely applied. For example, precise solutions of Hunter-Saxton equation and compressible Euler equation were obtained. This method was derived from the Lie symmetry analysis, and it is related to the conditional Lie-Bäcklund symmetry method   and the differential constraint method; its key step is to transform a nonlinear PDE into ODE systems based on the invariant subspaces . Then, we use the invariant subspaces we have obtained to construct the exact solutions of the nonlinearly PDE equations. So it is a dynamical system method by nature .
The outstanding feature of the invariant subspace method is its wide application range; it is also an algorithm that can construct more solutions and similar solutions for nonlinear partial differential equations.
However, there are still many problems in this method, such as how to construct more exact solutions according to different nonlinear equations. Furthermore, how to improve the efficiency of solving ordinary differential equations, such as with the help of Maple program, should be further studied.
In 1993, Camassa and Holm obtained a new class of completely integrable nonlinear shallow water wave equations Camassa-Holm (CH) equation
Equation (1.1) has a double Hamilton structure and is completely integrable. When , it has smooth isolated wave solutions, when , it has a peaked solitary wave solution the form of (where c is the velocity).
Dullin, Gottwald and Holm derived the DGH equation from a class of shallow water wave equations of non-local asymptotic form 
where is the linear dispersive term, when , Equation (1.2) is the Camassa-Holm (CH) equation, when is the dispersive term , Equation (1.2) becomes the generalized strongly dispersive DGH equation
where , is an unknown function, and are arbitrary constants. Equation (1.3) describes the unidirectional propagation of surface waves in shallow water.
In recent years, many scholars have also done a lot of research on DGH equation. In 2007, Tian Lixin studied the global existence of the solution and Blow-up phenomenon and Hamilton structure of DGH equation. By using the qualitative analysis of planar autonomous system method, Guo Bolin and Liu Zhengrong obtained the spike isolated wave solutions of DGH equation. Compared with other methods, the invariant subspace method can obtain many new exact solutions of the generalized strongly dispersive DGH equation, including rational function solutions, which are different from symmetric reduction solutions and soliton solutions    . The purpose of the current paper is to provide more exact solutions to the generalized strongly dispersive DGH Equation (1.3) by using the invariant subspace method. Firstly, we transform nonlinear partial differential Equation (PDE) into ordinary differential Equation (ODE) systems by using the invariant subspace method. Secondly, combining with the dynamical system method, we use the invariant subspaces which have obtained to construct the exact solutions of the generalized strongly dispersive DGH equation. Finally, the figures of the exact solutions are given.
This paper is organized as follows: In Section 2, the invariant subspace method is introduced briefly. In Section 3, all invariant subspaces and their basis functions are given. In Section 4, by using the results given in Section 3, the rational function solution, trigonometric function solution and exponential function solution of the Equation (1.3) are obtained. In Section 5, the conclusion and prospect are given.
In the following, we use G to represent the m-order differential operator, and by using the invariant subspace method, giving the nonlinear evolution equation determined by the differential operator and its corresponding invariant subspaces, and the basis functions in the invariant subspaces.
2. Invariant Subspace Method
Firstly, we consider a nonlinear evolution equation as follows:
where , G is a m-th order differential operator, if have p functions are linearly independent, and the n-dimensional linear space is  
If G satisfies , that is, for any constant , satisfy
The invariant subspace is allowed by the operator G, so the nonlinear Equation (2.1) have the solution as follows
where satisfy n-dimensional dynamical system  
Suppose the subspace composed of solutions to the following linear ODE
then the invariant condition with respect to G takes the form
Theorem 1. If the subspace is invariant under a nonlinear differential operator G of order m, then .
3. Invariant Subspace of Equation (1.3)
Firstly, Equation (1.3) can be written in the form of general evolution Equation (2.1)
the nonlinear operator G on the right side of Equation (3.1) is
then, by the invariant subspace theorem 1, we have , in what follows, the invariant subspaces of Equation (1.3) are given.
3.1. p = 2
In this case, the invariant condition of reads
where , substituting Equation (3.2) into Equation (3.3), we have a polynomial in
take the coefficients of the above polynomial to be zero, we obtain the overdetermined equations for solving
By solving the above equations, we can get three sets of solutions as follows:
According to the above three solutions, the two-dimensional invariant subspaces allowed by the differential operator of the Equation (3.2) are
3.2. p = 3
In this case, the invariant condition of reads
where , substituting the Equation (3.2) into the Equation (3.3), we have a polynomial in , and take its coefficients to be zero, we obtain the overdetermined equations for solving
By solving the above equations, we get a unique set of solutions as follows
Thus, we obtain the three-dimensional invariant subspace is
From the invariant subspace obtained when n = 2 and n = 3, we can see that the invariant subspace method is only related to the nonlinear terms of the equation. According to the similar calculation steps above, we can get the invariant subspaces when p = 4, 5, 6. But, when p = 7, there is no solution.
3.3. p = 4
The four-dimensional invariant subspace allowed by the differential operator of the Equation (3.2) is
3.4. p = 5
The five-dimensional invariant subspace allowed by the differential operator of the Equation (3.2) is
3.5. p = 6
The six-dimensional invariant subspace allowed by the differential operator of the Equation (3.2) is
4. Exact Solutions and Figures of Equation (1.3)
In this section, combining the invariant subspace method with the dynamical system method, by using the results obtained in Section 3, construing the exact solution of the Equation (1.3) in invariant subspaces. We can divide the solutions of generalized strongly dispersive DGH equation into three types: polynomial function solutions, exponential function solutions and trigonometric function solutions. Then, the figures of the exact solutions are given.
4.1. Polynomial Solution of the Equation (1.3)
When n = 2, for
where are arbitrary constants, from in Section 3, we get .
So, we can suppose that Equation (4.1) has the solution as follows:
Substituting (4.2) into Equation (4.1), we get the equation for , and take its coefficients to be zero, we get
By solving Equation (4.3), we get
The Polynomial solution of Equation (4.1) are obtained
Figure 1(a) is the rational solution of the Equation (4.1) when , , , and , .
Figure 1(b) is the rational solution of the Equation (4.1) when , and , .
4.2. Trigonometric Solution and Exponential Solution of Equation (1.3)
When n = 3, for
Figure 1. Polynomial solution of Equation (4.1).
where are arbitrary constants. Suppose that Equation (4.4) has exact solutions in the form of trigonometric and exponential functions:
substituting (4.5) into Equation (4.4), we get the equation for , and take its coefficients to be zero, we get
By solving Equation (4.3), we get
So, the solution in form of the trigonometric function of Equation (4.4) is
where are arbitrary constants, and .
Substituting (4.6) into Equation (4.4), we get the equation for , and take its coefficients to be zero, we get
By solving Equation (4.8), we get
Figure 2. Trigonometric solution of Equation (4.4).
Figure 3. Exponential solution of Equation (4.4).
So, the trigonometric solution of Equation (4.4) is
where are arbitrary constants, and , .
Firstly, by using the invariant subspace method, the invariant subspaces of the generalized strongly dispersive DGH equation are obtained. Secondly, we selected the partial invariant subspaces, and we get an exact solution other than the soliton of the equation, such as the solutions in form of polynomial function, trigonometric function and exponential function of the generalized strongly dispersive DGH equation. Meanwhile, the invariant subspace method derived from symmetry analysis, this method also can be used to solve the nonlinear PDEs. And symmetry analysis is an invariance analysis to some extent. The above theory provides a direction for future research. Moreover, the results in the present paper are verified by the maple procedure.
The authors would like to thank the Editor and Reviewers for their valuable comments.
 Ji, L.N. and Qu, C.Z. (2013) Conditional Lie-Bäcklund Symmetries and Invariant Subspaces to Nonlinear Diffusion Equations with Convection and Source. Studies in Applied Mathematics, 131, 266-301.
 Qu, C.Z. and Ji, L.N. (2013) Invariant Subspaces and Conditional Lie-Bäcklund Symmetries of Inhomogeneous Non-Linear Diffusion Equations. Studies in Applied Mathematics, 56, 2187-2203.
 Chang, L.N., Liu, H.Z. and Xin, X.P. (2020) Invariant Subspace Classification and Exact Explicit Solutions to a Class of Nonlinear Wave Equation. Qualitative Theory of Dynamical Systems, 19, Article No. 65.
 Galaktionov, V.A. (1995) Invariant Subspaces and New Explicit Solutions to Evolution Equations with Quadratic Nonlinearities. Proceedings of the Royal Society of Edinburgh Section A: Mathematics, 125, 225-246.
 Galaktionov, A. and Svirshchevskii, S.R. (2007) Exact Solutions and Invariant Subspaces of Nonlinear Partial Differential Equations in Mechanics and Physics. Chapman and Hall, London.
 Qu, C.Z. and Zhu, C.R. (2009) Classification of Coupled Systems with Two-Component Nonlinear Diffusion Equations by the Invariant Subspace Method. Journal of Physics A: Mathematical and Theoretical, 42, 1-27.
 Liu, H.Z. and Li, J.B. (2014) Symmetry Reductions, Dynamical Behavior and Exact Explicit Solutions to the Gordon Types of Equations. Journal of Computational and Applied Mathematics, 257, 144-156.
 Shen, S.F. and Qu, C.Z. (2012) Maximal Dimension of Invariant Subspaces to Systems of Nonlinear Evolution Equations. Chinese Annals of Mathematics, Series B, 33, 161-178.
 Guo, Z.H., Liu, X.C., Liu, X.X. and Qu, C.Z. (2019) Stability of Peakons for the Generalized Modified Camassa-Holm Equation. Journal of Differential Equations, 266, 7749-7779.
 Liu, H.Z. (2018) Invariant Subspace Classification and Exact Solutions to the Generalized Nonlinear D-C Equation. Applied Mathematics Letters, 83, 164-168.