斯特林公式：从历史到严格推导与应用

定义与核心形式

斯特林公式的严格数学表述包含两种常用形式：

对数形式适用于理论分析，其表达式为：

其中 表示当 时趋于 0 的余项。

指数形式适用于数值计算，其表达式为：

展开式中的系数与伯努利数相关，构成斯特林级数，可通过欧拉 - 麦克劳林公式推导得到。

推导方法一：积分近似与 Wallis 公式

第一步：对数阶乘的积分估计

对阶乘取自然对数，转化为求和问题：

将求和视为积分的近似，利用函数 的单调性，可得：

计算积分 ，代入得：

两端同除以 并取极限，可知 与 同阶，即 。

第二步：引入待定常数

假设 ，其中 为常数。取指数后得：

接下来需确定 。

第三步：利用 Wallis 公式定常数值

Wallis 公式表述为：

其中 ， 。将斯特林近似式代入分子分母：

代入 Wallis 公式并化简，最终解得 ，从而完成证明。

推导方法二：鞍点法与 Gamma 函数

第一步：Gamma 函数的积分表示

阶乘可通过 Gamma 函数推广至非整数情形： 。令 ，作变量替换得：

被积函数 在 处取最大值（鞍点），此时 ， 。

第二步：泰勒展开与高斯积分近似

在鞍点 附近展开 ：

当 时，积分主要贡献来自 邻域，可将积分限扩展至 ，并近似 。此时积分转化为高斯积分：

代入 Gamma 函数表达式，即得：

此方法展现了鞍点法在渐近分析中的核心思想：通过指数函数的集中性，将积分简化为极值点附近的二次近似。

推导方法三：Bohr-Mollerup 定理与函数方程

第一步：对数凸性与函数方程

Bohr-Mollerup 定理指出，满足以下条件的函数 唯一确定为 Gamma 函数： （函数方程）； ； 为凸函数。

构造函数 ，可证其单调递减且有下界，故极限 存在。结合函数方程与凸性条件，可唯一确定 ，从而导出斯特林公式。

误差分析与高精度展开

斯特林公式的余项可通过欧拉 - 麦克劳林公式精确表示：

其中 为伯努利数（如 ， ），余项 满足 。例如，取 时，公式精度提升为：

该展开在数值计算中具有重要应用，例如当 时，斯特林近似值为 ，与精确值 的相对误差仅为 ，加入 修正项后误差降至 。

应用领域与科学意义

斯特林公式的价值体现在理论与应用的双重层面：

在概率论与统计力学中，斯特林公式用于近似二项分布的概率质量函数，证明其渐近收敛于正态分布。例如，对二项分布 ，当 且 时，通过斯特林公式可推导出其密度函数趋于 。

在算法复杂度分析中，斯特林公式用于估计排列数、子集数等组合量的增长速度。例如，排序算法的时间复杂度 本质上源于阶乘的对数增长率。

在数值计算中，斯特林公式为计算器、编程语言中的阶乘函数提供底层实现。例如，Python 的 math.factorial 函数在大整数情形下即采用斯特林近似以平衡精度与效率。

在物理应用中，斯特林公式用于近似玻尔兹曼熵公式 ，其中 为系统微观状态数，通过阶乘的对数近似可简化熵的计算。

历史演变与学术影响

斯特林公式的发展折射了数学分析的历史进程：棣莫弗在 1733 年的《机会学说》中得到 ，但未能确定常数 ；斯特林通过双纽线积分的数值计算发现 ，并在 1730 年的著作中发表完整结果。欧拉后续将阶乘推广为 Gamma 函数，为公式的严格化奠定基础。20 世纪以来，随着渐近分析理论的发展，斯特林公式的余项估计、高维推广（如多元 Gamma 函数的渐近展开）成为研究热点，其思想也渗透到复分析、数论等领域。

从棣莫弗的初步探索到现代分析中的精细化研究，斯特林公式始终是 “近似” 与 “严格” 辩证统一的典范。它不仅是计算工具，更启发了鞍点法、欧拉 - 麦克劳林公式等一系列分析方法，成为连接离散与连续、有限与无穷的桥梁。在大数据与复杂系统研究日益重要的今天，斯特林公式所代表的渐近思维，仍将在科学前沿发挥不可替代的作用。