数值方法：斜率就是契约
引入| English | 释义 | 拼音 |
|---|---|---|
| Euler method | 欧拉法 | ōu lā fǎ |
| trapezoidal rule/ˈtræpɪzɔɪdl ruːl/ | 梯形法则 | tī xíng fǎ zé |
| Simpson's rule | 辛普森法则 | xīn pǔ sēn fǎ zé |
| round-off error | 舍入误差 | shě rù wù chā |
| truncation error | 截断误差 | jié duàn wù chā |
When no closed form exists
- The real rig has a nonlinear damper and a clipped floor input. Everything after "no closed form" is this unit: approximating honestly, and knowing the cost.
- Three instruments: the Euler method 欧拉法 (order 1), the trapezoidal rule 梯形法则 (order 2), Simpson's rule 辛普森法则 (order 4).
The slopes are the contract
- Measured on y' = -y and on the integral of e^x over [0,1]: Euler errors 4.0e-2, 1.9e-2, 9.4e-3 (halve h, halve the error); trapezoid 8.9e-3, 2.2e-3, 5.6e-4 (quarter it); Simpson 3.7e-5, 2.3e-6, 1.5e-7 (sixteenth it - measured 15.9).
Which method, which order?
Route each instrument to its measured slope.
Euler errors at h = 0.2, 0.1: 4.0e-2 and 1.9e-2. The error at h = 0.05 (two significant figures)?
难度 2/5 · 今日首次完成可获 140 金币及少量额外奖励
Order 1: halve h, halve the error.
Trapezoid 1.04e-2 at h = 1/4. At h = 1/8?
难度 2/5 · 今日首次完成可获 140 金币及少量额外奖励
Order 2: quarter it.
Match method to order and behaviour.
难度 2/5 · 今日首次完成可获 140 金币及少量额外奖励
Three slopes and one floor.
Choose h with evidence
- Budget 1e-4 over one second: Euler needs h ~ 1e-4 (10,000 steps); a fourth-order stepper needs h ~ 0.1 (ten steps).
- Push Simpson to a million panels and the error stops falling: round-off error 舍入误差 accumulates ~1e-16 per step - the U-curve's floor is real.
Simpson 2.4e-5 at h = 1/4. At h = 1/16 (two significant figures)?
难度 2/5 · 今日首次完成可获 140 金币及少量额外奖励
Two halvings: 2.4e-5 / 256.
Pushing Simpson to a million panels, the error stops falling because:
难度 1/5 · 今日首次完成可获 100 金币及少量额外奖励
The U-curve floor: truncation falls, round-off rises.
Three numerics slips
- Reporting a result without its error estimate.
- Assuming smaller h always helps. It helps until round-off wins.
- Validating against the same physics you encode. Validate on a closed form first.
A 1e-4 budget over one second needs which h for Euler (error about 0.19h)?
难度 1/5 · 今日首次完成可获 100 金币及少量额外奖励
The budget and the slope fix h; a fourth-order method gets there in ten.
Validate a solver against a closed-form case before trusting it on the real physics.
难度 1/5 · 今日首次完成可获 100 金币及少量额外奖励
e^-t and pi are this unit’s witnesses.
A numerical result without an error estimate is finished work.
难度 1/5 · 今日首次完成可获 100 金币及少量额外奖励
It is a number; the estimate makes it a result.
Recap
- Euler O(h); trapezoid O(h^2); Simpson O(h^4) - measured, not memorised.
- Read h off the budget and the slope; say so in the report.
- Truncation error 截断误差 falls, round-off rises: the U-curve is the whole story.