比例方程定步长方法的稳定性.pdfVIP

比例方程定步长方法的稳定性 杨占文 哈尔滨工业大学数学系 摘要 本文讨论了线性常系数比例延迟微分方程定步长方法的数值稳定性。应用离 散形式的Razumikhin定理，给出了两种口一方法的数值稳定区域，并证明了z． Jackiewicz的猜想是正确的。 1 5l。蓄 本文考虑比例延迟微分方程 佗o P璺2儿删，工(蛳 (1．1)’ 。 卜(0)=xo 比例方程系统可以用来模拟很多现象，例如：电子力学，非线性动力系统等等。 本文主要考虑线性常系数比例方程 (12) {工黧2删地(qtx垃。 。 。 【J(o)=‰ 其中a，b，X。是实常数且q∈(0,1)。 比例方程(1．2)的数值解方面最大的困难是有限的计算机内存。解决存储问题 的方法有两种。一是将(1．2)变换成常时滞变系数方程，另一种方法是直接对(1．2) 应用几何步长方法。应用几何步长数值方法的另一个原因是当应用定步长数值方 法时求解(1．2)时，得到的差分方程不是定阶的。最近一些年里，只有少数几位学 者讨论过定步长数值方法用于比例方程的渐近稳定性【lj。 但事实上，这两种方法可以看成是[2]的两种特殊的步长选择策略。于是， 这两种避免存储问题的方法的基本思想是应用一种变步长方法。注意到集合 s(t)={L：n∈z+使得f。∈【qt，以}中的点在计算过程中是必须存储的，雨为了在 很大的时间f保持计算精度，这些点的数量随着t-4m是趋于无穷的，因丽存储 问题仍然是个问题。这样，从本质上讲，这两种方法都没有根本地解决比倒力枉 数值计算的存储问题。这样本文将讨论比例方程定步长方法的渐近稳定性。 2 离敬系统的Razumikhin型定理 考虑下面的离散系统 {x(n十1)=厂(H，z(玎)，x(qcl(H)，t～，qc，()))，n≥0，。n 协)=‰ u‘¨ 5l 法可以定义(2．1)的平凡解的稳定性和渐近稳定性。 ”(1x1)蔓V(n，功≤v(Ixl)，n∈Z+ 使得 (2．2) ≤一w(1x(以)I) 若V(m，x(m))p(V(n，x(n))对所有的m∈以={【qH】，…，M)成立。 证明：(i)对任意的占0，取叩使得V(叼)≤“(占)。令1‰降玎，那么 v(o，x(O))s v(m，x(埘))≤v(，7)当0≤ms五 V(n+l，x(n+1))V(77) 于是△矿(-)0。另一方面当ln∈乇时有 x(n)1)≤0，矛盾。 由(2．2)得△矿(-)≤～w(I x(n)巨H且 (ii)设H0，吼O满足v(玑)=“(日)。由(i)得当1‰峰吼时有I 足。岛可芸嚣笔的整数。定义怫=高]，碍=[-面三豸亭丽] 52 』，={州，+1，…，／t。一1。则可以归纳得出 V(n