不完全贝塔函数与正则化不完全贝塔函数表计算器

输入计算

a（第一形状参数）

b（第二形状参数）

x 的初始值

步长

重复次数

结果

完全贝塔函数 B(a,b)

0.333333

last row x = 1: Bx = 0.333333, Ix = 1

x	Bx(a,b)	Ix(a,b)
0	0	0
0.02	0.01960267	0.058808
0.04	0.03842133	0.115264
0.06	0.056472	0.169416
0.08	0.07377067	0.221312
0.1	0.09033333	0.271
0.12	0.106176	0.318528
0.14	0.12131467	0.363944
0.16	0.13576533	0.407296
0.18	0.149544	0.448632
0.2	0.16266667	0.488
0.22	0.17514933	0.525448
0.24	0.187008	0.561024
0.26	0.19825867	0.594776
0.28	0.20891733	0.626752
0.3	0.219	0.657
0.32	0.22852267	0.685568
0.34	0.23750133	0.712504
0.36	0.245952	0.737856
0.38	0.25389067	0.761672
0.4	0.26133333	0.784
0.42	0.268296	0.804888
0.44	0.27479467	0.824384
0.46	0.28084533	0.842536
0.48	0.286464	0.859392
0.5	0.29166667	0.875
0.52	0.29646933	0.889408
0.54	0.300888	0.902664
0.56	0.30493867	0.914816
0.58	0.30863733	0.925912
0.6	0.312	0.936
0.62	0.31504267	0.945128
0.64	0.31778133	0.953344
0.66	0.320232	0.960696
0.68	0.32241067	0.967232
0.7	0.32433333	0.973
0.72	0.326016	0.978048
0.74	0.32747467	0.982424
0.76	0.32872533	0.986176
0.78	0.329784	0.989352
0.8	0.33066667	0.992
0.82	0.33138933	0.994168
0.84	0.331968	0.995904
0.86	0.33241867	0.997256
0.88	0.33275733	0.998272
0.9	0.333	0.999
0.92	0.33316267	0.999488
0.94	0.33326133	0.999784
0.96	0.333312	0.999936
0.98	0.33333067	0.999992
1	0.33333333	1

这个计算器能做什么

本工具会在给定形状参数 a 和 b 的情况下，针对一段 x 的取值区间，同时列出两个密切相关的特殊函数。下不完全贝塔函数 Bx(a,b) 是贝塔核函数从 0 到 x 的积分；正则化不完全贝塔函数 Ix(a,b) 则是这个积分再除以完全贝塔函数 B(a,b) 所得的结果。Ix(a,b) 同时也是 Beta(a,b) 分布的累积分布函数（CDF），并且是二项分布、学生 t 分布、F 分布等众多分布 CDF 背后的核心函数。

使用方法

输入形状参数 a 和 b（两者都必须为正数）。设置起始 x 值（即 x 的初始值）、步长以及重复次数（即表格的行数）。第 i 行对应的 x = 初始 x + i × 步长。所有取值都被限制在 [0, 1] 区间内；一旦 x 达到 1，表格便停止生成。顶部信息框会显示完全贝塔函数 B(a,b) 以及最后一行的结果，下方可滚动的表格则逐一列出每个 x 对应的 Bx(a,b) 与 Ix(a,b)。

公式详解

完全贝塔函数采用数值上更稳定的方式计算：B(a,b) = exp(lgamma(a) + lgamma(b) − lgamma(a+b))，其中对数伽马函数使用 Lanczos 近似法求得。正则化值 Ix(a,b) 则借助经典的《数值方法》（Numerical Recipes）连分式（即 Lentz 算法）来求解：令 bt = exp(lgamma(a+b) − lgamma(a) − lgamma(b) + a·ln x + b·ln(1−x))，当 x < (a+1)/(a+b+2) 时取 Ix = bt·betacf(a,b,x)/a，否则取 Ix = 1 − bt·betacf(b,a,1−x)/b。最后由 Bx(a,b) = Ix(a,b) × B(a,b) 得到结果。端点 I0 = 0 与 I1 = 1 会被精确处理，以避免出现 log(0)。

Advertisement

正则化不完全贝塔函数 I_x(a,b) 随着 x 从 0 变到 1 单调地从 0 增加到 1。

下不完全贝塔函数是被积函数下方从 0 到 x 的阴影面积。

计算实例

取 a = 1、b = 3：B(1,3) = Γ(1)Γ(3)/Γ(4) = (1·2)/6 = 1/3 ≈ 0.333333。在 x = 0.5 处，闭式公式给出 Bx = (1 − (1−x)^3)/3 = (1 − 0.125)/3 = 0.291667，于是 Ix = 0.291667 / 0.333333 = 0.875。再用对称性进行验算：Ix(1,3) = 1 − (1−x)^3 = 1 − 0.125 = 0.875，结果完全吻合。

常见问题

Bx 和 Ix 有什么区别？Bx(a,b) 是未经归一化的原始积分；Ix(a,b) = Bx(a,b)/B(a,b) 经过归一化处理，取值始终介于 0 和 1 之间。

为什么 a 和 b 必须为正数？只有当 a > 0 且 b > 0 时，相应的积分和伽马函数才能收敛。

如果步长让 x 超过了 1 怎么办？每个 x 都会被限制为最大 1，表格会在 x = 1 处停止；此时 Bx(a,b) = B(a,b)，Ix(a,b) = 1。