跳到正文
计算物理

Euler法与局部误差

65 分钟

Euler法与局部误差

本节属于“常微分方程初值问题”。目标不是记住一串数值分析名词,而是能够从问题定义出发,独立写出计算对象、推导公式、完成一个带数字的基准、实现可复查算法,并识别会让结论失效的边界。学完后应能在不看正文的情况下解释 Eulerlocal truncationglobal errortime gridfirst order 之间的关系。 对《Euler法与局部误差》,本段必须回到这条专属证据复核:显式Euler以左端斜率推进,是理解局部与全局误差的基线,不能因简单就默认适合所有动力学。

1. 数学对象与物理含义

显式Euler以左端斜率推进,是理解局部与全局误差的基线,不能因简单就默认适合所有动力学。

先固定本节的输入、输出和不变量。输入必须写明单位、尺度、数据类型、边界或初值;输出不仅是最终数字,还包括能证明计算可信的残差、守恒量、收敛阶、置信区间或状态码。Euler 是本节建模的起点,first order 是判断结果能否外推的关键边界。若题目没有给齐这些条件,应先补规格,不得默认成最方便的情形。 对《Euler法与局部误差》,本段必须回到这条专属证据复核:由Taylor展开得 y_{n+1}=y_n+h f(t_n,y_n),单步缺陷O(h²),固定终点累积后全局误差O(h)。

一个合格的问题规格至少回答四件事:求什么量、允许多大误差、在哪个参数范围内成立、失败时怎样显式返回。这样才能区分“程序跑完”“离散方程算对”和“物理结论可信”这三件不同的事。

2. 公式链与推导

由Taylor展开得 y_{n+1}=y_n+h f(t_n,y_n),单步缺陷O(h²),固定终点累积后全局误差O(h)。

推导时按以下顺序书写,而不是直接抄最后一式:

  1. 从上一节定义的变量或算子写出第一行等式,并标出求和范围、固定变量和边界条件。
  2. 代入离散表示、基函数、概率分布或时间推进式,保留决定主误差的第一项。
  3. 整理得到可计算形式,检查矩阵维数、量纲、符号与极限。
  4. 用零输入、解析极限或守恒关系反查;反查不过时,不能进入大规模计算。 对《Euler法与局部误差》,本段必须回到这条专属证据复核:y′=-y,y0=1,h=.1,三步为.9,.81,.729;真值e^-.3=.740818,终点误差-.0118,h减半时误差约减半。

上述公式只在本节给出的模型条件下成立。若使用近似,报告中必须显式写出控制近似的小参数;若使用范数或概率,必须说明归一化约定。公式看似相同但单位、内积、边界或随机解释不同,得到的是不同问题。

3. 已算完的基准例题

y′=-y,y0=1,h=.1,三步为.9,.81,.729;真值e^-.3=.740818,终点误差-.0118,h减半时误差约减半。

复算时先预测符号和数量级,再逐项代入。答案需保留关键中间量,不能只给最后一位小数。随后把控制参数扩大或缩小两倍,比较结果究竟按线性、平方、平方根、指数还是其他规律变化;这一步能迅速暴露漏系数、错单位和错误归一化。

本例也是实现的最小回归测试。手算值、程序值、二者差异和容差来源应同时保存。若只能通过放宽容差让测试变绿,应先解释误差预算,而不是修改期待值迎合程序。

4. 可执行算法或实验

写纯一步函数;固定T比较h,h/2,h/4;与解析解算最大误差;拟合阶数;显式处理不能整除T的最后一步。

把这段过程实现成六个可观察阶段:输入校验、对象构造、核心更新、停止判据、独立验证、结果归档。每一阶段输出一个能定位首个错误的证据。随机算法还要记录生成器、种子、热化和有效样本量;迭代算法要同时记录残差和迭代状态;网格算法要保存至少三层离散尺度;动力学算法要保存守恒量或概率随时间的漂移。 对《Euler法与局部误差》,本段必须回到这条专属证据复核:写纯一步函数;固定T比较h,h/2,h/4;与解析解算最大误差;拟合阶数;显式处理不能整除T的最后一步。

验收不能只看一张平滑图。至少运行正常、边界、退化、错误四类输入,并让错误输入得到明确失败状态。性能比较必须在相同问题、相同精度目标和相同硬件口径下进行。

5. 误差预算与适用边界

局部缺陷不等于每步真实误差,已有误差还会被稳定性放大;必须连同终点与范数报告。

误差表至少分开输入/测量误差、模型误差、离散或截断误差、迭代误差、舍入误差和抽样误差。只报告总误差会掩盖主导来源,也无法判断增加网格、样本或迭代是否值得。若两类误差相关,还应说明协方差或共同随机数的处理方式。

结论必须附适用范围:参数区间、稳定条件、边界类型、光滑性、独立性或遍历性假设。离开这些条件后,应重新验证,而不是沿用本节图表。

6. 最短失败案例

浮点循环可能多走一步;用整数步数并输出实际时间,避免把时间错位当数值误差。

请把这个失败案例真的实现出来,记录错误首次出现在哪一步:输入规格、算子组装、时间推进、停止条件还是统计解释。修复后同时保留失败测试,防止以后回归。仅写“可能不准确”不算分析;必须指出会多算或少算什么、哪个等式失效、残差或守恒量怎样变化,以及应该换用什么算法、边界或表示。 对《Euler法与局部误差》,本段必须回到这条专属证据复核:局部缺陷不等于每步真实误差,已有误差还会被稳定性放大;必须连同终点与范数报告。

7. 闭卷自检

  1. 用一句话区分 Eulerfirst order 在本节中的职责。
  2. 从定义重建上面的公式链,并在每个等号旁写成立条件。
  3. 不看答案重算数值例题,解释数量级和参数缩放。
  4. 写出正常、边界、退化、错误四组测试的预期状态。
  5. 说明误差预算中哪一项当前占主导,以及怎样用一次额外实验验证判断。 对《Euler法与局部误差》,本段必须回到这条专属证据复核:浮点循环可能多走一步;用整数步数并输出实际时间,避免把时间错位当数值误差。

五项中任何一项答不出,都应回到相应章节重算;继续浏览并不能替代真正掌握。

Practice

本课练习

6

先独立作答再提交;编程题会在隔离沙箱中真实编译、运行并对拍。

1单选:Euler法与局部误差 4

Euler法与局部误差计算实验中哪条工作流可信?

登录 后答题可以领小红花
2多选:Euler法与局部误差 4

复核Euler法与局部误差需要哪些证据?

多选题:必须选全正确项,漏选或多选均不得分。

登录 后答题可以领小红花
3计算:Euler法与局部误差 4

Euler法与局部误差基准包含8个网格块,每块8个自由度,总自由度是多少?

登录 后答题可以领小红花
4实验设计:Euler法与局部误差 4

为Euler法与局部误差写出模型、算法、基准、误差预算、四类测试和可复现记录。

【评分量表】模型2分;算法3分;验证3分;边界与复现2分。

登录 后答题可以领小红花
5u05独立题03:Euler法与局部误差 5

Euler法与局部误差网格有2块、每块7自由度,总自由度?

登录 后答题可以领小红花
6midterm独立题07:Euler法与局部误差 5

Euler法与局部误差网格有6块、每块11自由度,总自由度?

登录 后答题可以领小红花