基于辛euler法的数值算法研究.docxVIP

基于辛euler法的数值算法研究

0改进euter算法和改进euter算法组合而成算法

采用埃尔伯算法法研究了一阶微分方程的数值算法、改进的埃尔伯算法、隐藏的算法、gauss中等算法和二次方程的数值规律，结果表明，在一组方程算法中，用于计算二次微分方程的数值算法的增强算法由eil伯算法和隐子算法组成。因此，在同一研究中的新算法和正如所示的算法方面，作者提出了一个新算法，该算法由eil算法和gauss中点算法或eil算法组成，这与当前的欧算法进行了比较，并发现了该算法的缺点。

14常微分方程初值问题

设f(t,u)在区域G:0≤t≤T,|u|＜∞0≤t≤Τ,|u|＜∞上连续,求u=u(t)满足

其中u0是给定的初值,这就是一阶常微分方程的初值问题.为使问题(1)的解存在、唯一且连续依赖初值即初值u0问题(1)适定,还必须对右端f(t,u)加适当限制,通常要求f关于u满足Lipschitz条件:即存在常数L,使

对所有t∈[0,T]和u1,u2∈(-∞,∞)成立.在常微分方程数值解法中,对于常微分方程初值问题的数值解法都要满足上述的Lipschitz条件.

1.1t0的调值

最简单的数值解法是Euler法.将区间[0,T]作N等分,小区间的长度h=T/N称为步长,点列tn=nh(n=0,1,2,…,N)称为节点,t0=0由已知初值u(t0)=u0,可算出u(t)在t=t0的导数值u′(t0)=f(t0,u(t0))=f(t0,u0).利用Taylor展式

其中ζ∈(t0,t1),并略去二阶小量R0,得u1=u0+hf(t0,u0).u1就是u(t1)的近似值.利用u1又可算出u2,如此下去可算出u在所有节点上的值,一般递推公式为

un+1=un+hf(tn,un),n=0,1,…,N-1.(5)

这就是Euler法.

1.2改进euter法式

下面对Euler法运用数值积分法推导.将问题(1)写成等价的积分形式

u(t)=u0+∫0tf(τ,u(τ))dτ,t0=0,(6)

特别地

u(t1)=u0+∫0t1f(τ,u(τ))dτ,t0=0.

用左矩形公式近似右端积分,并用u1替代u(t1),即得u1=u0+hf(t0,u0),这就是Euler法(式(5)).我们也可用梯形公式近似上述积分,仍用u1替代u(t1),得

u1=u0+h2[f(t0,u0)+f(t1,u1)].u1=u0+h2[f(t0,u0)+f(t1,u1)].

一般而言,

un+1=un+h2[f(tn,un)+f(tn+1,un+1)],n=0,1,?,N?1un+1=un+h2[f(tn,un)+f(tn+1,un+1)],n=0,1,?,Ν-1(7)

称之为改进Euler法.

1.3有数值计算的迭代格式

隐式Euler法同Euler法的数值解法相似,将区间[0,T]作N等分,小区间的长度h=T/N称为步长,点列tn=nh(n=0,1,2,…,N)称为节点,t0=0是已知初值u(t0)=u0.在Euler法中,数值计算的迭代格式中运用初始值点(t0,u0)进行格式迭代计算.在改进Euler法中,数值计算的迭代格式中运用初始值点(t0,u0)和所求值点(t1,u1)进行格式迭代计算.根据上面的Euler法和改进Euler法的格式推导过程可以进行相似的演变.将Euler法中迭代点(t0,u0)用所求值点(t1,u1)代替,因此有数值计算的迭代格式

un+1=un+hf(tn+1,un+1),n=0,1,…,N-1.(8)

1.4相邻2点的中点法gauss中点法

对于问题(1)上面3种数值算法都可求解,但是,这样的求解结果并没有得到微分方程数值解研究者的满足.我们需要得到更好的数值算法.Gauss在对上述的3种数值算法反复研究后,再次改造Euler算法的算式格式.将运用迭代格式中相邻2点的中点作为算式的迭代值,即取t=0初始点(t0,u0)与后一点(t1,u1),将这2点的中点((t0+t1)/2,(u0+u1)/2)作为Euler法格式中的迭代点.即,这个中点迭代格式为

u1=u0+hf(t0+t12,u0+u12).u1=u0+hf(t0+t12,u0+u12).

在一般情况下,将区间[0,T]作N等分,小区间的长度h=T/N称为步长,点列tn=nh(n=0,1,2,…,N)称为节点,迭代格式为

un+1=un+hf(tn+tn+12,un+un+12),n=0,1,?,N?1un+1=un+