阶乘函数的解析延拓——伽马函数

数学小牛马

已于 2022-03-23 22:59:48 修改

阅读量3.4k

点赞数 1

分类专栏：数学文章标签：数学

于 2022-03-23 22:58:37 首次发布

本文链接：https://blog.csdn.net/qq_43409560/article/details/123699021

版权

数学专栏收录该内容

39 篇文章 3 订阅

订阅专栏

在这里插入图片描述

解析延拓

解析延拓可以看作是扩大原函数定义域的一种数学方式，例如 $f$ 函数扩张定义域至 $g$ 函数，且 $f$ 和 $g$ 函数再同一定义域下的定义（值）不变。

使用代数集合语言表达为：设 $A_0$ 为 $A$ 的子集， $f$ 为 $A_0$ 到 $B$ 的映射， $g$ 为 $A$ 到 $B$ 的映射， $\forall x\in A_0$ 有 $f (x) = g (x)$ ，则 $g$ 为 $f$ 的开拓，记为 $f=g|_{A_0}$ 。

很熟悉的例子是：
$\begin{aligned} S_n&=1+r+r^2+...+r^n+...\\ S_n&=\frac{1}{r-1}\quad(|r|<1)\\ \sum&_{i=0}^{+\infty}r^i=\frac{1}{r-1},(|r|<1) \end{aligned}$
所得结论的右式就是左式的解析延拓。

伽马函数

$\begin{aligned} \Gamma(s)=&\int_0^{+\infty}x^{s-1}e^{-x}dx\\=&\int_0^{1}x^{s-1}e^{-x}dx+\int_1^{+\infty}x^{s-1}e^{-x}dx\\ =&\Phi(s)+\Psi(s) \end{aligned}$

方程递推式

$\begin{aligned} \Gamma(s+1)=&\lim_{A\rightarrow+\infty}\int_0^{A}x^{s}e^{-x}dx\\=&-x^{s}e^{-x}|_{0}^{A}+\int_0^{A}e^{-x}dx^s\\ =&-A^{s}e^{-A}+s\int_0^{A}x^{s-1}e^{-x}dx\\ =&s\Gamma(s) \end{aligned}$

若参数是正整数，那么：
$\begin{aligned} \Gamma(s+1)=&s\Gamma(s)\\ =&(s+1)s*...*2*1\Gamma(1)\\ =&s!\int_0^{+\infty}e^{-x}dx\\ =&s! \end{aligned}$
可以看出对于非负整数， $\Gamma(n+1)=n!$ 。
当参数取零时：
$\begin{aligned} \Gamma(0)&=\int_0^{+\infty}x^{-1}e^{-x}dx\\ \lim_{\epsilon\rightarrow0}\int_\epsilon^{1}x^{-1}e^{-x}dx&\geq\frac{1}{e}\lim_{\epsilon\rightarrow0}\int_\epsilon^{1}x^{-1}dx\\ &=\frac{1}{e}\lim_{\epsilon\rightarrow0}lnx|_{\epsilon}^{1}\\ &=+\infty \end{aligned}$

所以得到伽马函数对负整数是没有定义的。
当参数取非整数时（以正数为例）：
$\begin{aligned} n < s&\le n+1\\ \Gamma(s+1)&=s\Gamma(s)\\ &=s(s-1)(s-n)\Gamma(s-n)\\ \end{aligned}$
如果已知 $\Gamma(s-n)$ 就可以算出 $\Gamma(s)$ 的值。

综上来看，伽马函数是阶乘函数的推广。

可视化

使用mathematica可视化实数平面图像与复平面图像，其中复平面取值为绝对值所以图像在实部小于零的情况下不绝对严谨。

在这里插入图片描述

f[x_] := Gamma[x];
Plot[f[x], {x, -3, 5}]

在这里插入图片描述

g[x_, y_] := Abs[Gamma[1.0*x + 1.0*y*I]];
Plot3D[g[x, y], {x, -3, 5}, {y, -2, 2}]

证明 $\Gamma(\frac{1}{2})=\sqrt{\pi}$

证明方式多种多样，这里暂时先给出一个例子

正态分布相关联：

正态分布是一类重要的分布，概率密度函数为：
$f_{\mu,\sigma}(x)=\frac{1}{\sqrt{2\pi \sigma ^2}}e^{-\frac{(x-\mu)^2}{2\sigma ^2}}dx$
为了更进一步地观察，将概率密度函数求积分：
$\int_{-\infty}^{+\infty}\frac{1}{\sqrt{2\pi \sigma ^2}}e^{-\frac{(x-\mu)^2}{2\sigma ^2}}dx=1$
接下来取 $\mu=0$ ， $\sigma=1$ 并使得积分对称得到：
$2\int_{0}^{+\infty}e^{-\frac{x^2}{2}}dx=\sqrt{2\pi}$
之后利用换元法 $u=\frac{x^2}{2}$ 发现 $dx=\frac{du}{\sqrt {2u}}$ ，带入上式：
$\sqrt2\int_{0}^{+\infty}u^{-\frac{1}{2}}e^{-u}dx=\sqrt{2\pi}$
最终得到：
$\int_{0}^{+\infty}u^{\frac{1}{2}-1}e^{-u}dx=\sqrt{\pi}$
可以发现等式左边的积分项为 $\Gamma(\frac{1}{2})$ 。