UA MATH564 概率论 高阶矩的计算:有限差算子方法1
UA MATH564 概率論 高階矩的計算:有限差算子方法1
- 有限差算子及其性質
- 0的有限差與Stirling數之間的關系
- 計算離散型隨機變量的n階矩
有限差算子(finite difference),也被稱為向前差分(forward difference)是一種非常古老的數學工具了,但它在離散型隨機變量的期望方面具有非常強大的實力。
有限差算子及其性質
假設f(x)f(x)f(x)是定義在R\mathbb{R}R上的實函數,定義有限差算子為
Δf(x)=f(x+1)?f(x)\Delta f(x) = f(x+1) - f(x)Δf(x)=f(x+1)?f(x)
定義位移算子為Mf(x)=f(x+1)Mf(x) = f(x+1)Mf(x)=f(x+1)
根據定義可以直接得到如下性質:
這里的1是恒等算子。性質1非常簡單,性質2也只是簡單套用定義:
Δf(x)=f(x+1)?f(x)=Mf(x)?1f(x)=(M?1)f(x)\Delta f(x) = f(x+ 1) - f(x) = Mf(x) - 1f(x) = (M-1)f(x)Δf(x)=f(x+1)?f(x)=Mf(x)?1f(x)=(M?1)f(x)
根據性質2可以知道有限差算子與位移算子之間的關系是:
Δ=M?1M=Δ+1\Delta = M - 1 \\ M = \Delta + 1Δ=M?1M=Δ+1
基于他們之間的關系可以進一步推得下面兩個非常有用的公式:
Δmf(x)=(M?1)mf(x)=∑j=0m(?1)jCmjMm?jf(x)=∑j=0m(?1)jCmjf(x+m?j)Mmf(x)=(1+Δ)mf(x)=∑j=0mCmjΔjf(x)\Delta^m f(x) = (M-1)^m f(x) = \sum_{j=0}^m (-1)^jC_m^jM^{m-j}f(x) =\sum_{j=0}^m (-1)^jC_m^jf(x+m-j) \\ M^mf(x) = (1+\Delta )^m f(x) = \sum_{j=0}^m C_m^j \Delta^j f(x)Δmf(x)=(M?1)mf(x)=j=0∑m?(?1)jCmj?Mm?jf(x)=j=0∑m?(?1)jCmj?f(x+m?j)Mmf(x)=(1+Δ)mf(x)=j=0∑m?Cmj?Δjf(x)
下面我們利用這兩個公式證明一個看似不可思議的定理。
0的有限差與Stirling數之間的關系
定理1
Δm0n=m!S2(n,m),n≥m\Delta^m 0^n = m! S_2(n,m),n \ge mΔm0n=m!S2?(n,m),n≥m
這個關系的含義是:0的nnn次冪的mmm階有限差等于mmm的階乘乘以第二類Stirling數。0的nnn次冪的mmm階有限差就是
Δm0n=(M?1)m0n=∑j=0m(?1)m?jCmjMj0n=∑j=0m(?1)m?jCmjjn\Delta^m 0^n = (M-1)^m 0^n = \sum_{j=0}^m (-1)^{m-j}C_m^jM^j0^n = \sum_{j=0}^m (-1)^{m-j} C_m^jj^nΔm0n=(M?1)m0n=j=0∑m?(?1)m?jCmj?Mj0n=j=0∑m?(?1)m?jCmj?jn
這個定理在計算離散型隨機變量的期望時有非常大的用處。
證明
第二類Stirling數有如下性質:
xn=∑j=0nS2(n,j)(x)j?Δmxn=∑j=0nS2(n,j)Δm(x)jx^n = \sum_{j=0}^n S_2(n,j)(x)_j \\ \Rightarrow \Delta^m x^n = \sum_{j=0}^n S_2(n,j)\Delta^m(x)_j xn=j=0∑n?S2?(n,j)(x)j??Δmxn=j=0∑n?S2?(n,j)Δm(x)j?
下面的計算需要用到引理1,先默認引理1是正確的,推導結束后再來證明引理1。
RHS=∑j=0nS2(n,j)(j)m(x)j?mI(m≤j)=∑j=mnS2(n,j)(j)m(x)j?mRHS = \sum_{j=0}^n S_2(n,j) (j)_m(x)_{j-m}I(m \le j) = \sum_{j=m}^n S_2(n,j) (j)_m(x)_{j-m}RHS=j=0∑n?S2?(n,j)(j)m?(x)j?m?I(m≤j)=j=m∑n?S2?(n,j)(j)m?(x)j?m?
令x=0x=0x=0,則當且僅當j=mj=mj=m時,(x)j?m=0!=1(x)_{j-m} = 0! = 1(x)j?m?=0!=1,否則(x)j?m=0(x)_{j-m} = 0(x)j?m?=0,所以
Δm0n=m!S2(n,m)\Delta^m 0^n = m!S_2(n,m)Δm0n=m!S2?(n,m)
引理1
Δm(x)j=(j)m(x)j?mI(m≤j)\Delta^m(x)_j = (j)_m(x)_{j-m}I(m \le j)Δm(x)j?=(j)m?(x)j?m?I(m≤j)
用數學歸納法,當m=1m=1m=1時,
Δ(x)j=(x+1)j?(xj)=(x+1)x(x?1)?(x+2?j)?x(x?1)?(x+2?j)(x+1?j)=[(x+1)?(x+1?j)](x)j?1=j(x)j?1\Delta (x)_j = (x+1)_j - (x_j) \\= (x+1)x(x-1)\cdots (x+2-j) - x(x-1) \cdots (x+2-j)(x+1-j) \\ = [(x+1)-(x+1-j)](x)_{j-1} = j(x)_{j-1}Δ(x)j?=(x+1)j??(xj?)=(x+1)x(x?1)?(x+2?j)?x(x?1)?(x+2?j)(x+1?j)=[(x+1)?(x+1?j)](x)j?1?=j(x)j?1?
顯然是服從引理1的;假設引理1對m=km=km=k成立,則
Δk+1(x)j=Δ(Δk(x)j)=(j)k(x+1)j?kI(k≤j)?(j)k(x)j?kI(k≤j)=(j)kI(k≤j)[(x+1)?(x+1?(j?k))](x)j?k?1=(j)kI(k≤j)(j?k)I(k+1≤j)(x)j?k?1=(j)k+1(x)j?k?1I(k+1≤j)\Delta^{k+1} (x)_j = \Delta \left( \Delta^{k}(x)_j \right)= (j)_k(x+1)_{j-k}I(k \le j) -(j)_k(x)_{j-k}I(k \le j) \\ = (j)_kI(k \le j)[(x+1)-(x+1-(j-k))](x)_{j-k-1} \\ = (j)_kI(k \le j)(j-k)I(k +1\le j)(x)_{j-k-1} \\ = (j)_{k+1}(x)_{j-k-1}I(k+1 \le j)Δk+1(x)j?=Δ(Δk(x)j?)=(j)k?(x+1)j?k?I(k≤j)?(j)k?(x)j?k?I(k≤j)=(j)k?I(k≤j)[(x+1)?(x+1?(j?k))](x)j?k?1?=(j)k?I(k≤j)(j?k)I(k+1≤j)(x)j?k?1?=(j)k+1?(x)j?k?1?I(k+1≤j)
即引理1對m=k+1m=k+1m=k+1也成立。
證畢
計算離散型隨機變量的n階矩
假設X=0,1,?X = 0,1,\cdotsX=0,1,?,并且P(X=j)=pjP(X = j) = p_jP(X=j)=pj?,則它的nnn階原點矩為
EXn=∑j=0∞jnpjEX^n = \sum_{j=0}^{\infty} j^n p_jEXn=j=0∑∞?jnpj?
替換掉jnj^njn,考慮
Mjxn=(x+j)n?Mj0n=jnM^jx^n = (x+j)^n \Rightarrow M^j0^n = j^nMjxn=(x+j)n?Mj0n=jn
因此
EXn=∑j=0∞pjMj0n=∑j=0∞pj(Δ+1)j0nEX^n = \sum_{j=0}^{\infty} p_jM^j0^n = \sum_{j=0}^{\infty} p_j(\Delta + 1)^j0^nEXn=j=0∑∞?pj?Mj0n=j=0∑∞?pj?(Δ+1)j0n
下一講介紹怎么用這個公式計算二項分布、Poisson分布、超幾何分布和負二項分布的nnn階矩。
總結
以上是生活随笔為你收集整理的UA MATH564 概率论 高阶矩的计算:有限差算子方法1的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH571B 试验设计VI 随
- 下一篇: UA MATH566 统计理论 截断数据