非线性Volterra积分方程
含causal算子分数阶非线性微分方程的拟线性方法

含causal算子分数阶非线性微分方程的拟线性方法2012正第32卷第1期河北大学(自然科学版)JournalofHebeiUniversity(NaturalScienceEdition)2O12V o1.32NO.1含causal算子分数阶非线性微分方程的拟线性方法王培光,李志芳(1.河北大学电子信息工程学院,河北保定071002;2.河北大学数学与计算机学院,河北保定071002)摘要:采用拟线性化方法讨论了含causal算子的分数阶非线性微分方程初值问题,通过构造2个单调迭代序列,证明了它们一致且平方收敛于给出问题的解.关键词:拟线性方法;causal算子;分数阶微分方程;平方收敛中图分类号:O175.1文献标志码:A文章编号:1000—1565(2012)01—0001—06 Quasilinearizati0nforsolutionofnonlinearcausalfractionaldifferentialequationsWANGPei—guang.LIZhi-fang.(1.CollegeofElectronicandInformationEngineering,HebeiUniversity,Baoding071002,C hina;2.CollegeofMathematicsandComputerScience,HebeiUniversity,Baoding071002,China )Abstract:Byusingthequasilinearizationmethodforcausalfractionaldifferentialequations, theau—thorsconstructtWOmonotonesequences,thenprovethattheybothconvergeuniformlyandq uadraticallytothesolutionofthegivenproblem.Keywords:quasilinearizationmethod;causaloperator;fractionaldifferentialequations;qu adraticcon—vergenceMSC2010:34A34在非线性微分方程解的定性问题的研究中,拟线性化方法得到了广泛的使用口].由于含causal算子微分方程系统模型可描述现实世界的许多问题,因而引起了人们的广泛关注.文献[2]利用上下解结合单调迭代方法,给出了一类一致收敛于含causal算子微分方程解的迭代序列.近年来分数阶微分方程引起了人们的广泛关注[3].然而关于用拟线性化方法研究含causal算子分数阶微分方程初值解的结果并未见到.本文将利用拟线性化方法对含causal算子的分数阶非线性微分系统两项和的初值问题(简称IVP)D(£)一(Qu)()+(Pu)(≠),(O)一0(1)进行研究,得到解的一致且平方收敛的结果.这里Q,P:E—E是连续causal算子,Dq 是Caputo分数阶导数,0<q<1,E—C(J×R,R)和t∈J—Eo,T].收稿日期:2011一O9~21基金项目:国家自然科学基金资助项目(10971045);河北省自然科学基金资助项目(A2009000151)第一作者:王培光(1963--),男,黑龙江哈尔滨人,河北大学教授,博士生导师,主要从事微分方程与控制理论方面的研究E-mail:pgwang@?2?河北大学(自然科学版)第32卷1预备知识利用下面定义和引理证明主要定理.定义1如果对于E—cEJ×R,R]中每对元素(z,),使z(s)一(s),有()(£)一()(£),其中0≤S≤£,t<T,T是任意正实数,则称Q:E—E是Causal算子.式(1)等价的V olterra分数阶积分方程咖为(f)一乱.+J.(一s)(()(+(P)(),其中I1是Gamma函数.定义2考虑初值问题式(1),若a,∈cq[J,R],满足Da(£)≤(()(£)+(j9)(£),a(O)≤o,D(£)≥()(£)+(Pa)(),卢(0)≥o,则称a,I9分别是式(1)的耦合下解和耦和上解.引理1E.若,训∈[.厂,R]分别是式(1)的耦合下解和耦合上解,且(f)≤叫(),t∈J和Q,P ∈[Q,R],其中Q一[(£,z):()≤z≤叫(£),t∈刀,则存在式(1)的唯一解z()满足()≤z(£)≤(t),t∈J.引理2E设,叫∈[.厂,R],Q∈cEJ×R.,R],若Dqv(t)≤Q(t,硼,叫),D(£)≥Q(t,,口),Q(t,zl,Y1)一Q(t,z2,Y2)≥一L[(l—2)+(1一Y2)],L≥0,其中Q是causal算子,zl≥.Tg2,Y1≥Y2,且(O)≤叫(O),则有()≤(),t∈J.证明设一+£E(3),一一£E(3),其中E口(3Lt)表示的0<q<1阶导数,£>0是任意小的实数,则有叫.>叫,<和叫.(O)>硼(0)>(O).由已知条件可得D0(£)=D(f)一3Eq(3Lt)≤Q(£,硼,∞)一3Eg(3Lt)≤Q(£,,硼.)+E.(3)一3kE(3.)≤Q(£,0,0)+2LsE(3Lt)一3kE(3Lt)<Q(f,0,硼o).同理可得D.(£)一D.叫(£)+3E(3)≥Q(,,)+3工£E(3Lt)≥Q(f,,)一LeE口(3)+3LsE(3Lt)≥Q(t,0,)一2LeE(3Lt)+3E(3Lt.)>Q(£,,0).下面证明(t)<硼.(£),t∈J.假设不然,则存在to∈(O,明使0(.)一Wo(),7-)0(£)<.(£),0≤t<to成立,由此可得D.(.)≥D训.(£.),则进一步有Q(£,叫o(o),硼0())>D0()≥D硼0(£o)>Q(,(),V o()),此不等式与Q(,训o(£o),o(£o))一Q(£,V o(),0(o))矛盾,其中()一砌.(.).因此(£)<叫.(£),t∈J成立.在()一eE(3Lt)一0(£)<叫0(£)一硼(£)+eE(3Lt)中令e一0,则有(£)≤硼(),t∈.,.证毕.引理3对于Caputo线性分数阶微分方程D":,lu+(Qu)(t),"(0)一0,第1期王培光等:含Causal算子分数阶非线性微分方程的拟线性方法?3?其中Q∈CqEJ,R],且对q为H~lder连续,其唯一解(£)一"oEq()+I(£一s)q-Eq,((一s)q)(Qu)(s)ds,t∈J,其中Eq(∑k=l,㈤一分别是含1个参数和2个参数的Mittag—Leffler方程.2主要结果定理1假设下列条件成立:(A)口.,∈Cq[J,R],a.(£)≤(£),t∈J,满足Dao(£)≤(邸o)(£)+(邱o)(),D‰(£)≥(Qo)(£)+(Pao)();(A2)Q,P∈C[Q,R],Frechet导数Qf,P,Q和P存在,连续且满足()(£)≥0,(P)(£)≤0,(£,)∈Q.(A3)()()≤0和(P)(£)≤0,(£,)∈n,则存在2个单调序列{a}和{)一致且平方收敛于式(1)的唯一解.证明由(QU)(t)≥0,(PU)(£)≤0,对U≥有下列不等式成立(Qu)()≤(Q)()+(Q)("一)(£),(2)(P")()≤()(£)+(P)("一)(f).(3)对于任意的口.(£)≤".≤U≤(£),t∈J;Q,P满足L(u1一U2)≥(1)(£)一(Qu2)()≥一L(u1一z),L>0,(4)L(u1一U2)≥(P1)(£)一(P2)(£)≥一L(u1一U2).(5)考虑下面初值问题fDqu(£)一F(,ao,;)一{(qo)(£)+(Qo)(一)(f)+(P卢o)()+(Pao)(一)(£),(6)I【"(0)一U0;rD(£)一G(t,a0,;)一{(.)(£)+(Q.)(一口.)()+(Pa0)(£)+(Pa.)("一ao)(£),(7)I【v(0)一"o,其中口.(O)≤M.≤(0).由不等式(2),(3)和(A)可知D口o(£)≤(0)(£)+(邱o)()三F(t,a0,;);Do()≥(Qo)(£)+(Pao)()≥(邸o)()+(Q.po)(ao一)(f)+(邱o)(£)+(Pao)(口.一)(£)兰F(t,a,;口o);D.ao(£)≤(q9o)(£)+(o)(£)≤(o)(£)+(Qo)(J岛一口o)(£)+(o)()+(Qao)(一口o)(£)三G(t,Olo,;)Do(£)≥(.)(£)+(Pao)(£)三G(t,口o,;a0).由于(Q)()≤0和(P)()≤0,可得F(t,口.,;)和G(t,a.,;")分别对于和是非增的.由引理1可知式(6)和(7)存在唯一解(n,),满足ao≤a,≤,t∈J.?4?河北大学(自然科学版)第32卷即Da1(£)一F(,口o,;);DI91(£)一G(,口o,;口1).由不等式(2)和(3)可得D(a)(£)一(邸o)(£)+(Qo)(J91一)(z)+(邱o)(£)+(PO/o)(p?一po)≤(邸1)()+(Qo)(一)()+(Qo)(o)(£)+(邱)(£)+(P卢)(一卢)(£)+(Pao)(p一)(£)一(.)(£)+(邱)(£)+[(P)()一(Pao)()]('一)(£)≤(Q81)()+(邱1)();D1()=(Q.)()+(Q.)(a1一a.)()+(.)()+(Pa0)(口1一O/o)(£)≥(1)(£)+(Qa1)(a0一a1)(£)+(Qo)(a1一ao)(£)十(1)()+(P口o)(口.一a1)(£)+(Pao)(口1一ao)()一()(£)+()(£)+[(Qj9.)()一(Qa1)(£)](a.一a.)(t)≥(Q1)(£)+(1)().因为(Q")(£)关于"是非减的,(P)(£)关于乱是非增的,因此应用引理2,可得a(£)≤(£),t∈J即d.(£)≤a()≤()≤(£),t∈J.(8)考虑下面一组初值问题D(£)一F(£,a1,;),(0)一0,(9)D(£)一G(£,口l,;),口(O)一o,(10)可推出下列不等式D1(£)≤(1)()+(邱1)(£)三F(t,a1,;);D()≥()()+()()≥(邸)()+(Q1)(al一)(£)+(邱1)(£)+(P口1)(口1一)(£)三F(t,口1,;口1);D()≤()(£)+(邱)(£)≤(邸)()+(Q-)(一口t)(£)+(1)(£)+(Pa)(一a1)()三G(t,口1,;);D1(£)≥(1)(£)+(1)(£)兰G(t,口1,;口1).由引理1可知式(9)和式(10)存在唯一解az,使a.≤a.,』92≤,t∈J成立.同样由于D.a2()≤(q}92)(£)+(邱2)();D2()≥(2)(£)+(2)().则应用引理2可得口.()≤(£),t∈J.综上可知O/0≤a≤口2≤≤≤,如此继续下去,可得ao≤O/1≤z≤…≤a≤≤…≤≤≤,(11)其中单调序列{a(£)),{(£))是下列线性方程初值问题Da井1()一F(t,口,;1),O~rH-1(O)=Uo,(12)D井1()=G(t,a,;口计1),+l(0)=o(13)的解.综上所述很容易得知{a(f)),{(£))序列一致收敛于式(1)的唯一解.下面证明收敛速度是2次的.为了证明收敛速度是2次,设P(£):()一a(£),(£)一()一(£),其中"()是式(1)的唯一解.利用a,的定义,中值定理以及(A),有第1期王培光等:含causal算子分数阶非线性微分方程的拟线性方法?5?DP()=D()一Da(£)一(Q)()+(Pu)()一r-(邸,广)()+(Q,r)(一)(f)+(,r)()+(P口)(一)(£)]一一(Q)q一1()一(Pa)q~-l(£)一(Q1)[q一q1]()一(Pa1)Eq一q,r1](£)≤[(Q口1)(£)一(Q甜)(£)]g,广(£)+[(尸d,r)()一(P)()]q,r()+响(£)=(Q1)q1()一(P1)[q,广1+1]g,r1()+Mq(£)≤(N+3N)ql(£)+(£)+Mq(t),即cDqp.(£)≤Mq()+(N1+导N2)gl()+N2z(£),(14)其中<,<,a<<g-,U<<,I(Q)(£)I≤M,I(P)(f)I≤M2,l(Q)()I≤N,l(PU)(£)I≤N2和M—M1+M2.同理,Dq(£)一D(t)一DU()=(Q1)(£)+(Q,r1)(口一a1)(f)+(1)(f)+(Pa,r1)(口一a1)()一(Q)(£)一(Pu)()一(Q)(口,r1一)()+(P)(口1一)()+(Q1)(口一a1)(£)+(Pa1)(口一Olin-1)(£)一一(Q)1(£)一(Pa)P~-i()+(Q,r1)P~-I(£)一(Q1)(£)+(Pa,r1)P~-I(£)一(P口,r1)(£)一[一(Q)()+(Q)(£)](£)+[一(P)(£)+(P口1)()]p1()一(Q1)p()一(P口1)()≤[(Q口1)(£)一(Qa)()]乡(£)+E(P口)()一(P)(£)]p,广()一[(Q.)()+(P口,r)()]户()一(Q1)(1一口1)p,r1(£)一(P)户l(£)一[(Q)(£)+(Pa)()](),其中a,广1<,<,a,r1<1<l,a1<口1<.但(Q1)(一1一口,r1)(£)一(Q)[q,rl+夕]p(£)≤N1E-}p~1()+专q,r2(£)].因此cDqq.(£)≤Mp()+(N.+3N)户()+q,r2(£),(15)其中M,N,N就是前面定义的常数.可以把不等式(14)和不等式(15)写成矢量的形式Dr≤Arn+其一(),A一(;7o111/1,B一,gn,,1N2N1+3NzN2+3N1N,应用引理3,把Br看作强迫项给出≤(£)≤r(£一s)r.r.(A(t0E—s)).Br1(s)ds≤≤r(£)≤l(£一s)?.(A(一s))?Br1(s)≤J0 .6.河北大学(自然科学版)第32卷令L一~----BE,(AT1"BE(AT)可以得到令一口,)可以得到定理证毕.参考文献:T.qBma,xIrl()IEq,(ATa),£∈.,.口JmaxIr(£)I≤LmaxIrn--1(£)I.JI-1-1LAKSHMIKANTHAMV,V ATSALAAS.Generalizedquasmnearizationfornonline arproblems[M].Dordrecht:Klu—werAcademicPublishers,1998.1-23LAKSHMIKANTHAMV,LEELAS,DRICIZahiaetc.Theoryofcausaldifferentialeq uations[M].WorldScientific:AtlantisPress,2009.[3]GENGFengjie.Differentialequationsinvolvingcausaloperatorswithnonlinearperiodic boundaryconditions[J].Mathe—maticalandComputerModelling,2008,48,859—866.[4]JANKOWSKIT.Boundaryvalueproblemswithcausaloperators[J].NonlinearAnalysis ,2008,68:3625—3632.Es]DIETHELMK,FORDNJ.Analysisoffractionaldifferentialequations[J].AnalAppl,20 02,265(2):229—248.[63CAPUTOM.LinearmodelsofdissipationwhoseQisalmostindependentⅡ[J].GeophysJRAstron,1967,13(5):529—539.[7]V ASNUDHARADEVIJ,MCREAFA,DRICIZ.Generaizedquasilinearizationforfracti onaldifferentialequations[J].CompMathAppl,2010,59(3):1057—1062.[83V ASNUDHARADEVIJ,SUSEELACH.Quasilinearizationforfractionaldifferentiale quations[J].CommunicationsinAppliedAnalysis,2008,12(4):407—418.[9]V ASNUDHARADEVIJ.Generalizedmonotonemethodforperiodicboundaryvaluepr oblemsofCaputofractionaldiffer—entialequations[J].CommunicationsinAppliedAnalysis,2008,12(4),399—406.[1O]王培光,高玮.集值微分方程初值问题的拟线性化方法[J].河北大学:自然科学版,2011,31(1):1—6.WANGPeiguang,GAOWei.Quasilinearizationofinitialvalueproblemforsetdifferentiale quati0ns[J]JournalofHebeiUniversity:NaturalScienceEdition,2011,31(1):1—6.(责任编辑:王兰英)。
Banach空间中n阶非线性微分积分方程初值问题解的存在与唯一性

Banach空间中n阶非线性微分积分方程初值问题解的存在
与唯一性
柴国庆;李必文
【期刊名称】《数学杂志》
【年(卷),期】2001(21)2
【摘要】用压缩映象原理研究Banach空间中定义在无界域上的n阶非线性Volterra型微分积分方程初值问题,获得了解的存在与唯一性及其迭代逼近新结果.【总页数】6页(P189-194)
【作者】柴国庆;李必文
【作者单位】湖北师范学院数学系;湖北师范学院数学系
【正文语种】中文
【中图分类】O175.15
【相关文献】
1.Banach空间半直线上奇异脉冲积分-微分方程初值问题解的存在性和唯一性 [J], 赵成龙;刘衍胜
2.Banach空间一阶脉冲积分-微分方程初值问题解的存在性 [J], 姜燕君;刘立山
3.Banach空间n阶非线性积分-微分方程初值问题解的存在性 [J], 柴国庆
4.Banach空间二阶微分方程初值问题解的存在唯一性 [J], 陈芳启;陈予恕
5.Banach空间中二阶常微分方程初值问题解的存在唯一性 [J], 洪世煌;胡适耕因版权原因,仅展示原文概要,查看原文内容请购买。
Chebyshev配置点法解Volterra型积分微分方程

Chebyshev配置点法解Volterra型积分微分方程吴华;张珏【摘要】采用Chebyshev配点法求解Voherra型积分微分方程,首先将Volterra 型积分微分方程重新写成一个第二类的线性积分方程组,然后将方程组中的被积函数用Lagrange基函数展开,再将Lagrange基函数用Chebyshev多项式展开,在L<,∞>范数下作误差分析,最后用数值算例来证明该方法的可行性.%A Chebyshev-collocation spectral method is developed for Volterra type integro-differential equations. The Volterra type integro-differential equation as two linear integral equations of the second kind is rewrited, and the integrand with Lagrange basis functions and the Lagrange basis functions in terms of the Chebyshev polynomials are expanded. An error analysis is conducted based on the L∞ -norm.Numerical results are presented to demonstrate effectiveness of the proposed method.【期刊名称】《上海大学学报(自然科学版)》【年(卷),期】2011(017)002【总页数】7页(P182-188)【关键词】Chebyshev配置点法;积分微分方程;Lagrange基函数;Chebyshev权【作者】吴华;张珏【作者单位】上海大学,理学院,上海,200444;上海大学,理学院,上海,200444【正文语种】中文【中图分类】O175.6由于 Volterra型积分微分方程带有记忆性质,与传统的常微分方程和偏微分方程相比,有着本质区别,因此,该方程的数值求解也更为困难.如何快速高效地求解这类方程,一直受到许多学者的关注.由于当解充分光滑时,用谱方法求出的近似解具有很高的精度,所以谱方法经常被运用到近似计算中.因此,本研究采用 Chebyshev配点法解下述 Volterra型积分微分方程:式中,函数 f(t),g(t)和核函数 R(t,s)都是充分光滑的.为了分析方便,将式 (1)和 (2)转换为定义在[-1,1]上的等价问题.通过变量替换,将式(1)和(2)转化为Jiang等[1]把给定的积分微分方程转换成第二类积分方程的方程组,因此,解积分微分方程的主要工作就转变成求解这个积分方程组.Atkinson[2]提出了多种求解第二类积分方程的数值方法,并且作了误差分析.Tang等[3]提出了用 Legendre配点法求解Volterra型积分方程的方法.Chen等[4-5]分别用Chebyshev配点法和Jacobi配点法求解了带有奇异核的第二类 Volterra型积分方程.Elnagar等[6]用Chebyshev谱方法求解了非线性 Volterra-Hammerstein积分方程.Tian[7]使用谱元法求解Volterra型积分方程,并将谱元法与梯形公式和Simpson公式作了比较,得到谱元法具有更高的收敛阶.Tang等[8]提出了谱处理技术在提高数值解的精度上的应用.另外,其他研究者还提出了获得带有奇异核 Volterra型积分方程和积分微分方程的高阶收敛性质的方法[9-11].本研究主要是利用 Chebyshev配点法求出Volterra型积分微分方程 (3)和 (4)的数值解.虽然Legendre点在计算中具有很高的稳定性,但是Legendre多项式的权是隐式的,不便于计算.然而,Chebyshev点和权能够直接使用,计算方便.而且Chebyshev多项式能够利用快速傅里叶变换 (fast Fourier transform,FFT),因此,Chebyshev谱方法更常用于工程计算中.一般情况下,Chebyshev谱方法多用于计算带有奇异核的积分方程,而不带奇异核的积分方程经常用 Legendre谱方法来计算.本研究采用 Chebyshev谱方法计算不带奇异核的积分方程,并对该方法作了严格的误差分析,用数值算例证明了其收敛阶.本研究中,C代表一个与N无关的正常数,但是依赖于所给方程 a(x),b(x),K(x,s)的边界.类似文献 [1]中的方法,将式 (3)和 (4)写为一种新的形式.记z=u′,则式 (3)等价于一个关于 z的第二类线性 Volterra型积分方程,即出,此方法能达到预期的谱收敛精度.之后,可用 Gauss积分公式来求解上面的方程组中的积分部分.【相关文献】[1] JIANG Y J.On spectral methods for Volterra type integro-differentialequations[J].JComput App lMath,2009,230(2):333-340.[2] ATKINSON K E. The numerical solution of intgral equations of the second kind [M]. Cambridge:Cambridge University Press,1997.[3] TANG T,XU X,CHENG J.On spectral methods for Volterra type integral equations and the convergence analysis[J].J ComputMath,2009,26:825-837.[4] CHEN Y P,TANG T.Convergence analysis for the Cheb yshev collocation methods to Volterra integral equations with aweakly singular kernel[J].SIAM JNumerAnal,2009,233:938-950.[5] CHEN Y P,TANG T.Convergence analysisof the Jacobi spectral-collocation methods for Volterra integral equations with a weakly singular kernel[J].Math Comput,2010,79:147-167.[6] ELNAGAR GN,KAZEM IM.Chebyshev spectral solution of nonlinear Volterra-Hammerstein integral equations[J].J Comput App Math,1996,76:147-158.[7] TIAN H C. Spectral methods for volterrra integral equations[D].Burnaby:Simon Fraser University,1995.[8] TANG T,XU X.Accuracy enhancement using spectral postprocessing for differential equations and integral equations[J].Commun Comput Phys,2009,5:779-792.[9] BRUNNER H.Polynomial spline collocation methods for Volterra integro-differential equationswith weakly singular kernels[J].IMA JNumer Anal,1986,6(2):221-139.[10] DIOGO T,MCKEE S,TANG T.Collocation methods for second-kind Volterra integral equations with weakly singular kernels[C]∥Proceedingsof The Royal Society of Edinburgh.1994,124A:199-210.[11] TANG T. Superconvergence of numerical solutions to weakly singular Volterra integro-differential equations[J].Numer Math,1992,61:373-382.。
Volterra型微分G积分方程的数值格式构造及理论分析

第36卷第4期2022年7月兰州文理学院学报(自然科学版)J o u r n a l o fL a n z h o uU n i v e r s i t y ofA r t s a n dS c i e n c e (N a t u r a l S c i e n c e s )V o l .36N o .4J u l .2022收稿日期:2022G03G13基金项目:国家自然科学基金(11861068);新疆维吾尔自治区自然科学基金杰出青年基金项目(2022D 01E 13)作者简介:罗紫洋(1996G),男,新疆昌吉人,在读硕士,研究方向:微分方程理论及数值模拟.E Gm a i l :466949841@q q .c o m.∗通讯作者:张新东(1981G),男,江苏徐州人,教授,硕士生导师,研究方向:微分方程理论及数值模拟.E Gm a i l :L i Ga o yu a n 1126@163.c o m ㊀㊀文章编号:2095G6991(2022)04G0010G05V o l t e r r a 型微分G积分方程的数值格式构造及理论分析罗紫洋,安文静,张新东(新疆师范大学数学科学学院,新疆乌鲁木齐830017)摘要:本文研究了一类V o l t e r r a 型微分G积分方程的数值格式构造及其理论分析.格式构造方面,利用有限差分方法进行时间和空间离散,对于积分项采用复合梯形求积公式进行处理.最后,给出了数值格式的稳定性分析和误差估计,其误差的收敛阶为Ο(τ+h 4),其中τ为时间步长,h 为空间步长.关键词:V o l t e r r a 型微分G积分方程;有限差分;复合梯形求积公式;稳定性;误差分析中图分类号:O 241.82㊀㊀㊀文献标志码:AN u m e r i c a l S c h e m eC o n s t r u c t i o na n dT h e o r e t i c a lA n a l y s i s f o r I n t e g r a l GD i f f e r e n t i a l E q u a t i o n s o fV o l t e r r aT y pe L U OZ i Gy a n g ,A N W e n Gj i n g ,Z HA N G X i n Gd o n g(S c h o o l o fM a t h e m a t i c a l S c i e n c e s ,X i n j i a n g N o r m a lU n i v e r s i t y ,U r u m qi 830017,C h i n a )A b s t r a c t :I n t h i s p a p e r ,an u m e r i c a l s c h e m ec o n s t r u c t i o na n dt h e o r e t i c a l a n a l y s i sw a sm a i n l yd e v e l o p e d f o rs o l v i n g t h e V o l t e r r ai n t e g r a l Gd i f f e r e n t i a le q u a t i o n s .I nt e r m so fs c h e m ec o n Gs t r u c t i o n ,t h e f i n i t ed i f f e r e n c ea p p r o x i m a t i o ni s s u e df o rt i m ea n ds pa c ed e r i v a t i v e ,a n dt h e c o m p o u n d t r a p e z o i d a l q u a d r a t u r e f o r m u l a i s s u e d f o r i n t e g r a l t e r m.F i n a l l y,t h eu n c o n d i t i o n a l s t a b i l i t y a n d c o n v e r g e n c ew e r e g i v e n .I na d d i t i o n ,t h ee r r o r a n a l y s i sw a s c a r r i e do u t ,w h i c h s h o w e d t h a t t h e c o n v e r g e n c e o r d e rw a s O (ι+h 4),w h e r e ιw a s t h e t i m e s t e p ,a n d h w a s t h e s p a c e s t e p.K e y wo r d s :V o l t e r r ai n t e g r a l Gd i f f e r e n t i a le q u a t i o n ;f i n i t ed i f f e r e n c e ;c o m p o u n dt r a p e z o i d a l q u a d r a t u r e f o r m u l a ;s t a b i l i t y ;e r r o r a n a l ys i s 0㊀引言V o l t e r r a 型微分G积分方程(V I D E s )源于19世纪末20世纪初,是在竞争或人口增长模型中所引入的一类具有积分内核的方程.V I D E s 中的积分项具有记忆性质[1],这一性质是其在物理领域有着广泛应用的重要原因.例如:粘弹性方程[2]㊁具有记忆性的热传导方程[3G4],以及核反应堆中的热交换过程等.由于求解此类方程的解析解较为困难,因此研究此类方程的数值解具有重要的理论价值和实际意义.处理微分G积分方程的主要方法为差分法.1897年,V o l t e r r a [5]在其著作中对该类方程的数值解法已有所研究.1974年,B r u n n e r [6]采用隐式R u n g e GK u t t a 方法求解V o l t e r r a 型微分G积分方程,并分析了此类方程数值解的稳定性.1992年,陈传淼等[7]利用内积近似一类微分G积分方程的积分项,并对空间项采用有限元方法,得到了相应的误差估计.1993年,汤涛[8]运用梯形求积技巧构造了一类微分G积分方程中积分项的数值格式,使其时间误差收敛阶可以达到Ο(τ3/2),并给出具体的理论分析.2005年,S h a h m o r a d [9]基于T a u方法分析了线性F r e d h o l m GV o l t e r r a 型微分G积分方程的有效误差.2006年,A m i r a l i ye v 等[10]在均匀网格中构建V o l t e r r a 型微分G积分方程的数值格式,并证明该格式的时间收敛阶为一阶均匀收敛.2008年,徐大[11]在希尔伯特空间中分析验证了有限差分方法求解线性V o l t e r r a 型方程的稳定性;同年,T a r i 和S h a h m o r a d [12]基于勒让德多项式研究二维线性F r e d h o l m 积分方程,并给出了数值格式的误差边界.2010年,W a z w a z [13]提出结合L a p l a c e 变换和A d o m i a n 分解法,构造非线性V o l t e r r a 型微分G积分方程的数值格式,并通过实验证明该方法的有效性;与此同时,Z a r e b Gn i a [14]利用S i n c 方法求解V o l t e r r a 型微分G积分方程,并验证该方法所具有的优良性,不仅收敛速度较快,且不存在当使用其他数值方法时常见的不稳定问题.2021年,C i m e n 和C a k i r [15]基于正交规则和积分恒等式构造了一类F r e d h o l m 型微分G积分方程一种新的差分格式,并分析证明该格式的稳定性和收敛性.2022年,S a n t r a 和M o h a Gpa t r a [16]利用复合梯形公式逼近V o l t e r r a 型微分G积分方程中的积分项,使其误差时间收敛阶为Ο(τ).本文考虑如下V o l t e r r a 型微分G积分方程(V I D E ):∂u (x ,t )∂t -∂2u (x ,t )∂x 2=ʏt0K (x ,t -s )u (x ,s )d s +f (x ,t ),(x ,t )ɪΩ;u (x ,0)=φ(x ),∀x ɪ[0,L ];u (0,t )=g 1(t ),u (L ,t )=g 2(t ),∀t ɪ(0,T ],ìîíïïïïïïïïïï(1)其中:Ω=[0,L ]ˑ(0,T ],空间L >0,且时间T >0;积分核K (x ,t -s )在Ω上光滑且有界;f ,φ,g 1和g 2为给定的光滑函数.1㊀紧致差分格式推导为方便起见,本文中不同地方的常数C 可以代表不同的数值.给定两个正整数M 和N ,定义h=L/M 为空间变量x 的步长,τ=T /N 为时间变量t 的步长,所以x j =j h ,j =0,1, ,M ,t n =n τ,n =0,1, ,N .用u n j 和U nj 分别表示函数u 在点(x j ,t n )处的精确解和数值解.为了方便表达及书写,给出如下算子定义,δ2xU j =U j -1-2U j +U j +1h 2,H U j =112(U j -1+10U j +U j +1)=(1+h 212δ2x )U j ,㊀㊀j =1,2, ,M -1;U j ,j =0,M .ìîíïïïï(2)由复合梯形求积公式可得ʏbag (t )d t =ðn -1i =0ʏt i +1t ig (t )d t =τ2ðn -1i =0[g (t i )+g (t i +1)]+R ,(3)其中:R 为整个区间的截断误差.由T n =τ2ðn -1i =0[g (t i )+g (t i +1)]和中值定理可得R =ʏbag (t )d t -T n =-b -a 12τ2g ᵡ(ξ),ξɪ[a ,b ].(4)引理1[17]㊀假设s (x )ɪC 6[a ,b ],则有112(s x x (x j +1)+10s x x (x j )+s x x (x j -1))-1h2(s (x j +1)-2s (x j )+s (x j -1))=h 4240s (6)(ξj ),其中:ξj ɪ(x j -1,x j +1),j =1,2, ,M -1.接下来对方程(1)进行逐项分析.对于方程(1)的时间项来说,利用T a y l o r 展开,可得∂u (x j ,t n )∂t n=u (x j ,t n +1)-u (x j ,t n )τ+O (τ).(5)对于方程(1)的空间项来说,H ∂2u (x j ,t n +1)∂x 2j=δ2x u (x j ,t n +1)+O (h 4).(6)对于方程(1)的积分项而言,由式(3)和(4)可得11第4期罗紫洋等:V o l t e r r a 型微分G积分方程的数值格式构造及理论分析ʏt n 0K (x j ,t n -s )u (x j ,s )d s =ðn -1i =0ʏt i +1t iK (x j,t n-s )u (x j,s )d s =τ2ðn -1i =0[K (x j ,t n -t i )u (x j ,t i )+K (x j ,t n -t i +1)u (x j ,t i +1)]+O (τ2).(7)如果对方程(1)作用H 算子,可得H ∂u (x j ,t n )∂t n-H ∂2u (x j ,t n +1)∂x 2j =Hʏt n 0K (x j ,t n -s )u (x j ,s )d s +H f (x j ,t n ),进而由式(5)㊁(6)和式(7)可得H u n +1j-τδ2x u n +1j =H u nj+H τ22ðn -1i =0(K n -i -1j u i +1j +K n -i j u i j )éëêêùûúú+H τf n j +R nj ,(8)其中:f nj =f (x j ,t n ),R n j ɤC (τ+h 4).省去式(8)中的截断误差,当j =1,2, ,M-1且n =1,2, ,N 时得到方程(1)的离散格式如下:H U n +1j -τδ2x U n +1j =H U nj +㊀H τ22ðn -1i =0(K n -i -1j U i +1j +K n -i j U i j )éëêêùûúú+H τf n j ,U 0j =φ(x ),U n +10=g 1(t n +1),U n +1M =g 2(t n +1).ìîíïïïïïï(9)2㊀稳定性分析本节将利用能量不等式的方法讨论数值格式的稳定性,并证明数值格式为无条件稳定.首先给出一些符号说明和引理.定义V h 为空间网格函数,V h ={V |V =(V 0,V 1, ,V M ),V 0=g 1(t n +1),V M =g 2(t n +1)}.以下内积和范数在证明中将用到,(U ,V )=h ðM -1j =1U jV j , U 2=(U ,U ).引理1[18]㊀假设U ɪV h ,则有-(δ2xU n H U n )ȡ23δx U n 2.定理1㊀紧致差分格式(9)是无条件稳定的,并且满足如下不等式:H U n +1 ɤC H U 0 +m a x 0ɤp ɤnH f p (),其中:C 为常数.证明㊀设m a x (x ,t )ɪΩK (x ,t ) =K x t .由式(9)可得H U n +1-τδ2xU n +1=H U n +H τ22ðn -1i =0(K n -i -1U i +1+K n -i U i )éëêêùûúú+H τf n ,(10)对式(10)两边同乘H U n +1,并在Ω上积分可得(H U n +1,H U n +1)-τ(δ2xU n +1,H U n +1)=(H U n ,H U n +1)+τ22ðn -1i =0(H (K n -i -1U i +1),H U n +1)+τ22ðn -1i =0(H (K n -i U i ),H U n +1)+τ(H f n ,H U n +1).(11)由引理1可知-(δ2xU n H U n )ȡ0.由式(11)可得如下不等式,(H U n +1,H U n +1)ɤ(H U n ,H U n +1)+τ22ðn -1i =0(H (K n -i -1U i +1),H U n +1)+τ22ðn -1i =0(H (K n -i U i ,H U n +1))+τ(H f n ,H U n +1).(12)接下来将用归纳法证明.由式(12)可知,当n =0时,(H U 1,H U 1)ɤ(H U 0,H U 1)+τ(H f 0,H U 1).再由柯西G施瓦茨不等式可得H U 1 ɤ H U 0 +τ H f 0 ɤC H U 0 + H f 0 ().假设对任意k ,k ɤn 时结论成立,即 H U k ɤC H U 0 +m a x 0ɤp ɤnH f p ().(13)由式(12),当k =n +1时,有(H U n +1,H U n +1)ɤ(H U n ,H U n +1)+τ22ðn -1i =0(H (K n -i -1U i +1),H U n +1)+τ22ðn -1i =0(H (K n -i U i ),H U n +1)+τ(H f n ,H U n +1).由上述不等式,结合式(13)可得 H U n +1 ɤ21㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀兰州文理学院学报(自然科学版)㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第36卷H U n+τ2K x t2ðn -1i =0H U i +1 +τ2K x t2ðn -1i =0H U i +τ H f n ɤC H U 0+m a x 0ɤp ɤn -1H f p()+τK x tT C H U 0+m a x 0ɤp ɤn -1H f p ()+τ H f nɤ(C +τK x tT C ) H U 0 +(C +τK x tT C +τ)m a x 0ɤp ɤnH f p ɤC H U 0+m a x 0ɤp ɤnH f p(),其中,C 为常数.定理1证明完毕.3㊀误差估计本节将讨论数值格式的收敛性,并推导出数值格式的误差收敛阶为O (τ+h 4).记e nj =u n j -U n j ,j =0,1, ,M 且n =0,1, ,N .定理2㊀假设u nj 为方程(1)的精确解,U n j 为离散格式(9)的数值解,则如下不等式成立, H e n +1jɤC (τ+h 4),其中,C 为常数.证明㊀设m a x (x ,t )ɪΩK (x ,t ) =K x t .用式(8)减去式(9),再由e nj 的定义可得H e n +1j -τδ2x e n +1j=H e n j+H [τ22ðn -1i =0(K n -i -1j e i +1j +K n -i j e i j )]+R n j .(14)将式(14)两边同乘H e n +1j,并在Ω上积分可得(H e n +1j ,H e n +1j )-τ(δ2x e n +1j ,H e n +1j)=(H e n j ,H e n +1j)+τ22ðn -1i =0(H (K n -i -1j e i +1j ),H e n +1j )+τ22ðn -1i =0(H (K n -i j e i j ),H e n +1j )+(R n j ,H e n +1j).(15)由引理1和式(15),有(H e n +1j ,H e n +1j )ɤ(H e n j ,H e n +1j)+τ22ðn -1i =0(H (K n -i -1j e i +1j ),H e n +1j )+τ22ðn -1i =0(H (K n -i j e i j ),H e n +1j )+(R n j ,H e n +1j).(16)类似定理1的证明,将采用归纳法.由式(16)可知,当n =0时,可得(H e 1j ,H e 1j )ɤ(H e 0j ,H e 1j )+(R 0j ,H e 1j ),再由柯西G施瓦茨不等式及R 0j ɤC (τ+h 4)可得 H e 1j ɤ H e 0j +R 0j ɤC (τ+h 4).假设对任意k ,k ɤn 时,结论成立,即 H e kj ɤC (τ+h 4).(17)由式(16),当k =n +1时,有(H e n +1j ,H e n +1j )ɤ(H e n j ,H e n +1j)+τ22ðn -1i =0(H (K n -i -1j e i +1j ),H e n +1j )+τ22ðn -1i =0(H (K n -i j e i j ),H e n +1j )+(R n j ,H e n +1j).由上述不等式,结合式(17)可得 H e n +1jɤ H e nj+τ2K x t 2ðn -1i =0H e i +1j +τ2K x t 2ðn -1i =0H e ij + R n j ɤC (τ+h 4)+τK x tT (C (τ+h 4))+C (τ+h 4)ɤm a x {C ,τK x tT C ,C }(τ+h 4)ɤC (τ+h 4),其中,C 为常数.定理2证明完毕.4㊀结论本文通过有限差分方法构造了一类V o l t e r r a型微分G积分方程的数值格式,分析了格式的稳定性和误差估计,得到其误差收敛阶为O (τ+h 4).相较于解决该问题原有的数值格式,在空间收敛阶上有进一步的提升.在今后的研究工作中,将以构建收敛速度更快,误差更低的数值格式为目标,并将此方法应用到分数阶微分G积分方程的数值求解.参考文献:[1]MO L F E R S D O R F L V.O ni d e n t i f i c a t i o no fm e m o r yk e r n e l s i n l i n e r t h e o r y of h e a t c o n d u c t i o n [J ].M a t h e Gm a t i c a lM e t h o d s i n A p pl i e dS c i e n c e s ,1994,17:919G932.[2]R C N A R D Y M.M a t h m e a t i c a l a n a l ys i so fv i s c o e l a s t i c f l o w s [J ].A n n u a lR e v i e w o fF l u i d M e c h a n i c ,1989,21:21G36.[3]G U R T I N M E ,P I P K I N AC .A g e n e r a l t h e o r y of h e a t c o n d u c t i o nw i t h f i n i t ew a v e s pe e d [J ].A r c h i v ef o rR a Gt i o n a lM e c h a n i c s a n dA n a l y s i s ,1968,31:113G126.[4]M I L L E R R K.A ni n t eg r o Gd i f f e r e n t i a le qu a t i o nf o r g i r dh e a tc o n d u c t o r s w i t h m e m o r y [J ].J o u r n a lo f 31第4期罗紫洋等:V o l t e r r a 型微分G积分方程的数值格式构造及理论分析M a t h e m a t i c a l A n a l y s i sa n d A p p l i c a t i o n s,1978,66:313G332.[5]V O L T E R R A V.S o p r aa l c u n e q u e s t i o n i d i i n v e r s i o n e d i i n t e g r a l i d e f i n i t i[J].A n n a l i d iM a t e m a t i c aP u r ae d A p p l i c a t a,1897,25:139G178.[6]B R U N N E R H.I m p l i c i t R u n g eGK u t t am e t h o d s o f o p t iGm a lo r d e rf o r V o l t e r r ai n t e g r oGd i f f e r e n t i a le q u a t i o n s [J].M a t h e m a t i c s o fC o m p u t a t i o n,1984,42:95G109.[7]C H E NC M,T HOM E E V,WA H L B I NLB.F i n t e e lGe m e n t a p p r o x i m a t i o n o f a p a r a b o l i c i n t e g r oGd i f f e r e n t iGa l e q u a t i o nw i t haw e a k l y s i n g u l a rk e r n e l[J].M a t h eGm a t i c a l o fC o m p u t a t i o n,1992,58:587G602.[8]T A N G T.Af i n i t ed i f f e r e n c es h e m ef o r p a r t i a l i n t eGg r oGd i f f e r e n t i a l e q u a t i o n sw i t haw e a k l y s i n g u l a rk e rGn e l[J].A p p l i e d N u m b e r i c a l M a t h e m a t i c s,1993,11:309G319.[9]S HA HMO R A DS.N u m e r i c a l s o l u t i o no f t h e g e n e r a l f o r ml i n e a rF r e d h o l mGV o l t e r r a i n t e g r oGd i f f e r e n t i a l eGq u a t i o n sb y t h eT a um e t h o dw i t ha ne r r o r e s t i m a t i o n [J].A p p l i e d M a t h e m a t i c sa n d C o m p u t a t i o n,2005,167:1418G1429.[10]AM I R A L I Y E V G M,S E V G I N S.U n i f o r m d i f f e rGe n c em e t h o df o r s i ng u l a r l yp e r t u r b e dV o l t e r r a i n t eGg r oGd i f f e r e n t i a l e q u a t i o n s[J].A p p l i e d M a t h e m a t i c sa n dC o m p u t a t i o n,2006,179:731G741.[11]X U D.S t a b i l i t y o f t h ed i f f e r e n c et y p e m e t h o d sf o r l i n e a rV o l t e r r ae q u a t i o ni n H i l b e r ts p a c e s[J].N uGm e r i s c h eM a t h e m a t i k,2008,109:571G595.[12]T A R IA,S HA HMO R A DS.Ac o m p u t a t i o n a lm e t hGo d f o r s o l v i n g t w oGd i m e n s i o n a l l i n e a rF r e d h o l mi nGt e g r a le q u a t i o n so fs e c o n d k i n d[J].T h e A n z i a mJ o u r n a l,2008,49:543G549.[13]WA Z WA Z A M.T h ec o m b i n e dL a p l a c et r a n s f o r mGA d o m i a nd e c o m p o s i t i o n m e t h o df o rh a n d l i n g n o nGl i n e a r V o l t e r r ai n t e g r oGd i f f e r e n t i a le q u a t i o n s[J].A p p l i e d M a t h e m a t i c sa n dC o m p u t a t i o n,2010,216:1304G1309.[14]Z A R E B N I A M.S i n c n u m e r i c a l s o l u t i o n f o r t h eV o l tGe r r ai n t e g r oGd if f e r e n t i a le q u a t i o n[J].C o mm u n i c aGt i o n s i n N o n l i n e a rS c i e n c ea n d N u m e r i c a lS i m u l aGt i o n,2010,15:700G706.[15]C I M E N E,C A K I R M.Au n i f o r m n u m e r i c a lm e t h o df o rs o l v i ng s i n g u l a r l yp e r t u r b e d F r e dh o l mi n t e g r oGd i f fe r e n t i a l p r o b l e m[J].C o m p u t a t i o n a l a n dA p p l i e dM a t h e m a t i c s,2021,42:1G14.[16]S A N T R AS,MOHA P A T R AJ.An o v e l f i n i t e d i f f e rGe n c e t e c h n i q u ew i t h e r r o r e s t i m a t ef o r t i m e f r a c t i o nGa l p a r t i a li n t e g r oGd i f f e r e n t i a le q u a t i o n o f V o l t e r r at y p e[J].J o u r n a l o f C o m p u t a t i o n a la n d A p p l i e dM a t h e m a t i c s,2022,400:113746.[17]L U O M,X U D,L I L M.A c o m p a c td i f f e r e n c e s c h e m ef o r a p a r t i a li n t e g r oGd i f f e r e n t i a le q u a t i o nw i t ha w e a k l y s i n g u l a rk e r n e l[J].A p p l i e d M a t h eGm a t i c a lM o d e l l i n g,2015,39:947G954.[18]A K B A R M.C o m p a c t f i n i t e d i f f e r e n c e s c h e m e f o r t h e s o l u t i o n f o r a t i m e f a r a c t i o n a l p a r t i a l i n t e g r oGd i f f e rGe n t i a l e q u a t i o n w i t ha w e a k l y s i n g u l a rk e r n e l[J].M a t h e m a t i c a l M e t h o d si n t h e A p p l i e d S c i e n c e s,2017,40:7627G7639.[责任编辑:赵慧霞]41㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀兰州文理学院学报(自然科学版)㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第36卷。
第二类Volterra积分方程的一种特殊解法

0 引 言
首 先考 虑第 二类 Vot , 者需要储 存 网格点 的所 前
有值 而 后者 只需 固定 数 目网格点 的信 息. 于方 程 对 () 1 的这 种储存 , 加可 以接 受应 用谱方法 的全 局基 更
函数. () 1
假 定 () 的 解 充分 光 滑 , 此 情况 下 , 必 要 1式 在 有 利用 高次数 值方 法 例如谱 方法来 解 方程 ()对 于方 1.
程()现已有很多数值方法例如配制法 、 1, 内积积分
法 , B u n r . 而 , 少 有 人 提 及 利 用 谱 逼 近 见 r n eC 然 ¨ 很
王 金 婵
( 州学 院数 学 系 , 东德 州 23 2 ) 德 山 5 0 3
擒 薹 : 究 了求 解 Vo er 积 分 方 程 的 方 法 , 点 介 绍 了 基 于 谱 方 法 解 决 V l ra型 积 分 方 程 的 一 种 新 研 l ra型 t 重 ot r e 的 数值 解 法 , g n r 配 制 法 得 到 充 分 的 应 用 , 进 行 了严 格 的误 差 分 析 , 明 在 核 函 数 和 原 函 数 充 分 光 滑 时 , 值 l e de e 并 表 数
得 到高 次精度 的 主要 困难 在 于解 ( ) 中 的积 3式
分项 , 别 的 , 于 充分 小 的值 , () 特 对 “ s 有很 少 的信
息可 以利 用 , 为此 , 把积 分 区间 [ 1 ]转换 到 [ 一 , 一
1 1上 , ,] 然后 选择 恰 当的求 积 规 则. 实 上 , 事 首先 作
一
个 线 性变换
收稿 日期 :2 1 0 — O 0 1— 7 2
作 1 介 :王金蝉( 90一 , , I l 18 )女 山东庆云人 , 讲师 , 硕士 , 究方向 : 研 信息 与计算科学.
第二类Volterra积分方程的一种特殊解法

第二类Volterra积分方程的一种特殊解法王金婵【摘要】In this paper,we study the method of solving the second kind of Volterra integral equations,propose a new numerical method for solving Volterra integral equations based on spectral method,Legendre preparation method is fully applied and a rigorous error analysis is done.The results show that the numerical error is the index fell.When kernel function and the original function is sufficiently smooth%研究了求解Volterra型积分方程的方法,重点介绍了基于谱方法解决Volterra型积分方程的一种新的数值解法,legendre配制法得到充分的应用,并进行了严格的误差分析,表明在核函数和原函数充分光滑时,数值误差是指数下降的.【期刊名称】《德州学院学报》【年(卷),期】2012(028)004【总页数】3页(P11-13)【关键词】legendre谱方法;Volterra积分方程;收敛性分析【作者】王金婵【作者单位】德州学院数学系,山东德州253023【正文语种】中文【中图分类】O175.5首先考虑第二类Volterra积分方程其中k(x,s,u(s))是核函数.假定(1)式的解充分光滑,在此情况下,有必要利用高次数值方法例如谱方法来解方程(1).对于方程(1),现已有很多数值方法例如配制法、内积积分法,见Brunner[1].然而,很少有人提及利用谱逼近法,在文献[2]中,作者用chebyshev谱方法在多重代数精度下求解第一类,但这些方法理论上不能获得高次精度.这种方法一直未能得到很好的应用.众所周知,Fredholm方程类似边界值问题见文献[5],因此,许多边界值问题有效的数值方法例如谱方法能直接用来解决Fredholm方程,而方程(1)类似初值问题,因此用谱方法解是非常困难的.主要原因在于方程(1)是一个局部方程而谱方法用到全部基函数.主要困难在于如何执行这种算法使其能最终得到精确谱.此外,方程(1)的数值解法有可能不同于标准初值问题.从这个意义上说,前者需要储存网格点的所有值而后者只需固定数目网格点的信息.对于方程(1)的这种储存,更加可以接受应用谱方法的全局基函数.不失一般性,假定解的区域为[-1,1],第二类线性积分方程一维形式如下设定N+1个配置点作为Legendre Gauss点集假定方程(2)在xi处成立,则有得到高次精度的主要困难在于解(3)式中的积分项,特别的,对于xi充分小的值,u(s)有很少的信息可以利用,为此,把积分区间[-1,xi]转换到[-1,1]上,然后选择恰当的求积规则.事实上,首先作一个线性变换则(3)式转换为然后利用N+1个点的Gauss积分法,以及Legendre权重{ωi},则有其中其中{θj}j=0,…N,与配置点相一致.用ui,0≤i≤N代替u(s(xj,θj)),用lagrange插值多项式表示u即.其中Fj是第j个lagrange基函数,代入(7)式,得从(8)式可以看出,为了计算u(xi)的近似值,需要的完全解信息和的半局部信息.其中-1≤s(xi,θj)≤xi,这不同于配制法或内积积分法,原因在于它们用到和的半局部信息.谱配制算法的执行令得到方程的矩阵形式其中矩阵A中元素如下给出下面讨论Fj(s(xi,θp)的计算效率,由于αp,j是Fj的离散的多项式系数,它的递推关系如下(见文献[5])其中由(10)式和(11)式可得结合LP(s)的递推公式,能有效的得出Fj(s(xi,θp)),这种情况也可以推广到非线性方程及二维情况(见文献[3]).下面从数值方面对Volterra方程进行收敛性分析,目的在于表明它的收敛率是指数型的,既谱精度可以从以提出的谱逼近中得到.引理1[5]假设N+1个点Gauss求积公式,及lagrange权重应用积分内积uφ,其中u∈Hm(I),I=(-1,1)m≥1,φ∈PN,则存在不依赖于N的常数C,使得其中引理2[5]假定u∈Hm(I),INu表示与它的N+1个Gauss点相关的插值多项式,即则引理3[5]假定Fj(x)是第N个Gauss点相关的Lagrange插值多项式,则其中是一个有界常量.引理4(Gronwall不等式),如果非负积分函数E(t)满足其中 G(t)是可积函数,则定理1 令u是Volterra方程(2)的精确解,假定其中uj由(8)给出,Fj(x)是同高斯点相关的第j个Lagange基函数,若,则对于这里N充分大,s(xi,θ)由(6)式给出,C是不依赖于N的常数(证明见文献[3]).由定理1知,收敛速度似乎是不可以选择的,应为Ο(Nm),而不应为(11)式给出的,如果引理3的估计能够进一步改进,这个结果应该是正确的,一个可能的改进是证明,假设这是正确的,则在定理1中应用‖J1‖L1(I)=O(N-m),在这种情况下,收敛阶O(N-m)能够得到.本文总结了Volterra积分方程的解法,并对他们的优缺点进行分析.同时给出了不同解法的收敛情况.重点介绍了基于谱方法解决二维第一类型Volterra积分方程的数值解法,借助离散Gronwall不等式,给出了一系列漂亮的结果.【相关文献】[1]H.Brunner.Couocation Methods for Volerra Integral and Related Functional Equations Methods[M].Cambridge University Press,2004.[2]H.Fujiwara.High-accurate Numerical Methed for Integral Equations of the First Kind under Muttipleprecision Arithmetic[M].Preprint RIMS,kyoto,University,2006. [3]T.Tang,X.Xu,J.Cheng.On Spectral Methods for Volerra Type Integral Equations and the Convergence Analysis[J]pwl Math,2007.[4]H.C.Tian.Spectral Methods for Volerra Integral Equations[M].Msc Thesis,Simon Fraser University,1995.[5]C.Caunto,M.Y.Hussaini,A.Quarteroni,etal.Spectral Methodes Fundamentals inSingle Domains[M].Springer-Verlag,2006.。
高阶分数阶微分方程系统的解的注记

1
引言
分数阶导数这一概念的的提出可追溯到 1695 年,其理论在 19 世纪得到了很大的发
展.有两种形式的分数阶导数,分别是 Caputo 型和 Riemann-Liouville 型,对于相当广 泛的一类函数,这两种形式的导数是相同的[1] .比如在阶数 α 为负实数和正整数时,它 们是等价的.分数阶微积分在数学、物理、化学、控制等领域都有广泛的应用,已经吸 引了许多科学家和工程学家的关注和研究[2-7] .关于单个分数阶微分方程解的存在性与唯 一性已有大量的研究[8-13] .但对分数阶微分方程系统 (分数阶微分方程组) 解的探讨还不 是很多[14] .文献 [15] 探讨了阶数为 q ∈ (2, 3] 的分数阶微分方程解的存在性与唯一性.文 献 [16] 探讨了初值为柯西型的阶数为 q ∈ (2, 3] 的分数阶微分方程系统初值问题解的存在 性和幂型估计.本文将在文献 [15,16] 的基础上,继续探讨分数阶微分方程系统解的存在 性与唯一性. 令 Mn 表示实数域或复数域上的 n 阶矩阵.若
高阶分数阶微分方程系统的解的注记∗
郑艳萍1 , 李胜利2
(1- 太原师范学院数学系,太原 030012; 2- 太原理工大学数学学院,太原 030024)
摘 要: 分数阶导数在描述不同物质的记忆与遗传性质方面提供了有力的工具.在科学和工程 的不同领域,都用分数阶微分方程组来描述动力系统.本文主要探讨分数阶微分方程 系统初值问题局部解的存在性与唯一性.对于线性系统,运用 Schur 分解定理,给出 其局部解的存在性与唯一性,并通过举例说明该方法是有效的.对于非线性系统,利 用 Schauder 不动点定理,给出了解的存在性;运用 Banach 不动点定理,给出了解的 唯一性. 关键词: Riemann-Liouville 型导数;局部解;三角形 分类号: AMS(2000) 35Q51; 35Q58 中图分类号: O175.14 文献标识码: A
一类Fredholm-Volterra型积分方程的数值求解

应用泛 函分析学报
第1 7 卷
数获取 了高 阶线性 F r e d h o l m— V o l t e r r a积分微分 方程的近似解 . Ya l s i n b a s 在文献 f 3 ] 中进 一步的 解决了非线性 F r e d h o l m— V o l t e r r a积分方程的近似解算法. Y u s u f o g l u和 E r b a s [ 】 基于插值的方 法求 出了线性 F r e d h o l m— V o l t e r r a积分 方程 的数值 解. Ma j e e d和 0mr a n 【 5 】 利用复化梯形 和复化 辛普森离散 的方法 求得 了线性 F r e d h o l m— V o l t e r r a积分方程数值 解. 本文应用 G a l e r k i n法离散 F r e d h o l m— V o l t e r r a 积分方程, 并以契贝晓夫多项式作为最佳平方 逼 近元, 解决 了这类方 程的数值解 问题 .
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
编辑版word 一类第二种非线性Volterra积分方程积分数值解方法 1前言 微分方程和积分方程都是描述物理问题的重要数学工具,各有优点.相对于某种情况来说,对于某种物理数学问题,积分方程对于问题的解决比微分方程更加有优势,使对问题的研究更加趋于简单化,在数学上,利用积分形式讨论存在性、唯一性往往比较方便,结果也比较完美,所以研究积分方程便得越来越有用,日益受到重视. 积分方程的发展,始终是与数学物理问题的研究息息相关.一般认为,从积分发展的源头可以追溯到国外的数学家克莱茵的著作《古今数学思想》,该书是被认为第一个清醒的认为应用积分方程求解的是Abel.Abel分别于1833年和1826年发表了两篇有关积分方程的文章,但其正式的名称却是由数学家du Bois-Raymond首次提出的,把该问题的研究正式命名为积分方程。所以最早研究积分方程的是Abel,他在1823年从力学问题时首先引出了积分方程,并用两种方法求出了它的解,第一的积分方程便是以Abel命名的方程.该方程的形式为:
baaxfdttxt)(
)(
)(,该方程称为广义Abel方程,式中a的值在(0,1)之间.当a=21时,
该式子便成为)()(xfdttxxxa.在此之前,Laplace于1782年所提出的求Laplace反变换问题,当时这个问题就要求解一个积分方程.但是Fourier其实已经求出了一类积分方程的反变换,这就说明在早些时候积分方程就已经在专业性很针对的情况下得到了研究,实际上也说明了Fourier在研究反变换问题是就相当于解出编辑版word
了一类积分方程.积分方程的形成基础是有两位数学家Fredholm和Volterra奠定的,积分方程主要是研究两类相关的方程,由于这两位数学家的突出贡献,所以这两个方程被命名为Fredholm方程和Volterra方程。后来又有德国数学家D.Hilbert进行了重要的研究,并作出了突出的贡献,由于D.Hilbert领头科学家的研究,所以掀起了一阵研究积分方程的热潮,并出现了很多重要的成果,后来该理论又推广到非线性部分。我国在60年代前,积分方程这部分的理论介绍和相关书本主要靠翻译苏联的相关书籍,那时研究的积分方程基本是一种模式,即用古典的方法来研究相关的积分方程问题,这样使得问题的研究变得繁琐、复杂,在内容方面比较单一、狭隘,甚至有些理论故意把积分方程的研究趋向于复杂化。随着数学研究的高速发展,特别是积分方程近年来的丰富发展,如此单一、刻板的解法已经不能跟上数学研究时代的步伐。在九十年代我国的数学专家路见可、钟寿国出版了《积分方程论》,该书选择2
L空间来讨论古典积分方程,并结合泛
函分析的算子理论来分析积分方程的相关问题。最近出版的比较适合一般读者阅览的积分方程的书有李星出版的《积分方程》,该书从最简单的方法分析研究积分方程的理论问题,并给日后打算研究泛函的读者提供了基本的实例。由于现代的计算机技术高速发展,对于一些比较复杂,难以求解的非线性积分方程逐渐采用了比较有效的数值解,常用的方法有逐次逼近法、Adomian分解法、配置法、haar小波方法、小波-Galerkin方法、泰勒展开等等一些方法. 现在积分方程的应用广泛,很多问题都可以引出积分方程,并可用积分方程来解决.像弹性弦问题、线性系统响应问题、人口增长问题、等时曲线问题等等.还有在空气动力学中研究分子运动,对于非均匀流体中悬浮晶粒的布朗位移,导致了以柯尔莫哥洛夫命名的一类非线性积分方程.在确定飞机机翼方面的研究编辑版word
中,对于气流、升力等问题的计算中也引出了积分方程的研究.现在很多积分方程方面的研究都取得了不错的进展. 2、预备知识 积分方程是一个在积分号下出现待求的函数的方程,称作积分方程。含一个未知函数的线性积分方程的一般形式为:
baxfdyyyxkx)()(),()(
bax,
[1] 我们把积分号中的上下限为常数的积分方程称为Fredholm方程。其中)(xf、),(yxk是已知函数,)(y是未知所要求的函数。一般称)(xf为自由项 ,),(yxk称为积分方程的核,而是积分方程的一个参数,方程的解与相关 而与Fredholm对应的是Volterra积分方程,与Fredholm的一个不同点是Volterra积分号下的上下限中的上限是一个变量而不是一个常数。本论文我们主要研究的是非线性的Volterra积分方程。非线性 Volterra积分方程的形式为xaxfdyyFyxkx)()((),()()
bax, 。同样,)(y为所求。
积分方程的解法: 一、 Fredholm积分方程一般的解法有:有限差分逼近法、逐次逼近法及解核、泛函修正平均法、Fredholm积分方程退化核解法、退化核近似代替法、待定系数法。 二、 Volterra积分方程的常用解法:有限差分逼近法、逐次逼近法、转化为常微分方程的初值问题、第二类Volterra积分方程的数值积分解法。 第二类线性volterra积分方程与第二类线性Fredholm积分方程的一个很大的差别是volterra方程的解不依赖于参数的值,即对于任何的参数该方程都有解,而且解具有唯一性。具体的证明过程参见[],所以得出定理:如果核),(yxk与自编辑版word
由函数 xf在Volterra积分方程的定域是连续函数,那么无论参数取何值该方程都有解,而且解具有唯一性。
3、Volterra积分方程 3.1第一类Volterra积分方程 3.1.1 第一类线性Volterra积分方程 形如:xaxfdyyyxk)()(),(
其中函数),(yxk、)(xf为已知的,)(x为所要求的未知函数,这样的方程叫第一类线性Volterra积分方程.一般来说第一类积分方程由于其不适定性,研究其解跟第二类积分方程有很大的不同,而且比较复杂,在此主要简单介绍一下. 3.1.2 第一类Volteraa积分方程的一种解法 在某种情况下第一类Volterra积分方程通常可以化为第二类volterra积分方程的解,一般对方程两边求导,当方程的),(yxk、)(xf可微,且0),(yxk,就把第
一类该方程化为第二类。化为第二类的形式为:),()()(),(),()(yxkxfdyyyxk
yxk
xxax
,
这样就可以用第二类积分方程的解法来求解.对于一种特殊的第一类Volterra积分方程:Abel方程,Abel方程是Volterra积分方程的一种特殊情况,其形式为: )()()(xfdyyxyxaa
其中当yx时,该方程出现弱奇性。其解可根据定理:假设Abel积分方程)()()(xfdyyxyxaa
的自由项)(xf是连续可微的,而且0)(af,则它有唯一的
解.即 dyyxyfdxdaxxaa1)(
)(sin
)(
3.2第二类volterra积分方程 编辑版word
3.2.1第二类线性Volterra积分方程 第二类方程的形式:
xaxfdyyyxkx)()(),()(
其中)(x是所要求的未知函数,是已知或是需要讨论的参数,跟Fredholm方程一样).(yxk是已知的函数,叫Volterra方程的核,当).(yxk=0时,Volterra方程可以看成特殊形式的Fredholm方程,而且Fredholm方程理论适合于Volterra方程。Volterra方程有自己的特点,例如,Volterr方程没人特征值,对于任意的自由项它都有解。对于第一类的Voltrra方程在某种条件下可以转化为第二类Volterra方程。与第二类线性Fredholm积分方程一样,第二类线性Volterra积分方程也有自己的迭核、解核,其迭核、解核的引出方法跟第二类线性Fredholm一样. 迭核:假设dtytktxkyxkxy),(),(),(2,则
xadttfyxkx)(),()(
22
...... xanndyyfyxkx)(),()( ,我们称),(yxkn为Volterra方程的迭核。
解核:我们称11),();,(nnnyxkyxR为解核
只要知道方程的迭核,就能求得方程的解核,从而求得方程的解。
3.2.2第二类线性Volerra积分方程的解法: 1、逐次逼近法 假设方程具有这样一个形式的解
niiinnxxxx010)()()()(
如果对于逐次方法来说该方程有解,解次方程一般令 )(0xf 编辑版word
xadyyxkxf01),()(
xadyyxkxf12),()(
…… xanndyyxkxf1),()(
那么,对于上述的级数一定收敛,即对级数niiix0)(收敛,可以证明对于任意
的参数方程都有解,依据定理:如果核),(yxk及自由项)(xf是连续的实函数.那么第二类线性Volterra方程 )()(),()(xfdyyyxkxxa
对于任意的参数存在一个唯一的连续解,而且解可以用逐次逼进法求出. 迭核:假设dtytktxkyxkxy),(),(),(2,则
xadttfyxkx)(),()(
22
...... xanndyyfyxkx)(),()( ,我们称),(yxkn为Volterra方程的迭核。
解核:我们称11),();,(nnnyxkyxR为解核, 只要知道方程的迭核,就能
求得方程的解核,从而求得方程的解。 3.2.3非线性第二类Volterra积分方程 非线性第二类volterra积分方程的形式如:xaxfdyyFyxkx)()((),()()
未知函数为)(x,而)(xf、),(yxk、)(xF都是已知的。当方程满足一定条件时,可用逐次逼近法求解。对于第一类非线性Volterra积分方程可以通过转化成第二类非线性Volterra积分方程求解。具体转化的过程参见[]。 非线性Volterra积分方程的数值解
3.3 卷积型Volterra方程的解法 3.3.1第二类卷积型Volterra积分方程的解