《微分方程数值解法》复习练习题试卷及答案.docxVIP

《微分方程数值解法》复习练习题试卷及答案

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

一、选择题

1.对于初值问题y=f(t,y),y(t0)=y0，数值解法本质上是将连续的微分方程问题转化为一系列关于自变量t的离散点tk上的（）问题。

A.代数方程

B.微分方程

C.积分方程

D.函数逼近问题

2.欧拉法(EulersMethod)的局部截断误差阶为（）。

A.O(1)

B.O(h)

C.O(h2)

D.O(h3)

3.改进欧拉法（HeunsMethod）也称为（）。

A.中点法

B.梯形法

C.改进的中点法

D.四阶龙格-库塔法

4.龙格-库塔法（Runge-KuttaMethods）的主要优点是（）。

A.计算量小

B.易于编程实现

C.对于刚性问题稳定性好

D.可以通过增加节点数来提高精度，无需修改公式

5.数值方法求得的解与精确解之间的差值，主要由（）引起。

A.舍入误差

B.截断误差

C.系统误差

D.随机误差

二、填空题

6.数值解法的一致性是指当步长h趋于零时，数值解y*_(n+1)收敛于精确解y(n+1)的极限性质。

7.数值方法的收敛阶是指局部截断误差的阶数。

8.求解常微分方程边值问题常用的方法有打靶法和有限差分法。

9.Adams外推法（Adams-Bashforth）是k步线性多步法，其局部截断误差为O(h^(k+1))。

10.若数值方法在每一步计算中仅使用前一步的值，则称为单步法。

三、计算题

11.(10分)用欧拉法求解初值问题y=t+y,y(0)=1，取步长h=0.1，计算y(0.3)的近似值，并计算该步的局部截断误差（用泰勒展开法推导）。

12.(10分)用改进欧拉法（HeunsMethod）求解初值问题y=t-y,y(0)=1，取步长h=0.1，计算y(0.1)的近似值。

13.(15分)写出求解初值问题y=y+t,y(0)=1的四阶龙格-库塔法（RK4）公式。

四、分析题

14.(10分)简述预测-校正方法（如Adams-Bashforth-Moulton）的基本思想和主要步骤，并说明其中预测公式和校正公式的作用。

15.(10分)解释什么是数值方法的稳定性。对于欧拉法，简要分析其稳定性条件（例如，对于y=λy的试验方程）。

16.(10分)比较欧拉法、改进欧拉法和四阶龙格-库塔法在求解初值问题时，在精度、计算量和实现复杂度方面的主要异同点。

试卷答案

一、选择题

1.A

2.B

3.B

4.D

5.B

二、填空题

6.一致性

7.收敛阶

8.打靶法，有限差分法

9.k+1

10.单步法

三、计算题

11.解：

*欧拉公式：y_(n+1)=y_n+h*f(t_n,y_n)

*本题f(t,y)=t+y,y(0)=1,h=0.1

*t_0=0,y_0=1

*t_1=t_0+h=0.1,y_1=y_0+h*f(t_0,y_0)=1+0.1*(0+1)=1.1

*t_2=t_1+h=0.2,y_2=y_1+h*f(t_1,y_1)=1.1+0.1*(0.1+1.1)=1.1+0.1*1.2=1.3

*t_3=t_2+h=0.3,y_3=y_2+h*f(t_2,y_2)=1.3+0.1*(0.2+1.3)=1.3+0.1*1.5=1.5

*所以y(0.3)的近似值为1.5。

*局部截断误差分析：

*y=t+y=y^(n+1)=y^n+h(t^n+y^n)+O(h2)

*精确解满足y^(n+1)=y^n+h(t^n+y^n)+y^(n+1)-y^n-h(t^n+y^n)=y^(n+1)-y^n-h(t^n+y^n)

*局部截断误差ε^(n+1)=y^(n+1)-y_(n+1)=y^(n+1)