数值分析实验六 常微分方程的数值解2.docxVIP

PAGE

PAGE1

数值分析实验六 常微分方程数值解

组号班级 学号 姓名 分数

一：实验目的

比较欧拉公式,隐式欧拉公式和改进欧拉公式在解微分方程时的结果，观察其精度的高低。

二：实验内容及基本知识介绍

科学技术中常常需要求解常微分方程的定解问题，本次将研究一些解一阶方程的初值问题的方法。

{y?f(x,y)

对于一阶方程 y(x

0

)?y

0

（1.1）

我们知道，在xy平面上，微分方程（1.1）的解y?y(x)称作它的积分曲线。积分曲线上

一点(x,y)的切线斜率等于函数f(x,y)的值。如果按函数f(x,y)在xy平面上建立一个方向场，那么，积分曲线上的每一点的切线方向与方向场在该点的方向相一致。

基于上述几何解释，我们从初始点 p

0

(x,y

0 0

)出发，先依方向场在该点的方向推进到

x?x

1

上一点p

1

，然后再从p

1

依方向场的方向推进到x?x

2

上一点p

2

，循此前进做出一条

折线p

0

pp...。

1 2

—般的，设已做出该折线的顶点p

n

，过p

n

(x,y

n n

)依方向场的方向再推进到

p

n?1

(x

n?1

,y

n?1

)，显然两个顶点p

n

，p

n?1

的坐标有关系

y ?y

n?1 n?f(x,y)，

x ?x n nn?1 n

即 y ?y ?hf(x,y) （2.1）

n?1 n n n

这就是著名的欧拉公式。

为了分析计算公式的精度，通常可用泰勒展开将y(x

n?1

)在x

n

处展开，则有

y(x

n?1

)?y(x

n

?y(x

n

h)

)?y(x

n

)h?h2y(?

2

) ，?

n

n?1

?(x

n

,x )。

n?1

在y ?y(x

n n

)的前提下，

f(x,y

n n

)?f(x

n

,y(x

n

))?y(x)。

n

于是可得欧拉法（2.1）的公式误差

y(x

n?1

)?y

?h2

n?1 2

y(?)?h2

n 2

y(x

n

)， （2.3）

称为此方法的局部截断误差。

如果对方程（1.1）从x

n

到x

n?1

积分，得

y(x

n?1

)?yx(

n

?)?xn?1x

n

f t(y,t (dt))

（2.4）

右端积分用左矩形公式hf(x

n

,y(x

n

))近似，再以y

n

代替y(x

n

)，y

n?1

代替y(x

n?1

)也得到

（2.1），局部截断误差也是（2.3）。

如果在（2.4）中右端积分用右矩形公式hf(x ,y(x ))近似，

n?1 n?1

则得到另一个公式

y ?y ?hf(x ,y )， （2.5）

n?1 n n?1 n?1

称为隐式欧拉公式。

隐式欧拉公式与欧拉公式有着本质的区别，后者是关于y 的一个直接计算公式，这类

n?1

公式称为显式的；然而公式（2.5）的右端含有未知的y ，它实际上是关于y 的一个函

n?1 n?1

数方程，这类公式称为隐式的。

相比较上述方法较复杂，因此，我们先用欧拉公式求得一个初步的近似值y ，称之为

n?1

预测值，预测值y 的精度可能很差，再用梯形公式将它校正一次，这个结果称为校正值，

n?1

而这样建立的预测-校正系统通常称为改进欧拉公式：

_nyn?1?y

_

n

?hf?x,y?

n n

h? _ ?

y ?y? ?f(x,y

)?f(x ,y

)? （2.13）

n?1 n 2? n n

n?1

n?1?

三：实验问题及方法、步骤

问题：取h?0.1，分别用欧拉公式,隐式欧拉公式和改进欧拉公式解微分方程

1?2x{y?y?2

1?2x

y(0)?1 0?x?1，并比较精度（解析解y? ）。

解：h?0.1，x

0

?0，x

1

?0.1，x

2

?0.2，x

3

?0.3，x

4

?0.4，x

5

?0.5,x

6

?0.6,

x ?0.7, x

7 8

?0.8,x

9

?0.9,x

10

?1根据式（2.1）得欧拉公式

y

n?1

?y ?h?y

n n

2x

n

/y?, n=0,1,2,3…

n

根据式（2.5）得隐式欧拉公式

y

n?1

?y ?