王林君, 吳 燕, 劉夢雪, 孟義平, 曹明鈺

(1. 江蘇大學(xué) 理學(xué)院, 江蘇 鎮(zhèn)江 212013； 2. 江蘇科技大學(xué) 理學(xué)院, 江蘇 鎮(zhèn)江 212003；3. 吉林大學(xué) 數(shù)學(xué)學(xué)院, 長春 130012)

分?jǐn)?shù)階微分方程廣泛應(yīng)用于自然界中各種現(xiàn)象的建模和分析中[1-2]. 不同于整數(shù)階微分算子, 分?jǐn)?shù)階微分算子是非局部的[3], 因此, 分?jǐn)?shù)階微分方程能更好地模擬自然界的物理現(xiàn)象和動態(tài)過程. 例如, 溶質(zhì)在天然多孔或斷裂介質(zhì)中的平流和分散[4], 在外部影響下模擬黏彈性材料的行為[5], 在數(shù)學(xué)心理學(xué)領(lǐng)域使用分?jǐn)?shù)階系統(tǒng)模擬人類的行為[6]等. 分?jǐn)?shù)階微分方程終值問題目前已有很多有效的數(shù)值方法： Ford等[7]通過將分?jǐn)?shù)階微分方程邊值問題轉(zhuǎn)化為Fredholm積分方程, 提出了非多項式配置法； Pedas等[8]利用樣條配置法求解了一類分?jǐn)?shù)階微分方程邊值問題. 本文針對一類分?jǐn)?shù)階微分方程終值問題, 給出一種混合配置算法.

考慮如下分?jǐn)?shù)階微分方程終值問題：

方程(3)-(4)等價于具弱奇核的Volterra積分方程[7]

(5)

在此基礎(chǔ)上, 文獻(xiàn)[9]結(jié)合不動點技巧給出了求解方程(3)-(4)的混合配置格式. 當(dāng)a>0時, 方程(1)-(2)是終值問題, Ford等[7]證明了其可轉(zhuǎn)化為等價的Fredholm積分方程:

(6)

并給出了其解的存在性定理與唯一性定理.

本文假設(shè)方程(1)-(2)滿足文獻(xiàn)[7,9]的條件, 給出方程(1)-(2)的混合配置解法.

首先, 基于打靶法將方程(1)-(2)轉(zhuǎn)化為初值問題(3)-(4), 步驟如下：

1) 假設(shè)初值條件(4)中的初值y0已知, 則用求解分?jǐn)?shù)階初值問題的混合配置法[9]進(jìn)行計算便可得到一個含有y0的解y(t,y0)；

2) 根據(jù)初值問題的解必經(jīng)過點(a,ya)從而有y(a,y0)=ya, 故可選用適當(dāng)?shù)牡ㄇ蠼獾玫統(tǒng)0；

3) 根據(jù)得到的y0, 將方程(1)-(2)轉(zhuǎn)化為初值問題(3)-(4).

其次, 在求解初值問題(3)-(4)過程中, 先將方程(3)-(4)轉(zhuǎn)化為方程(5), 然后采用文獻(xiàn)[9]中的混合配置法. 令N0表示非負(fù)整數(shù)集, 對于0<α<1和正整數(shù)m, 定義C[0,1]上有限維非多項式子空間[10]

Wα,m={i+jα:i,j∈N0,i+jα

t0j=t0+h0cj,j=1,2,…,l,

其中cj是配置參數(shù), 且滿足0≤c1

tik=ti+hidk,k=1,2,…,m,

其中0≤d1

由于方程(5)的解在原點附近具有奇性, 因此混合配置法就是在區(qū)間σ0上使用非多項式基底, 而在其他區(qū)間使用多項式基底. 定義一個有限維空間

令

z(t)∶=f(t,y(t)),

(7)

將式(7)代入方程(5)得

(8)

將式(8)代入式(7)可得z滿足非線性方程

(9)

最后把結(jié)果代入式(8)即得方程(5)的解.

注1[9]方程(10)-(11)的性質(zhì)由f(t,y(t))確定, 當(dāng)f(t,y(t))關(guān)于t和y(t)是線性時, 則得到的方程也是線性的.

例1考慮如下線性分?jǐn)?shù)階微分方程終值問題[7]：

其真解為

基于打靶法, 首先將其轉(zhuǎn)化為初值問題：

其等價于Volterra積分方程：

表1 例1的數(shù)值結(jié)果

例2考慮如下非線性分?jǐn)?shù)階微分方程終值問題:

其真解是y(t)=1-t2.

同理, 首先將其轉(zhuǎn)化為初值問題

(12)

令

(13)

則有

(14)

將式(14)代入式(13)可得

(15)

使用混合配置法求解方程(15), 將得到的解代入式(14), 所得終值問題的數(shù)值結(jié)果列于表2.

表2 例2的數(shù)值結(jié)果

表1和表2列出了當(dāng)m=2, 選取不同的N值時得到的初值條件y0值、誤差及收斂階. 由表1和表2可見, 用混合配置方法求解分?jǐn)?shù)階微分方程終值問題更有效.