求解二阶常微分方程的并行块方法.pdfVIP

上海师范大学硕士学位论文 第一章引言 第一章弟一早 亏I石引言 1．1数值方法及收敛性理论 考虑二阶常微分方程初值问题 可Ⅳ(z)=，(z，秒)，zXo， (1．1．1) v(：ro)=Vo，v(zo)=％， 其中，是连续的实值函数。二阶常微分方程在天体力学、理论物理与化学以及电子学等诸 多领域具有十分广泛的应用。近几十年来，关于求解二阶周期性常微分方程初值问题引 起了许多研究者的极大兴趣，已经构造了许多求解二阶常微分方程的数值方法。一般来 说，二阶常微分方程初值问题的数值方法可以分为两类：一类是通过引进新的变量将其 转化为一阶常微分方程组再进行求解，但是方程维数是原来的2倍；另一类是构造直接计 算的数值格式，保持方程维数不变。 求解二阶常微分方程初值问题(1．1．1)的线性多步法可以表示为 后 七 (1．1．2) ∑哟‰句=h2∑岛，(zn吖，蜘钾)，n=0，1，2，…， j=o j=o 解。定义线性多步法(1．1．2)的第一和第二特征多项式为 = = “O ％ ．r 州O 房 膏∑间 七∑细 假定方法(1．1．2)满足下列条件： l岛l≠o； (1)哟，岛∈R，并且Q知≠0，I,10l+l岛J≠o，∑冬o (2)p(()和矿(()是不可约的； (3)相容条件： p(1)=∥(1)=0，∥(1)=2盯(1)； (4)稳定性条件：设白是p(()的根，则l白I≤1，且若f白I=1，则其重数不超过2． 第一章引言 上海师范大学硕士学位论文 Henrici【22]i正明了满足上述四个条件的线性多步法(1．1．2)是收敛的。 典型的单步方法是如下形式的s-级Runge-Kutta-NystrSm方法： i=1，2，…，s， K=Yn+九c‘如+h2∑；：1aof(xn+c／h，K)， b,f(xn+c／h，K)， 必+1=必+九∑：l Yn+1=‰+嗽+h2∑：l瓦，(％+c‘^，K)． Kutta-NystrSm方法的阶条件。 就我们所知，求解二阶常微分方程的块方法的研究工作并不多。Chakravarti 和Wodand【5】首先提出了一类自启动的块方法，并得到k维块方法具有k+1阶精度。 T。， 有关二阶常微分方程计算方法的其它研究结果可以参见综述性论文【42】。 1．2数值方法的p稳定性 考虑试验问题 矿(z)=一A2y(z)，入∈R， (1．2．1) 其通解为 I+Id2I≠0． (1．2．2) y(x)=die‘蛔+dae一‘蛔，Idl 由于方程(1．2．1)的解具有周期性，因而希望数值