Wiener 过程可以通过几种等价的方式来刻画,上面给出的定义是其中最常见的一种。另一种常见的定义方式是利用我们在动机部分中提到的白噪声。在这种定义下,Wiener 过程是高斯白噪声的定积分,换句话说,高斯白噪声是 Wiener 过程的导数:
W(t)=∫0tη(s)ds(2.29)
dtdW(t)=η(s)(2.30)
为了理解这个关系为何成立,我们先从 [4] 中定义的随机过程的导数说起:
一个随机过程 X(t),t∈R,如果在二次均值意义下可导,其导数为 X′(t),那么它满足:
hX(t+h)−X(t)→X′(t)E[(hX(t+h)−X(t)−X′(t))2]→0(2.31)
其中 h→0 。
我们可以看到,这个定义与普通微积分中导数的定义非常相似,只是这里我们要求的是期望趋近于 0,并且是一个更弱的平方收敛。我们将在下一节中再次看到这种收敛方式。
从这个定义出发,我们可以计算 W(t) 导数的期望:
E[dtdW(t)]=E[h→0limhW(t+h)−W(t)]=h→0limhE[W(t+h)]−E[W(t)]=h→0limh0−0=0(2.32)
类似地,我们可以对随机过程导数的时间相关性给出一个通用性质的推导:
CX′(t1,t2)=E[k→0limkX(t1+k)−X(t1)⋅h→0limhX(t2+h)−X(t2)]=h→0limh1k→0limk1E[(X(t1+k)−X(t1))(X(t2+h)−X(t2))]=h→0limh1k→0lim(kE[X(t1+k)X(t2+h)]−E[X(t1+k)X(t2)]−kE[X(t1)X(t2+h)]+E[X(t1)X(t2)])=h→0limh1k→0lim(kCX(t1+k,t2+h)−CX(t1,t2+h)−kCX(t1+k,t2)−CX(t1,t2))=h→0limh1(∂t1∂CX(t1,t2+h)−∂t1∂CX(t1,t2))=∂t1∂t2∂2CX(t1,t2)(2.33)
因此我们可以得出结论,随机过程导数的时间相关性是协方差函数的二阶偏导。接下来我们只需要将其应用到 Wiener 过程上即可。
首先,假设 t1<t2,Wiener 过程的时间相关性为(详见这个 StackExchange 答案):
0CW(t1,t2)=E[W(t1)(W(t2)−W(t1))]=E[W(t1)W(t2)]−E[(W(t1))2]=E[W(t1)W(t2)]−t1=E[W(t1)W(t2)]=t1=min(t1,t2)(增量独立)(Var(W(t1))=t1)(2.34)
我们若假设 t2<t1,也会得到相同的结果,因此:
CW(t1,t2)=min(t1,t2)
接下来我们要计算它的二阶偏导数。只要 t1=t2,一阶偏导是很容易的(参考这个 StackExchange 答案):
∂t1∂min(t1,t2)={10if t1<t2if t1>t2=H(t2−t1)(除了 t1=t2)(2.35)
其中 H(x) 是 Heaviside 阶跃函数。而我们知道,该阶跃函数的导数就是 Dirac δ 函数(即使在间断点处未定义也无妨),因此:
CW′(t1,t2)=∂t1∂t2∂2min(t1,t2)=∂t2∂H(t2−t1)=δ(t2−t1)(2.36)
从公式 2.32 和 2.36 可见,这和我们在动机部分的公式 1.4 中对白噪声的定义具有相同的统计特性。由于均值为零,协方差就等于时间相关性,即
Cov(W′(t1),W′(t2))=CW′(t1,t2)。
现在我们只需要再证明其服从正态分布。根据上文定义,Wiener 随机过程的导数为:
dtdW(t)=h→0limhW(t+h)−W(t)(2.37)
由于 Wiener 过程的每个增量都服从正态分布(且相互独立),所以上述导数也是正态分布的,因为两个独立正态变量的差仍然服从正态分布。
因此我们可以得出结论:Wiener 过程的导数是一个均值为零、时间相关性为 Dirac δ 函数的高斯过程,这正是高斯白噪声的标准定义。于是我们便证明了公式 2.29 / 2.30 中所给的关系。