In this paper, an extended functional transformation is given to solve some nonlinear evolution equations. This function, in fact, is a solution of the famous KdV equation, so this transformation gives a transformatio...In this paper, an extended functional transformation is given to solve some nonlinear evolution equations. This function, in fact, is a solution of the famous KdV equation, so this transformation gives a transformation between KdV equation and other soliton equations. Then many new exact solutions can be given by virtue of the solutions of KdV equation.展开更多
In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination...In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.展开更多
In this paper, we extend the mapping transformation method through introducing variable coefficients.By means of the extended mapping transformation method, many explicit and exact general solutions with arbitrary fun...In this paper, we extend the mapping transformation method through introducing variable coefficients.By means of the extended mapping transformation method, many explicit and exact general solutions with arbitrary functions for some nonlinear partial differential equations, which contain solitary wave solutions, trigonometric function solutions, and rational solutions, are obtained.展开更多
This paper is to investigate the extended(2+1)-dimensional Konopelchenko-Dubrovsky equations,which can be applied to describing certain phenomena in the stratified shear flow,the internal and shallow-water waves, plas...This paper is to investigate the extended(2+1)-dimensional Konopelchenko-Dubrovsky equations,which can be applied to describing certain phenomena in the stratified shear flow,the internal and shallow-water waves, plasmas and other fields.Painleve analysis is passed through via symbolic computation.Bilinear-form equations are constructed and soliton solutions are derived.Soliton solutions and interactions are illustrated.Bilinear-form Backlund transformation and a type of solutions are obtained.展开更多
In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Ca...In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Calogero-Degasperis (VCCD) equation. New soliton and periodic solutions of many types are obtained. Furthermore, the soliton propagation is discussed under the effect of the variable coefficients.展开更多
In this paper,we propose an extended hybrid carrier system based on the weighted fractional Fourier transform to ensure the reliability of wireless communication.The proposed scheme improves the dispersion and compens...In this paper,we propose an extended hybrid carrier system based on the weighted fractional Fourier transform to ensure the reliability of wireless communication.The proposed scheme improves the dispersion and compensation capabilities of the hybrid carrier system for channel fading through the design of the signal power distribution,which has greatly reduced the probability of high-power distortion of the signal and improved the bit error rate performance as a result.Theoretical analysis has shown the superiority of the extended hybrid carrier system.With a lower cost of computational complexity increment,the proposed scheme obtains a performance improvement without occupying additional time-frequency physical resources.Compared with the existing hybrid carrier scheme,numerical simulation results have shown that the proposed extended hybrid carrier scheme has better anti-fading performance under the doubly-selective channel and improves the reliability of the wireless communication system effectively.展开更多
Jordan's lemma can be used for a wider range than the original one. The extended Jordan's lemma can be described as follows. Let f(z) be analytic in the upper half of the z plane (Imz≥0), with the exception o...Jordan's lemma can be used for a wider range than the original one. The extended Jordan's lemma can be described as follows. Let f(z) be analytic in the upper half of the z plane (Imz≥0), with the exception of a finite number of isolated singularities, and for P>o, if then where z=Rei and CR is the open semicircle in the upper half of the z plane.With the extended Jordan's lemma one can find that Laplace transform and Fourier transform are a pair of integral transforms which relate to each other.展开更多
Wavelet transform is used to analyze the scaling rule of temperature data (passive scalar) in Rayleigh Bénard convection flow from two aspects. By utilizing the method of extended self similarity (ESS), one can f...Wavelet transform is used to analyze the scaling rule of temperature data (passive scalar) in Rayleigh Bénard convection flow from two aspects. By utilizing the method of extended self similarity (ESS), one can find the obtained scaling exponent agrees well with the one obtained from the temperature data in a experiment of wind tunnel. And then we propose a newly defined formula based on wavelet transform, and can determine the scaling exponent ξ(q) of temperature data. The obtained results demonstrate that we can correctly extract ξ(q) by using the method which is named as wavelet transform maximum modulus (WTMM).展开更多
Wavetet transform was used to analyze the scaling law of temperature data (passive scalar) in Rayleigh-Bénard convection flow from two aspects. The first one was to utilize the method of extended self similarity,...Wavetet transform was used to analyze the scaling law of temperature data (passive scalar) in Rayleigh-Bénard convection flow from two aspects. The first one was to utilize the method of extended self similarity, presented first by Benzi et al., to study the scaling exponent of temperature data. The obtained results show that the inertial range is much wider than that one determined directly from the conventional structure function, and find the obtained scaling exponent agrees well with the one obtained from the temperature data in an experiment of wind tunnel. The second one was that, by extending the formula which was proposed by A. Arneodo et al. for extracting the scaling exponent ζ(q) of velocity data to temperature data, a newly defined formula which is also based on wavelet transform, and can determine the scaling exponent ξ(q) of temperature data was proposed. The obtained results demonstrate that by using the method which is named as WTMM (wavelet transform maximum modulus) ξ(q) correctly can be extracted.展开更多
Nonlinear response is an important factor affecting the accuracy of three-dimensional image measurement based on the fringe structured light method.A phase compensation algorithm combined with a Hilbert transform is p...Nonlinear response is an important factor affecting the accuracy of three-dimensional image measurement based on the fringe structured light method.A phase compensation algorithm combined with a Hilbert transform is proposed to reduce the phase error caused by the nonlinear response of a digital projector in the three-dimensional measurement system of fringe structured light.According to the analysis of the influence of Gamma distortion on the phase calculation,the algorithm establishes the relationship model between phase error and harmonic coefficient,introduces phase shift to the signal,and keeps the signal amplitude constant while filtering out the DC component.The phase error is converted to the transform domain,and compared with the numeric value in the space domain.The algorithm is combined with a spiral phase function to optimize the Hilbert transform,so as to eliminate external noise,enhance the image quality,and get an accurate phase value.Experimental results show that the proposed method can effectively improve the accuracy and speed of phase measurement.By performing phase error compensation for free-form surface objects,the phase error is reduced by about 26%,and about 27%of the image reconstruction time is saved,which further demonstrates the feasibility and effectiveness of the method.展开更多
On the basis of Hartmann Shack sensor imaging analysis, a new method is presented with which the wavefront slope can be obtained when the object is incoherent and extended. This method, which is demonstrated by both ...On the basis of Hartmann Shack sensor imaging analysis, a new method is presented with which the wavefront slope can be obtained when the object is incoherent and extended. This method, which is demonstrated by both theoretical interpreting and computer simulation, explains how to measure the wavefront slope difference between two sub apertures through the determination of image displacements on detector plane. It includes a fast and accurate digital algorithm for detecting wavefront disturbance, which is much suitable for realization in such electrical hardwares as digital signal processors.展开更多
A new method of unscented extended Kalman filter (UEKF) for nonlinear system is presented. This new method is a combination of the unscented transformation and the extended Kalman filter (EKF). The extended Kalman...A new method of unscented extended Kalman filter (UEKF) for nonlinear system is presented. This new method is a combination of the unscented transformation and the extended Kalman filter (EKF). The extended Kalman filter is similar to that in a conventional EKF. However, in every running step of the EKF the unscented transformation is running, the deterministic sample is caught by unscented transformation, then posterior mean of non- lineadty is caught by propagating, but the posterior covariance of nonlinearity is caught by linearizing. The accuracy of new method is a little better than that of the unscented Kalman filter (UKF), however, the computational time of the UEKF is much less than that of the UKF.展开更多
The adjacency matrix operations,which connect with configuration transformation correspondingly,can be used for analysis of configuration transformation of metamorphic mechanisms and the corresponding algorithm can ea...The adjacency matrix operations,which connect with configuration transformation correspondingly,can be used for analysis of configuration transformation of metamorphic mechanisms and the corresponding algorithm can easily be simulated by computer.But the adjacency matrix based on monochrome topological graph is not suitable for the topological representation of mechanisms with multiple joints.The method of adjacency matrix operations has its own limitations for analysis of configuration transformation of metamorphic mechanisms because it can only be used in the topological representation of mechanisms with single joints.In order to overcome the drawback of the adjacency matrix,a kind of new matrix named as extended adjacency matrix is proposed to express topological structures of all mechanisms.The extended adjacency matrix is not only suitable for the topological representation of mechanisms with single joints,but also can be used in that of mechanisms with multiple joints.On this basis,a method of matrix operations based on the extended adjacency matrix is proposed to analyze the configuration transformation of metamorphic mechanisms.The method is not only suitable for configuration analysis of metamorphic mechanisms with single joints as well as metamorphic mechanisms with multiple joints.The method is evaluated by calculating two examples representing metamorphic mechanisms with single joint and multiple joints respectively.It can be concluded that the method is effective and correct for analysis of configuration transformation of all metamorphic mechanisms.The proposed method is simple and easy to be achieved by computer programming.It provides a basis for structural synthesis of all metamorphic mechanisms.展开更多
The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, ma...The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.展开更多
J-integral has served as a powerful tool in characterizing crack tip status. The main feature, i.e. path- independence, makes it one of the foremost fracture parameters. In order to remain the path- independence for f...J-integral has served as a powerful tool in characterizing crack tip status. The main feature, i.e. path- independence, makes it one of the foremost fracture parameters. In order to remain the path- independence for fluid-driven cracks, J-integral is revised. In this paper, we present an extended J-in- tegral explicitly for fluid-driven cracks, e.g. hydraulically induced fractures in petroleum reservoirs, for three-dimensional (3D) problems. Particularly, point-wise 3D extended J-integral is proposed to char- acterize the state of a point along crack front. Besides, applications of the extended J-integral to porous media and thermally induced stress conditions are explored. Numerical results show that the extended J- integral is indeed path-independent, and they are in good agreement with those of equivalent domain integral under linear elastic and elastoplastic conditions. In addition, two distance-independent circular integrals in the K-dominance zone are established, which can be used to calculate the stress intensity factor (SIF).展开更多
Abstract By improving the extended homogeneous balance method, a general method is suggested to derive a new auto-Bgcklund transformation (BT) for (3-k l)-Dimensional Jimbo-Miwa (JM) equation. The auto-BT obtain...Abstract By improving the extended homogeneous balance method, a general method is suggested to derive a new auto-Bgcklund transformation (BT) for (3-k l)-Dimensional Jimbo-Miwa (JM) equation. The auto-BT obtained by using our method only involves one quadratic homogeneity equation written as a bilinear equation. Based on the auto-BT, two-soliton solution of the (3+1)-Dimensional JM equation is obtained.展开更多
By the application of the extended homogeneous balance method, we derive anauto-Backlund transformation (BT) for (2+1)-dimensional variable coefficient generalized KPequations. Based on the BT, in which there are two ...By the application of the extended homogeneous balance method, we derive anauto-Backlund transformation (BT) for (2+1)-dimensional variable coefficient generalized KPequations. Based on the BT, in which there are two homogeneity equations to be solved, we obtainsome exact solutions containing single solitary waves.展开更多
文摘In this paper, an extended functional transformation is given to solve some nonlinear evolution equations. This function, in fact, is a solution of the famous KdV equation, so this transformation gives a transformation between KdV equation and other soliton equations. Then many new exact solutions can be given by virtue of the solutions of KdV equation.
基金Project(2017YFC0404802)supported by the National Key R&D Program of ChinaProjects(U1965206,51979143)supported by the National Natural Science Foundation of China。
文摘In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.
文摘In this paper, we extend the mapping transformation method through introducing variable coefficients.By means of the extended mapping transformation method, many explicit and exact general solutions with arbitrary functions for some nonlinear partial differential equations, which contain solitary wave solutions, trigonometric function solutions, and rational solutions, are obtained.
基金Supported by the National Natural Science Foundation of China under Grant No.60772023the Open Fund under Grant No.SKLSDE-2011KF-03+2 种基金Supported project under Grant No.SKLSDE-2010ZX-07 of the State Key Laboratory of Software Development Environment,Beijing University of Aeronautics and Astronauticsthe National High Technology Research and Development Program of China(863 Program) under Grant No.2009AA043303the Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.200800130006,Chinese Ministry of Education
文摘This paper is to investigate the extended(2+1)-dimensional Konopelchenko-Dubrovsky equations,which can be applied to describing certain phenomena in the stratified shear flow,the internal and shallow-water waves, plasmas and other fields.Painleve analysis is passed through via symbolic computation.Bilinear-form equations are constructed and soliton solutions are derived.Soliton solutions and interactions are illustrated.Bilinear-form Backlund transformation and a type of solutions are obtained.
文摘In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Calogero-Degasperis (VCCD) equation. New soliton and periodic solutions of many types are obtained. Furthermore, the soliton propagation is discussed under the effect of the variable coefficients.
基金supported in part by the National Natural Science Foundation of China under Grant 61901140,in part by the National Natural Science Foundation of China under Grant 62171151in part by the Science and Technology on Communication Networks Laboratory under Grant 6142104190203in part by the Fundamental Research Funds for the Central Universities under Grant HIT.OCEF.2021012。
文摘In this paper,we propose an extended hybrid carrier system based on the weighted fractional Fourier transform to ensure the reliability of wireless communication.The proposed scheme improves the dispersion and compensation capabilities of the hybrid carrier system for channel fading through the design of the signal power distribution,which has greatly reduced the probability of high-power distortion of the signal and improved the bit error rate performance as a result.Theoretical analysis has shown the superiority of the extended hybrid carrier system.With a lower cost of computational complexity increment,the proposed scheme obtains a performance improvement without occupying additional time-frequency physical resources.Compared with the existing hybrid carrier scheme,numerical simulation results have shown that the proposed extended hybrid carrier scheme has better anti-fading performance under the doubly-selective channel and improves the reliability of the wireless communication system effectively.
文摘Jordan's lemma can be used for a wider range than the original one. The extended Jordan's lemma can be described as follows. Let f(z) be analytic in the upper half of the z plane (Imz≥0), with the exception of a finite number of isolated singularities, and for P>o, if then where z=Rei and CR is the open semicircle in the upper half of the z plane.With the extended Jordan's lemma one can find that Laplace transform and Fourier transform are a pair of integral transforms which relate to each other.
文摘Wavelet transform is used to analyze the scaling rule of temperature data (passive scalar) in Rayleigh Bénard convection flow from two aspects. By utilizing the method of extended self similarity (ESS), one can find the obtained scaling exponent agrees well with the one obtained from the temperature data in a experiment of wind tunnel. And then we propose a newly defined formula based on wavelet transform, and can determine the scaling exponent ξ(q) of temperature data. The obtained results demonstrate that we can correctly extract ξ(q) by using the method which is named as wavelet transform maximum modulus (WTMM).
文摘Wavetet transform was used to analyze the scaling law of temperature data (passive scalar) in Rayleigh-Bénard convection flow from two aspects. The first one was to utilize the method of extended self similarity, presented first by Benzi et al., to study the scaling exponent of temperature data. The obtained results show that the inertial range is much wider than that one determined directly from the conventional structure function, and find the obtained scaling exponent agrees well with the one obtained from the temperature data in an experiment of wind tunnel. The second one was that, by extending the formula which was proposed by A. Arneodo et al. for extracting the scaling exponent ζ(q) of velocity data to temperature data, a newly defined formula which is also based on wavelet transform, and can determine the scaling exponent ξ(q) of temperature data was proposed. The obtained results demonstrate that by using the method which is named as WTMM (wavelet transform maximum modulus) ξ(q) correctly can be extracted.
基金This work is funded by the Scientific and Technological Projects of Henan Province under Grant 152102210115.
文摘Nonlinear response is an important factor affecting the accuracy of three-dimensional image measurement based on the fringe structured light method.A phase compensation algorithm combined with a Hilbert transform is proposed to reduce the phase error caused by the nonlinear response of a digital projector in the three-dimensional measurement system of fringe structured light.According to the analysis of the influence of Gamma distortion on the phase calculation,the algorithm establishes the relationship model between phase error and harmonic coefficient,introduces phase shift to the signal,and keeps the signal amplitude constant while filtering out the DC component.The phase error is converted to the transform domain,and compared with the numeric value in the space domain.The algorithm is combined with a spiral phase function to optimize the Hilbert transform,so as to eliminate external noise,enhance the image quality,and get an accurate phase value.Experimental results show that the proposed method can effectively improve the accuracy and speed of phase measurement.By performing phase error compensation for free-form surface objects,the phase error is reduced by about 26%,and about 27%of the image reconstruction time is saved,which further demonstrates the feasibility and effectiveness of the method.
文摘On the basis of Hartmann Shack sensor imaging analysis, a new method is presented with which the wavefront slope can be obtained when the object is incoherent and extended. This method, which is demonstrated by both theoretical interpreting and computer simulation, explains how to measure the wavefront slope difference between two sub apertures through the determination of image displacements on detector plane. It includes a fast and accurate digital algorithm for detecting wavefront disturbance, which is much suitable for realization in such electrical hardwares as digital signal processors.
文摘A new method of unscented extended Kalman filter (UEKF) for nonlinear system is presented. This new method is a combination of the unscented transformation and the extended Kalman filter (EKF). The extended Kalman filter is similar to that in a conventional EKF. However, in every running step of the EKF the unscented transformation is running, the deterministic sample is caught by unscented transformation, then posterior mean of non- lineadty is caught by propagating, but the posterior covariance of nonlinearity is caught by linearizing. The accuracy of new method is a little better than that of the unscented Kalman filter (UKF), however, the computational time of the UEKF is much less than that of the UKF.
基金supported by National Natural Science Foundation of China (Grant No. 51075039, Grant No. 50705010)Beijing Municipal Natural Science Foundation of China (Grant No. 3082014, Grant No.3053017)Fundamental Research Funds for the Central Universities of China (Grant No. 2009CZ08)
文摘The adjacency matrix operations,which connect with configuration transformation correspondingly,can be used for analysis of configuration transformation of metamorphic mechanisms and the corresponding algorithm can easily be simulated by computer.But the adjacency matrix based on monochrome topological graph is not suitable for the topological representation of mechanisms with multiple joints.The method of adjacency matrix operations has its own limitations for analysis of configuration transformation of metamorphic mechanisms because it can only be used in the topological representation of mechanisms with single joints.In order to overcome the drawback of the adjacency matrix,a kind of new matrix named as extended adjacency matrix is proposed to express topological structures of all mechanisms.The extended adjacency matrix is not only suitable for the topological representation of mechanisms with single joints,but also can be used in that of mechanisms with multiple joints.On this basis,a method of matrix operations based on the extended adjacency matrix is proposed to analyze the configuration transformation of metamorphic mechanisms.The method is not only suitable for configuration analysis of metamorphic mechanisms with single joints as well as metamorphic mechanisms with multiple joints.The method is evaluated by calculating two examples representing metamorphic mechanisms with single joint and multiple joints respectively.It can be concluded that the method is effective and correct for analysis of configuration transformation of all metamorphic mechanisms.The proposed method is simple and easy to be achieved by computer programming.It provides a basis for structural synthesis of all metamorphic mechanisms.
基金The Natural Science Foundation of Tianjin, China under contract No.08JCZDZT00200
文摘The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.
文摘J-integral has served as a powerful tool in characterizing crack tip status. The main feature, i.e. path- independence, makes it one of the foremost fracture parameters. In order to remain the path- independence for fluid-driven cracks, J-integral is revised. In this paper, we present an extended J-in- tegral explicitly for fluid-driven cracks, e.g. hydraulically induced fractures in petroleum reservoirs, for three-dimensional (3D) problems. Particularly, point-wise 3D extended J-integral is proposed to char- acterize the state of a point along crack front. Besides, applications of the extended J-integral to porous media and thermally induced stress conditions are explored. Numerical results show that the extended J- integral is indeed path-independent, and they are in good agreement with those of equivalent domain integral under linear elastic and elastoplastic conditions. In addition, two distance-independent circular integrals in the K-dominance zone are established, which can be used to calculate the stress intensity factor (SIF).
基金Supported by National Natural Science Foundation of China under Grant No.11071209 the Natural Science Foundation of the Higer Education Institutions of Jiangsu Province under Grant No.10KJB110011
文摘Abstract By improving the extended homogeneous balance method, a general method is suggested to derive a new auto-Bgcklund transformation (BT) for (3-k l)-Dimensional Jimbo-Miwa (JM) equation. The auto-BT obtained by using our method only involves one quadratic homogeneity equation written as a bilinear equation. Based on the auto-BT, two-soliton solution of the (3+1)-Dimensional JM equation is obtained.
文摘By the application of the extended homogeneous balance method, we derive anauto-Backlund transformation (BT) for (2+1)-dimensional variable coefficient generalized KPequations. Based on the BT, in which there are two homogeneity equations to be solved, we obtainsome exact solutions containing single solitary waves.
基金Supported by National Natural Science Foundation of China (61135001, 61075029, 61074179, 61074155) and the Postdoctoral Science Foundation of China (20110491692)