一類非線性R-L分數(shù)階積分微分方程的數(shù)值解法

張盼盼,任正杰

(1 寧夏大學數(shù)學計算機學院,寧夏 銀川，750021；2 甘肅省永昌縣第一高級中學)

?

一類非線性R-L分數(shù)階積分微分方程的數(shù)值解法

張盼盼1,2,任正杰2

(1 寧夏大學數(shù)學計算機學院,寧夏 銀川，750021；2 甘肅省永昌縣第一高級中學)

利用Adomian多項式將分數(shù)階積分微分方程中的積分項離散化,進而得到原方程解的級數(shù)表達形式,數(shù)值算例驗證了該分解方法的有效性。

Adomian多項式;分數(shù)階積分微分方程;數(shù)值解

分數(shù)階微積分是一古老而新鮮的概念。早在17世紀末,就有像L’Hospital,Leibniz等數(shù)學家開始考慮分數(shù)階微積分,但由于缺乏應(yīng)用等眾多原因的支撐,使得分數(shù)階微積分的發(fā)展相對滯后。近年來,越來越多的分數(shù)階微積分方程廣泛出現(xiàn)在各大工程領(lǐng)域。例如,PID控制理論、粘彈性材料及混沌現(xiàn)象等反常問題[1,2],其中有些問題建模后得到的方程大多數(shù)都是非線性分數(shù)階積分微分方程。由于分數(shù)階微積分復(fù)雜的定義,增加了分數(shù)階方程的求解難度,分數(shù)階積分微分方程的求解對眾多學者來說更是挑戰(zhàn)。目前關(guān)于分數(shù)階積分微分方程數(shù)值解的文獻已有一些,主要是用小波的方法求解 Fredholm或Volterra型的分數(shù)階積分微分方程[3～6]。LI Huang等[7]采用Taylor級數(shù)的方式求解了線性分數(shù)階積分微分方程。然而,關(guān)于非線性分數(shù)階積分微分方程數(shù)值解的文獻甚少。

利用Adomian分解法求解代數(shù)方程具有強大的優(yōu)越性，其解為級數(shù)形式、收斂快、計算方便且容易在計算機上運算。它的核心思想是：把方程拆成n個部分，把解拆成n個項，同時非線性項用一種特殊的多項式An進行替代，然后由低階解分量逐漸向高階解分量逐一解出，從而得到方程的近似解析解，也可以得到精確解。該方法無需進行任何變換，很好的保持了原方程的物理性質(zhì)。在實際生活中不僅會遇到線性問題和確定問題,而且經(jīng)常會遇到非線性問題甚至是隨機問題。對于這些非線性問題,通常的辦法是將其線性化近似處理或采用擾動技術(shù)。然而,這些方法在求解時無形中改變了原問題,得到的解往往不能滿足實際需要。Adomian分解法不會存在這樣的問題,計算量很小，而且很快就能收斂到真解。所以,Adomian分解法是求解非線性方程的有力工具[8,9],筆者利用Adomian多項式將R-L分數(shù)階積分微分方程離散化,得到原方程解的級數(shù)表達式,并使用數(shù)值算例證明該算法的有效性。

1 預(yù)備知識

定義[1]設(shè)f(x)∈L[a,b],α>0。則稱

(1)

為Riemann-Liouville分數(shù)階積分,其中t∈[a,b], Γ(α)為Gamma函數(shù)。

2 主要結(jié)果

設(shè)

(2)

(3)

由Adomian分解法知,設(shè)

(4)

其中

將(4)式代入(3)式,整理得

(5)

(6)

(7)

3 解的收斂性

唯一性定理 方程(2)有唯一解?0<γ<1,其中γ=λML/Γ(α+2)。

證明:設(shè)

則有

由方程(2)的條件得

=γ‖x-y‖

易見,當γ∈(0,1)時,H即為壓縮映射。由引理知,H中存在唯一的不動點,即方程(2)有唯一解。

收斂性定理 方程(2)的解收斂?g(t)有界且0<γ<1,其中γ=λML/Γ(α+2)。

=γ‖Sp-1-Sq-1‖

所以

‖Sp-Sq‖≤γ‖Sp-1-Sq-1‖≤γ2‖Sp-2-Sq-2‖≤…≤γq‖Sp-q-S0‖

令p=q+1,則有

‖Sp-Sq‖≤γq‖S1-S0‖

‖Sp-Sq‖ =‖Sq+1-Sq+Sq+2-Sq+1+…+Sp-Sp-1‖

≤‖Sq+1-Sq‖+‖Sq+2-Sq+1‖+…+‖Sp-Sp-1‖

≤[γq+γq+1+…+γn-1]‖S1-S0‖

因為g(t)有界,由(5),(6)式知,x1(t)有界,因此當γ<1,q→∞時,‖Sp-Sq‖→0,即函數(shù)列{Sp}是Banach空間B上的柯西列。

4 數(shù)值算例

例1 求解分數(shù)階積分微分方程

其初始條件為x(0)=1。這類方程的解析解一般不易求得,但由文獻知,當α=1時,此方程的解析解為x(t)=et。因此,以α=1為例檢驗筆者算法的有效性。

圖1 非線性R-L分數(shù)階積分微分方程的數(shù)值解與精確解的比較

表1 非線性R-L分數(shù)階積分微分方程的數(shù)值結(jié)果及誤差(α=1)

例2 求解分數(shù)階積分微分方程

圖2 分數(shù)階積分微分方程的數(shù)值解與精確解的比較

表2 α=1時分數(shù)階積分微分方程的當數(shù)值解與精確解的比較

[1] 陳文,孫洪廣,李西成.力學與工程問題的分數(shù)階導(dǎo)數(shù)建模[M].北京：科學出版社,2010.

[2] 張亞鵬,高峰.分數(shù)階粘彈性積分本構(gòu)模型[J].濟南大學學報:自然科學版,2012,26(1):102-105.

[3] 尹建華,任建婭,儀明旭.Legendre小波求解非線性分數(shù)階Fredholm積分微分方程[J].遼寧工程技術(shù)大學學報:自然科學版,2012,31(3):405-408.

[4] Xiaohua Ma,Chengming Huang.Numerical solution of fractional integro-differential equations by a hybrid collocation method[J].Applied Mathematics and Computation,2013,219(12):6 750-6 760.

[5] Mustafa Gulsu.Numerical approach for solving fractional Fredholm integro-differential Equation[J].International Journal of Computer Mathematics,2013,29(16):1-22.

[6] Zhu Li,Fan Qibin.Numerical solution of nonlinear fractional-order Volterra integro-differential equations by SCW[J].Commun nonlinear Sci Numer Simulat,2013,18(5):1 203-1 213.

[7] LI Huang,Xianfang Li.Approximate solution of Fracional integro-differential equations by Taylor expansion method[J].Computer Math Appl,2011,24(62):1 127-1 134.

[8] G Adomian.Random Volterra integra equations[J].Math Comput Model,1995,22(8):101-102.

[9] G Adomian,R Rach.Modified Adomian polynomials[J].Math Comput Model,1996,24(11):39-46.

[10] 張恭慶,林源渠.泛函分析講義[M].北京:北京大學出版社,2011.

(責任編輯：朱寶昌)

Numerical Solution of a Class of Nonlinear R-L Fractional Integro-differential Equation

ZHANG Pan-pan1，2,REN Zheng-jie2

(1 Department of Mathematics and Computer Science,Ningxia University, Yinchuan, 750021;2 The First High School of Yongchang, Yongchang Gansu;China)

In this paper, Adomian polynomial is used to discrete fractional integro-differential equation to obtain solution of the original equation, numerical examples show that this method is effective to approximate the numerical result.

Adomian polynomials;fractional integro-differential equations;numerical solution

10.3969/J.ISSN.1672-7983.2015.02.010

2015-04-02; 修改稿收到日期: 2015-05-24

O241.8

A

1672-7983(2015)02-0047-05

張盼盼(1986-),男,碩士研究生。主要研究方向：復(fù)分析在力學中的應(yīng)用。