AAM  >> Vol. 7 No. 7 (July 2018)

    基于首次积分法求解两个非线性薛定谔方程的精确解
    The First Integral Method for Solving Exact Solutions of Two Nonlinear Schrodinger Equations

  • 全文下载: PDF(625KB) HTML   XML   PP.962-969   DOI: 10.12677/AAM.2018.77113  
  • 下载量: 189  浏览量: 1,807  

作者:  

张清梅,熊梅,陈龙伟:云南财经大学,统计与数学学院,云南 昆明

关键词:
首次积分方法广义非线性薛定谔方程高阶色散非线性薛定谔方程行波解The First Integral Method The Generalized Nonlinear Schrodinger Equation The High Order Dispersion Nonlinear Schrodinger Equation Travelling Wave Solutions

摘要:

首次积分方法是由冯第一个提出的对非线性偏微分方程进行可靠处理的积分方法,该方法是基于交换代数环的理论。本文将利用首次积分法对广义非线性薛定谔方程和高阶色散非线性薛定谔方程的精确行波解进行研究。即先通过引入恰当的行波变换,将非线性薛定谔方程化为常微分方程,再根据多项式除法原理,得到两个非线性薛定谔方程的精确行波解。

The first integral method proposed by Feng is very reliable integral method for solving nonlinear partial differential equations, which is based on the ring theory of commutative algebra. In this paper, exact travelling wave solutions of the generalized nonlinear Schrodinger equation and the high order dispersion nonlinear Schrodinger equation are studied by using the first integral me-thod. By introducing the travelling wave transformations, two nonlinear Schrodinger equations have been transformed into ordinary differential equations. Then according to the division theorem of polynomial, exact travelling wave solutions of two nonlinear Schrodinger equations are obtained.

1. 引言

非线性现象出现在各种各样的科学应用中,如等离子体物理、固态物理、流体动力学等。为了更好地理解这些非线性现象,许多数学家和物理学家都在努力寻找其精确解,运用这些精确解可以预测和解释一些重要的物理现象。所以非线性偏微分方程(PDE)在物理现象研究中起着重要的作用。近年来,许多方法被用来寻找非线性PDE的解,其中一种方法被称为首次积分法,首次积分法是基于交换代数环的理论,由Feng [1] 第一个提出的对非线性PDE进行可靠处理,并由他自己进一步发展了这个理论 [2] [3] [4] 。这一方法已被许多学者应用于解决在科学和工程中遇到的不同的偏微分方程模型 [5] [6] [7] 。为了得到非线性PDE更丰富的精确解,前人提出了如下这些方法,如tanh-sech函数法 [8] ,雅可比椭圆函数展开法 [9] ,F展开法 [10] ,指数函数法 [11] ,Hirota双线性法 [12] ,动力系统分岔理论 [13] 等等。

在本文中,我们将考虑如下广义非线性薛定谔方程 [14] 和高阶色散非线性薛定谔方程 [15] :

i u t r 2 u x x + c 3 | u | 2 u = i [ ( s 0 + s 2 | u | 2 ) u ] x c 5 | u | 4 u , (1)

i q t a q x x b q x x x x + c ( | q | 2 + d | q | 4 ) q = 0. (2)

这里 r 2 , c 3 , c 5 , s 0 , s 2 , a , b , c , d 为实常数, u , q 为复函数。

2. 首次积分法的基本介绍

考虑一个如下典型的非线性PDE

P ( u , u t , u x , u x x , u t t , u x t , u x x x , ) = 0. (3)

为了将方程(3)化为常微分方程(ODE),我们介绍这个变换

u ( x , t ) = u ( ξ ) , ξ = x c t . (4)

其中c是常数。

t ( . ) = c ξ ( . ) , x ( . ) = ξ ( . ) , 2 t 2 ( . ) = c 2 2 ξ 2 ( . ) , (5)

步骤1:通过上述变换可将PDE转换为如下ODE形式

Q ( u , u ξ , u ξ ξ , ) = 0. (6)

步骤2:假设这里的ODE(6)能被写成

u ( x , t ) = u ( ξ ) . (7)

步骤3:现在我们引入一个新的独立变量,即

X ( ξ ) = u ( ξ ) , Y ( ξ ) = u ξ ( ξ ) . (8)

并给出了ODE的系统

{ X ξ ( ξ ) = Y ( ξ ) , Y ξ ( ξ ) = F ( X ( ξ ) , Y ( ξ ) ) . (9)

步骤4:根据常微分方程的定性理论 [16] 。如果我们能在相同条件下找到方程(9)的积分,那么方程(9)的解就可以直接得到。但是,目前还没有一个系统的理论告诉我们如何去找到首次积分,因此我们将运用除法原理去寻找方程(9)的首次积分,它将会减少常微分方程(6)的一阶积分。通过解这个方程可以得到方程(3)的精确解。现在,我们来介绍除法原理。

定理(除法原理):设 P ( x , y ) Q ( x , y ) 为两个变量x和y在域 C ( x , y ) 中的多项式,并令 P ( x , y ) C ( x , y ) 中不可约。如果 Q ( x , y ) 在所有零点 G ( x , y ) 上可约,那么在 C ( x , y ) 中存在一个多项式 G ( x , y ) ,使得 Q ( x , y ) = P ( x , y ) G ( x , y )

3. 应用

在这一节中,我们将利用首次积分法来详细求解两个非线性薛定谔方程的精确解。

首先我们考虑方程(1)。假设方程(1)的解满足如下形式

u ( x , t ) = φ ( ξ ) e i ( x ω t ) . (10)

其中 ξ = x λ t ,这里 λ , ω 是实参数。

将方程(10)代入(1)中,并令虚部和实部分别为零,整理可得

( 2 r 2 λ s 0 ) φ 3 s 2 φ 2 φ = 0 , (11)

r 2 φ + ( ω + s 0 r 2 ) φ ( c 3 + s 2 ) φ 3 + c 5 φ 5 = 0. (12)

将方程(11)积分一次代入(12)得

φ + ω + r 2 λ r 2 φ + c 3 r 2 φ 3 + c 5 r 2 φ 5 = 0. (13)

为了计算简便,我们定义: k 1 = ω + r 2 λ r 2 k 2 = c 3 r 2 k 3 = c 5 r 2

根据方程(13)和(8)我们得

{ X ( ξ ) = Y ( ξ ) , Y ( ξ ) = k 1 X ( ξ ) + k 2 X ( ξ ) 3 + k 3 X ( ξ ) 5 . (14)

其次,我们考虑方程(2)。假设方程(2)的解满足如下形式

q ( x , t ) = ϕ ( ξ ) e i ( k x ω t ) . (15)

其中 ξ = β x λ t ,这里 β , λ , k , ω 是实参数。

将方程(15)代入(2)中,并令虚部和实部为零,整理可得

4 b k β 3 ϕ + ( λ + 2 a β k + 4 b β k 3 ) ϕ = 0. (16)

b β 4 ϕ ( 4 ) + ( a β 2 + b β 2 k 2 ) ϕ + ( ω a k 2 b k 4 ) ϕ + c ϕ 3 + d c ϕ 5 = 0. (17)

对(16)微分一次代入(17)中,令 p 1 = λ + 2 a β k + 4 b β k 3 , p 2 = a β 2 + b β 2 k 2 , p 3 = ω a k 2 b k 4 。则

ϕ 4 b k p 3 β 3 p 1 b β 4 4 p 2 b k β 3 ϕ 4 b c k β 3 p 1 b β 4 4 p 2 b k β 3 ϕ 3 4 b c d k β 3 p 1 b β 4 4 p 2 b k β 3 ϕ 5 = 0. (18)

根据方程(18)和(8)我们得

{ X ( ξ ) = Y ( ξ ) , Y ( ξ ) = 4 b k p 3 β 3 p 1 b β 4 4 p 2 b k β 3 X ( ξ ) + 4 b c k β 3 p 1 b β 4 4 p 2 b k β 3 X ( ξ ) 3 + 4 b c d k β 3 p 1 b β 4 4 p 2 b k β 3 X ( ξ ) 5 . (19)

这里如果我们仍然定义

k 1 = 4 b k p 3 β 3 p 1 b β 4 4 p 2 b k β 3 , k 2 = 4 b c k β 3 p 1 b β 4 4 p 2 b k β 3 , k 3 = 4 b c d k β 3 p 1 b β 4 4 p 2 b k β 3 方程(19)和(14)有相同的首次积分,要找到它们的精确解,我们仅需讨论其中一个方程即可,下面我们将给出对方程(14)的求解过程。

现在,根据除法原理,我们假设 X = X ( ξ ) Y = Y ( ξ ) 是方程(14)的非平凡解。且

P ( X , Y ) = i = 0 m a i ( X ) Y i . (20)

在复数域 C ( X , Y ) 中是一个不可约多项式,即

P ( X ( ξ ) , Y ( ξ ) ) = i = 0 m a i ( X ( ξ ) ) Y ( ξ ) i = 0. (21)

这里 a i ( X ) ( i = 1 , 2 , , m ) 是X的多项式以及 a m ( X ) 0 。方程(21)称为(14)的首次积分。定义 P ( X ( ξ ) , Y ( ξ ) ) 关于X和Y的多项式,且 d P d ξ | ( 20 ) = 0 。根据除法定理,我们可以看到在复数域 C ( X , Y ) 中存在一个多项式 H ( X , Y ) = h ( X ) + g ( X ) Y ,使得

d P d ξ | ( 20 ) = ( d P d X d X d ξ + d P d Y d Y d ξ ) | ( 20 ) = ( h ( X ) + g ( X ) Y ) ( i = 0 m a i ( X ) Y i ) . (22)

3.1. 情形一

在方程(21)中我们假设 m = 1 ,即方程(22)可以被写成

i = 0 1 a i ( X ) Y i + 1 + i = 0 1 i a i ( X ) Y i 1 ( Y ( ξ ) ) = ( h ( X ) + g ( X ) Y ) ( i = 0 1 a i ( X ) Y i ) . (23)

然后,在方程(23)中通过比较两边 Y i ( i = 2 , 1 , 0 ) 的系数,我们得到

a 1 ( X ) = g ( X ) a 1 ( X ) , (24)

a 0 ( X ) = h ( X ) a 1 ( X ) + g ( X ) a 0 ( X ) , (25)

a 1 ( X ) ( k 1 X + k 2 X 3 + k 5 X 5 ) = h ( X ) a 0 ( X ) . (26)

因为 a i ( x ) ( i = 0 , 1 ) 是X的一个多项式,从方程(24)中我们得出 a 1 ( X ) 是一个常数,且 g ( X ) = 0 。为了便于计算,我们令 a 1 ( X ) = 1 。平衡 h ( X ) a 0 ( X ) 的阶数,我们得出 deg ( h ( x ) ) = 2 。假设 h ( X ) = A X 2 + B X + C ,这里的 A 0 ,那么

a 0 ( X ) = 1 3 A X 3 + 1 2 B X 2 + C X + D . (27)

其中D是任意的积分常数。将 a 0 ( X ) , a 1 ( X ) , h ( X ) 代入方程(26)中,通过比较系数得到如下代数方程组

{ 1 3 A 2 = k 3 5 6 A B = 0 4 3 A C + 1 2 B 2 = k 2 A D + 3 2 B C = 0 B D + C 2 = k 1 C D = 0 (28)

用Maple求解它们,可得到

D = 0 , B = 0 , A = ± 3 k 3 , C = ± k 1 , (29)

这里 k 1 , k 3 < 0 且满足约束条件

k 1 k 3 = 3 16 k 2 2 . (30)

把条件(29)代入(21)中,得

Y = 1 3 A X 3 C X . (31)

再由方程(8)有

φ ( ξ ) = 1 3 A φ ( ξ ) 3 C φ ( ξ ) . (32)

为了求解(32)的解,我们引入伯努利方程 [17] 有

Z = a Z + b Z p , (33)

这里 a , b , p R , a b 0 , p 1 。则方程(33)有如下形式的解

Z ( ξ ) = [ a / b ξ 0 e a ( 1 p ) + 1 ] 1 p 1 , (34)

= { { a 2 b [ 1 + tanh ( a ( p 1 ) 2 ξ ln ξ 0 2 ) ] } 1 p 1 if ξ 0 > 0 , { a 2 b [ 1 + coth ( a ( p 1 ) 2 ξ ln ( ξ 0 ) 2 ) ] } 1 p 1 if ξ 0 < 0 , { a b } 1 p 1 if ξ 0 = 0. (35)

这里 ξ 0 为任意常数。令 a = C , b = 1 3 A , p = 3 代入方程(35)中有

u ( ξ ) = { { ± 1 2 3 k 1 k 3 [ 1 + tanh ( k 1 ξ ln ξ 0 2 ) ] } 1 2 if ξ 0 > 0 , { ± 1 2 3 k 1 k 3 [ 1 + coth ( k 1 ξ ln ( ξ 0 ) 2 ) ] } 1 2 if ξ 0 < 0 , { ± 3 k 1 k 3 } 1 2 if ξ 0 = 0. (36)

则广义非线性薛定谔方程(1)的行波解可写为

u ( x , t ) = { { ± 1 2 3 k 1 k 3 [ 1 + tanh ( k 1 ( x λ t ) ln ξ 0 2 ) ] } 1 2 e ( x ω t ) if ξ 0 > 0 , { ± 1 2 3 k 1 k 3 [ 1 + coth ( k 1 ( x λ t ) ln ( ξ 0 ) 2 ) ] } 1 2 e ( x ω t ) if ξ 0 < 0 , { ± 3 k 1 k 3 } 1 2 e ( x ω t ) if ξ 0 = 0. (37)

3.2. 情形二

在方程(21)中我们假设 m = 2 ,通过比较两边 Y i ( i = 2 , 1 , 0 ) 系数,我们得到

a 2 ( X ) = g ( X ) a 2 ( X ) , (38)

a 1 ( X ) = h ( X ) a 2 ( X ) + g ( X ) a 1 ( X ) , (39)

a 0 ( X ) + 2 a 2 ( X ) ( k 1 X + k 2 X 3 + k 5 X 5 ) = h ( X ) a 1 ( X ) + g ( X ) a 0 ( X ) , (40)

a 1 ( X ) ( k 1 X + k 2 X 3 + k 5 X 5 ) = h ( X ) a 0 ( X ) . (41)

因为 a i ( x ) ( i = 0 , 1 , 2 ) 是X的一个多项式,从方程(38)中我们得出 a 2 ( X ) 是一个常数,且 g ( X ) = 0 。为了便于计算,我们令 a 2 ( X ) = 1 。平衡 h ( X ) , a 1 ( X ) a 0 ( X ) 的阶数,我们得出 deg ( h ( x ) ) = 2 deg ( a 1 ( x ) ) = 3 deg ( a 0 ( x ) ) = 6 。假设 h ( X ) = A X 2 + B X + C ,这里的 A 0 ,那么

a 1 ( X ) = 1 3 A X 3 + 1 2 B X 2 + C X + D . (42)

a 0 ( X ) = ( 1 18 A 2 1 3 k 3 ) X 6 + ( 1 6 A B ) X 5 + ( 1 3 A C + 1 8 B 2 + 1 2 k 2 ) X 4 + ( 1 3 A D + 1 2 B C ) X 3 + ( 1 2 B D + 1 2 C 2 k 1 ) X 2 + D C X + E . (43)

其中D是任意的积分常数。将 a 0 ( X ) , a 1 ( X ) , a 2 ( X ) , h ( X ) 代入方程(41)中,通过比较系数得到如下代数方程组

{ 1 18 A 3 = 2 3 A k 3 2 9 A 2 B = 5 6 B k 3 7 18 A 2 C + 7 24 A B 2 = 5 6 A k 2 + 4 3 C k 3 1 3 A 2 D + A B C + 1 8 B 2 = B k 2 + D k 3 7 6 A B D + 5 6 A C 2 + 1 8 B 2 C = 4 3 A k 1 + 3 2 C k 2 4 3 A C D + 1 2 D B 2 + B C 2 = 3 2 B k 1 + 2 D k 2 3 2 B C D + 1 2 C 2 + A E = 2 C k 1 D C 2 + B E = D k 1 E C = 0 , (44)

用Maple求解它们,可得到

D = 0 , B = 0 , A = ± 2 3 k 3 , C = ± 2 k 1 , (45)

这里 k 1 , k 3 < 0 且满足约束条件

k 1 k 3 = 3 16 k 2 2 . (46)

把条件(37)代入(21)中,得

Y = 1 6 A X 3 1 2 C X . (47)

再由方程(8)有

φ ( ξ ) = 1 6 A φ ( ξ ) 3 1 2 C φ ( ξ ) . (48)

同理,令 a = 1 2 C , b = 1 6 A , p = 3 代入方程(34)和(35)中,则方程(48)和(32)有相同的精确解。

4. 总结

在本文中,首次积分法成功地求解了广义非线性薛定谔方程和高阶色散非线性薛定谔方程,丰富了其解空间。与其他方法不同的是,它是求解精确行波解的一种有效的方法,即我们可以轻松地借助计算机软件完成其复杂的代数计算,从而得到了更精确的行波解,因此,该方法可以推广到求解更多非线性偏微分方程的精确解。

文章引用:
张清梅, 熊梅, 陈龙伟. 基于首次积分法求解两个非线性薛定谔方程的精确解[J]. 应用数学进展, 2018, 7(7): 962-969. https://doi.org/10.12677/AAM.2018.77113

参考文献

[1] Feng, Z.S. (2002) The First Integral Method to Study the Burgers Korteweg-de Vries Equation. Physics Letters A, 35, 343-349.
https://doi.org/10.1088/0305-4470/35/2/312
[2] Feng, Z.S. (2002) On Explicit Exact Solutions to the Compound Burgers Korteweg-de Vries Equation. Physics Letters A, 293, 57-66.
https://doi.org/10.1016/S0375-9601(01)00825-8
[3] Feng, Z.S. (2008) Travelling Wave Behavior for a generalized Fisher Equation. Chaos, Solitons & Fractals, 38, 481-488.
https://doi.org/10.1016/j.chaos.2006.11.031
[4] Feng, Z.S. and Knobel, R. (2007) Travelling Waves to a Burgers Korteweg-de Vries Equation with Higher Order Nonlinearities. Journal of Mathematical Analysis and Applications, 328, 1435-1450.
https://doi.org/10.1016/j.jmaa.2006.05.085
[5] Abdoon, M.A. (2015) Programming First Integral Method General Formula for the Solving Linear and Nonlinear Equations. Applied Mathematics, 6, 568-575.
https://doi.org/10.4236/am.2015.63051
[6] Seadawy, A. and Sayed, A. (2017) Soliton Solutions of Cubic-Quintic Nonlinear Schrodinger and Variant Boussinesq Equations by the First Integral Method. Filomat, 214, 4199-4208.
https://doi.org/10.2298/FIL1713199S
[7] Ibrahim, S. and El-Ganaini, A. (2007) The First Integral Method to the Nonlinear Schrodinger Equations in Higher Dimensions. Abstract and Applied Analysis, 18, 1187-1197.
[8] Ma, W.X. (1993) Travelling Wave Solutions to a Seventh Order Generalized KdV Equation. Physics Letters A, 180, 221-224.
https://doi.org/10.1016/0375-9601(93)90699-Z
[9] Liu, S., Fu, Z., Liu, S.D. and Zhao, Q. (2001) Jacobi Elliptic Function Expansion Method and Periodic Wave Solutions of Nonlinear Wave Equations. Physics Letters A, 289, 69-74.
https://doi.org/10.1016/S0375-9601(01)00580-1
[10] Zhang, S. (2006) The Periodic Wave Solutions for the (2 + 1) Dimensional Konopelchenko Dubrovsky Equations. Chaos, Solitons & Fractals, 30, 1213-1220.
https://doi.org/10.1016/j.chaos.2005.08.201
[11] He, J.H. and Zhang, L.N. (2008) Generalized Solitary Solution and Compac-ton-Like Solution of the Jaulent-Miodek Equations Using the Exp-Function Method. Physics Letters A, 372, 1044-1047.
https://doi.org/10.1016/j.physleta.2007.08.059
[12] Wang, D.S. (2009) A systematic Method to Construct Hirotas Transforma-tions of Continuous Soliton Equations and Its Applications. Computers and Mathematics with Applications, 58, 146-153.
https://doi.org/10.1016/j.camwa.2009.03.077
[13] Li, Y., Shan, W.R., Shuai, T.P. and Rao, K. (2015) Bifurcation Analysis and Solutions of a Higher Order Nonlinear Schrodinger Equation. Mathematical Problems in Engineering, 3, 1-10.
https://doi.org/10.1155/2015/408586
[14] Geng, Y.X. and Li, J.B. (2010) Exact Explicit Traveling Wave Solutions for Two Nonlinear Schrodinger Type Equations. Applied Mathematics and Computation, 217, 1509-1521.
https://doi.org/10.1016/j.amc.2009.06.031
[15] Xu, G.Q. (2011) New Types of Exact Solutions for the Fourth-Order Disperesive Cubic-Quintic Nonlinear Schrodinger Equation. Applied Mathematics and Computation, 217, 5967-5971.
https://doi.org/10.1016/j.amc.2010.12.008
[16] Ding, T.R. and Li, C.Z. (1996) Ordinary Differential Equations. Peking Univer-sity Press, Peking.
[17] Ma, W.X. and Fuchsstciner, B. (1996) A Explicit and Exact Solutions to Kolmogorov-Petrovskli-Piskunov Equation. International Journal of Non-Linear Mechanics, 31, 329-338.
https://doi.org/10.1016/0020-7462(95)00064-X