跳转至

怎样对 Gamma、Beta 积分作带状态的可靠近似?

先备知识

熟悉端点比较界、复合 Simpson 误差界与状态式算法接口。

学习目标

能读取误差字段与三种状态,不把数值稳定误写成数学证书。

牵引问题

截断误差和有限区间求积误差怎样合成一个可靠上界?

探索与猜想

先解析控制被删去的端点,再由已证明的四阶导数界认证有限主部。

概念与理论

误差合同

程序把反常端点界记录为 endpoint_error_bound,把 Simpson 界记录为 quadrature_error_bound,两者相加为 total_error_bound。外部提供的四阶导数界是数学输入,程序不从采样中证明它。

from mathbook_examples.parametric_integrals import beta_integral, gamma_integral

met = gamma_integral(2.0, 1e-5, 4096, fourth_derivative_bound=10.0)
exhausted = gamma_integral(2.0, 1e-12, 2, fourth_derivative_bound=10.0)
diagnostic = beta_integral(2.0, 3.0, 1e-5, 512)

三种状态

  • target_met:端点界与条件式求积界之和不超过目标。
  • budget_exhausted:仍有有限有效总界,但面板预算不足。
  • uncertified:缺少四阶导数界,只报告诊断值,不返回总误差证书。

uncertified 不等于数值错误,budget_exhausted 也不等于积分发散。

例题与迁移

例 1:证书达标

只有状态达标且总界不超过容差时才可声称满足目标。

例 2:缺失数学输入

诊断值可用于探索,但网格最大值不能冒充全区间导数界。

即时检验与回望

即时检验 1

谁证明四阶导数界?

答案

调用者;函数只把它作为条件式证书输入。

即时检验 2

预算耗尽时还有界吗?

答案

有,返回有限总界,但该界大于目标容差。

常见误区与后续

  • 采样最大值不是全区间导数上界证书。
  • 浮点舍入和外部函数求值误差不包含在解析求积界内。

习题与答案

习题 1

端点误差字段是什么?

答案

endpoint_error_bound

习题 2

总误差字段是什么?

答案

total_error_bound

习题 3

达标状态是什么?

答案

target_met

习题 4

预算不足状态是什么?

答案

budget_exhausted

习题 5

缺少导数界状态是什么?

答案

uncertified

习题 6

未认证是否表示积分发散?

答案

不表示,只是缺少求积误差证书。

习题 7

预算耗尽是否没有任何界?

答案

不是,仍返回当前有限总界。

习题 8

四阶导数界可由网格猜测吗?

答案

不可作为证书。

习题 9

总界怎样组成?

答案

端点界加有限主部求积界。

习题 10

可靠报告至少包含什么?

答案

数值、截断区间、误差字段、评估次数和状态。

答案

还应说明导数界来自外部数学证明。