θ-方法求解非线性延迟微分方程的渐近稳定性分析.pdfVIP

第 27卷 哈尔滨师范大学 自然科学学报 Vo1．27，No．32011 第3期 NATURALSCIENCESJOURNALOFHARBINNORMALUNIVERSITY l5}一方法求解非线性延迟微分方程的 渐近稳定性分析 蔡承文 ，陈 勇 ，刘 兴。 (1．南京工业职业技术学院；2．黑龙江省特种设备检验研究院；3．哈尔滨工业大学) 【摘要】 主要研究非线性延迟微分方程的数值稳定性．文中给出一个充分条 件，使得在该条件下，由 一方法所求的数值解可以保持解析解的渐近稳定性． 关键词：延迟微分方程 ；0一方法；渐近稳定性 1 龙阁库塔方法求解非线性延迟 稳分方程 0 弓l言 应用 Runge—Kutta方法于方程(1)，可得： 主要考虑如下形式的方程 ： r )： ，y()，y(一 ))，≥to, (1) n =Y+h∑ t+ ^，，)，i=1， 』 1 ty(t)= (t)，t≤to … - 其中7是一个正常数，(t)是一个连续函数 ， Y川=Y+h∑bf(t+cih，，)，(3) ： [t。，+∞)×R ×R 一R是一个给定的函数 且满足 其中 表示Y(t+cih—下)的一种近似．取h= l一 2√’(t，1，Y)一／(t，2，Y)≤ ， 其中m≥1为一给定的整数．则 可以被定 Otl『l一 2I】2，Vt≥t0，1，2，Y∈RⅣ， 义成 因此 (3)变为 lIf(t，，Y。)一／t，，Y2)ll≤卢llY 一Y2lJ， Vt≥ to，Yl，Y2， ∈R ， =Y+h∑0 t+cjh，，一)，i=1， J= L 其中 ·，·表示R 上的一个内积，lI·ll表 … ， ， 示相应的范数 ，Ol和 均为实数．为了研究(1)的 稳定性，同时考虑下面的方程： Y =Y+h∑bf(t十cih，，一)．(4) )=八t,z()，(一 ))，≥to, (2) 同理可得，方程(2)对应的数值解为： ty()= (￡)，t≤t0， z：= +h∑0 t+cjh，，一)，i=1， 其中lfr(t)是一个连续函数． J 1 … 定理 1 (参见文献 [2])若 ≤一0c，则 ， ， )一z(训 }