隨機過程學習筆記——概論
隨機過程學習筆記——概論
參考教材:陸大jin《隨機過程及其應用》
1. 隨機過程
1.1 基本概念
隨機過程是這樣一個過程,它不能用一個時間t的確定性函式表示,它在每個時刻的狀態是隨機的。
- 例子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去逐步測量,它得到的總是一條正弦曲線。
- 例子2:
某個餐廳,一個時刻到達的人數是隨機的,也就是隨機變數。如果研究一個上午各個時間的總人數,那麼這就是一個隨機過程。
1.2 描述隨機過程的方法
- 方法一:固定時間
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)
- 方法二:不固定時間,一次特定的實驗,隨著時間去逐步測量整個過程,因此這將得到一個確定的樣本函式,記為 ξ ( 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種型別。
-
離散引數離散型隨機過程。引數集(也就是時間)是離散的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)也是離散的。
比如伯努利過程:一個質點每過一秒就會等概率地在數軸上向左移動或者向右移動。 -
連續引數離散型隨機過程。引數集(也就是時間)是連續的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)是離散的。
比如脈衝數字通訊系統,傳送的訊號是脈寬為 T 0 T_0 T0的脈衝訊號,但是脈衝的幅度是隨機變數。 -
連續引數連續型隨機過程。引數集(也就是時間)是連續的,固定時間 t k t_k tk,隨機變數 ξ ( t k ) \xi(t_k) ξ(tk)也是連續的。
比如正弦波過程(例子1)和電晶體噪聲。 -
離散引數連續型隨機過程。引數集(也就是時間)是離散的,固定時間 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))=x1∗P(ξ(t1)=x1)+x2∗P(ξ(t1)=x2)+...+xn∗P(ξ(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))=i∑j∑xixjP(ξ(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)])=i∑j∑(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. 兩個或兩個以上隨機過程的聯合分佈和數字特徵
注意,這一節考慮的是兩個不同的隨機過程(注意和一個隨機過程的兩次不同的實現是有區別的)
相關文章
- 隨機過程複習筆記隨機筆記
- 李巨集毅機器學習課程筆記-1.機器學習概論機器學習筆記
- 機率論 學習筆記筆記
- 機器學習課程筆記機器學習筆記
- 期望 與 機率論 學習筆記筆記
- 學習C過程中的筆記系列-2筆記
- 群論學習筆記筆記
- 隨機過程(高斯隨機過程、譜分析、白噪聲)隨機
- 「馬爾可夫決策過程」學習筆記馬爾可夫筆記
- SAP Commerce Cloud 的構建過程學習筆記Cloud筆記
- 【方法論】機器學習演算法概覽機器學習演算法
- lua課程學習筆記筆記
- Python零基礎學習筆記(九)——隨機數Python筆記隨機
- Android學習過程的Cursor遊標填坑筆記Android筆記
- 吳恩達機器學習筆記 —— 1 緒論:初識機器學習吳恩達機器學習筆記
- 大資料 機器學習 演算法概論大資料機器學習演算法
- 資訊理論理論學習筆記筆記
- 林軒田機器學習基石課程學習筆記14 — Regularization機器學習筆記
- 林軒田機器學習基石課程學習筆記15 — Validation機器學習筆記
- 機器學習學習筆記機器學習筆記
- 學習筆記:數論分塊筆記
- 數論學習筆記 (1):整除筆記
- 達內課程學習筆記筆記
- CUDA學習筆記-1: CUDA程式設計概覽筆記程式設計
- PHP學習筆記系列:PHP生成器概覽PHP筆記
- 谷歌機器學習課程筆記(4)——降低損失谷歌機器學習筆記
- go lang學習筆記——channel機理及排程理解Go筆記
- 林軒田機器學習基石課程學習筆記12 — Nonlinear Transformation機器學習筆記ORM
- 林軒田機器學習基石課程學習筆記4 — Feasibility of Learning機器學習筆記
- 林軒田機器學習基石課程學習筆記9 — Linear Regression機器學習筆記
- 軒田機器學習基石課程學習筆記11 — Linear Models for Classification機器學習筆記
- 林軒田機器學習基石課程學習筆記13 — Hazard of Overfitting機器學習筆記
- 林軒田機器學習基石課程學習筆記6 — Theory of Generalization機器學習筆記
- 林軒田機器學習基石課程學習筆記7 — The VC Dimension機器學習筆記
- 林軒田機器學習基石課程學習筆記10 — Logistic Regression機器學習筆記
- 林軒田機器學習基石課程學習筆記3 — Types of Learning機器學習筆記
- 林軒田機器學習基石課程學習筆記8 — Noise and Error機器學習筆記Error
- 林軒田機器學習技法課程學習筆記9 — Decision Tree機器學習筆記