隨機過程學習筆記——概論

RealEchooo發表於2020-10-06


參考教材:陸大jin《隨機過程及其應用》

1. 隨機過程

1.1 基本概念

隨機過程是這樣一個過程,它不能用一個時間t的確定性函式表示,它在每個時刻的狀態是隨機的。

  1. 例子1:
  • y = s i n ( w t + φ ) y=sin(wt+\varphi) y=sin(wt+φ),其中 φ = π / 2 \varphi=\pi/2 φ=π/2。顯然這是一個確定性過程,這個過程可以用時間t的確定性函式表示。
  • y = s i n ( w t + φ ) y=sin(wt+\varphi) y=sin(wt+φ),其中 φ \varphi φ服從 [ − π , π ] [-\pi,\pi] [π,π]的均勻分佈。顯然這是一個隨機過程,這個過程不可以用時間t的確定性函式表示,每個時刻 t i t_i ti y y y的取值是隨機的, y y y它又是隨時間變化的。這個函式可以有一個現實的意義,比如生產的一批振盪器,他們的初始相位顯然是隨機的,但是選中一臺振盪器後,它的相位就是固定的,也就是說,這個隨機過程,如果隨著時間t去逐步測量,它得到的總是一條正弦曲線。
  1. 例子2:
    某個餐廳,一個時刻到達的人數是隨機的,也就是隨機變數。如果研究一個上午各個時間的總人數,那麼這就是一個隨機過程

1.2 描述隨機過程的方法

  1. 方法一:固定時間 t t t,則隨機過程 ξ ( t ) \xi(t) ξ(t)在此時刻是隨機變數,即: ξ ( t i ) \xi(t_{i}) ξ(ti)是隨機變數, t i t_i ti是具體的時刻。
    既然是隨機變數,那麼每個時刻都可以用分佈函式來描述。
    一維分佈函式 F t 1 ( x 1 ) = P ( ξ ( t 1 ) < x 1 ) F_{t_1}(x_{1})=P(\xi(t_{1})<x_1) Ft1(x1)=P(ξ(t1)<x1),這個只能描述某一時刻,隨機過程 ξ ( t ) \xi(t) ξ(t)的分佈。
    n維分佈函式 F t 1 , t 2 , . . . , t n ( x 1 , x 2 , . . . , x n ) = P ( ξ ( t 1 ) < x 1 , ξ ( t 2 ) < x 2 , . . . , ξ ( t n ) < x n ) F_{t_1,t_2,...,t_n}(x_{1},x_{2},...,x_{n})=P(\xi(t_{1})<x_1,\xi(t_{2})<x_2, ... ,\xi(t_{n})<x_n) Ft1,t2,...,tn(x1,x2,...,xn)=P(ξ(t1)<x1,ξ(t2)<x2,...,ξ(tn)<xn),這個能描述隨機過程 ξ ( t ) \xi(t) ξ(t),在不同時刻的相互關係,因此這也被稱為一族有限維分佈函式。可以看到這一族有限維分佈函式不僅描述了某一時刻 ξ ( t ) \xi(t) ξ(t)的統計規律,還描述了不同時刻 ξ ( t i ) \xi(t_i) ξ(ti), ξ ( t j ) \xi(t_j) ξ(tj)之間的相互關係。

注:

  • ξ ( t ) \xi(t) ξ(t) ξ ( w , t ) \xi(w,t) ξ(w,t)的含義相同
  • ξ ( t ) \xi(t) ξ(t) t = t k t=t_k t=tk的具體取值稱為,隨機過程 ξ ( t ) \xi(t) ξ(t) t = t k t=t_k t=tk的狀態,記為 x ( t k ) x(t_k) x(tk)
  1. 方法二:不固定時間,一次特定的實驗,隨著時間去逐步測量整個過程,因此這將得到一個確定的樣本函式,記為 ξ ( k ) ( t ) \xi^{(k)}(t) ξ(k)(t),有時為了避免混淆,也用 x k ( t ) x_k(t) xk(t)表示。在例子1中的隨機過程,每個 x k ( t ) x_k(t) xk(t)都是正弦函式。

2. 隨機過程的分類和舉例

根據時間的連續與否,每個時刻隨機變數取值的連續與否,可以分為4種型別。

  1. 離散引數離散型隨機過程。引數集(也就是時間)是離散的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)也是離散的。
    比如伯努利過程:一個質點每過一秒就會等概率地在數軸上向左移動或者向右移動。

  2. 連續引數離散型隨機過程。引數集(也就是時間)是連續的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)是離散的。
    比如脈衝數字通訊系統,傳送的訊號是脈寬為 T 0 T_0 T0的脈衝訊號,但是脈衝的幅度是隨機變數。

  3. 連續引數連續型隨機過程。引數集(也就是時間)是連續的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)也是連續的。
    比如正弦波過程(例子1)和電晶體噪聲。

  4. 離散引數連續型隨機過程。引數集(也就是時間)是離散的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)也是連續的。
    比如每隔單位時間就對電晶體噪聲進行抽樣。

3. 隨機過程的數字特徵

從以上的分析來看,隨機過程的概率密度分佈函式族可以完善地描述隨機過程的統計特徵,但是要確定這兩個東西,並加以分析,往往比較困難。為了便於分析,需要引入隨機過程的數字特徵。

3.1 均值(數學期望)

對於一隨機過程 ξ ( t ) \xi(t) ξ(t),其某一時刻 t 1 t_1 t1的取值 ξ ( t 1 ) \xi(t_1) ξ(t1)是隨機變數(這一點要時刻記住,不然下面的公式就會理解不了),對於這個隨機變數,可以用概率論中的均值來描述,即:
μ ξ ( t 1 ) = E ( ξ ( t 1 ) ) \mu_{\xi}(t_1)=E(\xi(t_1)) μξ(t1)=E(ξ(t1))
對於連續型隨機過程,假設其分佈函式為 F ξ t 1 ( x ) = P ( ξ ( t 1 ) ≤ x ) F_{\xi_{t_1}}(x)=P(\xi(t_1)\leq x) Fξt1(x)=P(ξ(t1)x),概率密度為 f ξ t 1 ( x ) f_{\xi_{t_1}}(x) fξt1(x),那麼均值:
μ ξ ( t 1 ) = E ( ξ ( t 1 ) ) = ∫ − ∞ + ∞ x f ξ t 1 ( x ) d x \mu_{\xi}(t_1)=E(\xi(t_1)) =\int_{-\infty}^{+\infty}xf_{\xi_{t_1}}(x)dx μξ(t1)=E(ξ(t1))=+xfξt1(x)dx
對於離散型隨機過程,假設其分佈函式為 F ξ t 1 ( x ) = P ( ξ ( t 1 ) ≤ x ) F_{\xi_{t_1}}(x)=P(\xi(t_1)\leq x) Fξt1(x)=P(ξ(t1)x),那麼均值:
μ ξ ( t 1 ) = E ( ξ ( t 1 ) ) = x 1 ∗ P ( ξ ( t 1 ) = x 1 ) + x 2 ∗ P ( ξ ( t 1 ) = x 2 ) + . . . + x n ∗ P ( ξ ( t 1 ) = x n ) \mu_{\xi}(t_1) = E(\xi(t_1)) \\ = x_1*P(\xi(t_1)=x_1)+x_2*P(\xi(t_1)=x_2)+...+x_n*P(\xi(t_1)=x_n) μξ(t1)=E(ξ(t1))=x1P(ξ(t1)=x1)+x2P(ξ(t1)=x2)+...+xnP(ξ(t1)=xn)

技巧
有時候並不需要真的去求概率密度為 f ξ t 1 ( x ) f_{\xi_{t_1}}(x) fξt1(x)才能算均值。
比如: y = s i n ( w t + φ ) y=sin(wt+\varphi) y=sin(wt+φ),其中 φ \varphi φ服從 [ − π , π ] [-\pi,\pi] [π,π]的均勻分佈,求 μ ξ ( t 1 ) \mu_{\xi}(t_1) μξ(t1)。這個問題就變成了概率論中經典的已知X的分佈,且Y=g(X),求E(Y)。所以:
μ ξ ( t 1 ) = ∫ − ∞ + ∞ x f ξ t 1 ( x ) d x = ∫ − π π s i n ( w t 1 + φ ) 1 2 π d φ \mu_{\xi}(t_1)= \int_{-\infty}^{+\infty}xf_{\xi_{t_1}}(x)dx= \int_{-\pi}^{\pi}sin(wt_1+\varphi)\frac{1}{2\pi} d\varphi μξ(t1)=+xfξt1(x)dx=ππsin(wt1+φ)2π1dφ
按照上式計算就會簡單很多。

3.2 方差(二階中心矩)

這個方差的定義與概率論中隨機變數的方差定義完全一致。
σ ξ 2 ( t 1 ) = D ( ξ ( t 1 ) ) = E ( [ ξ ( t 1 ) − μ ξ ( t 1 ) ] 2 ) = E ( [ ξ ( t 1 ) − E ( ξ ( t 1 ) ) ] 2 ) \sigma_{\xi}^2(t_1)=D(\xi(t_1))=E([\xi(t_1)-\mu_{\xi}(t_1)]^2) \\ =E([\xi(t_1)-E(\xi(t_1))]^2) σξ2(t1)=D(ξ(t1))=E([ξ(t1)μξ(t1)]2)=E([ξ(t1)E(ξ(t1))]2)

3.3 自相關函式(簡稱:相關函式)

自相關函式是為了描述兩個不同的時刻 t 1 , t 2 t_1,t_2 t1,t2的隨機變數 ξ ( t 1 ) , ξ ( t 2 ) \xi(t_1),\xi(t_2) ξ(t1),ξ(t2)之間的關係,記為 R ξ ξ ( t 1 , t 2 ) R_{\xi\xi}(t_1,t_2) Rξξ(t1,t2),這是一個二階混合原點矩。注意這兩個隨機變數並不獨立
R ξ ξ ( t 1 , t 2 ) = E ( ξ ( t 1 ) ξ ( t 2 ) ) = ∑ i ∑ j x i x j P ( ξ ( t 1 ) = x i , ξ ( t 2 ) = x j ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ x 1 x 2 f ξ t 1 ξ t 2 ( x 1 , x 2 ) d x 1 d x 2 R_{\xi\xi}(t_1,t_2)=E(\xi(t_1)\xi(t_2)) \\ = \sum_{i}\sum_{j}x_ix_jP(\xi(t_1)=x_i,\xi(t_2)=x_j) \\ =\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}x_1x_2f_{\xi_{t_1}\xi_{t_2}}(x_1, x_2)dx_1dx_2 Rξξ(t1,t2)=E(ξ(t1)ξ(t2))=ijxixjP(ξ(t1)=xi,ξ(t2)=xj)=++x1x2fξt1ξt2(x1,x2)dx1dx2
f ξ t 1 ξ t 2 ( x 1 , x 2 ) f_{\xi_{t_1}\xi_{t_2}}(x_1, x_2) fξt1ξt2(x1,x2)是隨機變數 ξ ( t 1 ) , ξ ( t 2 ) \xi(t_1),\xi(t_2) ξ(t1),ξ(t2)的聯合概率密度。

3.4 自協方差函式(簡稱:協方差函式)

類似地,可以寫出兩個變數的二階混合中心矩,記為 C ξ ξ ( t 1 , t 2 ) C_{\xi\xi}(t_1,t_2) Cξξ(t1,t2)
C ξ ξ ( t 1 , t 2 ) = E ( [ ξ ( t 1 ) − μ ξ ( t 1 ) ] [ ξ ( t 2 ) − μ ξ ( t 2 ) ] ) = ∑ i ∑ j ( x i − μ ξ ( t 1 ) ) ( x j − μ ξ ( t 2 ) ) P ( ξ ( t 1 ) = x i , ξ ( t 2 ) = x j ) = ∫ − ∞ + ∞ ∫ − ∞ + ∞ ( x 1 − μ ξ ( t 1 ) ) ( x 2 − μ ξ ( t 2 ) ) f ξ t 1 ξ t 2 ( x 1 , x 2 ) d x 1 d x 2 C_{\xi\xi}(t_1,t_2)=E([\xi(t_1)-\mu_{\xi}(t_1)][\xi(t_2)-\mu_{\xi}(t_2)]) \\ = \sum_{i}\sum_{j}(x_i-\mu_{\xi}(t_1))(x_j-\mu_{\xi}(t_2))P(\xi(t_1)=x_i,\xi(t_2)=x_j) \\ =\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}(x_1-\mu_{\xi}(t_1))(x_2-\mu_{\xi}(t_2))f_{\xi_{t_1}\xi_{t_2}}(x_1, x_2)dx_1dx_2 Cξξ(t1,t2)=E([ξ(t1)μξ(t1)][ξ(t2)μξ(t2)])=ij(xiμξ(t1))(xjμξ(t2))P(ξ(t1)=xi,ξ(t2)=xj)=++(x1μξ(t1))(x2μξ(t2))fξt1ξt2(x1,x2)dx1dx2
f ξ t 1 ξ t 2 ( x 1 , x 2 ) f_{\xi_{t_1}\xi_{t_2}}(x_1, x_2) fξt1ξt2(x1,x2)是隨機變數 ξ ( t 1 ) , ξ ( t 2 ) \xi(t_1),\xi(t_2) ξ(t1),ξ(t2)的聯合概率密度。
上式只是定義式,通過化簡可以得到:
C ξ ξ ( t 1 , t 2 ) = R ξ ξ ( t 1 , t 2 ) − μ ξ ( t 1 ) μ ξ ( t 2 ) C_{\xi\xi}(t_1,t_2)=R_{\xi\xi}(t_1,t_2)-\mu_{\xi}(t_1)\mu_{\xi}(t_2) Cξξ(t1,t2)=Rξξ(t1,t2)μξ(t1)μξ(t2)
t 1 = t 2 t_1=t_2 t1=t2時:
σ ξ 2 ( t 1 ) = C ξ ξ ( t 1 , t 1 ) = R ξ ξ ( t 1 , t 1 ) − μ ξ ( t 1 ) μ ξ ( t 1 ) \sigma_{\xi}^2(t_1)=C_{\xi\xi}(t_1,t_1)=R_{\xi\xi}(t_1,t_1)-\mu_{\xi}(t_1)\mu_{\xi}(t_1) σξ2(t1)=Cξξ(t1,t1)=Rξξ(t1,t1)μξ(t1)μξ(t1)
這就類似於 D ( X ) = E ( X 2 ) − ( E ( X ) ) 2 D(X)=E(X^2)-(E(X))^2 D(X)=E(X2)(E(X))2

4. 兩個或兩個以上隨機過程的聯合分佈和數字特徵

注意,這一節考慮的是兩個不同的隨機過程(注意和一個隨機過程的兩次不同的實現是有區別的)

相關文章