阶乘函数的解析延拓——伽马函数
解析延拓可以看作是扩大原函数定义域的一种数学方式
解析延拓
解析延拓可以看作是扩大原函数定义域的一种数学方式,例如fff函数扩张定义域至ggg函数,且fff和ggg函数再同一定义域下的定义(值)不变。
使用代数集合语言表达为:设A0A_0A0为AAA的子集,fff为A0A_0A0到BBB的映射,ggg为AAA到BBB的映射,∀x∈A0\forall x\in A_0∀x∈A0有f(x)=g(x)f(x)=g(x)f(x)=g(x),则ggg为fff的开拓,记为f=g∣A0f=g|_{A_0}f=g∣A0。
很熟悉的例子是:
Sn=1+r+r2+...+rn+...Sn=1r−1(∣r∣<1)∑i=0+∞ri=1r−1,(∣r∣<1) \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} SnSn∑=1+r+r2+...+rn+...=r−11(∣r∣<1)i=0+∞ri=r−11,(∣r∣<1)
所得结论的右式就是左式的解析延拓。
伽马函数
Γ(s)=∫0+∞xs−1e−xdx=∫01xs−1e−xdx+∫1+∞xs−1e−xdx=Φ(s)+Ψ(s) \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} Γ(s)===∫0+∞xs−1e−xdx∫01xs−1e−xdx+∫1+∞xs−1e−xdxΦ(s)+Ψ(s)
方程递推式
Γ(s+1)=limA→+∞∫0Axse−xdx=−xse−x∣0A+∫0Ae−xdxs=−Ase−A+s∫0Axs−1e−xdx=sΓ(s) \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} Γ(s+1)====A→+∞lim∫0Axse−xdx−xse−x∣0A+∫0Ae−xdxs−Ase−A+s∫0Axs−1e−xdxsΓ(s)
-
若参数是正整数,那么:
Γ(s+1)=sΓ(s)=(s+1)s∗...∗2∗1Γ(1)=s!∫0+∞e−xdx=s! \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} Γ(s+1)====sΓ(s)(s+1)s∗...∗2∗1Γ(1)s!∫0+∞e−xdxs!
可以看出对于非负整数,Γ(n+1)=n!\Gamma(n+1)=n!Γ(n+1)=n!。 -
当参数取零时:
Γ(0)=∫0+∞x−1e−xdxlimϵ→0∫ϵ1x−1e−xdx≥1elimϵ→0∫ϵ1x−1dx=1elimϵ→0lnx∣ϵ1=+∞ \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} Γ(0)ϵ→0lim∫ϵ1x−1e−xdx=∫0+∞x−1e−xdx≥e1ϵ→0lim∫ϵ1x−1dx=e1ϵ→0limlnx∣ϵ1=+∞所以得到伽马函数对负整数是没有定义的。
-
当参数取非整数时(以正数为例):
n<s≤n+1Γ(s+1)=sΓ(s)=s(s−1)(s−n)Γ(s−n) \begin{aligned} n < s&\le n+1\\ \Gamma(s+1)&=s\Gamma(s)\\ &=s(s-1)(s-n)\Gamma(s-n)\\ \end{aligned} n<sΓ(s+1)≤n+1=sΓ(s)=s(s−1)(s−n)Γ(s−n)
如果已知Γ(s−n)\Gamma(s-n)Γ(s−n)就可以算出Γ(s)\Gamma(s)Γ(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}]
证明Γ(12)=π\Gamma(\frac{1}{2})=\sqrt{\pi}Γ(21)=π
证明方式多种多样,这里暂时先给出一个例子
正态分布相关联:
正态分布是一类重要的分布,概率密度函数为:
fμ,σ(x)=12πσ2e−(x−μ)22σ2dx f_{\mu,\sigma}(x)=\frac{1}{\sqrt{2\pi \sigma ^2}}e^{-\frac{(x-\mu)^2}{2\sigma ^2}}dx fμ,σ(x)=2πσ21e−2σ2(x−μ)2dx
为了更进一步地观察,将概率密度函数求积分:
∫−∞+∞12πσ2e−(x−μ)22σ2dx=1 \int_{-\infty}^{+\infty}\frac{1}{\sqrt{2\pi \sigma ^2}}e^{-\frac{(x-\mu)^2}{2\sigma ^2}}dx=1 ∫−∞+∞2πσ21e−2σ2(x−μ)2dx=1
接下来取μ=0\mu=0μ=0,σ=1\sigma=1σ=1并使得积分对称得到:
2∫0+∞e−x22dx=2π 2\int_{0}^{+\infty}e^{-\frac{x^2}{2}}dx=\sqrt{2\pi} 2∫0+∞e−2x2dx=2π
之后利用换元法u=x22u=\frac{x^2}{2}u=2x2发现dx=du2udx=\frac{du}{\sqrt {2u}}dx=2udu,带入上式:
2∫0+∞u−12e−udx=2π \sqrt2\int_{0}^{+\infty}u^{-\frac{1}{2}}e^{-u}dx=\sqrt{2\pi} 2∫0+∞u−21e−udx=2π
最终得到:
∫0+∞u12−1e−udx=π \int_{0}^{+\infty}u^{\frac{1}{2}-1}e^{-u}dx=\sqrt{\pi} ∫0+∞u21−1e−udx=π
可以发现等式左边的积分项为Γ(12)\Gamma(\frac{1}{2})Γ(21)。
更多推荐
所有评论(0)