日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

UA MATH564 概率论 高阶矩的计算:有限差算子方法1

發布時間:2025/4/14 编程问答 34 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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)

根據定義可以直接得到如下性質:

  • Mmf(x)=f(x+m)M^mf(x) = f(x+m)Mmf(x)=f(x+m)
  • Δf(x)=(M?1)f(x)\Delta f(x) = (M-1)f(x)Δf(x)=(M?1)f(x)
  • 這里的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=0m?(?1)jCmj?Mm?jf(x)=j=0m?(?1)jCmj?f(x+m?j)Mmf(x)=(1+Δ)mf(x)=j=0m?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),nm

    這個關系的含義是: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=0m?(?1)m?jCmj?Mj0n=j=0m?(?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=0n?S2?(n,j)(x)j??Δmxn=j=0n?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=0n?S2?(n,j)(j)m?(x)j?m?I(mj)=j=mn?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(mj)

    用數學歸納法,當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(kj)?(j)k?(x)j?k?I(kj)=(j)k?I(kj)[(x+1)?(x+1?(j?k))](x)j?k?1?=(j)k?I(kj)(j?k)I(k+1j)(x)j?k?1?=(j)k+1?(x)j?k?1?I(k+1j)

    即引理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的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。