微分方程数值解习题试卷及答案.docxVIP

微分方程数值解习题试卷及答案

考试时间：______分钟总分：______分姓名：______

一、

简述微分方程数值解的必要性和局限性。

二、

已知初值问题：$y=-y+1$,$y(0)=0$。

1.用欧拉法（步长$h=0.1$）计算$y(0.3)$的近似值，保留四位小数。

2.写出欧拉法求该初值问题的近似解的一般公式。

3.估计用步长$h=0.1$的欧拉法计算$y(0.3)$时，局部截断误差的主项。

三、

比较欧拉法、改进欧拉法（中点法）和四阶龙格-库塔法（RK4）在求解初值问题$y=f(t,y)$时，关于收敛速度和计算量的优缺点。

四、

判断下列初值问题在所给步长下，用欧拉法求解是否稳定：

1.$y=y$,$y(0)=1$,$h=0.5$

2.$y=-2y$,$y(0)=1$,$h=0.1$

五、

已知初值问题：$y=t+y$,$y(0)=1$。

1.用四阶龙格-库塔法（RK4）求$y(0.5)$的近似值，取步长$h=0.1$，保留四位小数。

2.写出该初值问题精确解的表达式，并计算$y(0.5)$的精确值，与RK4的结果比较。

3.估计用步长$h=0.1$的RK4方法计算$y(0.5)$时的全局截断误差主项。

六、

解释什么是数值方法的收敛性？证明欧拉法是收敛的。

七、

简述线性多步法（如Adams方法）与单步法（如欧拉法、龙格-库塔法）在计算效率和稳定性方面的主要区别。

八、

对于初值问题$y=f(t,y)$，解释如何通过稳定性区域的概念来选择合适的数值方法和步长。

九、

设计一个简单的算法（用伪代码或流程图描述，文字说明辅助），描述如何使用四阶龙格-库塔法（RK4）求解一个给定的二阶常微分方程初值问题。例如，可以将二阶方程$y=g(t,y,y)$化为两个一阶方程的一阶组，然后说明如何用RK4求解这个一阶组。

试卷答案

一、

微分方程在实际问题中往往难以找到解析解，或者解析解过于复杂。数值解方法可以将微分方程转化为离散的代数方程组，通过计算机进行迭代计算，得到近似解。这使得我们能够处理更复杂、更实际的模型，获得可计算的、有意义的结果。然而，数值解依赖于离散化的过程，引入了截断误差和舍入误差，且解的精度受步长限制，计算量也通常比解析解大。此外，数值解的稳定性问题也需要考虑。

二、

1.欧拉法公式：$y_{n+1}=y_n+hf(t_n,y_n)$。

令$f(t,y)=-y+1$，$t_0=0$,$y_0=0$,$h=0.1$。

$y_1=y_0+0.1f(0,0)=0+0.1(-0+1)=0.1$。

$y_2=y_1+0.1f(0.1,y_1)=0.1+0.1(-0.1+1)=0.1+0.1\times0.9=0.19$。

$y_3=y_2+0.1f(0.2,y_2)=0.19+0.1(-0.19+1)=0.19+0.1\times0.81=0.19+0.081=0.271$。

$y_4=y_3+0.1f(0.3,y_3)=0.271+0.1(-0.271+1)=0.271+0.1\times0.729=0.271+0.0729=0.3439$。

近似值$y(0.3)\approx0.3439$。

2.一般公式：$y_{n+1}=y_n+h(-y_n+1)=y_n(1-h)+h$。

3.欧拉法的局部截断误差主项为$O(h^2)$。根据$y=-y+1$，其精确解为$y(t)=1-e^{-t}$。代入欧拉公式：

$y(t+h)\approxy(t)+h(-y(t)+1)=y(t)-hy(t)+h=(1-h)y(t)+h$。

精确解：$y(t+h)=1-e^{-(t+h)}=1-e^{-t}e^{-h}=y(t)e^{-h}$。

误差：$E_{n+1}=y(t+h)-y_{n+1}=y(t)e^{-h}-(1-h)y(t)-h$。

$E_{n+1}=y(t)(e^{-h}-1+h)-h$。

当$h\to0$，$e^{-h}\approx1-h+\frac{h^2}{2}$。

$E_{n+1}\app