各态历经性-随机过程

💡解答
要分析 X={Xt,t∈T}X = \{X_t, t \in T\}X={Xt,t∈T} 和 Y={Yt,t∈T}Y = \{Y_t, t \in T\}Y={Yt,t∈T}的联合平稳性及相关性,需从联合宽平稳过程的定义(均值函数为常数、自相关函数仅与时间差有关、互相关函数仅与时间差有关)入手。
一、分析均值函数
-
对Xt=cos(ωt+Φ)X_t = \cos(\omega t + \Phi)Xt=cos(ωt+Φ),由于 Φ\PhiΦ 服从 [0,2π)[0, 2\pi)[0,2π)上的均匀分布,即fΦ(ϕ)=12π,ϕ∈[0,2π)f_\Phi(\phi) = \frac{1}{2\pi}, \phi \in [0, 2\pi)fΦ(ϕ)=2π1,ϕ∈[0,2π),则:
μX(t)=E[Xt]=∫02πcos(ωt+ϕ)⋅12πdϕ \mu_X(t) = E[X_t] = \int_{0}^{2\pi} \cos(\omega t + \phi) \cdot \frac{1}{2\pi} d\phi μX(t)=E[Xt]=∫02πcos(ωt+ϕ)⋅2π1dϕ
令u=ωt+ϕu = \omega t + \phiu=ωt+ϕ,则 du=dϕdu = d\phidu=dϕ,积分区间变为 [ωt,ωt+2π)[\omega t, \omega t + 2\pi)[ωt,ωt+2π),积分结果为:
μX(t)=12π∫ωtωt+2πcosu du=12π[sinu]ωtωt+2π=0 \mu_X(t) = \frac{1}{2\pi} \int_{\omega t}^{\omega t + 2\pi} \cos u \, du = \frac{1}{2\pi} \left[ \sin u \right]_{\omega t}^{\omega t + 2\pi} = 0 μX(t)=2π1∫ωtωt+2πcosudu=2π1[sinu]ωtωt+2π=0 -
对 Yt=sin(ωt+Φ)Y_t = \sin(\omega t + \Phi)Yt=sin(ωt+Φ),同理:
μY(t)=E[Yt]=∫02πsin(ωt+ϕ)⋅12πdϕ \mu_Y(t) = E[Y_t] = \int_{0}^{2\pi} \sin(\omega t + \phi) \cdot \frac{1}{2\pi} d\phi μY(t)=E[Yt]=∫02πsin(ωt+ϕ)⋅2π1dϕ
令 u=ωt+ϕu = \omega t + \phiu=ωt+ϕ,积分结果为:
μY(t)=12π∫ωtωt+2πsinu du=12π[−cosu]ωtωt+2π=0 \mu_Y(t) = \frac{1}{2\pi} \int_{\omega t}^{\omega t + 2\pi} \sin u \, du = \frac{1}{2\pi} \left[ -\cos u \right]_{\omega t}^{\omega t + 2\pi} = 0 μY(t)=2π1∫ωtωt+2πsinudu=2π1[−cosu]ωtωt+2π=0
因此,μX(t)=0\mu_X(t) = 0μX(t)=0,μY(t)=0\mu_Y(t) = 0μY(t)=0,均为与 ttt 无关的常数。
二、分析自相关函数
(1)XXX的自相关函数RX(t,t+τ)R_X(t, t + \tau)RX(t,t+τ)
RX(t,t+τ)=E[XtXt+τ]=E[cos(ωt+Φ)cos(ω(t+τ)+Φ)]
\begin{align*}
R_X(t, t + \tau) &= E[X_t X_{t + \tau}] \\
&= E\left[ \cos(\omega t + \Phi) \cos(\omega (t + \tau) + \Phi) \right]
\end{align*}
RX(t,t+τ)=E[XtXt+τ]=E[cos(ωt+Φ)cos(ω(t+τ)+Φ)]
利用三角恒等式 cosAcosB=12[cos(A−B)+cos(A+B)]\cos A \cos B = \frac{1}{2}[\cos(A - B) + \cos(A + B)]cosAcosB=21[cos(A−B)+cos(A+B)],令 A=ωt+ΦA = \omega t + \PhiA=ωt+Φ,B=ω(t+τ)+ΦB = \omega (t + \tau) + \PhiB=ω(t+τ)+Φ,则 A−B=−ωτA - B = -\omega \tauA−B=−ωτ,A+B=2ωt+ωτ+2ΦA + B = 2\omega t + \omega \tau + 2\PhiA+B=2ωt+ωτ+2Φ,代入得:
RX(t,t+τ)=E[12cos(−ωτ)+12cos(2ωt+ωτ+2Φ)]=12cos(ωτ)+12E[cos(2ωt+ωτ+2Φ)]
\begin{align*}
R_X(t, t + \tau) &= E\left[ \frac{1}{2} \cos(-\omega \tau) + \frac{1}{2} \cos(2\omega t + \omega \tau + 2\Phi) \right] \\
&= \frac{1}{2} \cos(\omega \tau) + \frac{1}{2} E\left[ \cos(2\omega t + \omega \tau + 2\Phi) \right]
\end{align*}
RX(t,t+τ)=E[21cos(−ωτ)+21cos(2ωt+ωτ+2Φ)]=21cos(ωτ)+21E[cos(2ωt+ωτ+2Φ)]
对第二项,由于 ( \Phi ) 均匀分布,( 2\Phi ) 也服从 ( [0, 4\pi) ) 上的均匀分布,积分
∫02πcos(2ωt+ωτ+2ϕ)⋅12πdϕ=0
\int_{0}^{2\pi} \cos(2\omega t + \omega \tau + 2\phi) \cdot \frac{1}{2\pi} d\phi = 0
∫02πcos(2ωt+ωτ+2ϕ)⋅2π1dϕ=0
因此:
RX(t,t+τ)=12cos(ωτ)
R_X(t, t + \tau) = \frac{1}{2} \cos(\omega \tau)
RX(t,t+τ)=21cos(ωτ)
可见,RXR_XRX仅与时间差 τ\tauτ有关。
(2)YYY的自相关函数 RY(t,t+τ)R_Y(t, t + \tau)RY(t,t+τ)
类似地,利用三角恒等式sinAsinB=12[cos(A−B)−cos(A+B)]\sin A \sin B = \frac{1}{2}[\cos(A - B) - \cos(A + B)]sinAsinB=21[cos(A−B)−cos(A+B)],计算 E[YtYt+τ]E[Y_t Y_{t + \tau}]E[YtYt+τ]:
RY(t,t+τ)=E[sin(ωt+Φ)sin(ω(t+τ)+Φ)]=E[12cos(−ωτ)−12cos(2ωt+ωτ+2Φ)]=12cos(ωτ)
\begin{align*}
R_Y(t, t + \tau) &= E\left[ \sin(\omega t + \Phi) \sin(\omega (t + \tau) + \Phi) \right] \\
&= E\left[ \frac{1}{2} \cos(-\omega \tau) - \frac{1}{2} \cos(2\omega t + \omega \tau + 2\Phi) \right] \\
&= \frac{1}{2} \cos(\omega \tau)
\end{align*}
RY(t,t+τ)=E[sin(ωt+Φ)sin(ω(t+τ)+Φ)]=E[21cos(−ωτ)−21cos(2ωt+ωτ+2Φ)]=21cos(ωτ)
(第二项积分结果为 000,推导同XXX 的自相关)
因此,RYR_YRY仅与时间差 τ\tauτ有关**。
三、分析互相关函数 ( R_{XY}(t, t + \tau) )
RXY(t,t+τ)=E[XtYt+τ]=E[cos(ωt+Φ)sin(ω(t+τ)+Φ)]
\begin{align*}
R_{XY}(t, t + \tau) &= E[X_t Y_{t + \tau}] \\
&= E\left[ \cos(\omega t + \Phi) \sin(\omega (t + \tau) + \Phi) \right]
\end{align*}
RXY(t,t+τ)=E[XtYt+τ]=E[cos(ωt+Φ)sin(ω(t+τ)+Φ)]
利用三角恒等式 cosAsinB=12[sin(A+B)+sin(B−A)]\cos A \sin B = \frac{1}{2}[\sin(A + B) + \sin(B - A)]cosAsinB=21[sin(A+B)+sin(B−A)],令 A=ωt+ΦA = \omega t + \PhiA=ωt+Φ,B=ω(t+τ)+ΦB = \omega (t + \tau) + \PhiB=ω(t+τ)+Φ,则
A+B=2ωt+ωτ+2ΦA + B = 2\omega t + \omega \tau + 2\PhiA+B=2ωt+ωτ+2Φ,B−A=ωτB- A=\omega\tauB−A=ωτ,代入得:
RXY(t,t+τ)=E[12sin(ωτ)+12sin(2ωt+ωτ+2Φ)]=12sin(ωτ)+12E[sin(2ωt+ωτ+2Φ)]
\begin{align*}
R_{XY}(t, t + \tau) &= E\left[ \frac{1}{2} \sin(\omega \tau) + \frac{1}{2} \sin(2\omega t + \omega \tau + 2\Phi) \right] \\
&= \frac{1}{2} \sin(\omega \tau) + \frac{1}{2} E\left[ \sin(2\omega t + \omega \tau + 2\Phi) \right]
\end{align*}
RXY(t,t+τ)=E[21sin(ωτ)+21sin(2ωt+ωτ+2Φ)]=21sin(ωτ)+21E[sin(2ωt+ωτ+2Φ)]
对第二项,同理积分结果为 ( 0 ),因此:
RXY(t,t+τ)=12sin(ωτ)
R_{XY}(t, t + \tau) = \frac{1}{2} \sin(\omega \tau)
RXY(t,t+τ)=21sin(ωτ)
可见,RXYR_{XY}RXY仅与时间差 τ\tauτ有关。
四、结论
- 联合平稳性:由于 ( X ) 和 ( Y ) 的均值函数为常数,且自相关函数、互相关函数均仅与时间差 τ\tauτ有关,因此 XXX 和YYY 联合宽平稳。
- 相关性:互相关函数RXY(τ)=12sin(ωτ)R_{XY}(\tau) = \frac{1}{2} \sin(\omega \tau)RXY(τ)=21sin(ωτ)不恒为 000,因此XXX 和 YYY 相关。
(注:若进一步分析“不相关”需看互协方差
CXY(τ)=RXY(τ)−μXμYC_{XY}(\tau) = R_{XY}(\tau) - \mu_X \mu_YCXY(τ)=RXY(τ)−μXμY,
由于 μX=μY=0\mu_X = \mu_Y = 0μX=μY=0,故 (CXY(τ)=RXY(τ)( C_{XY}(\tau) = R_{XY}(\tau)(CXY(τ)=RXY(τ),不恒为 ( 0 ),因此相关。)

解答:各态历经性分析
各态历经性需验证时间平均是否以概率 1 等于统计平均(包括均值和自相关函数)。
步骤1:判断是否为宽平稳过程
首先验证 {X(t)}\{X(t)\}{X(t)} 是宽平稳过程:
-
均值函数:
μX(t)=E[X(t)]=E[X]=1⋅13+2⋅13+3⋅13=2 \mu_X(t) = E[X(t)] = E[X] = 1 \cdot \frac{1}{3} + 2 \cdot \frac{1}{3} + 3 \cdot \frac{1}{3} = 2 μX(t)=E[X(t)]=E[X]=1⋅31+2⋅31+3⋅31=2
(与 ttt 无关) -
自相关函数:
RX(t1,t2)=E[X(t1)X(t2)]=E[X2]=12⋅13+22⋅13+32⋅13=143 R_X(t_1, t_2) = E[X(t_1)X(t_2)] = E[X^2] = 1^2 \cdot \frac{1}{3} + 2^2 \cdot \frac{1}{3} + 3^2 \cdot \frac{1}{3} = \frac{14}{3} RX(t1,t2)=E[X(t1)X(t2)]=E[X2]=12⋅31+22⋅31+32⋅31=314
(仅与时间差 τ=t2−t1\tau = t_2 - t_1τ=t2−t1 有关,实际上与 τ\tauτ 无关)
因此,{X(t)}\{X(t)\}{X(t)} 是宽平稳过程。
步骤2:计算时间平均
对于随机过程 X(t)=XX(t) = XX(t)=X(不随 ttt 变化),其时间平均(均值和自相关)为:
-
时间均值:
⟨X(t)⟩=limT→∞12T∫−TTX(t) dt=limT→∞12T∫−TTX dt=X \langle X(t) \rangle = \lim_{T \to \infty} \frac{1}{2T} \int_{-T}^{T} X(t) \, dt = \lim_{T \to \infty} \frac{1}{2T} \int_{-T}^{T} X \, dt = X ⟨X(t)⟩=T→∞lim2T1∫−TTX(t)dt=T→∞lim2T1∫−TTXdt=X -
时间自相关:
⟨X(t)X(t+τ)⟩=limT→∞12T∫−TTX2 dt=X2 \langle X(t)X(t+\tau) \rangle = \lim_{T \to \infty} \frac{1}{2T} \int_{-T}^{T} X^2 \, dt = X^2 ⟨X(t)X(t+τ)⟩=T→∞lim2T1∫−TTX2dt=X2
步骤3:验证各态历经性
各态历经性要求“时间平均以概率 1 等于统计平均”:
-
均值的各态历经性:
统计均值为 μX=2\mu_X = 2μX=2,但时间均值 ⟨X(t)⟩=X\langle X(t) \rangle = X⟨X(t)⟩=X(XXX 可取 1、2、3)。
显然 XXX 不恒等于 2(例如 P(X=1)=13P(X=1)=\frac{1}{3}P(X=1)=31),因此 ⟨X(t)⟩≠μX\langle X(t) \rangle \neq \mu_X⟨X(t)⟩=μX(以正概率不相等)。 -
自相关的各态历经性:
统计自相关为 RX(τ)=143R_X(\tau) = \frac{14}{3}RX(τ)=314,但时间自相关 ⟨X(t)X(t+τ)⟩=X2\langle X(t)X(t+\tau) \rangle = X^2⟨X(t)X(t+τ)⟩=X2(X2X^2X2 可取 1、4、9)。
X2X^2X2 不恒等于 143\frac{14}{3}314,因此自相关也不满足各态历经性。
结论
随机过程 {X(t),−∞<t<+∞}\{X(t), -\infty < t < +\infty\}{X(t),−∞<t<+∞}
不具有均值的各态历经性,也不具有自相关函数的各态历经性。
