学过微积分的人,肯定都接触过Euler积分,按教科书上的说法,这是两种含有参变量的定积分,但其实没那么玄乎,它们只是两个函数.其中第一型Euler积分叫\(B\)-函数,第二型Euler积分叫\(\Gamma\)-函数,这两个函数的定义如下:\begin{align} \label{eq: beta} B (m, n) & = \int_0^1 x^{m-1} (1-x)^{n-1} \text{d} x \\ \label{eq: gamma} \Gamma (n) & = \int_0…

http://cos.name/2013/01/lda-math-gamma-function/ 1. 神奇的Gamma函数1.1 Gamma 函数诞生记学高等数学的时候,我们都学习过如下一个长相有点奇特的Gamma函数 Γ(x)=∫∞0tx−1e−tdt 通过分部积分的方法,可以推导出这个函数有如下的递归性质 Γ(x+1)=xΓ(x) 于是很容易证明,Γ(x) 函数可以当成是阶乘在实数集上的延拓,具有如下性质 Γ(n)=(n−1)! 学习了Gamma 函数之后,多年以来我一直有两个疑问: 这个…

Gamma函数 当n为正整数时,n的阶乘定义如下:n! = n * (n - 1) * (n - 2) * … * 2 * 1. 当n不是整数时,n!为多少?我们先给出答案. 容易证明,Γ(x + 1) = x * Γ(x),当n为正整数时,显然有Γ(n) = (n – 1)!. 计算(1/2)! 先给一个神奇的公式,证明不详述. (1) 定义如下函数 令上式p = 1,q = 1/2,同时根据对称性原理,有 (2) 同时容易证明 (3) 令p = 1/2,结合(2)(3)式,有 由于B关于q递…

神奇的gamma函数(上) 神奇的gamma函数(下) gamma函数的定义及重要性质 \[\Gamma(x)=\int_0^{\infty}t^{x-1}e^{-t}dt\] \[\Gamma(x+1) = x \Gamma(x)\] \[\Gamma(n) = (n-1)! \] \[\Gamma(0) = 1\] \[\Gamma({1\over 2}) = 2\int_0^{+\infty}e^{-u^2}du = \sqrt\pi\] gamma函数的图像 在matlib中,我们可以方…

Gamma 函数 Time Limit:1000MS Memory Limit:65536KBTotal Submit:237 Accepted:138 Description  Pollux最近在复习概率论与数理统计,他发现里面有很多有意思的积分,今天Pollux碰到了一个Gamma函数,定义如下: Input  第一行为一个整数T,表示测数数据的组数.接下去T行,每行一个整数n, (1<= n <=1000). Output  每组测试数据输出一行,对于每个n, 输出T(n)%1999 S…

1. Γ(⋅) 函数定义 Γ(α)=∫∞0tα−1e−tdt 可知以下基本性质: Γ(α+1)=αΓ(α)(分部积分法) Γ(1)=1 ⇒ Γ(n+1)=n! Γ(12)=π√ 2. 常见变形 对于 a>0,p>0: ⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪∫∞0xp−1e−axdx=a−pΓ(p)∫∞0x−p−1e−ax−1dx=a−pΓ(p)∫∞0xp−1e−ax2dx=12a−p2Γ(p2) 对其进行简单证明,∫∞0xp−1e−axdx=a−1a−(p−1)∫∞0(ax)p−1e−axda…

1. Γ(⋅) 函数 Γ(α)=∫∞0tα−1e−tdt 可知以下基本性质: Γ(α+1)=αΓ(α) Γ(1)=1 ⇒ Γ(n+1)=n! Γ(12)=π√ 2. 指数幂分布(exponential power distribution) f(x)=12q+1qΓ(q+1q)σexp(−12∣∣x−μσ∣∣2) 之所以说,指数幂分布是一种对正态分布的推广, q=2 ⇒ 正态分布 q=1 ⇒ 拉普拉斯分布…

Technorati Tags: LDA主题模型…

基础知识部分 为了方便理解,首先会对(Luminance)的相关概念做一个简单介绍.如果已经了解就跳到后面吧. 我们用Radiant energy(辐射能量)来描述光照的能量,单位是焦耳(J),因为光实际是以一定速度在传播的电磁波,所以把单位时间内的传播的Radiant energy(能量)称作radiant  flux(辐射通量),用来描述他的能量表现,单位瓦特(Watt). Radiant intensity(辐射强度)用来指定radiant  flux(辐射通量)的方向,正式的来说,他是用…

ndims(A)返回A的维数size(A)返回A各个维的最大元素个数length(A)返回max(size(A))[m,n]=size(A)如果A是二维数组,返回行数和列数nnz(A)返回A中非0元素的个数 MATLAB的取整函数:fix(x), floor(x) :,ceil(x) , round(x)(1)fix(x) : 截尾取整. >> fix( [3.12 -3.12]) ans = 3 -3(2)floor(x):不超过x 的最大整数.(高斯取整) >> floor(…