杂志的下一篇文章
第三类退化Daehee数
下一篇特刊文章
关于一个单价函数的长度问题
期刊上的上一篇文章
定向代数与Hochschild上同调群
 
 
订购文章重印
字体类型:
宋体 佐治亚州 宋体,Verdana
字体大小:
澳大利亚 澳大利亚 澳大利亚
行距:
列宽:
背景:
第条

用Legendre小波求解分数阶微分方程组的一个新的分数阶导数运算矩阵

通过
艾丁·塞瑟
1,*
塞尔维·阿尔顿
2
1
土耳其伊斯坦布尔尤尔迪斯技术大学数学工程系,34200
2
土耳其伊斯坦布尔尤尔迪斯技术大学,34200
*
信件应寄给的作者。
数学 2018,6(11), 238;https://doi.org/10.3390/math6110238
收到的提交文件:2018年10月2日/修订日期:2018年10月29日/接受日期:2018年10月31日/发布日期:2018年11月5日

摘要

:
本文介绍了一种用勒让德小波运算矩阵法(LWOMM)求解分数阶微分方程组(FDE)的新的数值方法。我们首先利用勒让德小波和移位勒让德多项式的一些显著特征,在一些特殊条件下,建立了分数阶导数的运算矩阵。然后,利用这些运算矩阵将分数阶微分方程组转换为代数方程组。在本文的最后,通过几个例子说明了该方法的有效性和正确性。将该方法与几种公认的方法进行了比较,结果表明,勒让德小波运算矩阵方法的优点在于其计算的准确性和可理解性。

1.简介

微分算子和积分算子是数学模型的基础,它们也被用来理解自然和人工系统的工作原理。因此,微分方程和积分方程在理论和实践上都具有重要意义。这样的方程有着广泛的应用,包括在物理科学(如物理和工程)以及社会科学中。微分方程组,如微分方程,经常用于弹性理论、动力学、流体力学、振荡和量子动力学等问题。
对微分算子和积分算子的兴趣导致了对分数微分算子和整数算子的探索,进一步深入研究了这些问题。由于一个问题,分数微积分的起源出现在1695年莱布尼茨给L'Hópital的一封信中。近年来,分数微积分因其能够简化许多物理、工程和经济现象而备受关注,例如流体动力学交通模型、阻尼定律、连续介质和统计力学、扩散过程、固体力学、控制理论、有色噪声、粘弹性、电化学、,电磁等。
由于分数阶微分方程(FDE)的各种解无法解析地找到,因此需要数值和近似方法。许多研究人员在数值求解FDE和此类方程组时研究了许多技术。其中一些技术是Song和Wang提出的Adomian分解方法[1]配置法,改进的运算矩阵法[2,,4],∑enol和Dolapçonf引入的摄动迭代法[5],计算矩阵方法,如Khader等人[6]Ertürk和Momani证明的微分变换方法[7]Gupta等人给出的变分迭代法和拉普拉斯变换法[8]以及Ghazanfari和Ghazanfari研究的分数复变换方法[9]等等。Kilbas等人[10]包括研究分数阶微分方程和分数阶积分微分方程。此外,利用勒让德多项式运算矩阵法给出了有限差分方程及其方程组的数值解[11],伯恩斯坦运算矩阵法[12],Genocchi运算矩阵法[13]雅可比运算矩阵法[14]切比雪夫小波运算矩阵法[15]多项式最小二乘法(PLSM)[16],勒让德类小波运算矩阵法(LWPT)[17]和Genocchi类小波运算矩阵方法[18].
本文利用勒让德小波运算矩阵方法对分数阶微分方程组进行了数值分析。所提出的方法最重要的优点是,它提供了一个可以理解的过程,可以将FDE和此类方程组简化为代数方程组。首先,我们首先介绍一些基本定义和基本关系第2节第3节分别是。然后,在第4节,分数阶导数的运算矩阵自然地表示为分数阶微分方程的线性和非线性系统。第5节给出了用所介绍的方法进行测试的五个示例。最后,最后一部分包括结论。

2.基本定义

Liouville_Caputo分数阶导数、移位勒让德多项式和勒让德小波的定义如下[19,20].
定义 1
Liouville_Caputo分数导数 u个 定义为[19]
D类 α u个 ( t吨 ) = 1 Γ ( n个 α ) 0 t吨 u个 ( n个 ) ( ξ ) ( t吨 ξ ) α + 1 n个 d日 ξ , n个 1 < α n个 ,   n个 N个 .
Liouville_Caputo分数导数的一些特征如下:
D类 α C类 = 0 ,
哪里 C类 是一个常量。此外,还有
D类 α t吨 β = { 0 , β N个 0   n个 d日   β < α Γ ( β + 1 ) Γ ( β + 1 α ) x β α , β N个 0   n个 d日   β α   o个 第页   β N个   n个 d日   β > α ,
在哪儿 α α 分别表示最大整数小于或等于 α ,且最小整数大于或等于 α .
Liouville_Caputo分数阶导数是整数阶导数的线性运算
D类 α ( η u个 ( t吨 ) + ζ v(v) ( t吨 ) ) = η D类 α u个 ( t吨 ) + ζ D类 α v(v) ( t吨 ) ,
哪里 η ζ 都是常量。
定义 2
b条 分别是一个称为母小波的单一函数的伸缩和平移参数。如果 b条 连续变换,得到以下连续小波族[21,22]:
ψ b条 ( t吨 ) = | | 1 / 2 ψ ( t吨 b条 ) , , b条 R(右) , 0 .
定义 三。
P(P) ( t吨 ) 隐含移位的勒让德多项式 .然后 P(P) ( t吨 ) 可以表述为[21]
P(P) ( t吨 ) = k个 = 0 ( 1 ) + k个 ( + k个 ) ! ( k个 ) ! t吨 k个 ( k个 ! ) 2 ,
正交条件为
0 1 P(P) ( t吨 ) P(P) n个 ( t吨 ) d日 t吨 = { 1 2 + 1 , (f) o个 第页   = n个 0 , (f) o个 第页   n个 .
定义 4
n个 k个 为任意正整数, 是移位的勒让德多项式的阶,以及 t吨 是标准化时间。然后是勒让德小波 ψ n个 ( t吨 ) = ψ ( k个 , n个 , , t吨 ) 在间隔上定义 [ 0 , 1 ] 由[21,22].
ψ n个 ( t吨 ) = { 2 k个 + 1 2 + 1 2 P(P) ( 2 k个 t吨 n个 ) , n个 2 k个 t吨 n个 + 1 2 k个 0 , o个 t吨 小时 e(电子) 第页 w个 e(电子) ,
哪里 = 0 , 1 , , M(M) ;   n个 = 0 , 1 , , ( 2 k个 1 ) .系数 + 1 2 表示正交。
定义 5
u个 ( t吨 ) v(v) ( t吨 ) 是定义的函数 [ 0 , 1 ] 然后根据勒让德小波展开为[21,22]
u个 ( t吨 ) = n个 = 0 = 0 c n个 ψ n个 ( t吨 ) ,
哪里 c n个 = ( u个 ( t吨 ) , ψ n个 ( t吨 ) ) ,其中 ( . , . ) 意味着内部产品。如果等式(9)中的无穷级数被截断,那么它可以表示为
u个 ( t吨 ) n个 = 0 2 k个 1 = 0 M(M) c n个 ψ n个 ( t吨 ) = C类 T型 ψ ( t吨 ) ,
哪里 C类 ψ ( t吨 ) 是矩阵,如下所示
C类 = [ c 0 , 0 , c 0 , 1 , , c 0 , M(M) , , c 2 , M(M) , , c 2 k个 1 , 0 , c 2 k个 1 , 1 , , c 2 k个 1 , M(M) ] T型 ψ = [ ψ 0 , 0 , ψ 0 , 1 , , ψ 0 , M(M) , , ψ 2 , M(M) , , ψ 2 k个 1 , 0 , ψ 2 k个 1 , 1 , , ψ 2 k个 1 , M(M) ] T型 .

3.基本关系

萨达曼迪和德汉[11]利用移位勒让德多项式导出分数阶导数的运算矩阵。在本节中,我们展示了如何通过借鉴Saadatmandi和Dehghan,在某些特殊条件下导出分数阶导数的Legendre小波运算矩阵[11]. 此外,与Mohammadi所示导数的Legendre小波运算矩阵相关的定理和推论[21]此处引用如下。
定理 1
ψ ( t吨 ) 是方程(8)中引入的勒让德小波向量。然后 ψ ( t吨 ) 表示为[21,22]
d日 ψ ( t吨 ) d日 t吨 = D类 ψ ( t吨 ) ,
哪里 D类 2 k个 ( M(M) + 1 ) 导数的运算矩阵,可以表示为
D类 = [ U型 O(运行) O(运行) O(运行) U型 O(运行) O(运行) O(运行) U型 ] ,
哪里 U型 是一个 ( M(M) + 1 ) ( M(M) + 1 ) 矩阵及其 ( 第页 , ) 第个 元素写入为
U型 第页 , = { 2 k个 + 1 ( 2 第页 1 ) ( 2 1 ) , 第页 = 2 , , ( M(M) + 1 ) ,   = 1 , , 第页 1     ( 第页 + )   古怪的 0 , 否则 .
推论 1
使用方程(12) n个 第个 导数可以表述为[21]
d日 n个 ψ ( t吨 ) d日 t吨 n个 = D类 n个 ψ ( t吨 ) ,
哪里 D类 n个 n个 第个 矩阵的幂 D类 .
引理 1
ψ ( t吨 ) 是方程(8)中引入的勒让德小波向量。假设 k个 = 0 ,然后
D类 α ψ 第页 ( t吨 ) = 0 ,   第页 = 0 , 1 , , α 1 ,   α > 0 .
证明。 
通过使用方程式(8)中的方程式(2)和(4),可以获得所需结果。□
定理 2
ψ ( t吨 ) 是方程(8)中引入的勒让德小波向量。假设 k个 = 0 α > 0 ,然后
D类 α ψ ( t吨 ) D类 ( α ) ψ ( t吨 ) ,
哪里 D类 ( α ) ( M(M) + 1 ) x ( M(M) + 1 ) 分数阶导数的运算矩阵 α > 0 ,   N个 1 < α N个 在Liouville_Caputo意义上,可以说
D类 ( α ) = ( 0 0 0 0 0 0 小时 = α α ξ α , 0 , 小时 小时 = α α ξ α , 1 , 小时 小时 = α α ξ α , , 小时 小时 = α 第页 ξ 第页 , 0 , 小时 小时 = α 第页 ξ 第页 , 1 , 小时 小时 = α 第页 ξ 第页 , , 小时 小时 = α ξ , 0 , 小时 小时 = α ξ , 1 , 小时 小时 = α ξ , , 小时 ) ,
哪里 ξ 第页 , , 小时 写为
ξ 第页 , , 小时 = 2 第页 + 1 2 + 1 = 0 ( 1 ) 第页 + + 小时 + ( 第页 + 小时 ) ! ( + ) ! ( 第页 小时 ) ! 小时 ! Γ ( 小时 α + 1 ) ( ) ! ( ! ) 2 ( 小时 + α + 1 ) .
在中考虑 D类 ( α ) 第一个 α 行都为零.
证明。 
假设 ψ 第页 ( t吨 ) 第页 第个 向量的元素 ψ ( t吨 ) 在方程式(11)中引入,其中 第页 = n个 M(M) + ( + 1 ) ,   = 0 , 1 , , M(M) ,   n个 = 0 , 1 , , ( 2 k个 1 ) .然后 ψ 第页 ( t吨 ) 可以表示为
ψ 第页 ( t吨 ) = 2 k个 + 1 2 第页 + 1 2 P(P) 第页 ( 2 k个 t吨 n个 ) χ [ n个 2 k个 , n个 + 1 2 k个 ] .
接受这一点 k个 = 0 ,通过使用移位的勒让德多项式,我们得到
ψ 第页 ( t吨 ) = 2 第页 + 1 2 小时 = 0 第页 ( 1 ) 第页 + 小时 ( 第页 + 小时 ) ! ( 第页 小时 ) ! ( 小时 ! ) 2 t吨 小时 χ [ 0 , 1 ] .
如果我们使用方程(3)、(4)和(21),那么我们有
D类 ( α ) ψ 第页 ( t吨 ) = 2 第页 + 1 2 小时 = 0 第页 ( 1 ) 第页 + 小时 ( 第页 + 小时 ) ! ( 第页 小时 ) ! ( 小时 ! ) 2 D类 α ( t吨 小时 ) χ [ 0 , 1 ] = 2 第页 + 1 小时 = α 第页 ( 1 ) 第页 + 小时 ( 第页 + 小时 ) ! ( 第页 小时 ) ! ( 小时 ! ) Γ ( 小时 α + 1 ) t吨 小时 α χ [ 0 , 1 ] , 第页 = α , , .
近似值 t吨 小时 α 通过 ( + 1 ) Legendre小波的项,然后我们得到
t吨 小时 α b条 小时 , ψ ( t吨 ) ,
哪里
b条 小时 , = 0 1 t吨 小时 α ψ ( t吨 ) d日 t吨 = 2 + 1 2 = 0 ( 1 ) + ( + ) ! ( ) ! ( ! ) 2 0 1 t吨 小时 + α d日 t吨 = 2 + 1 = 0 ( 1 ) + ( + ) ! ( ) ! ( ! ) 2 ( 小时 + α + 1 ) .
利用方程(22)和(24),我们得到
D类 α ψ 第页 ( t吨 ) 2 第页 + 1 小时 = α 第页 = 0 ( 1 ) 第页 + 小时 ( 第页 + 小时 ) ! ( 第页 小时 ) ! ( 小时 ! ) Γ ( 小时 α + 1 ) b条 小时 , ψ ( t吨 ) χ [ 0 , 1 ] = = 0 ( 小时 = α 第页 ξ 第页 , , 小时 ) ψ ( t吨 ) χ [ 0 , 1 ] , 第页 = α , , ,
在哪儿 ξ 第页 , , 小时 如方程式(19)所示。此外,如果我们使用引理1,那么我们可以写
D类 α ψ 第页 ( t吨 ) = 0 , 第页 = 0 , 1 , , α 1 , α > 0 .
结合方程(25)和(27),可以得到结果。□

4.分数阶微分方程的求解系统

在本节中,采用勒让德小波运算矩阵方法获得分数阶微分方程组的数值解。考虑一个分数阶微分方程系统,如下所示:
D类 η 1 u个 1 ( t吨 ) = U型 1 ( t吨 , u个 1 , u个 2 , , u个 ) , D类 η 2 u个 2 ( t吨 ) = U型 2 ( t吨 , u个 1 , u个 2 , , u个 ) , D类 η n个 u个 ( t吨 ) = U型 ( t吨 , u个 1 , u个 2 , , u个 ) ,
哪里 U型 是的线性/非线性函数 t吨 , u个 1 , u个 2 , , u个 , D类 η 是的衍生物 u个 按照的顺序 η 在Liouville–Caputo意义上 N个 1 η < N个 ,并受初始条件的约束
u个 1 ( t吨 0 ) = u个 10 , d日 u个 1 d日 t吨 ( t吨 0 ) = u个 11 , d日 2 u个 1 d日 t吨 2 ( t吨 0 ) = u个 12 , , d日 n个 1 u个 1 d日 t吨 n个 1 ( t吨 0 ) = u个 1 ( n个 1 ) u个 2 ( t吨 0 ) = u个 20 , d日 u个 2 d日 t吨 ( t吨 0 ) = u个 21 , d日 2 u个 2 d日 t吨 2 ( t吨 0 ) = u个 22 , , d日 n个 1 u个 2 d日 t吨 n个 1 ( t吨 0 ) = u个 2 ( n个 1 ) u个 ( t吨 0 ) = u个 0 , d日 u个 d日 t吨 ( t吨 0 ) = u个 1 , d日 2 u个 d日 t吨 2 ( t吨 0 ) = u个 2 , , d日 n个 1 u个 d日 t吨 n个 1 ( t吨 0 ) = u个 ( n个 1 ) .
首先,近似 u个 1 ( t吨 ) , u个 2 ( t吨 ) , , u个 ( t吨 ) D类 η 1 u个 1 ( t吨 ) , D类 η 2 u个 2 ( t吨 ) , , D类 η n个 u个 ( t吨 ) ,我们获得
u个 1 ( t吨 ) n个 = 0 2 k个 1 = 0 M(M) c 1 n个 , ψ n个 , = C类 1 T型 ψ ( t吨 ) u个 2 ( t吨 ) n个 = 0 2 k个 1 = 0 M(M) c 2 n个 , ψ n个 , = C类 2 T型 ψ ( t吨 ) u个 ( t吨 ) n个 = 0 2 k个 1 = 0 M(M) c n个 n个 , ψ n个 , = C类 T型 ψ ( t吨 ) ,
哪里 C类 ,   = 1 , 2 , , 是未知向量,并且 ψ ( t吨 ) 是方程式(8)中引入的向量。如果我们使用方程(17),那么我们有
D类 η 1 u个 1 ( t吨 ) C类 1 T型 D类 ( η 1 ) ψ ( t吨 ) D类 η 2 u个 2 ( t吨 ) C类 2 T型 D类 ( η 2 ) ψ ( t吨 ) D类 η n个 u个 ( t吨 ) C类 T型 D类 ( η n个 ) ψ ( t吨 ) .
将方程(29)和(30)代入方程(27),我们得到
R(右) 1 ( t吨 ) = C类 1 T型 D类 ( η 1 ) ψ ( t吨 ) U型 1 ( t吨 , C类 1 T型 ψ ( t吨 ) , C类 2 T型 ψ ( t吨 ) , , C类 T型 ψ ( t吨 ) ) R(右) 2 ( t吨 ) = C类 2 T型 D类 ( η 2 ) ψ ( t吨 ) U型 2 ( t吨 , C类 1 T型 ψ ( t吨 ) , C类 2 T型 ψ ( t吨 ) , , C类 T型 ψ ( t吨 ) ) R(右) ( t吨 ) = C类 T型 D类 ( η n个 ) ψ ( t吨 ) U型 ( t吨 , C类 1 T型 ψ ( t吨 ) , C类 2 T型 ψ ( t吨 ) , , C类 T型 ψ ( t吨 ) )
如果 U型 是的线性函数 t吨 , u个 1 , u个 2 , , u个 , 然后我们生产 2 k个 ( M(M) + 1 ) n个 通过实现线性方程
0 1 ψ j个 ( t吨 ) R(右) ( t吨 ) d日 t吨 = 0 ,   j个 = 1 , , 2 k个 ( M(M) + 1 ) n个 ,   = 1 , 2 , , .
此外,通过将方程(28)中的初始条件替换为方程(30),我们得到
u个 1 ( t吨 0 ) C类 1 T型 ψ ( t吨 0 ) = u个 10 , d日 u个 1 d日 t吨 ( t吨 0 ) C类 1 T型 D类 ψ ( t吨 0 ) = u个 11 , , d日 n个 1 u个 1 d日 t吨 n个 1 ( t吨 0 ) C类 1 T型 D类 n个 1 ψ ( t吨 0 ) = u个 1 ( n个 1 ) u个 2 ( t吨 0 ) C类 2 T型 ψ ( t吨 0 ) = u个 20 , d日 u个 2 d日 t吨 ( t吨 0 ) C类 2 T型 D类 ψ ( t吨 0 ) = u个 21 , , d日 n个 1 u个 2 d日 t吨 n个 1 ( t吨 0 ) C类 2 T型 D类 n个 1 ψ ( t吨 0 ) = u个 2 ( n个 1 ) u个 ( t吨 0 ) C类 T型 ψ ( t吨 0 ) = u个 0 , d日 u个 d日 t吨 ( t吨 0 ) C类 T型 D类 ψ ( t吨 0 ) = u个 1 , , d日 n个 1 u个 d日 t吨 n个 1 ( t吨 0 ) C类 T型 D类 n个 1 ψ ( t吨 0 ) = u个 ( n个 1 ) .
A类 2 k个 ( M(M) + 1 ) 通过组合方程(32)和(33)生成线性方程组。对于向量的未知系数,可以获得这些线性方程的解 C类 。因此, u个 1 ( t吨 ) , u个 2 ( t吨 ) , , u个 ( t吨 ) 可以计算公式(27)中引入的。
如果 U型 是的非线性函数 t吨 , u个 1 , u个 2 , , u个 , 然后我们首先计算 R(右) 1 ( t吨 ) , R(右) 2 ( t吨 ) , , R(右) ( t吨 ) 2 k个 ( M(M) + 1 ) n个 点,为了获得更好的结果,请使用第一个 2 k个 ( M(M) + 1 ) n个 移位勒让德的根 P(P) 2 k个 ( M(M) + 1 ) ( t吨 ) 然后,这些方程与方程(33)一起产生 2 k个 ( M(M) + 1 ) 非线性方程组。这些非线性方程的解可以用牛顿迭代法求得。因此, u个 1 ( t吨 ) , u个 2 ( t吨 ) , , u个 ( t吨 ) 可以计算公式(27)中引入的。

5.示例

在本节中,为了证明所介绍方法的适用性和强大性,我们给出了五个分数阶微分方程线性和非线性系统的解。
例子 1
我们首先考虑以下线性分数阶微分方程组[7,8]:
D类 α u个 ( t吨 ) = u个 ( t吨 ) + v(v) ( t吨 ) D类 α v(v) ( t吨 ) = u个 ( t吨 ) + v(v) ( t吨 ) ,
从属于
u个 ( 0 ) = 0 , v(v) ( 0 ) = 1 .
α = 1 众所周知
u个 ( t吨 ) = e(电子) t吨 t吨 , v(v) ( t吨 ) = e(电子) t吨 余弦 t吨 .
对该示例进行了检查 M(M) = 2 , k个 = 0 、和 α = 0.9 , 0.7 , 0.5 。当获得的结果与精确解匹配时 α = 1 ,如所示图1,我们可以清楚地观察到 α 接近1时,我们的结果接近精确解。我们还使用勒让德多项式运算矩阵法(LPOMM)解决了这个问题,并将结果与LWOMM进行了比较。的数值计算u个(t吨)和v(v)(t吨)何时 α = 0.9 在中显示表1表2.
例子 2
我们考虑了以下非线性分数阶微分方程组[13]:
D类 2 u个 ( t吨 ) = 8 u个 ( t吨 ) + v(v) 2 ( t吨 ) 4 t吨 6 + 4 t吨 + 8 t吨 2 π 1 D类 1 2 v(v) ( t吨 ) = t吨 2 D类 u个 ( t吨 ) + v(v) ( t吨 ) t吨 4 2 t吨 + 32 t吨 5 2 5 π 1
u个 ( 0 ) = 0 , v(v) ( 0 ) = 1 , u个 ( 1 ) = 1 , v(v) ( 1 ) = , u个 ( 0 ) = 0 , u个 ( 1 ) = .
这个系统的精确解是
u个 ( t吨 ) = t吨 , v(v) ( t吨 ) = 2 t吨 + 1
使用参数 M(M) = k个 = 0 ,我们将所提出的方法和LPOMM应用于解决此问题,并表明我们的方法更加有效和有用。我们的数值结果支持这样的观点,即我们的解比近似解LPOMM更接近精确解。近似解和精确解的比较见表3表4.
例子 三。
我们考虑了具有初始条件的非线性分数阶微分方程组[8]
D类 α u个 ( t吨 ) = u个 ( t吨 ) 2 D类 α v(v) ( t吨 ) = u个 2 ( t吨 ) + v(v) ( t吨 ) ,
u个 ( 0 ) = 1 , v(v) ( 0 ) = 0 .
α = 1 众所周知
u个 ( t吨 ) = e(电子) ( t吨 2 ) , v(v) ( t吨 ) = t吨 e(电子) t吨 .
参数 M(M) = 2 ,   k个 = 0 、和 α = 0.5 , 0.7 , 0.9 已被利用。我们的结果与精确解的比较 α = 1 显示在中图2数据支持了这一点 α 近似值为1时,我们的结果近似于精确解。我们还使用LPOMM解决了这个问题,并将结果与LWOMM进行了比较。最后,我们给出了u个(t吨)和v(v)(t吨)何时 α = 0.9 在里面表5表6.
例子 4
我们考虑了具有初始条件的非线性FDE系统[13]
D类 α u个 ( t吨 ) = 1002 u个 ( t吨 ) + 1000 v(v) 2 ( t吨 ) D类 α v(v) ( t吨 ) = u个 ( t吨 ) v(v) ( t吨 ) v(v) 2 ( t吨 )
u个 ( 0 ) = 1 , v(v) ( 0 ) = 1
α = 1 众所周知
u个 ( t吨 ) = e(电子) 2 t吨 , v(v) ( t吨 ) = e(电子) t吨
对该示例进行了分析 M(M) = , k个 = 0 、和 α = 0.9 , 0.7 , 0.5 当所获得的结果与精确解相匹配时 α = 1 ,如所示图3,我们可以清楚地观察到 α 接近1时,我们的结果接近精确解。我们还使用LPOMM解决了这个问题,并将结果与LWOMM进行了比较。的数值计算u个(t吨)和v(v)(t吨)何时 α = 0.99 也显示在表7表8.
例子 5
我们考虑了以下分数阶布鲁塞尔系统[16,17]:
D类 α u个 ( t吨 ) = 2 u个 ( t吨 ) + u个 2 ( t吨 ) v(v) ( t吨 ) D类 α v(v) ( t吨 ) = u个 ( t吨 ) u个 2 ( t吨 ) v(v) ( t吨 )
u个 ( 0 ) = 1 , v(v) ( 0 ) = 1 .
α = 1 α = 0.98 由Chang和Isah使用LWPT提交[17]Bota和Caruntu使用PLSM[16]. 这些解决方案 α = 98 由提供
u个 L(左) W公司 P(P) T型 ( t吨 ) = 1 1.0791 t吨 + 0.2711 t吨 2 0.0638 t吨 , v(v) L(左) W公司 P(P) T型 ( t吨 ) = 1 + 0.0151 t吨 + 0.4185 t吨 2 0.2624 t吨 u个 P(P) L(左) S公司 M(M) ( t吨 ) = 1 1.08655 t吨 + 0.311138 t吨 2 + 0.0243682 t吨 , v(v) P(P) L(左) S公司 M(M) ( t吨 ) = 1 + 0.0349127 t吨 + 0.333424 t吨 2 0.184414 t吨 .
参数 M(M) = 2 , k个 = 0 、和 α = 0.98 使用了。我们的结果与Bota和Caruntu引入的近似解的比较[16]还有Chang和Isah[17]何时 α = 0.98 显示在中图4最后,我们还提供了以下方面的数值计算u个(t吨)和v(v)(t吨)何时 α = 0.98 在里面表9表10.

6.结论

本文从分数阶导数的一个新的运算矩阵出发,在某些特殊条件下,对分数阶微分方程组进行了研究。为了在Maple中求解分数阶微分方程的线性和非线性系统,我们还系统化了一种非常实用的算法。Maple生成的所有数值结果和图形表示都证实了Legendre小波运算矩阵方法是非常有效和适用的。作为下一步,本文介绍的方法可以应用于分数阶偏微分方程及其系统、分数阶积分方程及其系统以及分数阶积分微分方程。这些方程至少与分数阶微分方程一样重要,在科学、工程和技术中具有非常重要的意义。

作者贡献

所有作者对本文的写作贡献均等。所有作者阅读并批准了最终手稿。

基金

这项研究没有得到外部资助。

利益冲突

作者声明没有利益冲突。

工具书类

  1. 宋,L。;Wang,W.一种新的改进Adomian分解方法及其在分数阶微分方程中的应用。申请。数学。模型。 2013,37, 1590–1598. [谷歌学者] [交叉参考]
  2. Balaji,S.Legendre小波运算矩阵法求解分数阶Riccati微分方程。J.埃及。数学。Soc公司。 2015,23, 263–270. [谷歌学者] [交叉参考]
  3. Rehman,M。;Khan,R.A.解分数阶微分方程的Legendre小波方法。Commun公司。非线性科学。数字。模拟。 2011,16, 4163–4173. [谷歌学者] [交叉参考]
  4. Mohammadi,F。;侯赛尼,M.M。;Mohyud-Din,S.T.勒让德小波的新运算矩阵及其在分数阶边值问题中的应用。国际期刊系统。科学。 2011,32, 7371–7378. [谷歌学者] [交叉参考]
  5. ∑烯醇,M。;伊利诺伊州多拉普松。关于分数阶微分方程的扰动迭代算法。沙特国王大学。 2016,28, 69–74. [谷歌学者] [交叉参考]
  6. 卡德尔,M.M。;达纳夫,T.S。;Hendy,A.S.一种求解高阶分数阶微分方程组的计算矩阵方法。申请。数学。模型。 2013,37, 4035–4050. [谷歌学者] [交叉参考]
  7. 埃蒂尔克,V.S。;Momani,S.使用微分变换方法求解分数阶微分方程组。J.计算。申请。数学。 2008,215, 142–151. [谷歌学者] [交叉参考]
  8. 古普塔,S。;Kumar博士。;辛格,J.通过拉普拉斯变换对分数阶微分方程组的数值研究。J.埃及。数学。Soc公司。 2015,23, 256–262. [谷歌学者] [交叉参考]
  9. 加扎法里,B。;Ghazanfari,A.G.用分数复变换法求解分数微分方程组。亚洲J.应用。科学。 2012,5, 438–444. [谷歌学者] [交叉参考]
  10. 基尔巴斯,A.A。;斯利瓦斯塔瓦,H.M。;J.J.特鲁希略。分数阶微分方程的理论与应用第1版。;爱思唯尔科学:荷兰阿姆斯特丹,2006年;第540页。[谷歌学者]
  11. Saadatmandi,A。;Dehghan,M.解分数阶微分方程的新运算矩阵。计算。数学。申请。 2010,59, 1326–1336. [谷歌学者] [交叉参考]
  12. 阿尔什布尔,M.H.T。;巴塔尼赫,A.S。;哈希姆,I。;基于新分数Bernstein函数的运算矩阵的分数阶微分方程的解。沙特国王大学。 2017,29, 1–18. [谷歌学者] [交叉参考]
  13. Isah,A。;Phang,C.通过Genocchi多项式求解非线性分数阶微分方程的导数的新运算矩阵。沙特国王大学。 2017. [谷歌学者] [交叉参考]
  14. 多哈,E.H。;Bhrawy,A.H。;Ezz-Eldien,S.S.一种新的雅可比运算矩阵:求解分数阶微分方程的应用。申请。数学。模型。 2012,36, 4931–4943. [谷歌学者] [交叉参考]
  15. Mohammadi,F.利用分数阶导数的Chebyshev小波运算矩阵对Bagley-Torvik方程进行数值求解。国际期刊Adv.Appl。数学。机械。 2014,2, 83–91. [谷歌学者]
  16. 博塔,C。;Caruntu,B.使用多项式最小二乘法的分数阶布鲁塞尔振子系统的近似分析解。公司高级差异。埃克。 2015. [谷歌学者] [交叉参考]
  17. Chang,P。;Isah,A.Legendre通过小波多项式变换的分数阶导数小波运算矩阵及其在求解分数阶布鲁塞尔振子系统中的应用。《物理学杂志》。Conf.序列号。 2016,693, 012001. [谷歌学者] [交叉参考]
  18. Isah,A。;Phang,C.Genocchi类小波运算矩阵及其在求解非线性分数阶微分方程中的应用。打开物理。 2016,14, 463–472. [谷歌学者] [交叉参考]
  19. Miller,K.S。;B.罗斯。分数阶微积分和分数阶微分方程简介第1版。;威利:霍博肯,新泽西州,美国,1993年;第384页。[谷歌学者]
  20. Oldham,K.B。;斯潘尼尔,J。分数微积分第1版。;学术出版社:纽约,纽约,美国,1974年;第234页。[谷歌学者]
  21. Mohammadi,F。;Hosseini,M.M.一种新的Legendre小波导数运算矩阵及其在求解奇异常微分方程中的应用。J.弗兰克尔。仪器。 2011,348, 1787–1796. [谷歌学者] [交叉参考]
  22. Mohammadi,F。;侯赛尼,M.M。;Mohyud-Din,S.T.Legendre小波Galerkin方法用于求解具有非解析解的常微分方程。国际期刊系统。科学。 2011,42, 579–585. [谷歌学者] [交叉参考]
图1。 α = 0.9 , 0.7 , 0.5 在示例1中:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
图1。 α = 0.9 , 0.7 , 0.5 在示例1中:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
数学06 00238 g001a数学06 00238 g001b
图2。 α = 0.9 , 0.7 , 0.5 例如3:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
图2。将我们的解与精确解进行比较 α = 0.9 , 0.7 , 0.5 例如3:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
数学06 00238 g002
图3。 α = 0.9 , 0.7 , 0.5 对于示例4:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
图3。将我们的解与精确解进行比较 α = 0.9 , 0.7 , 0.5 对于示例4:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
数学06 00238 g003
图4。 α = 0.98 例如5:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
图4。 α = 0.98 例如5:()我们的解决方案u个(t吨); 和(b条)我们的解决方案v(v)(t吨).
数学06 00238 g004
表1。的数值解u个(t吨)何时 α = 0.9 通过所介绍的方法和实施例1的LPOMM实现。
表1。的数值解u个(t吨)何时 α = 0.9 通过所介绍的方法和实施例1的LPOMM实现。
t吨 u个 L(左) W公司 O(运行) M(M) M(M) u个 L(左) P(P) O(运行) M(M) M(M) 绝对误差
00.3 × 10 9 00.3 × 10 9
0.10.14837843300.14837843250.12 × 10 9
0.20.32176452830.32176452770.633 × 10 9
0.30.52015828620.52015828550.65 × 10 9
0.40.74355970670.74355970590.83 × 10 9
0.50.99196878980.99196878900.8 × 10 9
0.61.2653855361.2653855350.53 × 10 9
0.71.5638099441.5638099430.95 × 10 9
0.81.8872420141.8872420140.733 × 10 9
0.92.2356817482.2356817480.62 × 10 9
12.6091291442.6091291440.3 × 10 9
表2。的数值解v(v)(t吨)何时 α = 0.9 通过所介绍的方法和实施例1的LPOMM实现。
表2。的数值解v(v)(t吨)何时 α = 0.9 通过引入的方法和实施例1的LPOMM获得。
t吨 v(v) L(左) W公司 O(运行) M(M) M(M) v(v) L(左) P(P) O(运行) M(M) M(M) 绝对误差
0110.2 × 10 9
0.11.1522708991.152270900−0.5 × 10 9
0.21.2748018581.2748018580.599 × 10 9
0.31.3675928771.367592878−0.67 × 10 9
0.41.4306439561.430643957−0.13 × 10 8
0.51.4639550941.463955094−0.1 × 10 8
0.61.4675262911.467526293−0.13 × 10 8
0.71.4413575491.441357551−0.171 × 10 8
0.81.3854488661.385448868−0.1461 × 10 8
0.91.2998002441.299800245−0.15 × 10 8
11.1844116801.184411682−0.18 × 10 8
表3。使用该方法获得的数值结果与近似解LPOMM和精确解的比较u个(t吨)例2中。
表3。与近似解LPOMM和精确解相比,使用引入的方法获得的数值结果u个(t吨)例2中。
t吨 精确解决方案 u个 L(左) W公司 O(运行) M(M) M(M) u个 L(左) P(P) O(运行) M(M) M(M)
00−0.12 × 10 9 0
0.10.0010.010000000050.001000000000
0.20.0080.020000000160.008000000000
0.30.0270.037500000210.027000000000
0.40.0640.070000000200.064000000000
0.50.1250.125000000010.125000000000
0.60.2160.210000000000.216000000000
0.70.3430.332499999980.343000000000
0.80.5120.499999999960.512000000000
0.90.7290.719999999930.729000000000
110.999999999891
表4。用该方法得到的数值结果与近似解LPOMM和精确解的比较v(v)(t吨)例2中。
表4。用该方法得到的数值结果与近似解LPOMM和精确解的比较v(v)(t吨)例2中。
t吨 精确解决方案 v(v) L(左) W公司 O(运行) M(M) M(M) v(v) L(左) P(P) O(运行) M(M) M(M)
0110.9999999998
0.11.0021.0348413671.165803114
0.21.0161.0575876281.283854751
0.31.0541.0865376681.374911456
0.41.1281.1399903701.459729770
0.51.2501.2362446181.559066242
0.61.4321.3935992961.693677414
0.71.6861.6303532901.884319831
0.82.0241.9648054822.151750038
0.92.4582.4152547582.516724579
1333
表5。我们的解决方案u个(t吨)何时 α = 0.9 通过所述方法和实施例3的LPOMM获得。
表5。我们的解决方案u个(t吨)何时 α = 0.9 通过所述方法和实施例3的LPOMM获得。
t吨 u个 L(左) W公司 O(运行) M(M) M(M) u个 L(左) P(P) O(运行) M(M) M(M) 绝对误差
011−0.37 × 10 10
0.11.0648163201.064816320−0.132 × 10 9
0.21.1302485881.130248589−0.1375 × 10 9
0.31.1962968071.1962968070.44 × 10 10
0.41.2629609741.262960975−0.808 × 10 9
0.51.3302410901.330241091−0.532 × 10 9
0.61.3981371561.398137156−0.168 × 10 9
0.71.4666491711.466649171−0.756 × 10 9
0.81.5357771341.535777134−0.1375 × 10 9
0.91.6055210461.6055210460.468 × 10 9
11.6758809081.6758809080.163 × 10 9
表6。我们的解决方案v(v)(t吨)何时 α = 0.9 通过所述方法和实施例3的LPOMM获得。
表6。我们的解决方案v(v)(t吨)何时 α = 0.9 通过所述方法和实施例3的LPOMM获得。
t吨 v(v) L(左) W公司 O(运行) M(M) M(M) v(v) L(左) P(P) O(运行) M(M) M(M) 绝对误差
0−0.1 × 10 9 0−0.1 × 10 9
0.10.13817626030.1381762609−0.7 × 10 9
0.20.31336033610.3133603363−0.2 × 10 9
0.30.52555222590.5255522263−0.37 × 10 9
0.40.77475193040.7747519309−0.5 × 10 9
0.51.0609594501.060959450−0.5 × 10 9
0.61.3841747831.384174784−0.8 × 10 9
0.71.7443979321.744397932−0.47 × 10 9
0.82.1416288942.141628895−0.7 × 10 9
0.92.5758676722.5758676720.3 × 10 9
13.0471142643.047114264−0.1 × 10 9
表7。的数值解u个(t吨)何时 α = 0.99 例如示例4。
表7。的数值解u个(t吨)何时 α = 0.99 例如示例4。
t吨 u个 L(左) W公司 O(运行) M(M) M(M) u个 L(左) P(P) O(运行) M(M) M(M) 绝对误差
011−0.4 × 10 10
0.10.81443515290.8144351528−0.21 × 10 10
0.20.66394252330.6639425233−0.15 × 10 9
0.30.54299472290.5429947230−0.34 × 10 9
0.40.44606436360.4460643636−0.42 × 10 9
0.50.36762405680.36762405680
0.60.30214641420.3021464147−0.68 × 10 9
0.70.24410404870.2441040481−0.24 × 10 9
0.80.18796956990.18796956900.53 × 10 9
0.90.12821559200.12821559050.461 × 10 9
10.05931472480.05931472220.144 × 10 8
表8。的数值解v(v)(t吨)何时 α = 0.99 例如示例4。
表8。的数值解v(v)(t吨)何时 α = 0.99 例如示例4。
t吨 v(v) L(左) W公司 O(运行) M(M) M(M) v(v) L(左) P(P) O(运行) M(M) M(M) 绝对误差
010.99999999990.79 × 10 10
0.10.90256018370.90256018370.1498 × 10 9
0.20.81526464870.8152646488−0.119 × 10 9
0.30.73711165770.7371116578−0.8 × 10 10
0.40.66709947370.6670994739−0.11 × 10 9
0.50.60422635970.6042263600−0.24 × 10 9
0.60.54749057850.5474905788−0.19 × 10 9
0.70.49589039320.4958903935−0.26 × 10 9
0.80.44842406640.4484240670−0.443 × 10 9
0.90.40408986140.4040898620−0.5898 × 10 9
10.36188604080.3618860417−0.719 × 10 9
表9。的数值解u个(t吨)何时 α = 0.98 通过引入的方法、LWPT和实施例5的PLSM获得。
表9。的数值解u个(t吨)何时 α = 0.98 通过引入的方法、LWPT和实施例5的PLSM获得。
t吨 u个 L(左) W公司 O(运行) M(M) M(M) u个 L(左) W公司 P(P) T型 u个 P(P) L(左) S公司 M(M)
0111
0.10.89420248260.89473720.8944807482
0.20.79506969160.79451360.7953304656
0.30.70260162680.69894640.7026953614
0.40.61679828830.60765280.6167216448
0.50.53765967610.52025000.5375555250
0.60.46518579020.43635520.4653432112
0.70.39937663060.35558560.4002309126
0.80.34023219730.27755840.3423648384
0.90.28775249020.20189080.2918911978
10.24193750950.12820000.2489562000
表10。的数值解v(v)(t吨)何时 α = 0.98 通过引入的方法、LWPT和实施例5的PLSM获得。
表10。的数值解v(v)(t吨)何时 α = 0.98 通过引入的方法、LWPT和实施例5的PLSM获得。
t吨 v(v) L(左) W公司 O(运行) M(M) M(M) v(v) L(左) W公司 P(P) T型 v(v) P(P) L(左) S公司 M(M)
0111
0.11.0080693071.00543261.006641096
0.21.0194799611.01766081.018844188
0.31.0342319591.03511021.035502792
0.41.0523253041.05620641.055510424
0.51.0737599951.07937501.077760600
0.61.0985360321.10304161.101146836
0.71.1266534151.12563181.124562648
0.81.1581121431.14557121.146901552
0.91.1929122171.16128541.167057064
11.2310536381.17120001.183922700

分享和引用

MDPI和ACS样式

Secer,A。;阿尔顿,南卡罗来纳州。通过勒让德小波求解分数微分方程组的分数导数新运算矩阵。数学 2018,6, 238.https://doi.org/10.3390/math6110238

AMA风格

A区,阿尔顿S区。通过勒让德小波求解分数微分方程组的分数导数新运算矩阵。数学. 2018; 6(11):238.https://doi.org/10.3390/math6110238

芝加哥/图拉宾风格

塞瑟、艾丁和塞尔维·阿尔顿。2018.“通过Legendre小波求解分数微分方程组的分数导数新运算矩阵”数学6,11号:238。https://doi.org/10.3390/math6110238(网址:https://doi.org/10.3390/math6110238)

请注意,从2016年第一期开始,该杂志使用文章编号而不是页码。查看更多详细信息在这里.

文章指标

返回页首顶部