劉 蕊，周永芳

（河北工業(yè)大學(xué) 理學(xué)院，天津 300401）

0 引言

分?jǐn)?shù)階微分方程邊值問題廣泛地應(yīng)用于數(shù)學(xué)、工程、化學(xué)等領(lǐng)域中。但是由于它的解析解很難求解，因此，其數(shù)值解研究尤為重要。學(xué)者們提出了各種數(shù)值方法求解分?jǐn)?shù)階微分方程的數(shù)值解，包括Adomian分解法[1]，有限差分法[2]，同倫分析法[3]，譜方法[4-6]等。近年來，再生核理論廣泛應(yīng)用于各種數(shù)學(xué)、物理和工程問題中。學(xué)者們提出了在再生核空間中求解不同分?jǐn)?shù)階微分方程的數(shù)值方法。Mohammed Al-Smadi等[7]求解了時間分?jǐn)?shù)階邊值問題，將再生核空間法和迭代法相結(jié)合，在再生核Hilbert空間中提供了分?jǐn)?shù)階邊值問題的近似解，并且用歸納法分析了這種方法的收斂性和穩(wěn)定性。Omar Abu Arqub[8-9]考慮了在ABC分?jǐn)?shù)階導(dǎo)數(shù)意義下的分?jǐn)?shù)階微分方程，文獻(xiàn)[8]利用再生核空間研究了分?jǐn)?shù)階積分微分方程的近似解，文獻(xiàn)[9]利用再生核空間法和算子不動點理論研究了Riccati分?jǐn)?shù)階微分方程和Bernoulli分?jǐn)?shù)階微分方程的近似解，計算速度快。Hossein Beyrami 等[10]在梯度網(wǎng)格上研究了帶有弱奇異核的分?jǐn)?shù)階積分微分方程的近似解，在再生核空間中給出了誤差估計，并討論了近似解的穩(wěn)定性。Esmail Babolian等[11]在再生核空間中考慮了利用再生核空間法求解線性和非線性微分方程數(shù)值解產(chǎn)生的誤差。ZhongChen等[12]在分?jǐn)?shù)階導(dǎo)數(shù)意義下的加權(quán)再生核空間中考慮帶有弱奇異核的分?jǐn)?shù)階微分積分方程近似解，收斂速度快，同時計算精度高。文獻(xiàn)[13-14]中考慮了下列分?jǐn)?shù)階微分方程邊值問題的近似解：

式中：Dα為Caputo分?jǐn)?shù)階導(dǎo)數(shù)，1＜α≤2；a0,a1,b0,b1為常數(shù)；b(t),c(t)為充分光滑函數(shù)；f(t)為給定的函數(shù)，且

作者利用打靶法將邊值問題轉(zhuǎn)化為初值問題，并將方程轉(zhuǎn)化為帶有弱奇異核的分?jǐn)?shù)階積分方程求解方程的近似解。本文將在再生核空間中討論這類分?jǐn)?shù)階邊值問題的近似解。

1 預(yù)備知識

定義1（Riemann-Liouville分?jǐn)?shù)階積分）對于區(qū)間[0，1]上的連續(xù)函數(shù)y(t)，有

式中，Γ(·)為Gamma函數(shù)。

引理1 若K(s,t)∈C(D)，那么存在且在[0，1]區(qū)間上連續(xù)，并且有

具體的證明過程參見文獻(xiàn)[15]中的定理2.3。

定義2（Caputo分?jǐn)?shù)階導(dǎo)數(shù)）對于區(qū)間[0，1]上的連續(xù)函數(shù)y(t)，有

式中，m=[α]+1，[α]為小于α的最大整數(shù)。由式（1）中1＜α≤2 知，m=2。

定義3（再生核空間）設(shè)H為Hilbert 空間，H中的元素是集合X上的復(fù)合函數(shù)，對于任意的y(t),z(t)∈H，定義內(nèi)積為

若存在一個關(guān)于t的函數(shù)Rs(t)，使得對于任意的函數(shù)y(t)∈H，有

則稱Hilbert空間H為再生核空間，其中函數(shù)Rs(t)為再生核空間H的再生核函數(shù)。

定義4 Hilbert空間W1[0,1]={y(t)|y(t)為絕對連續(xù)函數(shù)，y′(t)∈L2[0,1]}。

在W1[0,1]空間中，定義內(nèi)積為

式中，y(t),z(t)∈W1[0,1]。范數(shù)為

在文獻(xiàn)[16]中，得到W1[0,1]空間的再生核函數(shù)rs(t)的表達(dá)式為

定義5 Hilbert空間W2[0,1]={y(t)|y(i)(t)(i=0,1,2)為絕對連續(xù)函數(shù)，

在W2[0,1]空間中，定義內(nèi)積為

式中，y(t),z(t)∈W2[0,1]。范數(shù)為

Hilbert 空間W2[0,1]為再生核空間。設(shè)Rs(t)為W2[0,1]的再生核函數(shù)，那么根據(jù)W2[0,1]空間的內(nèi)積定義，可以得到如下等式：

利用分部積分，可以得到

由再生核函數(shù)Rs(t)∈W2[0,1]，Rs(t)滿足W2[0,1]空間的邊值條件

由y(t)∈W2[0,1]，y(t)同樣滿足W2[0,1]邊值條件，那么

于是有再生核空間W2[0,1]的再生核函數(shù)Rs(t)所滿足的條件：

再由函數(shù)的連續(xù)性以及跳躍條件，可以得到

W2[0,1]的再生核函數(shù)Rs(t)的一般形式為

式中，ci,di,i=0,1,…,4,5 由式（2）～式（4）共12個方程唯一確定。

下面引入算子L:W2[0,1]→W1[0,1]，由定義2

將方程（1）的邊值條件齊次化，則方程（1）可以表示為如下算子方程：

式中，y(t)∈W2[0,1]，當(dāng)y(t)∈W2[0,1]時，f(t)∈W1[0,1]。記L*為L的共軛算子。

定理1 (LRs)(i)(t)∈W2[0,1],i=0,1,2。

式中，M1,M2,M3為常數(shù)。那么就有下列不等式成立：

接下來證明(LRs)(t)滿足W2[0,1]的邊值條件。

綜合可得，(LRs)(t)∈W2[0,1]。

由文獻(xiàn)[16]中定理6.3知，若Rs(t)為W2[0,1]的再生核函數(shù)，那么無論是關(guān)于變量s還是關(guān)于變量t，總有函數(shù)

證明成立。

引理2 算子L:W2[0,1]→W1[0,1]為有界線性算子。

證明 顯然L為線性算子。下證算子的有界性。

由再生核的性質(zhì)，可以得到(Ly)(i)(t),i=0,1可以表示為

那么由Schwarz不等式可以得到

那么就有

證明成立。

那么?i=1,2,…，有

證明成立。

式中，γij為正交化系數(shù)。

2 方程的精確解和近似解

定理2 方程（5）的精確解y(t)可以表示為

證明 根據(jù)引理5，y(t)∈W2[0,1]，有

通過對方程級數(shù)形式的精確解（6）截斷，得到方程（5）的近似解為

證明成立。

3 收斂性分析和誤差估計

證明 根據(jù)引理5，y(t)∈W2[0,1]，當(dāng)i=0,1,2 時，有

證明成立。

式中：h=max|tj+1-tj|,j=1,2,…,n-1；a為確定的實數(shù)。

證明 取t∈[tj,tj+1]，其中j=1,2,…,n-1

那么

根據(jù)上式可以得出結(jié)論

同理可得，y(t)∈W2[0,1]，那么存在a2為常數(shù)，使得

由定理3可以得到，對于任意ε＞0，當(dāng)時，有

證明成立。

4 數(shù)值算例

在這部分內(nèi)容中，所有數(shù)值結(jié)果由Mathematic 11.0 計算所得。應(yīng)用本文提出的算法，在[0,1]上取=1,2,…,n。由式（7）知道，方程的近似解為

對n取不同的值，求解下列分?jǐn)?shù)階微分方程邊值問題的近似解。數(shù)值結(jié)果說明了算法的有效性。

例1 考慮下面分?jǐn)?shù)階微分方程邊值問題的數(shù)值解

式中：t∈[0,1],α=1.1；f(t) 是使得方程的解為y(t)=-0.3t3+0.1t+1 的函數(shù)。數(shù)值結(jié)果見表1，其中en=|yn(t)-y(t) |，分別取n=10,20。

表1 算例1 數(shù)值結(jié)果Tab.1 Numerical result for Example 1

例2 考慮下面分?jǐn)?shù)階微分方程邊值問題：

式中：t∈[0,1],α=1.2 ；f(t) 是使得方程的解為y(t)=-0.65t3+0.4t+3 的函數(shù)。數(shù)值結(jié)果見表2，其中en=|yn(t)-y(t) |，分別取n=10,25。

表2 算例2 數(shù)值結(jié)果Tab.2 Numerical result for Example 2

5 結(jié)論

本文主要介紹了在再生核空間中求解分?jǐn)?shù)階微分方程邊值問題近似解的數(shù)值方法，利用再生核空間的良好性質(zhì)獲得了這類方程級數(shù)形式的精確解，通過截斷方程級數(shù)形式的精確解獲得了方程的近似解，并在再生核空間中證明了所提方法的收斂性，給出了誤差估計。這種方法同樣可以應(yīng)用于其他線性分?jǐn)?shù)階微分方程數(shù)值解的求解。今后將進(jìn)一步在再生核空間中研究非線性分?jǐn)?shù)階微分方程數(shù)值解法，同時提高再生核空間數(shù)值解法的精確度。