非线性时间分数阶微分方程的exp(-Φ(ξ))解法

U (ξ ) = u ( x, y , z, t ) ; ξ = x + y + z − rt α , Γ (1 + α )
(9)
517
王亚东，张新东
其中 r 是一个任意的非零常数，通过在(7)中使用链式规则，我们可以得到
Dtα u = σ t
du α Dt ξ , dξ
σ t 称为分形指标，我们通常令 σ t = r ，其中 r 是一个常数。
The Exp(-Φ(ξ)) Method for the Nonlinear Time Fractional Differential Equations
Yadong Wang, Xindong Zhang*
College of Mathematics Sciences, Xinjiang Normal University, Urumqi Xinjiang Received: Jun. 29th, 2017; accepted: Jul. 17th, 2017; published: Jul. 20th, 2017
Advances in Applied Mathematics 应用数学进展, 2017, 6(4), 515-522 Published Online July 2017 in Hans. /journal/aam https:///10.12677/aam.2017.64062
2(λ (ξ + c ) + 2) ln − Φ 4 (ξ ) = , λ 2 (ξ + c )
(16)
⑤ 当 λ 2 − 4 µ = 0, µ = 0, λ = 0 时，
Φ 5 (ξ ) = ln (ξ + c ) ,
2. 准备工作
本节我们主要介绍一些文章中将用到的关于分数阶微积分的一些预备知识。 具体请参考文献[15] [16] [17] [18]。 定义 2.1 假设 φ ( x ) 是定义在包含 x0 和 x0 领域内的所有点的区间上的函数，如果对于任意的正常数
ε ，存在正常数 p 、 δ ，当 x − x0 < δ 时，有
Keywords
Time Fractional KdV-ZK Equation, Fractional Derivative, Exact Solutions, Exp(-Φ(x)) Method
非线性时间分数阶微分方程的exp(-Φ(ξ))解法
王亚东，张新东*
新疆师范大学数学科学学院，新疆 乌鲁木齐
dη Eη kxη dxη
( ) = kE
η
( kx ) ， k 为常数。
η
这些导数在研究高度不规则和不可微函数的分数阶可微性方面发挥着重要的作用。
3. Exp(-Φ(ξ))方法的基本思想
本节我们考虑非线性分数阶微分方程的一般形式：
F u, Dtα u, u x , u y , uz , u xx , u yy , uzz , u xy , u xz , = 0 ,
Abstract
In this paper, the exp(-Φ(ξ)) method is used to construct the approximate analytical solution of Korteweg-de Vries Zakharov Kuznetsov (KdV-ZK) equation via the (3 + 1)-dimensional. The results of numerical example show that the method is very useful in solving nonlinear time fractional differential equations.
利用方程式(9)将程式(8)转化为下面的非线性常微分方程
Q (U ,U ′,U ′′,) = 0 ,
(10)
其中 ' 表示 U 于 ξ 的导数。下面我们对方程式(10)进行逐项积分。 步骤 2：根据 exp ( −Φ (ξ ) ) 的方法，方程(10)的近似解析解可以用以下的形式表示：
U (ξ ) =
φ ( x ) − φ ( x0 ) < pε η , 0 < η < 1 ,
则称 φ ( x ) 在 x = x0 处局部分数连续，即满 ) = φ ( x0 ) ,
因此，如果
φ ( x ) ∈ Cη ( a, b ) ,
(2)
516
王亚东，张新东
其中 0 < η < 1 是分形维数，则称函数 φ ( x ) 在区间 ( a, b ) 上局部分数阶连续。 定义 2.2 函数 φ ( x ) : R → R, X φ ( X ) 是指数为 η 的不可微函数，并关于指数 η 满足 Hölder 方程，则 当 x, y ∈ X 时，我们有
(5)
其中 ∆η (φ ( x ) − φ ( x0 ) ) ≅ Γ (1 + η ) ∆ (φ ( x ) − φ ( x0 ) ) 。 注：如果 ω ( x ) = (ϕ u )( x ) ，则当 u ( x ) = h ( x ) 时，存在 φ (η ) ( h ( x ) ) ， h (1) ( x ) 和 φ (1) ( h ( x ) ) ， h (η ) ( x ) 使 得
λ 2 − 4µ − λ 2 − 4 µ tanh ξ + c) − λ ( 2 , Φ1 ( ξ ) = ln µ 2
(13)
② 当 λ 2 − 4 µ < 0, µ ≠ 0 时，
4µ − λ 2 − 4 µ − λ 2 tan ξ + c) − λ ( 2 , Φ 2 (ξ ) = ln 2µ
Copyright © 2017 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). /licenses/by/4.0/
(
)
(8)
其中 Dtα 是在关与 u 关于 x, y , z, t 的局部分数导数， F 是 u = u ( x, y, z , t ,) 和 u 的各阶导数的多项式，其中 包含非线性项和 u 的最高阶导数。寻求分数阶微分方程的近似解析解的主要步骤如下[19]： 步骤 1：考虑下面的非线性分数阶复变换，使用分数阶变量变换[20] [21] [22] [23]：
∑ n = 0 an ( exp ( −Φ (ξ ) ) )
N
n
,
(11)
这里 an ( an ≠ 0 ) 是待定常数， Φ ( x ) 满足下列辅助常微分方程：
Φ′ (ξ ) = exp ( −Φ (ξ ) ) + µ exp ( Φ (ξ ) ) + λ ,
(12)
因此辅助方程式(12)有如下不同形式的解： ① 当 λ 2 − 4 µ > 0, µ ≠ 0 时，
φ ( x ) − φ ( x0 ) =O ( x − x0 ) .
η
(
)
定义 2.4 假设 φ ( x ) ∈ Cη ( a, b ) ，则 φ ( x ) 在 x = x0 处的 η 阶局部分数导数定义如下：
= φ ( x0 ) ∆η (φ ( x ) − φ ( x0 ) ) dηφ ( x ) = lim , 0 <η <1, η dxη x = x x → x0 x − x ( ) 0 0
收稿日期：2017年6月29日；录用日期：2017年7月17日；发布日期：2017年7月20日
摘
要
本文主要研究(3+1)维Korteweg-de Vries Zakharov (KdV-ZK)方程的exp(-Φ(ξ))解法。利用exp(-Φ(ξ)) 方法获得所研究方程的近似解析解。数值算例表明，该方法在求解非线性分数阶微分方程的近似解析解 时非常有效。
(14)
③ 当 λ 2 − 4 µ > 0, µ = 0, λ ≠ 0 时，
λ Φ 3 (ξ ) = − ln , sinh ( λ (ξ + c ) ) + cosh ( λ (ξ + c ) ) − 1
(15)
④ 当 λ 2 − 4 µ > 0, µ ≠ 0, λ ≠ 0 时，
*
通讯作者。
文章引用: 王亚东, 张新东. 非线性时间分数阶微分方程的 exp(-Φ(ξ))解法[J]. 应用数学进展, 2017, 6(4): 515-522. https:///10.12677/aam.2017.64062
王亚东，张新东
关键词
时间分数阶KdV -ZK方程，分数阶导数，精确解，exp(-Φ(ξ))方法
φ ( x ) − φ ( x0 ) < C x − y , 0 < η < 1 .
定义 2.3 如果函数 φ ( x ) : R → R, X φ ( X ) 满足
(3)
φ ( x ) − φ ( x0 ) = O ( x − x0 ) , 0 < η < 1 ,
η
(4)
则称 φ ( x ) 是关于指数 η 局部分数阶连续的。 注：函数 φ ( x ) 在空间 Cη ( a, b ) 中，当且仅当对任意的 x0 ∈ [a, b] ， 0 < η < 1 ，满足
dηω ( x ) dx dx
注：
η
= φ(
( (h ( x ) h
η)
(1)
( x ) )) ( x ) ))
η
, .