### Python 实现二级欧拉法求解常微分方程 二级欧拉法是对经典欧拉法的一种改进，它通过引入更多的中间点来提高精度。这种方法也被称为改进欧拉法或梯形法则，在每一步中利用预测值和校正值的平均值作为最终的结果。 以下是基于 Python 的实现代码： ```python import numpy as np import matplotlib.pyplot as plt def improved_euler_method(f, y0, t_span, h): """ 使用改进欧拉法（二级欧拉法）求数值解 参数: f (function): 微分方程 dy/dt = f(t, y) y0 (float): 初始条件 y(t0) = y0 t_span (tuple): 时间区间 (t_start, t_end) h (float): 步长 返回: ts (list): 时间节点列表 ys (list): 对应时间节点上的数值解列表 """ t_start, t_end = t_span ts = np.arange(t_start, t_end + h, h) ys = [y0] for i in range(len(ts) - 1): t_i = ts[i] y_i = ys[-1] # 预测步 k1 = f(t_i, y_i) y_pred = y_i + h * k1 # 校正步 k2 = f(ts[i + 1], y_pred) y_new = y_i + h * (k1 + k2) / 2 ys.append(y_new) return ts.tolist(), ys # 定义微分方程 dy/dt = f(t, y) def f(t, y): return -2 * t * y # 示例：dy/dt = -2ty # 设置参数 y0 = 1 # 初始条件 y(0) = 1 t_span = (0, 2) # 时间范围 h = 0.1 # 步长 # 调用改进欧拉法 ts, ys = improved_euler_method(f, y0, t_span, h) # 绘图展示结果 plt.plot(ts, ys, label="Improved Euler Method", marker='o') plt.xlabel("Time (t)") plt.ylabel("Solution (y)") plt.title("Numerical Solution of ODE using Improved Euler Method") plt.legend() plt.grid(True) plt.show() ``` #### 解析 上述代码实现了改进欧拉法的核心逻辑[^1]。具体来说： - 函数 `improved_euler_method` 接受微分方程定义函数 `f`、初始条件 `y0`、时间跨度 `t_span` 和步长 `h`。 - 在每次迭代中，先执行 **预测步** 计算初步估计值 \( k_1 \)，再通过 **校正步** 结合下一步的时间点计算更精确的斜率 \( k_2 \)[^2]。 - 最终取两者的加权平均值作为当前步的新值[^3]。 此方法相较于传统欧拉法显著提高了准确性，但仍属于显式单步法范畴。 --- 相关问题