关于α→0且x+α非小情形下Γ(x+α)的近似方法与上下界的技术咨询

嗨，我来帮你梳理这个问题的解法和思路~ 当α趋近于0且x+α不是很小（也就是x本身是正实数，且α的绝对值足够小，保证x+α>0）时，我们可以用泰勒展开、对数凸性等工具来得到更精准的近似和可靠的上下界，下面具体展开：

一、更精准的近似方法 #

直接用$\Gamma(x+\alpha) \approx \Gamma(x)$效果差，是因为忽略了α带来的一阶及高阶修正项，我们可以从Γ函数的对数泰勒展开入手，这是最通用的方法：

1. 一阶近似（最实用）
   利用digamma函数（$\psi(z) = \frac{d}{dz} \ln\Gamma(z)$）的定义，$\ln\Gamma(x+\alpha)$在$\alpha=0$处的一阶泰勒展开为：

   lnΓ(x+α) = lnΓ(x) + αψ(x) + o(α)
   

   两边取指数后得到：

   Γ(x+α) ≈ Γ(x) · exp(αψ(x))
   

   其中digamma函数$\psi(x)$有明确的表达式和近似：

   • 对于$x>0$，$\psi(x) = -\gamma + \sum_{k=0}^\infty \left( \frac{1}{k+1} - \frac{1}{x+k} \right)$（$\gamma$是欧拉常数≈0.5772）
   • 当$x$较大时，$\psi(x) \approx \ln(x) - \frac{1}{2x} - \frac{1}{12x^2} + o\left(\frac{1}{x^{2}\right)$，这时候代入近似式就会和你提到的$x\to\infty$时的$\Gamma(x+\alpha)\approx\Gamma(x)x}\alpha$结果对齐（因为$\exp\left(\alpha\left(\ln x - \frac{1}{2x}\right)\right) = x^\alpha \cdot \exp\left(-\frac{\alpha}{2x}\right)$，当$x\to\infty$时$\exp\left(-\frac{\alpha}{2x}\right)\approx1$）。

2. 二阶近似（更高精度）
   如果需要更精准的结果，可以保留泰勒展开的二阶项（用到trigamma函数$\psi'(z)=\frac{d}{dz}\psi(z)$）：

   lnΓ(x+α) = lnΓ(x) + αψ(x) + \frac{\alpha^2}{2}\psi'(x) + o(\alpha^2)
   Γ(x+α) ≈ Γ(x) · exp\left(αψ(x) + \frac{\alpha^2}{2}\psi'(x)\right)
   

   trigamma函数$\psi'(x) = \sum_{k=0}^\infty \frac{1}{(x+k)^2}$，对于大$x$也有近似$\psi'(x)\approx\frac{1}{x} + \frac{1}{2x^2} + \frac{1}{6x^3} + o\left(\frac{1}{x^3}\right)$。

举个例子：$x=1$，$\alpha=0.1$，$\Gamma(1)=1$，$\psi(1)=-\gamma\approx-0.5772$，代入一阶近似得$1\cdot\exp(-0.1\times0.5772)\approx0.944$，而真实值$\Gamma(1.1)\approx0.9514$，误差仅约0.8%，比直接取$\Gamma(x)=1$的误差小很多。

二、上下界估计 #

我们可以利用$\ln\Gamma(z)$的凸性（$\ln\Gamma$在$z>0$时是严格凸函数）和digamma函数的单调性（$\psi(z)$在$z>0$时严格递增）来推导上下界：

情况1：α>0（小正数） #

根据凸函数的切线不等式（凸函数的图像在切线上方），有：

lnΓ(x+α) ≥ lnΓ(x) + αψ(x)

即：

Γ(x+α) ≥ Γ(x) · exp(αψ(x))

同时，利用中值定理，存在$0<\theta<1$使得：

lnΓ(x+α) - lnΓ(x) = αψ(x + θα)

因为$\psi(z)$递增，所以$\psi(x) < \psi(x+\theta\alpha) < \psi(x+\alpha)$，因此：

Γ(x) · exp(αψ(x)) < Γ(x+α) < Γ(x) · exp(αψ(x+α))

如果结合大$x$时$\psi(x)\approx\ln x - \frac{1}{2x}$，还能得到更直观的上下界：

Γ(x) · x^α · exp\left(-\frac{\alpha}{2x}\right) < Γ(x+α) < Γ(x) · x^α

情况2：α<0（小负数，令α=-β，β>0，保证x-β>0） #

类似地，凸性给出：

lnΓ(x-β) ≥ lnΓ(x) - βψ(x)

即：

Γ(x-β) ≥ Γ(x) · exp(-βψ(x))

再用中值定理，存在$0<\theta<1$使得：

lnΓ(x-β) - lnΓ(x) = -βψ(x - θβ)

因为$\psi(z)$递增，所以$\psi(x - \theta\beta) < \psi(x)$，因此：

Γ(x) · exp(-βψ(x)) < Γ(x-β) < Γ(x) · exp(-βψ(x - β))

三、补充说明 #

如果你需要更简洁的上下界，还可以利用Γ函数的单调性：当$x>0$时，$\Gamma(z)$在$z\in(0,1)$时递减，$z>1$时递增。不过这个性质得到的界通常不如基于digamma的界紧凑，更推荐用前面的凸性推导结果。

总之，基于对数泰勒展开的近似比直接取$\Gamma(x)$精准得多，而利用凸性和digamma单调性得到的上下界在α很小时是非常紧的。

备注：内容来源于stack exchange，提问作者mathseeker