1、随机过程
通俗的理解:随机变量是单个变量X是个随机值,比如扔硬币,正面还是反面是个随机值,它的数字特征有均值、方差。而过程则增加一个因素T,通常是时间,如时刻t抛一次硬币,在s时刻抛一次硬币,则可以研究两个时刻抛硬币的均值、方差及相关性。

从图上看,单个随机变量用一个二维平面表示,单个变量的随机过程需要三维平面表示。在每个时间t,它就是一个随机变量。
随机过程的均值函数:mX(t),此记法表示,随机变量X,在t对应的均值,因为t是个变量,所以不同t会对应不同的均值,所以称之为函数;m是mean的缩写。所以mX(t)=E[X(t)],对应于上图:在t时刻切片下的变量的图形均值。
随机过程的方差函数:DX(t),如均值函数一样,是随机变量X,在t对应的方差,因为t是个变量,所以不同t会对应不同的方差,所以称之为方差函数;DX(t)=E[(X(t)−EX(t))2],计算方法和随机变量的方差一样。
随机过程的协方差函数:BX(s,t),研究随机变量X,在s时间切片下X(s)和在t时间切片下X(t)的相互关系。如果把X(s)看为变量Y,X(t)看为变量Z,也就是相当于Cov(Y,Z),所以BX(s,t)=E[{X(s)−mX(s)}{X(t)−mX(t)}]=E[X(s)X(t)]−mX(s)mX(t)。像协方差公式展开后,把E(YZ)叫着相关矩,这里则是相关函数:RX(s,t)=E[X(s)X(t)]。
上面研究是X本身在不同时刻的相关性,也称为自协方差函数、自相关函数,对应如果是X在t时间,Y在s时间的相互关系:
互协方差函数:BXY(s,t)=E[{X(s)−mX(s)}{Y(t)−mY(t)}]
互相关函数:RXY(s,t)=E[X(s)Y(t)]
两函数和W(t)=X(t)+Y(t)的相关性:
他们的均值函数为mW(t)=E[X(t)+Y(t)]=mX(t)+mY(t)
相关函数:RW(t,s)=E[(X(t)+Y(t))(X(s)+Y(s)]=E[X(t)X(s)+X(t)Y(s)+Y(t)X(s)+Y(t)Y(s)]=RX(t,s)+RXY(t,s)+RYX(t,s)+RY(t,s)当X、Y不相关,则RW(t,s)=RX(t,s)+RY(t,s)
2、复随机过程:
{Xt,t∈T}{Yt,t∈T}是取实值的两个随机过程,则Zt=Xt+iYt为复随机过程。
均值函数:mZ(t)=EXt+iEYt
方差函数:如按实变量公式:DZ(t)=E[(Zt−EZt))2]=E{[(Xt+iYt)−E(Xt+iYt)]2},展开:=E{(Xt−EXt)2+2i(Xt−EXt)(Yt−EYt)−(Yt−EYt)2},可以看出里面公式有可能为负数,与方差的定义不符。在此可以知道,复数不能直接展开。因为方差表示变量的振幅,对于复数,则用模表示距离:DZ(t)=E{∣(Xt+iYt)−E(Xt+iYt)∣2},(此方式也是实变量的计算方法,只是实函数并不受影响而已,所以省了模的符号)。DZ(t)=E{[(Xt+iYt)−E(Xt+iYt)][(Xt+iYt)−E(Xt+iYt)]}=E{[(Xt−EXt)+i(Yt−EYt)][(Xt−EXt)+i(Yt−EYt)]}=E{[(Xt−EXt)+i(Yt−EYt)][(Xt−EXt)−i(Yt−EYt)]}=E{(Xt−EXt)2+(Yt−EYt)2}
协方差函数:BZ(s,t)=E{[Zs−E(Zs)][Zt−E(Zt)]},此方式保证当s=t时,即是复变量的方差。但可能还是会想,不同时刻,另一个为什么是取它的共轭?
首先复数可看成二维空间,由X(实部)、Y(虚部)轴共同确定。把复数转成Zs=∣Zs∣eiθs形式,,可看出相互关系不再仅仅是比大小,而是比方向(相位)的关系和模长的关系。如果BZ(s,t)=E{∣Zs∣eiθs∗∣Zt∣eiθt}=E[∣Zs∣∣Zt∣eiθs+iθt](假设均值都为0)。此式中模的积能表达模关系,但相位是绝对关系,如果都增加一个常数,协方差也会跟着改变;和协方差要减去均值去抵消绝对数值关系的初衷有点出入,而且在平稳的研究过程中,仅与τ=t−s有关。而如果是i(θs−θt)则其相对性,就具有很好的研究性质。因此∣Zs∣e−iθs=Zs(欧拉公式e−iθs=cos(−θ)+isin(−θ)=cos(θ)−isin(θ))。
相关函数:RZ(s,t)=E(ZsZˉt)。协方差函数展开:BZ(s,t)=E[ZsZˉt]−E[Zs]E[Zˉt],前一项即是相关函数。
同理上面的可称为自协方差函数和自相关函数,那对于BXY(s,t)=E{[Xs−E(Xs)][Yt−E(Yt)]},RXY(s,t)=E[XsYˉt]称为互协方差函数和互相关函数
3、计算正态随机过程的数字特征
例:正态随机过程X(t)=Y+Zt,(t>0),其中Y,Z是相互独立的N(0,1)随机变量,则X的一维概率密度:
EX(t)=E[Y+Zt]=EY+EZt=0
E[(X(t)−μX(t))2]=E[(Y+Zt)2]=EY2+E(2YZt)+E(Z2t)=1+t2
所以:
一维的概率密度:fX(t)(x)=2π(1+t2)1e−2(1+t2)x2
对于此例是正态随机过程,如果非正态,则需要通过分布来推导:P(X(t)≤x)=∫−∞∞P(Y=y)P(Z≤tx−y∣Y=y)dy
由Y,Z相互独立:fX(t)(x)=∫−∞∞fY(y)∫−∞tx−yfz(z)dzdy,对X求导fX(t)(x)=∫−∞∞fY(y)fz(tx−y)t1dy
此例代入Y,Z的概率密度计算:fX(t)(x)=∫−∞∞2π1e−2y22π1e−2t2(x−y)2t1dy,通过合并e的系数,然后对y进行配方fX(t)(x)=2πt1∫−∞∞exp[−2t21+t2(y−1+t2x)2−2(1+t2)x2]dy,再利用∫−∞∞e−2y2dy=2π去掉积分后,即可得fX(t)(x)=2π(1+t2)1e−2(1+t2)x2
二维的概率密度:
先求出协方差BX(s,t)==E[{X(s)−mX(s)}{X(t)−mX(t)}]=E[(Y+Zs)(Y+Zt)]=1+st
代入多维正态概率密度公式去求出:fX(s,t)(xs,xt)