例 平稳高斯过程

高斯过程由期望函数和协方差函数决定

由严格平稳,有:

μX(s+h)=μX(s)Cx(s,t)=Cx(s+h,t+h)\mu_X(s+h)=\mu_X(s) \\ C_x(s,t)=C_x(s+h,t+h)

tT\forall t \in T,有:

μx(t)=μx(0)T=[0,T0]   Cx(t,s)=Cx~(ts)其中Cx~是一元函数\mu_x(t)=\mu_x(0) \\ T=[0,T_0] \ \ \ C_x(t,s)=\tilde{C_x}(|t-s|) 其中\tilde{C_x}是一元函数

对高斯过程严格平稳,意味着:

期望函数为常数,协方差函数仅依赖于时间间隔ts|t-s|

一般来说,如果一个过程XX满足上述性质,就称之为**(宽)平稳过程二次平稳过程**。

平稳增量

有时过程本身不满足平稳的条件,但其增量是平稳的。令X=(Xt,tT)X=(X_t,t \in T)是一个随机过程,$T\subset R $是一个区间,如果对于所有的 t>s, t,sTt>s,\ t,s\in T 使得t+h,s+hTt+h,s+h\in Thh,有:

XtXs=dXt+hXs+hX_t-X_s \overset{d}{=} X_{t+h}-X_{s+h}

则称XX是具有平稳增量的。

独立增量

如果对于任意的tiTt_i\in T,且t1<t2<...<tn, n1t_1<t_2<...<t_n ,\ n\geq 1,有Xt2Xt1,Xt3Xt2,...,XtnXtn1X_{t_2}-X_{t_1},X_{t_3}-X_{t_2},...,X_{t_n}-X_{t_{n-1}} 是相互独立的随机变量,则称XX具有独立增量。即不重合的时间区间上的增量独立。

平稳独立增量

如果随机过程XX既具有平稳增量,又具有独立增量,则称XX具有平稳独立增量LevyLevy过程)

两个代表性的LevyLevy过程:齐次泊松过程布朗运动

齐次泊松过程

随机过程(Nt , t[0,)N_t\ ,\ t\in[0,\infty))称为齐次泊松过程或带有强度为λ\lambda或速率λ>0\lambda >0 的泊松过程,如果满足:

  1. 过程从00出发 (No=0)(N_o=0)

  2. 过程具有平稳独立增量

  3. t>0,NtP(λt)\forall t>0,N_t \sim P(\lambda t) (参数为λt\lambda t的泊松分布),即:

P(Nt=k)=(λt)kk!eλt   ,k=0,1,2... P(N_t=k)=\frac{(\lambda t)^k}{k!}e^{-\lambda t} \ \ \ ,k=0,1,2...

期望、方差

E(Nt)=λtVar(Nt)=λtE(N_t)=\lambda t \\ Var(N_t)=\lambda t

自相关系数

RN(t,s)=E[NtNs]R_N(t,s)=E[N_tN_s]

tst \geq s时:

RN(t,s)=E[(NtNs+Ns)Ns]=E[(NtNs)Ns]E[Ns2]=E[NtNs]E[Ns]E[Ns2]=λ(ts)(λs)+λs+(λs)2=λ2ts+λsR_N(t,s)=E[(N_t-N_s+N_s)N_s] \\ =E[(N_t-N_s)N_s]-E[N_s^2] \\ =E[N_t-N_s]E[N_s]-E[N_s^2] \\ = \lambda(t-s)*(\lambda s)+\lambda s+(\lambda s)^2 \\ = \lambda^2 ts+\lambda s

同理,当$t \leq s $时:

RN(t,s)=λ2ts+λtR_N(t,s)=\lambda^2 ts+\lambda t

因此:

RN(t,s)=λ2ts+λmin(t,s)R_N(t,s)=\lambda^2 ts+\lambda min(t,s)

协方差函数

CN(t,s)=RN(t,s)E[Nt][Ns]=λ min(t,s)C_N(t,s)=R_N(t,s)-E[N_t][N_s]=\lambda \ min(t,s)

齐次泊松过程等价定义

Nt=${n:Tnt},t0N_t= \$ \{n: T_n\leq t \} ,t \geq 0

其中$A\$ A表示给定集合AA中的元素个数。

Tn=Y1+Y2+...+YnT_n=Y_1+Y_2+...+Y_n

YiY_i 是相互独立的服从指数分布Exp(λ)Exp(\lambda)的独立同分布的序列。

YiY_i是第i1i-1次跳跃到第ii次跳跃的时间间隔,TnT_n是发生第nn次跳跃的时刻。

一个计数过程,如果事件发生的间隔时间是独立同分布的指数随机变量,这个计数过程就是泊松过程。

*练习

{Nt,t0}\{N_t,t\geq0 \}是强度为λ\lambda的泊松过程,令Mt=0tNsds,t0M_t=\int_{0}^t N_s ds ,t \geq 0,计算E[Mt]E[M_t]Var(Mt)Var(M_t)

E[Mt]=0tE[Ns]ds=0tλs ds=0.5λt2E[M_t]=\int_{0}^t E[N_s] ds =\int_{0}^t \lambda s \ ds =0.5\lambda t^2

Var(Mt)=E[Mt2]E2[Mt]Var(M_t)=E[M_t^2]-E^2[M_t]

E[Mt2]=E[(0tNsds)2]=E[(0tNsds)(0tNrdr)]=E[(0t0tNsNrdsdr)]=0t0tE[NsNr]dsdr=0t0t(λ2sr+λmin(s,r)) dsdr=14λ2t4+13λt3E[M_t^2]=E[(\int_{0}^t N_s ds)^2]=E[(\int_{0}^t N_s ds)(\int_{0}^t N_r dr)] \\ =E[(\int_{0}^t \int_{0}^t N_s N_r ds dr)] \\ =\int_{0}^t \int_{0}^t E[N_s N_r] ds dr \\=\int_{0}^t \int_{0}^t( \lambda^2 sr+\lambda min(s,r))\ ds dr \\ =\frac{1}{4}\lambda^2 t^4+\frac{1}{3}\lambda t^3