解二阶微分方程的Runge—Kutta—Nystriim 方法的正则性.pdfVIP

第19卷 第2期 黑 龙 江 大 学 自 然 科 学 学 报 vo1．19No．2 2002年 6月 JOURNAL OF NATURAL SCIENCE OF HEILONGJIANG UNIVERSITY June．2002 文章编号：1001—7011(2002)02—0007—04 解二阶微分方程的Runge—Kutta—Nystriim 方法的正则性 邹 巾英， 丁效华，金承 日 (哈尔滨工业大学 威海分校，山东 威海 264209) 摘 要：研究Runge—Kutta—Nystrt~m 方法平衡点与二阶常微分方程 系统平衡态一致 的条 件，该方法称为正则的。给 出了一些推断正则性和强正则性的准则并给出一个例子。 关键词：Runge—Kutta—NystrSm方法；二阶初值 问题 ；正则性 中图分类号：O241．81 文献标识码：A 1 引 言 本文的 目的是研究二阶常微分方程系统的平衡态与用 Runge—Kutta—NystriSm 方法 (RKN方法)求解 二阶常微分方程系统的平衡态相一致 的条件 。即研究用数值方法求解系统时，是否会出现伪平衡的条件。不可 能出现伪平衡 (即系统的平衡态与方法的平衡态相一致)的方法称为是正则的。对于一阶常微分方程系统，RK 方法、线性多步法和一般线形多步法的正则性 已被众多的文献 q所研究。 设 ． 厂是一连续函数，考虑求解二次初值 问题： 【 ， y()=Yo y(to)=Y。 的 Runge—Kutta—Nystr6m 方法 ： fiY=y+hy+h：U,s() ‘{ly ：)}+ (yJ) (1．2) lYf= +qhy+h2S f(y，) i=1，2，’。，‘． ≥3) 其 中h0，Y、Y 分别逼近于微分方程的真解在 处的值 y()及其导数值 y(t3，Y逼近于 y + )， 1，2，…，S．与Runge—Kutta方法类似 [2]，Runge—Kutta—Nystrbm方法-nI~A表示为下例更紧凑的格式： (1．3) 二阶常微分方程的平衡 点是指满足 (y)=0的点，方法 (1．2)达到平衡 是指计算步 (y，Y)一 (Y ， ) 有 Y+1=y，Y +1=Y =0。 收稿 日期 ：200l一09—27 基金项 目：国家 自然科学基金资助项 目；哈尔滨工业大学校 内基金资助项 目(2OOO一007)；哈尔滨工业大学威海分校校 内基金资助 项 目(2000-031) 作者简介：邹 巾英 (1966一)，女，山东威海人，哈尔滨工业大学威海分校讲师，主要研究方向：常微分方程数值解． 黑 龙 江 大 学 自 然 科 学 学 报 第 l9卷 考虑系统 f善 ()o { f(Yo (1．4) lYiFh2Z口．』()， l，2，…，s． 定义集合S={yeR： )=o)和集合S= ∈Rq：(3，， 一， )满足 (1．4)}。 若Y是微分方程的平衡点，令 多=y， 。=y， l，…，S，可知 (1．4)成立，即ScS．但是ScS是否成立， 即方法达到平衡时，微分方程是否也已达到平衡