关于伽马函数
1、
函数格式:GammaP(a,b,value)
a,b:服从Γ(a,b)分布的函数.
value:代表服从Γ(a,b)分布的函数在x≥value时的概率值.
函数功能:概率统计的gamma函数在x≥value时的概率值.
2、
函数格式:Gammln(a)
a:Γ(a)当中的值
当a≤1的时候,采用多项式逼近的方式来求得Gammln(a)
函数功能:函数成功返回ln(Γ(a))的值.
3、
函数格式:Gamma(a)
a:Γ(a)当中的值
函数原理:
Γ(a+1)=(a+1)*Γ(a)
Γ(a)*Γ(1-a)=π/sin(π*a)
Γ(a)=Exp(Gammln(a))
函数功能:函数成功返回Γ(a)的值.
4、
函数格式:Gamma(a,value)
a:Γ(a,x)当中的a,a>0
value:Γ(a,value)当中的x,value>0
注意:上面a、b参数的值必须都大于0,且都必须小于171,否则程序会崩溃
函数功能:返回不完全伽马函数Γ(a,value)的值.
例子:
问:已知某元件寿命服从X~Γ(2,0.5)分布[单位小时],随便取一个元件,求该元件寿命大于4小时的概率.
答:执行本函数GammaP(2,0.5,4)后返回0.40600584970982即此值即为其概率值.
源代码:
Public Function Gamma(ByVal a As Double, ByVal value As Double) As Double '注意a必须大于0 '不完全伽马函数求值P(a,value) If value = 0 Then Return 0 ElseIf value > 1.0E+35 Then Retu