文章目录
随机过程的微分和积分过程1. 极限2. 连续 (这个连续特指的是均方意义下的连续)3. 导数(均方可导) 4. 微分变换 5. 积分
随机过程通过连续时间系统的分析1. 冲激响应法(用卷积乘法进行计算) 2. 频谱法
随机过程的微分和积分过程
1. 极限
随机变量序列的极限:
随机变量序列
{
X
n
}
\{X_n\}
{Xn} 依概率收敛于随机变量
X
X
X
lim
n
→
∞
X
n
=
P
X
\lim_{n\rightarrow\infin}X_n \xlongequal{P} X
n→∞limXnP
X随机变量序列
{
X
n
}
\{X_n\}
{Xn} 依均方收敛于随机变量
X
X
X(一般主要用依均方收敛的概念,并且若已经依均方收敛则必然依概率收敛)
lim
n
→
∞
E
{
∣
X
n
−
X
∣
2
}
=
0
\lim_{n\rightarrow\infin}E\{|X_n-X|^2\} = 0
n→∞limE{∣Xn−X∣2}=0 随机过程的极限:
随机过程
X
(
t
)
X(t)
X(t) 依概率收敛于随机变量
X
X
X
lim
t
→
t
0
X
(
t
)
=
P
X
\lim_{t\rightarrow t_0}X(t) \xlongequal{P} X
t→t0limX(t)P
X随机过程
{
X
n
}
\{X_n\}
{Xn} 依均方收敛于随机变量
X
X
X
lim
t
→
t
0
E
{
∣
X
(
t
)
−
X
∣
2
}
=
0
\lim_{t\rightarrow t_0}E\{|X(t)-X|^2\} = 0
t→t0limE{∣X(t)−X∣2}=0
2. 连续 (这个连续特指的是均方意义下的连续)
lim
Δ
t
→
0
E
{
∣
X
(
t
+
Δ
t
)
−
X
(
t
)
∣
2
}
=
0
\lim_{\Delta t\rightarrow 0}E\{|X(t+\Delta t)-X(t)|^2\} = 0
Δt→0limE{∣X(t+Δt)−X(t)∣2}=0 记作:
均方连续的充要条件:
R
X
(
t
1
,
t
2
)
R_X(t_1,t_2)
RX(t1,t2) 在
t
1
=
t
2
t_1 = t_2
t1=t2 时连续。平稳过程的均方连续:
lim
τ
→
0
R
X
(
τ
)
=
R
X
(
0
)
\lim_{\tau \rightarrow 0}R_X(\tau)=R_X(0)
τ→0limRX(τ)=RX(0) 就是说平稳过程在零点的连续性决定了整个过程的连续性。在均方连续的条件下,求极限和求均值可以互换次序。
3. 导数(均方可导)
随机过程
X
(
t
)
X(t)
X(t) 均方可导的充要条件:相关函数
R
X
(
t
1
,
t
2
)
R_X(t_1,t_2)
RX(t1,t2) 在自变量
t
1
=
t
2
=
t
t_1=t_2=t
t1=t2=t 的时候存在二阶混合偏导数。平稳随机过程均方可导的条件:相关函数
R
X
(
τ
)
R_X(\tau)
RX(τ) 当自变量
τ
=
0
\tau=0
τ=0 时要存在二阶偏导数
R
X
′
′
(
0
)
R_X''(0)
RX′′(0)
4. 微分变换
随机过程导数的运算与均值的运算可以交换次序:
E
[
d
X
(
t
)
d
t
]
=
d
E
[
X
(
t
)
]
d
t
E[\dfrac{dX(t)}{dt}]=\dfrac{dE[X(t)]}{dt}
E[dtdX(t)]=dtdE[X(t)]随机过程与其均方导数的互相关函数:
R
X
X
˙
(
t
1
,
t
2
)
=
∂
R
X
(
t
1
,
t
2
)
∂
t
2
R
X
˙
X
(
t
1
,
t
2
)
=
∂
R
X
(
t
1
,
t
2
)
∂
t
1
\begin{matrix} R_{X\dot{X}}(t_1,t_2)=\dfrac{\partial R_X(t_1,t_2)}{\partial t_2}\\ ~~\\ R_{\dot{X}X}(t_1,t_2)=\dfrac{\partial R_X(t_1,t_2)}{\partial t_1} \end{matrix}
RXX˙(t1,t2)=∂t2∂RX(t1,t2) RX˙X(t1,t2)=∂t1∂RX(t1,t2) 实际上这里可以看成协方差的底对谁求导,协方差本身就对相应的时间进行求导。对平稳随机过程来说(记
τ
=
t
2
−
t
1
\tau = t_2-t_1
τ=t2−t1):
R
X
X
˙
(
τ
)
=
∂
R
X
(
τ
)
∂
τ
R
X
˙
X
(
τ
)
=
−
∂
R
X
(
τ
)
∂
τ
\begin{matrix} R_{X\dot{X}}(\tau)=\dfrac{\partial R_X(\tau)}{\partial \tau}\\ ~~\\ R_{\dot{X}X}(\tau)=-\dfrac{\partial R_X(\tau)}{\partial \tau} \end{matrix}
RXX˙(τ)=∂τ∂RX(τ) RX˙X(τ)=−∂τ∂RX(τ) 注意此时会有一个结论:平稳过程
X
(
t
)
X(t)
X(t) 在任一固定时刻
t
1
t_1
t1 的取值
X
(
t
1
)
X(t_1)
X(t1) 与其同一时刻的导数取值
X
˙
(
t
1
)
\dot{X}(t_1)
X˙(t1) 不相关。随机过程与其均方导数的自相关函数:
R
X
˙
(
t
1
,
t
2
)
=
∂
R
X
X
˙
(
t
1
,
t
2
)
∂
t
1
R
X
˙
(
t
1
,
t
2
)
=
∂
R
X
˙
X
(
t
1
,
t
2
)
∂
t
2
R
X
˙
(
t
1
,
t
2
)
=
∂
R
X
(
t
1
,
t
2
)
∂
t
1
∂
t
2
\begin{matrix} R_{\dot{X}}(t_1,t_2)&=\dfrac{\partial R_{X\dot{X}}(t_1,t_2)}{\partial t_1}\\ ~~\\ R_{\dot{X}}(t_1,t_2)&=\dfrac{\partial R_{\dot{X}X}(t_1,t_2)}{\partial t_2}\\ ~~\\ R_{\dot{X}}(t_1,t_2)&=\dfrac{\partial R_{X}(t_1,t_2)}{\partial t_1\partial t_2} \end{matrix}
RX˙(t1,t2) RX˙(t1,t2) RX˙(t1,t2)=∂t1∂RXX˙(t1,t2)=∂t2∂RX˙X(t1,t2)=∂t1∂t2∂RX(t1,t2)对平稳随机过程来说(记
τ
=
t
2
−
t
1
\tau = t_2-t_1
τ=t2−t1):
R
X
˙
(
τ
)
=
−
d
2
d
τ
2
R
X
(
τ
)
R_{\dot{X}}(\tau) = -\dfrac{d^2}{d\tau^2}R_X(\tau)
RX˙(τ)=−dτ2d2RX(τ)导数的功率谱密度:
G
X
˙
(
ω
)
=
ω
2
G
X
(
ω
)
G_{\dot{X}}(\omega)=\omega^2G_X(\omega)
GX˙(ω)=ω2GX(ω)
这里有一种简单记忆的形式
X
→
X
X
X\rightarrow XX
X→XX,
X
˙
→
X
˙
X
˙
\dot{X}\rightarrow \dot{X}\dot{X}
X˙→X˙X˙,然后再按照对底变换和对对应时间求导等价的原则进行变换即可。
5. 积分
这里要说明一下积分的定义,因为感觉这个和高数中对积分定义的差别还是比较大的。
lim
Δ
t
k
→
0
E
{
∣
∑
k
=
1
n
X
(
t
k
)
Δ
t
k
−
∫
a
b
X
(
t
)
d
t
∣
2
}
=
0
\lim_{\Delta t_k\rightarrow 0}E\{|\sum_{k=1}^{n}X(t_k)\Delta t_k - \int_a^bX(t)dt|^2 \} = 0
Δtk→0limE{∣k=1∑nX(tk)Δtk−∫abX(t)dt∣2}=0 这实际上是一个依均方收敛的定义,实际上是先对每一个样本函数进行横向的积分最后全部整合到纵向,然后再在纵向上使用依均方收敛的定义。积分变换与微分变换是十分相似的。
随机过程通过连续时间系统的分析
1. 冲激响应法(用卷积乘法进行计算)
对于一般的随机过程:
R
X
Y
(
t
1
,
t
2
)
=
R
X
(
t
1
,
t
2
)
⊗
h
(
t
2
)
R
Y
X
(
t
1
,
t
2
)
=
R
X
(
t
1
,
t
2
)
⊗
h
(
t
1
)
R
Y
(
t
1
,
t
2
)
=
R
X
(
t
1
,
t
2
)
⊗
h
(
t
1
)
⊗
h
(
t
2
)
\begin{matrix} R_{XY}(t_1,t_2)=R_X(t_1,t_2)\otimes h(t_2)\\ ~~\\ R_{YX}(t_1,t_2)=R_X(t_1,t_2)\otimes h(t_1)\\ ~~\\ R_{Y}(t_1,t_2)=R_X(t_1,t_2)\otimes h(t_1)\otimes h(t_2) \end{matrix}
RXY(t1,t2)=RX(t1,t2)⊗h(t2) RYX(t1,t2)=RX(t1,t2)⊗h(t1) RY(t1,t2)=RX(t1,t2)⊗h(t1)⊗h(t2)对于平稳的随机过程:
R
X
Y
(
τ
)
=
h
(
τ
)
⊗
h
(
τ
)
R
Y
X
(
τ
)
=
R
X
(
τ
)
⊗
h
(
−
τ
)
R
Y
(
τ
)
=
R
X
(
τ
)
⊗
h
(
τ
)
⊗
h
(
−
τ
)
\begin{matrix} R_{XY}(\tau)=h(\tau)\otimes h(\tau)\\ ~~\\ R_{YX}(\tau)=R_X(\tau)\otimes h(-\tau)\\ ~~\\ R_{Y}(\tau)=R_X(\tau)\otimes h(\tau)\otimes h(-\tau) \end{matrix}
RXY(τ)=h(τ)⊗h(τ) RYX(τ)=RX(τ)⊗h(−τ) RY(τ)=RX(τ)⊗h(τ)⊗h(−τ)实际上从记忆公式的角度来说,这个和上面求自相关函数的导数是十分相似的,将
X
→
X
X
X\rightarrow XX
X→XX,卷积上时间和变换对应的底是等价的。由于卷积的计算相对比较麻烦,因此一般对于平稳过程的求解都会采用频谱法进行计算而不会采用冲激响应法进行计算。
2. 频谱法
注意频谱法要求线性系统的输入
X
(
t
)
X(t)
X(t) 和
Y
(
t
)
Y(t)
Y(t) 均为平稳过程:一般的求解过程如下:
首先知道过程
X
(
t
)
X(t)
X(t) 的相关函数
R
X
(
τ
)
R_X(\tau)
RX(τ)。然后使用傅里叶变换求解该相关函数对应的频谱密度函数
G
X
(
ω
)
G_X(\omega)
GX(ω)。采用公式
G
Y
(
ω
)
=
∣
H
(
ω
)
∣
2
G
X
(
ω
)
G_Y(\omega) = |H(\omega)|^2G_X(\omega)
GY(ω)=∣H(ω)∣2GX(ω)
|