非线性延迟积分微分方程线性多步法的渐近稳定性.docVIP

——————————————————————————————————————————————— 非线性延迟积分微分方程线性多步法的渐近稳定性 应??用??数??学 MATHEMATICAAPPLICATA 2010,23(1):194??197 非线性延迟积分微分方程线性 多步法的渐近稳定性 余越昕,江春华 (湘潭大学数学系,湖南湘潭411105) 摘要:本文研究求解非线性延迟积分微分方程的线性多步法的渐近稳定性,其中积分部分采用复化梯形公式计算,结果表明:在问题真解渐近稳定的条件下,A??稳定的线性多步法也是渐近稳定的. 关键词:延迟积分微分方程;线性多步法;渐近稳定性 中图分类号:O241.8??????AMS(2000)主题分类:65L06 文献标识码:A??????文章编号:1001??9847(2010)01??0194??04 2002年,针对线性延迟积分微分方程 y??(t)=Ly(t)+My(t-??)+K y(??)d??,t-??t(1) 这里L,M,K!CN?N为常矩阵,Koto[1]讨论了方程理论解渐近稳定的充分条件并研究了Runge??Kutta方法的渐近稳定性.最近,针对非线性延迟积分微分方程初值问题 y??(t)=ft,y(t),y(t-??), y(t)= (t), N[2] g(t,??,y(??))d,t-??tt#0,-???t?0,(2)这里??gt;0是实常数,f%[0,+amp;)?CN?CN?CN?CN,g%[0,+amp;)?R?CN?CN和 %[-??,0]?C是给定的连续映射,张诚坚和Vandewalle研究了其理论解和数值方法 (Runge??Kutta方法)的稳定性和渐近稳定性.其后,文[3,4]分别研究了单支方法与一般线性方法关于问题(2)的数值稳定性.但作为常见的数值方法,线性多步法关于非线性延迟积分微分方程的数值稳定性研究尚未见于文献,因此本文着重研究求解问题(2)的线性多步法,获得了方法渐近稳定的条件. 首先我们给出问题(2)理论解的稳定性条件,为此引入其扰动问题 z??(t)=ft,z(t),z(t-??), z(t)=!(t), * g(t,??,z(??))d,t-??tt#0,-???t?0,(3)收稿日期:2009??04??16 基金项目:国家自然科学基金资助项目,湖南省教育厅重点项目(09A093),湖南省自然科学基金资助项目(08JJ6002),,,,,:. N195这里!%[-??,0]?C是给定的连续映射,并设问题(2)和(3)分别有唯一真解y(t)和z(t). 对问题(2)和(3),我们设映射f,g满足 R(u1-u2,f(t,u1,v,w)-f(t,u2,v,w))???u1-u2?2, t#0,u1,u2,v,w!CN, ?f(t,u,v1,w1)-f(t,u,v2,w2)??#1?v1-v2?+#2?w1-w2?, t#0,u,v1,v2,w1,w2!CN, ?g(t,?,s1)-g(t,?,s2)??%?s1-s2?, t#0,t-?????t,s1,s2!CN, 范数,并将满足条件(4)??(6)的所有初值问题(2)记为问题类D(?,#1,#2,%). 注意本文所研究的延迟积分微分方程属于文[5]研究的Volterra泛函微分方程,因此将文[5]中定理2.2应用于延迟积分微分方程初值问题(2)立得 定理1??若问题(2)属于问题类D(?,#1,#2,%),那么当?+#1+#2%??lt;0时成立 t?+amp;(4)(5)(6)这里?,#1,#2和%是适度大小的实常数,(+,+)和?+?分别是复空间CN中的内积及相应的lim?y(t)-z(t)?=0.(7) 上式表征问题(2)的渐近稳定性. 将求解常微分方程的线性多步法 amp;(E)yn=h?(E)f(tn,yn) 用于求解延迟积分微分方程初值问题(2),有 amp;(E)yn=h?(E)f(tn,yn,yn-m,Yn),n=0,1,2,,, t(8)(9)这里h=??/m,m是一任意给定的正整数,tn=nh,E是位移算子:Eyn=yn+1,yn是y(tn)的逼 近,当n?0时,yn= (tn),Yn是 m-1 g(t,??,y(??))d??的逼近,由于A??稳定的线性多步法的阶nnt-??n不高于2阶,因此我们采用复