对时间分数阶电报方程构造了一类普遍性差分方法，采用傅里叶方法分析该类差分方法的稳定性和收敛性；最后，通过数值试验验证本文方法求解分数阶电报方程的有效性。选取不同θ值进行比较分析，数值结果表明当θ取0.5附近时，数值解的精度较好；表明普遍性差分方法求解时间分数阶电报方程是有效的。 A kind of universal difference method is constructed for the time fractional telegraph equation. The stability and convergence of the difference method are analyzed by Fourier method. Finally, the effectiveness of the method for solving the fractional telegraph equation is verified by numerical experiments. The comparison of different θ values is carried out. The numerical results show that the numerical solution is better when θ is around 0.5. Therefore, it is effective to solve the time-fractional telegraph equation by the universal difference method.
吴立飞，杨晓忠*
华北电力大学数理学院，北京
收稿日期：2019年8月21日；录用日期：2019年9月5日；发布日期：2019年9月12日
对时间分数阶电报方程构造了一类普遍性差分方法，采用傅里叶方法分析该类差分方法的稳定性和收敛性；最后，通过数值试验验证本文方法求解分数阶电报方程的有效性。选取不同θ值进行比较分析，数值结果表明当θ取0.5附近时，数值解的精度较好；表明普遍性差分方法求解时间分数阶电报方程是有效的。
关键词 :时间分数阶电报方程，普遍性差分方法，傅里叶方法，稳定性，数值试验
Copyright © 2019 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
随着分数阶微积分在工程和其他方面的重要作用，最近几年，更多的学者研究分数阶微积分的求解方法，其中就有很多关于分数阶电报方程的研究 [
随着分数阶电报方程在各个领域中的应用不断增强，并且一般情况下解析解很难给出，因此分数阶电报方程的数值算法的研究也随之不断发展，当前已有许多学者和专家给出了不同的求解分数阶电报方程数值解的方法。综合来看，目前为止最有效、最常用的方法有差分方法、有限元方法、谱方法、同伦分析法等。Ford [
针对时间分数阶电报方程，本文构造了一类普遍性差分方法，采用傅里叶方法分析差分格式的稳定性和收敛性；通过数值算例验证理论分析，说明普遍性差分方法求解分数阶电报方程的有效性。
考虑如下时间分数阶电报方程：
∂ 2 α u ( x , t ) ∂ t 2 α + ∂ α u ( x , t ) ∂ t α = k ∂ 2 u ( x , t ) ∂ x 2 + f ( x , t ) , 0 < x < L , t > 0 , 1 2 < α ≤ 1 (1)
其初始条件为： u ( x , 0 ) = ϕ 0 ( x ) ， ∂ u ( x , 0 ) ∂ t = ϕ 1 ( x ) ，边界条件为： u ( 0 , t ) = 0 ， u ( L , t ) = 0.
为了得到分数阶电报方程的显–隐差分格式，首先对求解区域作一网格剖分：取空间步长 h = L / M 和时间步长 τ = T / N ，其中M，N都是自然数；有 x i = i h ( i = 1 , 2 , ⋯ , M ) , M h = L ； t k = k τ ( k = 1 , 2 , ⋯ , N ) , N τ = T 网格节点为 ( x i , t k ) .
时间分数阶导数可采用L1公式 [
∂ α u ( x i , t k + 1 ) ∂ t α = τ − α Γ ( 2 − α ) ∑ j = 0 k c j ∇ t u ( x i , t k − j ) + o ( τ 2 − α ) (2)
和
∂ 2 α u ( x i , t k + 1 ) ∂ t 2 α = τ − 2 α Γ ( 3 − 2 α ) ∑ j = 0 k - 1 d j δ t 2 u ( x i , t k − j ) + 2 d k ( u ( x i , t 1 ) − u ( x i , t 0 ) − τ u t ( x i ) ) + o ( τ 3 − 2 α ) (3)
当 j = 0 或者 j = k 时上式会出现项 u i − 1 ，用初始边界条件去逼近 u i − 1 ， u i − 1 = u i 1 − 2 k u 1 ( x i ) + O ( τ 3 ) ，误差忽略不计。
空间二阶导数采用中心差分离散：
∂ 2 u ( x i , t k + 1 ) ∂ x 2 = δ x 2 [ θ u ( x i , t k + 1 ) + ( 1 − θ ) u ( x i , t k ) ] + o ( h 2 ) (4)
其中 δ x 2 u ( x i , t k ) = u ( x i − 1 , t k ) − 2 u ( x i , t k ) + u ( x i + 1 , t k ) ， ∇ t u ( x i , t k ) = u ( x i , t k + 1 ) − u ( x i , t k ) ， δ t 2 u ( x i , t k ) = u ( x i , t k − 1 ) − 2 u ( x i , t k ) + u ( x i , t k + 1 ) ， c j = ( j + 1 ) 1 − α − j 1 − α ， d j = ( j + 1 ) 2 − 2 α − j 2 − 2 α ， j = 0 , 1 , 2 , ⋯ , k − 1.
把式(2, 3, 4)代入分数阶电报方程(1)得到其普遍性差分格式为：
τ − 2 α Γ ( 3 − 2 α ) ∑ j = 0 k - 1 c j δ t 2 u ( x i , t k − j ) + τ − α Γ ( 2 − α ) ∑ j = 0 k d j ∇ t u ( x i , t k − j ) = k h 2 δ x 2 [ θ u ( x i , t k + 1 ) + ( 1 − θ ) u ( x i , t k ) ] + f ( x i , t k + 1 ) (5)
其中， 0 ≤ θ ≤ 1 。令 u i k 表示函数 u ( x , t ) 在点 ( x i , t k ) 的值； f i k 表示函数 f ( x , t ) 在点 ( x i , t k ) 的值；设 a = τ − 2 α Γ ( 3 − 2 α ) ， b = τ − α Γ ( 2 − α ) ， r = k h 2 ，则普遍性差分格式可写为如下形式：
a ∑ j = 0 k − 1 d j δ t 2 u i k − j + 2 d k [ u ( x i , t 1 ) − u ( x i , t 0 ) − τ u t ( x i ) ] + b ∑ j = 0 k − 1 c j ∇ t u i k − j = r δ x 2 ( θ u i k + 1 + ( 1 − θ ) u i k ) + f i k + 1 .
普遍性差分格式改成矩阵形式：
{ [ ( 2 a + b ) E + r θ G ] U 1 = [ ( 2 a + b ) E − r ( 1 − θ ) G ] U 0 + 2 a τ H + F 1 , A U k + 1 = B U k + q 1 U k − 1 + ⋯ + q k − 2 U 2 + q k − 1 U 1 + q k U 0 + 2 a d k τ H + F k + 1 . (6)
其中， q 0 = a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) ， q j = a ( − d j − 1 + 2 d j − d j + 1 ) + b ( c j − c j + 1 ) , j = 1 , 2 , 3 , ⋯ , k − 2 ， q k − 1 = a ( − d k − 2 + 2 d k − 1 − 2 d k ) + b ( c k − 1 − c k ) ， q k = a ( 2 d k − d k − 1 ) + b c k ， U k = [ u 1 k , u 2 k , ⋯ , u m − 1 k ] ′ ， U 0 = [ φ 0 ( x 1 ) , φ 0 ( x 2 ) , ⋯ , φ 0 ( x m − 1 ) ] ′ ， H = [ φ 1 ( x 1 ) , φ 1 ( x 2 ) , ⋯ , φ 1 ( x m − 1 ) ] ′ ， F = [ f 1 , f 2 , ⋯ , f m − 1 ] ′ ， A = ( a d 0 + b c 0 ) E + r θ G ， B = ( a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) ) E − r ( 1 − θ ) G ，
G = [ 2 − 1 − 1 2 − 1 ⋱ ⋱ ⋱ ⋱ ⋱ − 1 − 1 2 ] .
由矩阵A定义可知，A是严格对角占优矩阵，所以分数阶电报方程的普遍性差分格式是唯一可解的。
假设 u ˜ i n 是普遍性差分格式的近似解，定义误差 ε i n = u i n − u ˜ i n ， ε 0 n = ε M n = 0 ， ε n = ( ε 1 n , ε 2 n , ⋯ , ε M − 1 n ) 。满足如下误差扰动方程：
{ ( 2 a + b + 2 r θ ) ε i 1 − r θ ( ε i + 1 1 + ε i − 1 1 ) = ( 2 a + b − 2 r ( 1 − θ ) ) ε i 0 + ( 1 − θ ) r ( ε i + 1 0 + ε i − 1 0 ) , ( a + b + 2 r θ ) ε i k + 1 − r θ ( ε i + 1 k + 1 + ε i − 1 k + 1 ) = [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 2 r ( 1 − θ ) ] ε i k + r ( 1 − θ ) ( ε i + 1 k + ε i − 1 k ) + ∑ j = 1 k − 1 q j ε i k − j + q k ε i 0 . (7)
下面使用Fourier方法分析普遍性差分格式的稳定性。在此，定义网格函数
ε n ( x ) = { ε l n , x ∈ ( x l − 1 / 2 , x l + 1 / 2 ] , l = 1 , 2 , 3 , ⋯ , M − 1 , 0 , x ∈ [ L 1 , h / 2 ] ∪ ( L 2 − h / 2 , L 2 ] .
对 ε n ( x ) 进行傅立叶展开：
ε n ( x ) = ∑ j = − ∞ + ∞ v j n e i 2 j π x L ( i = − 1 ) , n = 0 , 1 , 2 , ⋯ , N
其中 v j n = 1 L ∫ L 1 L 2 ε n e i 2 j π x L d x , j = ± 0 , ± 1 , ⋯ .
定义范数： ‖ ε n ( x ) ‖ 2 = ∑ j = 1 M − 1 h | ε j n | 2 ( ε 1 n , ε 2 n , ⋯ , ε M − 1 n ) 。由于 ε 0 n = ε M n = 0 ，有 ‖ ε n ‖ 2 = ∫ L 1 L 2 | ε n ( x ) | 2 d x = ‖ ε n ( x ) ‖ 2 .
从而利用Parseval定理 ‖ ε n ( x ) ‖ 2 = L ∑ j = − ∞ + ∞ | v j n | 2 ，得出
‖ ε n ( x ) ‖ 2 = ∑ j = 1 M − 1 h | ε j n | 2 = L ∑ j = − ∞ + ∞ | v j n | 2 , n = 0 , 1 , ⋯ , N
基于以上分析，可以假设方程误差解的形式为 ε l n = v n e i λ l h ( λ = 2 j π L ) .
当 k = 1 时，
( 2 a + b + 2 r θ − r θ e i λ h − r θ e − i λ h ) v 1 = ( 2 a + b − 2 r ( 1 − θ ) + ( 1 − θ ) r ( e i λ h + e − i λ h ) ) v 0 (8)
因 e ± i λ h = cos λ h ± i sin λ h ， sin 2 λ h 2 = − 1 4 ( e i λ h − 2 + e − i λ h ) ，可得
| v 1 | = | 2 a + b − 4 r ( 1 − θ ) sin 2 λ h 2 | | 2 a + b + 4 r θ sin 2 λ h | | v 0 | .
由 | 2 a + b − 4 r ( 1 − θ ) sin 2 λ h 2 | | 2 a + b + 4 r θ sin 2 λ h | ≤ 1 ，可得 2 r ( 1 − 2 θ ) ≤ 2 a + b ；因此，普遍性差分格式的稳定条件为：当 0 ≤ θ < 1 2 时， r ≤ 2 a + b 2 ( 1 − 2 θ ) ；当 1 2 ≤ θ ≤ 1 时，无条件稳定。
当 k > 1 时，
( a + b + 2 r θ − r θ ( e i λ h + e − i λ h ) ) v k + 1 = [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 2 r ( 1 − θ ) + r ( 1 − θ ) ( e i λ h + e − i λ h ) ] v k + ∑ j = 1 k − 1 q j v k − j + q k v 0 (9)
( a + b + 4 r θ sin 2 λ h ) | v k + 1 | = [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | v k | + ∑ j = 1 k − 1 q j | v k − j | + q k | v 0 | ≤ { | [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | + ∑ j = 1 k − 1 | q j | + | q k | } | v 0 | ≤ { | [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | + b c 1 + a ( d 0 − d 1 ) } | v 0 | ≤ | [ a + b − 4 r ( 1 − θ ) sin 2 λ h ] | | v 0 |
由 | a + b − 4 r ( 1 − θ ) sin 2 λ h 2 | | a + b + 4 r θ sin 2 λ h | ≤ 1 ，可得 2 r ( 1 − 2 θ ) ≤ 2 a + b ；此时普遍性差分格式的稳定条件为：当 0 ≤ θ < 1 2 时， r ≤ a + b 2 ( 1 − 2 θ ) ；当 1 2 ≤ θ ≤ 1 时，无条件稳定。
所以，有如下定理：
定理2：当 1 2 ≤ θ ≤ 1 时，分数阶电报方程的普遍性差分格式为无条件稳定的；当 0 ≤ θ < 1 2 且 r ≤ a + b 2 ( 1 − 2 θ ) 时，分数阶电报方程的普遍性差分格式是稳定的。
定义精确解和数值解之间的误差为
e j k = u ( x j , t k ) − u j k ( k = 1 , 2 , ⋯ , N ; j = 1 , 2 , ⋯ , M − 1 )
和
e k = [ e 1 k , e 2 k , ⋯ , e M − 1 k ] T ， R k = [ R 1 k , R 2 k , ⋯ , R M − 1 k ] T .
则初始条件和边界条件为
e i 0 = 0 , i = 0 , 1 , ⋯ , M , e 0 k = e M k = 0 , k = 0 , 1 , ⋯ , N
下面使用Fourier方法分析普遍性差分格式的收敛性。在此，定义网格函数
e k ( x ) = { 0 , 0 ≤ x ≤ h 2 , e j k , x j − h 2 < x ≤ x j + h 2 , j = 1 , 2 , ⋯ , M − 1 , 0 , L − h 2 < x ≤ L ( k = 1 , 2 , ⋯ , N )
和
R k ( x ) = { 0 , 0 ≤ x ≤ h 2 , R j k , x j − h 2 < x ≤ x j + h 2 , j = 1 , 2 , ⋯ , M − 1 , 0 , L − h 2 < x ≤ L , ( k = 1 , 2 , ⋯ , N )
e k ( x ) , R k ( x ) 的Fourier展开形式为：
e k ( x ) = ∑ l = − ∞ ∞ η k ( l ) e i 2 π l x / L ，其中系数 η k ( l ) = 1 L ∫ 0 L e k ( x ) e − i 2 π l x / L d x .
R k ( x ) = ∑ l = − ∞ ∞ ξ k ( l ) e i 2 π l x / L ，其中系数 ξ k ( l ) = 1 L ∫ 0 L R k ( x ) e − i 2 π l x / L d x .
同样，关于离散2模，有
‖ e k ‖ 2 2 = ∑ j = 1 M − 1 h | e j k | 2 = ∑ l = − ∞ ∞ | η k ( l ) | 2 , (10)
‖ R k ‖ 2 2 = ∑ j = 1 M − 1 h | R j k | 2 = ∑ l = − ∞ ∞ | ξ k ( l ) | 2 . (11)
基于上面的分析，假设 e j k = η k e i σ j h ， R j k = ξ k e i σ j h ，其中 σ = 2 π m / L .
{ ( 2 a + b + 2 r θ ) e i 1 − r θ ( e i + 1 1 + e i − 1 1 ) = ( 2 a + b − 2 r ( 1 − θ ) ) e i 0 + ( 1 − θ ) r ( e i + 1 0 + e i − 1 0 ) + R i 1 , ( a + b + 2 r θ ) e i k + 1 − r θ ( e i + 1 k + 1 + e i − 1 k + 1 ) = [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 2 r ( 1 − θ ) ] e i k + r ( 1 − θ ) ( e i + 1 k + e i − 1 k ) + ∑ j = 1 k − 1 q j e i k − j + q k e i 0 + R i k + 1 . (12)
其中 | R i k | ≤ K ( τ 2 − α + h 2 ) .
引理1：设 η k ( k = 1 , 2 , ⋯ , N ) 是方程(10)的解，存在一个正数 c 2 ，使得
| η k | ≤ c 2 a k − 1 − 1 ξ 1 , k = 1 , 2 , ⋯ , N .
证明由定义，可知 e 0 = 0 ，有 η 0 ≡ η 0 ( m ) = 0.
另外，有
‖ R k ‖ 2 ≤ c 1 L ( τ 2 − α + h 2 ) , (13)
(11)式等式右边的系数数列是收敛的，因此存在一个正数 c 2 使得
| ξ k | ≡ | ξ k ( m ) | ≤ c 2 ξ 1 ≡ c 2 | ξ 1 ( m ) | , k = 1 , 2 , ⋯ , N . (14)
因此有
| ξ k | ≤ c 2 ξ 1 , k = 1 , 2 , ⋯ , N .
参考稳定性证明，使用数学归纳法可得
当 k = 1 时，有
( 2 a + b + 2 r θ ) e i 1 − r θ ( e i + 1 1 + e i − 1 1 ) = ( 2 a + b − 2 r ( 1 − θ ) ) e i 0 + ( 1 − θ ) r ( e i + 1 0 + e i − 1 0 ) + R i 1 ,
| η 1 | = 1 | 2 a + b + 4 r θ sin 2 λ h | ξ 1 ≤ q 0 − 1 ξ 1 .
当 k > 1 时，
( a + b + 2 r θ ) e i k + 1 − r θ ( e i + 1 k + 1 + e i − 1 k + 1 ) = [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 2 r ( 1 − θ ) ] e i k + r ( 1 − θ ) ( e i + 1 k + e i − 1 k ) + ∑ j = 1 k − 1 q j e i k − j + q k e i 0 + R i k + 1
( a + b + 4 r θ sin 2 λ h ) | η k + 1 | = [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | η k | + ∑ j = 1 k − 1 q j | η k − j | + | R k | ≤ { | [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | q k − 1 + ∑ j = 1 k − 1 | q j | q k − j − 1 + 1 } c 2 | ξ 1 | ≤ { | [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | + ∑ j = 1 k − 1 | q j | + q k } c 2 q k − 1 | ξ 1 | ≤ { | [ a ( 2 d 0 − d 1 ) + b ( c 0 − c 1 ) − 4 r ( 1 − θ ) sin 2 λ h ] | + b c 1 + a ( d 0 − d 1 ) } c 2 q k − 1 | ξ 1 | ≤ | [ a + b − 4 r ( 1 − θ ) sin 2 λ h ] | c 2 q k − 1 | ξ 1 |
即
| η k + 1 | ≤ | a + b − 4 r ( 1 − θ ) sin 2 λ h | a + b + 4 r θ sin 2 λ h c 2 q k − 1 | ξ 1 | ≤ c 2 q k − 1 | ξ 1 | ,
lim n → ∞ c n − 1 n α = lim n → ∞ n − α n 1 − α − ( n − 1 ) 1 − α = lim n → ∞ n − 1 1 − ( 1 − 1 n ) 1 − α = 1 1 − α ,
由 q k = a ( 2 d k − d k − 1 ) + c k = τ − 2 α Γ ( 3 − 2 α ) ( 2 d k − d k − 1 ) + τ − α Γ ( 2 − α ) c k ，可知 q k − 1 n α ≤ c k − 1 n α ， lim n → ∞ q k − 1 n α ≤ lim n → ∞ c k − 1 n α = 1 1 − α .
则存在正常数 C 1 ，使得
| η k + 1 | ≤ c 2 q k − 1 | ξ 1 | ≤ ( k τ ) α Γ ( 2 − α ) c k k α c 2 | ξ 1 | ≤ T α ( 1 − α ) Γ ( 2 − α ) c 2 | ξ 1 | = C 1 | ξ 1 | , k = 1 , 2 , ⋯ , N .
其中 C 1 = T α ( 1 − α ) Γ ( 2 − α ) c 2 .
同稳定性分析，由 | a + b − 4 r ( 1 − θ ) sin 2 λ h | a + b + 4 r θ sin 2 λ h ≤ 1 可知，普遍性差分格式的收敛条件为：
当 0 ≤ θ < 1 2 时， r ≤ a + b 2 ( 1 − 2 θ ) ；当 1 2 ≤ θ ≤ 1 时，无条件收敛。
综上，可有
定理3：当 1 2 ≤ θ ≤ 1 时，分数阶电报方程的普遍性差分格式是无条件收敛的，收敛阶是 O ( τ 2 − α + h 2 ) ；当 0 ≤ θ < 1 2 时，分数阶电报方程的普遍性差分格式收敛的条件是 r ≤ a + b 2 ( 1 − 2 θ ) .
基于Intel Core i5-2400 CPU@3.10 GHz，在MatlabR2014a环境下进行数值试验；我们将通过数值试验验证前面的理论分析。
例1考虑如下时间分数阶电报方程 [
∂ 1.6 u ( x , t ) ∂ t 1.6 + ∂ 0.8 u ( x , t ) ∂ t 0.8 = ∂ 2 u ( x , t ) ∂ x 2 + f ( x , t ) , 0 < x < 1 , t > 0
其中 f ( x , t ) = x ( 1 − x ) ( 2 Γ ( 1.4 ) t 0.4 + 2 Γ ( 2.2 ) t 1.2 ) + 2 ( 1 + t 2 ) ，
初边值条件： u ( x , 0 ) = x ( 1 − x ) , u t ( x , 0 ) = 0 , 0 < x < 1 ； u ( 0 , t ) = u 1 ( t ) = 0 , u ( L , t ) = u 2 ( t ) = 0.
则该方程有精确解： u ( x , t ) = ( 1 + t 2 ) x ( 1 − x ) .
取 M = 200 ， N = 100 时，其解析解曲面见图1；此时 r = 1.0 e + 04 ，令 r c = a + b 2 ( 1 − 2 θ ) ，当 θ = 0. 2 时， r c = 4.57537 e + 03 < r ，不满足普遍性差分格式稳定条件，其数值解呈指数增长趋势，其数值解曲面见图2；当 θ = 0. 4 时， r c = 1.37261 e + 04 > r ，满足普遍性差分格式稳定条件，其数值解曲面见图3~图4可以很好逼近解析解曲面、且光滑；当 θ ≥ 0. 5 时，普遍性差分格式无条件稳定，其解析解曲面与解析解曲面相近，且光滑，可以很好逼近解析解；与前面普遍性差分格式的稳定性分析一致。
当 θ ≥ 0. 5 时，普遍性差分格式无条件稳定。考察 θ = 0. 5 和1.0时普遍性差分格式的稳定性，分别取 M = 200 ， N = 50 和 M = 400 ， N = 200 两种情况，由图5~图8可知，隐式格式和C-N格式无条件稳定，并且普遍性差分格式离散解较好地逼近了解析解。
图1. 精确解曲面
取M = 200，N = 100，差分格式解与解析解的误差情况见表1。由表1可知，当θ = 0.5时(C-N格式)的精度比θ = 1 (隐式格式)的精度要好，且数值解均小于解析解；当θ = 0.4时数值解的绝对误差与θ = 0.5时基本相同，但此时数值解大于解析解的值。因此当θ取0.4至0.5中间某值得时候精度可取的最好。如图9所以，当θ = 0.45时的误差要小于其他三种情况。由表1和图9可以看出，当θ取值在0.5附近时精度较高，对于例1而言，当θ = 0.45时普遍性差分格式的精度最好。对方程进行离散时网格划分越密集，所得的差分格式解越平滑越接近于解析解，采用普遍性差分格式，尤其是θ = 0.5附近值时普遍性差分方法求解时间分数阶电报方程是可行的。
图2. 当θ = 0.2时普遍性差分格式解曲面
图3. 当θ = 0.4时普遍性差分格式解曲面
图4. 当θ = 0.5时普遍性差分格式解曲面
| x | 精确解 | θ = 1 | 误差值 | θ = 0.5 | 误差值 | θ = 0.4 | 误差值 |
|---|---|---|---|---|---|---|---|
| 0.1 | 0.1800 | 0.179790 | 2.094338E−04 | 0.179981 | 1.841547E−05 | 0.180019 | −1.987445E−05 |
| 0.2 | 0.3200 | 0.319629 | 3.05689E−04 | 0.319965 | 3.484993E−05 | 0.320032 | −3.244954E−05 |
| 0.3 | 0.4200 | 0.419515 | 4.47175E−04 | 0.419952 | 4.772436E−05 | 0.420039 | −3.988088E−05 |
| 0.4 | 0.4800 | 0.479447 | 5.28257E−04 | 0.479944 | 5.590459E−05 | 0.480043 | −4.371741E−05 |
| 0.5 | 0.5000 | 0.499424 | 5.54648E−04 | 0.499941 | 5.870740E−05 | 0.500044 | −4.489234E−05 |
| 0.6 | 0.4800 | 0.479447 | 5.28257E−04 | 0.479944 | 5.590459E−05 | 0.480043 | −4.371741E−05 |
| 0.7 | 0.4200 | 0.419515 | 4.47175E−04 | 0.419952 | 4.772436E−05 | 0.420039 | −3.988088E−05 |
| 0.8 | 0.3200 | 0.319629 | 3.05689E−04 | 0.319965 | 3.484993E−05 | 0.320032 | −3.244954E−05 |
| 0.9 | 0.1800 | 0.179790 | 2.94338E−04 | 0.179981 | 1.841547E−05 | 0.180019 | −1.987445E−05 |
表1. 普遍性差分格式解的误差分析(取M = 200, N = 100)
图5. N = 200，M = 50时C-N (θ = 0.5)格式解
图6. N = 400，M = 200时C-N (θ = 0.5)格式解
图7. N = 200，M = 50时隐式(θ = 1)格式解
图8. N = 400，M = 200时隐式(θ = 1)格式解
图9. 当θ取不同值时普遍性差分格式的误差曲面
本文对时间分数阶电报方程构造了一类普遍性差分方法，分析了该种差分方法解的存在唯一性，稳定性和收敛性，最后进行了数值实验；对比解析解与数值解，结果表明该类方法尤其是当θ取0.5附近时，对数值求解时间分数阶电报方程是有效的、可行的。
国家科技重大专项子课题(2017ZX07101001-01)，中央高校基金科研业务专项资金资助(2018MS168)。
吴立飞,杨晓忠. 分数阶电报方程一类有效的普遍性差分方法A Kind of Efficient Universal Differential Methods for Fractional Telegraph Equations[J]. 应用数学进展, 2019, 08(09): 1544-1555. https://doi.org/10.12677/AAM.2019.89181
https://doi.org/10.1063/1.2890375
https://doi.org/10.1016/j.cam.2013.02.009
https://doi.org/10.1016/j.apm.2013.11.035