∑

数学知识体系

Observatory Archive of Mathematics
⌕2026/8/31
概念

误差估计

先验误差估计在求解前由网格尺寸与阶数给出误差界,后验误差估计由残差驱动自适应网格加密,把计算资源集中到误差大的区域。

所属主题:计算数学 ↗
阅读路径

参考可汗学院 Get ready 机制:先修概念 → 当前概念 → 进阶概念,✓ 表示已读。

当前概念

误差估计计算数学

01定义

有限元的误差理论分两层:先验估计在求解前承诺 ∣u−uh∣≤Chp|u-u_h|\leq Ch^p(网格越细、阶数越高,误差越小);后验估计在求解后由残差定位误差大的区域,驱动自适应网格加密——把算力精准投放到需要的地方。
h → 0 误差 ~ h^p(斜率 = p) 残差大 ⇒ 局部加密 先验:‖u−uₕ‖ ≤ C hᵖ |u|(解前可知收敛阶) 后验:η_K ≈ 单元 K 的误差贡献(解后计算)
先验给收敛率,后验指明加密位置

02核心要点

01

先验估计

分片 pp 次多项式、网格尺寸 hh:H1H^1 误差 O(hp)O(h^p)、L2L^2 误差 O(hp+1)O(h^{p+1})(对偶论证升一阶)——收敛阶只取决于离散阶数,与具体问题无关。

02

后验估计与自适应

残差型(单元残差 + 边跳跃)与恢复型(ZZ 梯度恢复)两大流派:可靠上界 + 效率下界保证指标可信;加密-求解-估计循环直至目标精度。

03

hp 自适应

光滑区升高阶数(pp 加密,指数收敛),奇异性区细化网格(hh 加密):hp 策略在裂纹、角点等奇性问题中比纯 h 或纯 p 高效一个量级。

03关键公式

∥u−uh∥H1≤Chp ∣u∣Hp+1,ηK≳∥u−uh∥K\|u-u_h\|_{H^1}\leq Ch^{p}\,|u|_{H^{p+1}},\quad \eta_K\gtrsim\|u-u_h\|_K

04历史沿革

Céa 1971 年奠基有限元逼近论;Ainsworth 与 Oden 1990 年代系统化后验估计理论;自适应有限元如今是 PDE 求解器的标准功能。

05应用与延伸

裂纹扩展仿真(应力奇异性自动加密)、气动外形的激波捕捉、电子器件热分析、心脏电生理的多尺度建模。

06交互演示

泰勒截断误差:交错级数逐项逼近增多项数:截断误差 ≤ 下一项的绝对值

07相关概念