Advances in Applied Mathematics
Vol. 08  No. 11 ( 2019 ), Article ID: 32822 , 11 pages
10.12677/AAM.2019.811198

Travelling Wave Solutions of a Class of Generalized Klein-Gordon Equations

Yuee Lv1, Qi Liu2

1College of Information Science and Engineering, Linyi University, Linyi Shandong

2College of Mathematics and Statistics, Linyi University, Linyi Shandong

Received: Oct. 10th, 2019; accepted: Oct. 29th, 2019; published: Nov. 5th, 2019

ABSTRACT

In this paper, the travelling wave solutions of a class of generalized Klein-Gordon equations are studied by Riccati equation method. Some new travelling wave solutions of special forms are obtained.

Keywords:Riccati Equation, Traveling Wave Solution

广义Klein-Gordon方程的行波解

吕月娥1,刘琪2

1临沂大学信息科学与工程学院,山东 临沂

2临沂大学数学与统计学院,山东 临沂

收稿日期:2019年10月10日;录用日期:2019年10月29日;发布日期:2019年11月5日

摘 要

本文运用黎卡提方程法解决Klein-Gordon方程的行波解问题,得到了几类特殊形式的行波解。

关键词 :黎卡提方程,行波解

Copyright © 2019 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 基础知识

形如

d y d t = A ( t ) y 2 + B ( t ) y + C (t)

的方程称为黎卡提方程,其中 A ( t ) , B ( t ) C ( t ) 在所讨论的区间内为连续函数。1841年法国数学家刘维尔证明了在一般情形下黎卡提方程的解不能运用函数的有限次积分以及有限次代数运算而得到,然而他得到了特殊类型的黎卡提方程可化为伯努利方程求解的方法。国内外在黎卡提方程的推广研究侧重点各有差异,比如代数微分方程,离散微分方程,KdV方程的无穷序列复合型类孤子解等。

1) 运用黎卡提方程,可以求解非线性波动方程

p ( u , u t , u x , u t t , u t x , u x x , ) = 0

的行波解,该方法主要步骤如下:

我们首先引入波动变量 ξ = x c t 使得 u ( x , t ) = u ( ξ ) ,因此偏微分方程化为关于 u = u ( ξ ) 的常微分方程

p ( u , c u , u , c 2 u , c u , u , ) = 0 (1)

2) 然后我们假设方程(1)的解可以表示为

u ( ξ ) = k = 0 N a k ( ξ ) y ( μ ξ ) k (2)

a k ( ξ ) ( k = 0 , 1 , , N )是关于 ξ 的函数,是需要进一步确定的常数, y ( t ) 是黎卡提方程的解。之后有

{ d d t y k ( t ) = k [ A ( t ) y k + 1 + B ( t ) y + C ( t ) y k 1 ] d 2 d t 2 y k ( t ) = k ( k + 1 ) A 2 ( t ) y k + 2 + k [ ( 2 k + 1 ) A ( t ) B ( t ) + A ( t ) ] y k + 1 + k [ 2 k A ( t ) C ( t ) + k B 2 ( t ) + B ( t ) ] y k + k [ ( 2 k 1 ) B ( t ) C ( t ) + C ( t ) ] y k 1 + k ( k 1 ) C 2 ( t ) y k 2 (3)

其他导数也是如此。因此, d d ξ u ( ξ ) d 2 d ξ 2 u ( ξ ) 分别是 N + 1 N + 2 次黎卡提微分方程的多项式。

3) 将(2)代入到方程(1)中,通过(3)得到关于 y ( u ξ ) 的代数方程。为了确定参数N,一般我们会平衡得到的方程中最高阶的线性项与最高阶的非线性项,并收集所有相同 y ( u ξ ) 的系数,把这些系数都消去,最后得到一个包括函数 a k ( ξ ) ( k = 0 , 1 , , N ),c和 μ 的方程组。

4) 求解步骤3的代数方程组得到了函数 a k ( ξ ) ( k = 0 , 1 , , N ),c和 μ 的所有可能值,确定了这些函数和参数,知道N是正整数,通过(2)得到封闭形式的解析解 u ( x , t )

注1:如果我们用黎卡提方程 y ( t ) = y 2 ( t ) + 1 的解 y ( t ) = tanh ( t ) ,当 a k ( ξ ) ( k = 0 , 1 , , N ),则黎卡提方程法可转化为tanh思想 [1]。

注2:若我们用黎卡提方程 y ( t ) = ( y 2 ( t ) + α y ( t ) + β ) 的解 ,当 ω ( t ) 是二阶线性微分方程 ω ( t ) + α ω ( t ) + β ω ( t ) = 0 时,取 μ = 1 ,则黎卡提方程法转化为 ( G G ) -展开法 [2] [3]。

2. 广义Klein-Gordon方程行波解

广义Klein-Gordon方程为

u t t a 2 u x x + b u + f u 2 k u 3 = 0 (4)

a , b , f 和k都为常数, u ( x , t ) 是未知函数。引入波动变量 ξ = x c t ,则方程(4)转化为

b u + f u 2 k u 3 + ( c 2 a 2 ) u ξ ξ = 0 (5)

则平衡 u ξ ξ u 3 后有 N + 2 = 3 N ,解得 N = 1 。因此,(2)可表示为

u ( ξ ) = a 1 ( ξ ) y ( μ ξ ) + a 0 ( ξ ) (6)

将(6)代入到方程(5)中,并计算 y ( μ ξ ) 的系数,就可得出关于 a 1 ( ξ ) a 0 ( ξ ) ,c和 μ 的方程组

{ y 3 : k a 1 3 ( ξ ) + 2 μ 2 ( c 2 a 2 ) a 1 ( ξ ) A 2 ( μ ξ ) = 0 y 2 : f a 1 2 ( ξ ) 3 k a 0 ( ξ ) a 1 2 ( ξ ) + ( c 2 a 2 ) [ 2 μ a 1 ( ξ ) A ( μ ξ ) + μ 2 a 1 ( ξ ) ( A ( μ ξ ) + 3 A ( μ ξ ) B ( μ ξ ) ) ] = 0 y 1 : b a 1 ( ξ ) + 2 f a 0 ( ξ ) a 1 ( ξ ) 3 k a 0 2 ( ξ ) a 1 ( ξ ) + ( c 2 a 2 ) [ a 1 ( ξ ) + 2 μ a 1 ( ξ ) B ( μ ξ ) + μ 2 a 1 ( ξ ) ( B ( μ ξ ) + 2 A ( μ ξ ) C ( μ ξ ) + B 2 ( μ ξ ) ) ] = 0 y 0 : b a 0 ( ξ ) + f a 0 2 ( ξ ) 3 k a 0 3 ( ξ ) + ( c 2 a 2 ) [ a 0 ( ξ ) + 2 μ a 1 ( ξ ) C ( μ ξ ) + μ 2 a 1 ( ξ ) ( C ( μ ξ ) + B ( μ ξ ) C ( μ ξ ) ) ] = 0 (7)

上述代数系统依赖于黎卡提方程中系数 A ( t ) , B ( t ) C ( t ) 的选择,下面我们选取系数 A ( t ) , B ( t ) C ( t ) 的一些特殊情况,基于这些选择和上述方程组我们构造了推广的Klein-Gordon方程(4)的特殊精确解 [4]。

1) 如果我们取 A ( t ) = 1 , B ( t ) = 0 , C ( t ) = 1 ,则可得到黎卡提方程 d d t y ( t ) = y 2 ( t ) + 1 ,解为

y ( t ) = tanh ( t ) 。这时解方程组(7)得

{ a 0 = f 3 k a 1 = ± 3 k b + f 2 3 k 2 μ 2 = 3 k b + f 2 6 k ( c 2 a 2 ) , ( c a ) (8)

{ a 0 = 2 f 3 k a 1 = ± b k μ 2 = b 2 ( c 2 a 2 ) , ( c a ) (9)

其中c为任意常数。根据(9)可把解(6)写为

u ( x , t ) = ± 3 k b + f 2 3 k 2 y ( μ ξ ) + f 3 k (10)

u ( x , t ) = ± b k y ( μ ξ ) + 2 f 3 k (11)

3 k b + f 2 6 k ( c 2 a 2 ) > 0 时,有纽结解

u 1 ( x , t ) = ± 3 k b + f 2 3 k 2 tanh ( 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) ) + f 3 k (12)

k , b , f , c 取1,a取0时,解(12)的图像为图1

Figure 1. Image of solution (12)

图1. 解(12)的图像

3 k b + f 2 6 k ( c 2 a 2 ) < 0 时,有复周期解

u 2 ( x , t ) = ± i 3 k b + f 2 3 k 2 tan ( 3 k b + f 2 6 k ( a 2 c 2 ) ( x c t ) ) + f 3 k (13)

b 2 ( c 2 a 2 ) > 0 时,有纽结解

u 3 ( x , t ) = ± b k tanh ( b 2 ( c 2 a 2 ) ( x c t ) ) + 2 f 3 k (14)

k , b , f , c 取1,a取0时,解(14)的图像为图2

Figure 2. Image of solution (14)

图2. 解(14)的图像

b 2 ( c 2 a 2 ) < 0 时,有复周期解

u 4 ( x , t ) = ± i b k tan ( b 2 ( a 2 c 2 ) ( x c t ) ) + 2 f 3 k (15)

如果我们用黎卡提方程 d d t y ( t ) = y 2 ( t ) + 1 的解 y ( t ) = coth ( t ) ,则

3 k b + f 2 6 k ( c 2 a 2 ) > 0 时,有纽结解

u 5 ( x , t ) = ± 3 k b + f 2 3 k 2 coth ( 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) ) + f 3 k (16)

k , b , f , c 取1,a取0时,解(16)的图像为图3

3 k b + f 2 6 k ( c 2 a 2 ) < 0 时,有复周期解

u 6 ( x , t ) = ± i 3 k b + f 2 3 k 2 cot ( 3 k b + f 2 6 k ( a 2 c 2 ) ( x c t ) ) + f 3 k (17)

Figure 3. Image of solution (16)

图3. 解(16)的图像

b 2 ( c 2 a 2 ) > 0 时,有纽结解

u 7 ( x , t ) = ± b k coth ( b 2 ( c 2 a 2 ) ( x c t ) ) + 2 f 3 k (18)

b 2 ( c 2 a 2 ) < 0 时,有复周期解

u 4 ( x , t ) = ± i b k cot ( b 2 ( a 2 c 2 ) ( x c t ) ) + 2 f 3 k (19)

当±都取+时,行波解与文 [5] 中采用tanh法得到的行波解相同。

2) 如果我们取 A ( t ) = 1 , B ( t ) = α , C ( t ) = β ,则可得到黎卡提方程 d d t y ( t ) = ( y 2 ( t ) + α y ( t ) + β ) ,解为 y ( t ) = ω ( t ) ω ( t ) ,而 ω ( t ) 是二阶线性微分方程 ω ( t ) + α ω ( t ) + β ω ( t ) = 0 的解(13)。这时解方程组(7)得

{ a 0 = f 3 k + α 2 a 1 a 1 = ± 4 ( 3 k b + f 2 ) 3 k 2 ( α 2 4 β ) μ 2 = 2 ( 3 k b + f 2 ) 3 k ( c 2 a 2 ) ( α 2 4 β ) (20)

其中 c , α β ( c a , α 2 4 β , k 0 )为任意常数。这将得到孤波解

如果,则我们可得

ω ( t ) = exp ( α 2 t ) ( C 1 cosh ( α 2 4 β 2 t ) + C 2 sinh ( α 2 4 β 2 t ) )

ξ = x c t 时有

u 1 ( x , t ) = ± 3 k b + f 2 3 k 2 [ C 1 sinh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) + C 2 cosh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) C 1 cosh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) + C 2 sinh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) ] + f 3 k (21)

C 1 取1, C 2 取2, k , b , f , c 取1,a取0时,解(21)的图像为图4

Figure 4. Image of solution (21)

图4. 解(21)的图像

u 2 ( x , t ) = ± 3 k b + f 2 3 k 2 [ C 1 sinh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) + C 2 cosh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) C 1 cosh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) C 2 sinh ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) ] + f 3 k (22)

C 1 取1, C 2 取2, k , b , f , c 取1,a取0时,解(22)的图像为图5

如果 α 2 4 β < 0 ,则我们可得

ω ( t ) = e α 2 t ( C 1 cos ( 4 β α 2 2 t ) + C 2 sin ( 4 β α 2 2 t ) )

ξ = x c t 时有

Figure 5. Image of solution (22)

图5. 解(22)的图像

u 3 ( x , t ) = ± 3 k b f 2 3 k 2 [ C 1 sin ( 3 k b + f 2 6 k ( a 2 c 2 ) ξ ) + C 2 cos ( 3 k b + f 2 6 k ( a 2 c 2 ) ξ ) C 1 cos ( 3 k b + f 2 6 k ( a 2 c 2 ) ξ ) + C 2 sin ( 3 k b + f 2 6 k ( a 2 c 2 ) ξ ) ] + f 3 k (23)

u 4 ( x , t ) = ± 3 k b f 2 3 k 2 [ C 1 sin ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) + C 2 cos ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) C 1 cos ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) C 2 sin ( 3 k b + f 2 6 k ( c 2 a 2 ) ξ ) ] + f 3 k (24)

其中 C 1 , C 2 为任意常数。如果取 C 1 = 1 , C 2 = 0 ,则解(3-22) (3-23)可写为纽结解

u 5 ( x , t ) = ± 3 k b + f 2 3 k 2 tanh ( 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) ) + f 3 k (25)

且解(23) (24)可写为周期解

u 6 ( x , t ) = ± 3 k b f 2 3 k 2 tan ( 3 k b + f 2 6 k ( a 2 c 2 ) ( x c t ) ) + f 3 k (26)

3) 如果我们取 A ( t ) = α , B ( t ) = 0 , C ( t ) = β ,则可得到黎卡提方程 d d t y ( t ) = α y 2 ( t ) + β

解为 y ( t ) = 1 α Q ( t ) Q ( t ) ,而 J 1 2 , Y 1 2 是取 1 2 时的贝塞尔函数 [6],其中 C1,C2 是任意常数。这时解方程组(7)得

{ a 0 = f 3 k a 1 = ± α ( 3 k b + f 2 ) 3 k 2 β μ 2 = 3 k b f 2 6 k α β ( a 2 c 2 ) (27)

(28)

其中 c , α ( c a )为任意常数。当解为(27)时,得到行波解

u 1 ( x , t ) = ± 1 α α ( 3 k b + f 2 ) 3 k 2 β Q ( 3 k b f 2 6 k α β ( a 2 c 2 ) ( x c t ) ) Q ( 3 k b f 2 6 k α β ( a 2 c 2 ) ( x c t ) ) + f 3 k (29)

u 2 ( x , t ) = ± 1 α α ( 3 k b + f 2 ) 3 k 2 β Q ( 3 k b f 2 6 k α β ( a 2 c 2 ) ( x c t ) ) Q ( 3 k b f 2 6 k α β ( a 2 c 2 ) ( x c t ) ) + f 3 k (30)

当解为(28)时,得到行波解

u 3 ( x , t ) = ± 1 α α b k β Q ( b 2 α β ( a 2 c 2 ) ( x c t ) ) Q ( b 2 α β ( a 2 c 2 ) ( x c t ) ) + 2 f 3 k (31)

u 4 ( x , t ) = ± 1 α α b k β Q ( b 2 α β ( a 2 c 2 ) ( x c t ) ) Q ( b 2 α β ( a 2 c 2 ) ( x c t ) ) + 2 f 3 k (32)

如果我们取 A ( t ) = 1 , B ( t ) = 0 , C ( t ) = α 2 ,则可得到黎卡提方程 d d t y ( t ) = y 2 ( t ) α 2 ,解为 y ( t ) = α + 2 α C 1 e 2 α t C 2 C 1 e 2 α t ,其中 C 1 , C 2 是任意常数 [7]。这时解方程组(7)得

{ a 0 = f 3 k a 1 = ± 3 k b + f 2 3 k 2 α 2 μ 2 = 3 k b + f 2 6 k α 2 ( c 2 a 2 ) (33)

(34)

其中c和 α ( c a )为任意常数。当解为(33)时,得到行波解

u 1 ( x , t ) = ± 3 k b + f 2 3 k 2 α 2 ( α + 2 α C 1 e 2 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) C 2 C 1 e 2 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) ) + f 3 k (35)

C 1 , α 取1, C 2 取2, k , b , f , c 取1,a取0时,解(35)的图像为图6

Figure 6. Image of solution (35)

图6. 解(35)的图像

u 2 ( x , t ) = ± 3 k b + f 2 3 k 2 α 2 ( α + 2 α C 1 e 2 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) C 2 C 1 e 2 3 k b + f 2 6 k ( c 2 a 2 ) ( x c t ) ) + f 3 k (36)

当解为(34)时,得到行波解

u 3 ( x , t ) = ± b k α 2 ( α + 2 α C 1 e 2 b 2 ( c 2 a 2 ) ( x c t ) C 2 C 1 e 2 b 2 ( c 2 a 2 ) ( x c t ) ) + 2 f 3 k (37)

u 4 ( x , t ) = ± b k α 2 ( α + 2 α C 1 e 2 b 2 ( c 2 a 2 ) ( x c t ) C 2 C 1 e 2 b 2 ( c 2 a 2 ) ( x c t ) ) + 2 f 3 k (38)

本文主要利用黎卡提方程作为辅助方程,研究了一类波动方程的行波解问题,得到了一些新的行波解的形式。相比原有的方法,该方法更加容易理解与计算。

文章引用

吕月娥,刘 琪. 广义Klein-Gordon方程的行波解
Travelling Wave Solutions of a Class of Generalized Klein-Gordon Equations[J]. 应用数学进展, 2019, 08(11): 1689-1699. https://doi.org/10.12677/AAM.2019.811198

参考文献

  1. 1. Wazwaz, A.-M. (2005) The Tanh Method for a Reliable Treatment of the K(n, n) and the KP(n, n) Equations and Its Variants. Applied Mathematics and Computation, 170, 361-379.
    https://doi.org/10.1016/j.amc.2004.12.008

  2. 2. Zhang, S., Tong, J. and Wang, W. (2008) Generalized G'/G-Expansion Method for the mKDV Equation with Variable Coefficients. Physics Letters A, 372, 2254-2257.
    https://doi.org/10.1016/j.physleta.2007.11.026

  3. 3. Wang, M.L., Li X.Z. and Zhang, J.L. (2008) The G'/G-Expansion Method and Travelling Wave Solutions of Nonlinear Evolution Equations in Mathematical Physics. Physics Letters A, 372, 417-423.
    https://doi.org/10.1016/j.physleta.2007.07.051

  4. 4. 陈向华, 白晓东. 常微分方程[M]. 内蒙古: 内蒙古大学出版社, 2002.

  5. 5. 张玮玮. 一类特殊类型的Riccati方程的求解[J]. 安庆师范学院学报(自然科学版), 2015, 21(2): 110-111.

  6. 6. 李晓琴. 二阶非线性微分方程Riccati方程的解法及应用[J]. 聊城大学学报, 2012, 25(2): 19-21.

  7. 7. 李迎娣. 二阶变系数线性齐次微分方程的一些解法[J]. 西北民族大学学报(自然科学版), 2015, 936(3): 1-2.

期刊菜单