第八章 Gauss过程与 Brown运动 81多维正态分布 设X是n维随机变量,称X服从n维正态分布,如果它的特征函数为 (1)=exp{t14--tl},并记为X~N(A,∑)。易知,EX=,Var(X)=∑。如 果≠0,则X的分布密度为f(x)= 定理8.11:多元正态分布随机变量的边际分布仍然是正态分布 定理8.12:X~N(μ,∑)为n维正态分布台对任意n维向量t tX~N(tu,t∑t 定理8.3:X~N(,∑)为n维正态分布,Y=AX,其中A。n为常值矩阵,则 Y~N(A,A∑A) 定理814:随机变量512的联合分布是正态分布,则512相互独立等价于51,2 不相关。 例81.1:若51253,4的联合分布为零均值的正态分布,则 E5925354=E5152E5354+E5153E5294+E5154E5253° 证明:设其特征函数为 q(1,12,1214)=Eexp(∑15k)=exp(-l2l) 设4=∑o1) 28120kA1)=exp(
第八章 Gauss 过程与 Brown 运动 8.1 多维正态分布 设 X 是 n 维随机变量,称 X 服从 n 维正态分布,如果它的特征函数为 } 2 1 (t) exp{ jt t t T T ϕ = µ − Σ ,并记为 X ~ N(µ,Σ) 。易知,EX = µ,Var(X ) = Σ。如 果 Σ ≠ 0,则 X 的分布密度为 ( ) ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ − − Σ − Σ = − ( ) ( ) 2 1 exp (2 ) 1 1 2 1 2 µ µ π f x x x T n 。 定理 8.1.1:多元正态分布随机变量的边际分布仍然是正态分布。 定理 8.1.2 : X ~ N(µ,Σ) 为 n 维正态分布 ⇔ 对任意 n 维向量 t , t X ~ N(t ,t t) 。 T T T µ Σ 定理 8.1.3: X ~ N(µ,Σ) 为n 维正态分布,Y = AX ,其中 为常值矩阵,则 。 Am×n ~ ( , ) T Y N Aµ AΣA 定理 8.1.4:随机变量 1 2 ξ ,ξ 的联合分布是正态分布,则 1 2 ξ ,ξ 相互独立等价于 1 2 ξ ,ξ 不相关。 例 8.1.1 : 若 1 2 3 4 ξ ,ξ ,ξ ,ξ 的联合分布为零均值的正态分布,则 Eξ 1 ξ 2ξ 3ξ 4 = Eξ 1 ξ 2Eξ 3ξ 4 + Eξ 1 ξ 3Eξ 2ξ 4 + Eξ 1 ξ 4Eξ 2ξ 3。 证明:设其特征函数为 ) 2 1 ) exp( 2 1 exp( ' ) 2 1 ( , , , ) exp( ) exp( 4 1 4 , 1 , 4 1 1 2 3 4 ∑ ∑ ∑ = = = = − = − = = − Σ k k k k l k k l l k k k t t t u t t t t E j t t t σ ϕ ξ (设 ∑ )。 = = 4 l 1 k kl l u σ t 1
=1,,,1)-(u4+∑1o4) q(1,22t1)-(a1+1)=-49(1t2,l24) a-o o(1,12t2,t4)+u1l2o(1,12t2t4) at,at, aL2 -σ12(t12t2,4)+l429p(1l2,l3n) = u,om p([2412)+24.2ot,,,)-ano,n,,) at, at, at, l3o2o(t123,t4)+(a1l2+o2l1)o(1,121,t4)-412u29(1t2,t3,t) a a,oa,.o%4,)-4口1D12,) +{c22+a31l2l14)-(a12+a)ou,,,) 0.0,,4)+24,, o(t,l2,t3,t4)-oi2ll4g(1,2,l2,t4) o{tn2,t3,t4)-( 4o(t1 (o1l2l43+a2u1l43+ol12)o(1,t2,,t1)+1l112g(12,1,t4) 从而E5525=9 12O34+13O24+23O14 at, at,at,at 4 oran,mn=-ou,1≤k,≤4 因此E1点2534=E12E34+E5153E524+E5253E5154 定理815:y=NCmD其中 >0.则給定Y2时 H1的条件分布是p维正态分布,并且条件期望和方差分别为 E(2)=1+V12V2(2-42),Vam(Y|2)=V1-V2V2V2
( ) ( ) ( ) 1 2 3 4 1 1 1 1 2 3 4 4 1 1 2 3 4 1 1 1 ( )] , , , 2 1 , , , [ ( )] 2 1 , , , [ t t t t u u u t t t t t t t t u t t k k k ϕ ϕ ϕ σ ϕ = − + = − = − + ∂ ∂ ∑= ( ) ( ) ( ) ( ) 12 1 2 3 4 1 2 1 2 3 4 1 2 3 4 1 2 1 2 3 4 2 1 2 1 2 , , , , , , , , , , , , t t t t u u t t t t t t t t u u t t t t t u t t σ ϕ ϕ ϕ ϕ ϕ = − + + ∂ ∂ = − ∂ ∂ ∂ ( ) ( ) ( ) ( ) , , , ( ) ( , , , ) ( , , , ) , , , , , , , , , 3 12 1 2 3 4 13 2 23 1 1 2 3 4 1 2 3 1 2 3 4 1 2 3 4 1 2 3 1 2 3 4 3 1 2 3 12 1 2 3 4 3 2 1 3 u t t t t u u t t t t u u u t t t t t t t t u u u t t t t t u u u t t t t t t t σ ϕ σ σ ϕ ϕ σ ϕ ϕ ϕ ϕ = + + − − ∂ ∂ = + ∂ ∂ ∂ ∂ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = ∂ ∂ kl l k t u σ 。 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 14 2 3 24 1 3 34 1 2 1 2 3 4 1 2 3 4 1 2 3 4 13 24 23 14 1 2 3 4 13 2 23 1 4 1 2 3 4 12 34 1 2 3 4 12 3 4 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 4 4 1 2 3 1 2 3 4 13 2 23 4 1 2 3 4 4 1 23 4 2 12 12 1 2 3 4 3 4 12 1 2 3 4 4 3 4 3 2 1 4 ( ) , , , , , , ( ) , , , ( ) , , , , , , , , , , , , , , , [ ] , , , ( ) , , , , , , , , , 1 u u u u u u t t t t u u u u t t t t t t t t u u u t t t t t t t t u u t t t t t t t t u u u u t t t t t u u u t t t t u u t t t t t u t u t t t t u u t t t t t u t t t t u σ σ σ ϕ ϕ σ σ σ σ ϕ σ σ ϕ σ σ ϕ σ ϕ ϕ ϕ σ σ ϕ σ σ ϕ σ ϕ σ ϕ ϕ − + + + + + − + = − + ∂ ∂ − − + ∂ ∂ + ∂ ∂ + − ∂ ∂ = ∂ ∂ ∂ ∂ ∂ 从而 0 12 34 13 24 23 14 1 2 3 4 4 1 2 3 4 | σ σ σ σ σ σ ϕ ξ ξ ξ ξ = + + ∂ ∂ ∂ ∂ ∂ = t= t t t t E t kl l k k l t t E σ ϕ ξ ξ = − ∂ ∂ ∂ = =0 2 | ,1 ≤ k,l ≤ 4 因此 Eξ 1 ξ 2ξ 3ξ 4 = Eξ 1 ξ 2Eξ 3ξ 4 + Eξ 1 ξ 3Eξ 2ξ 4 + Eξ 2ξ 3Eξ 1 ξ 4。 定理 8.1.5:设 。则给定 时 的条件分布是 ~ ( , ), , 0 21 22 11 12 2 1 2 1 >⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = V V V V N V V Y Y Y p q µ µ µ 其中µ Y2 Y1 p 维正态分布,并且条件期望和方差分别为 : ( ) ( ) 21 1 2 2 1 2 11 12 22 1 1 2 1 12 22 E Y | Y V V (Y ), Var Y | Y V V V V − − = µ + − µ = − 。 2
证明:定义 1222 B 1-22,则“是正态分布, 01人y2 Y2 且 1 2 BArly B V1-V22210 这表明 H1-V122Y2和Y2是独立的,因此给定Y2,H的条件分布是p维正态分布 E(x1|x2)=E(1-V2F2Y2+222|y2)=E(x1|y2)=A1+V2l2(x2 ar(x|2)=ar(x-F1F2y2+F222|2)=a(1-V2F2Y2|y2)=V1-V2F22 82Gau过程 定义82.1:随机过程X(),t∈T称为Gaus过程,如果对于任意的n, ,2…tn∈T,随机向量(X(1)X(2)…X(tn)为n维正态分布,即其特征函数 为 a,…,un)=exp∑u1(4)-∑∑ur(2) 其中少()=EX(),r(s,t)=cov(X(s),X() 引理82:x=(x1",…,x)是k维实正态随机序列,如果x均方收敛于 X,即X(均方收敛于X1≤i≤k,那么X是k维正态随机向量 证明:令1 ()…,),m=EX=(…,山) ()=E( )( ),=(G 叨人k×k E(X-)(X-) 由X均方收敛于X知lim1"=p,limo/=on,1≤i,l≤k,从而 u xpijuu-ou 为 imx=X, limp(u)=p(n),o{)是X的特征函数 故 o()=exp{n-nxn,从而x是k维正态随机向量
证明:定义 ,则 是正态分布, 且 ,这表明 是独立的,因此给定 , 的条件分布是 ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − = = ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − =⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − − 2 2 1 1 12 22 2 1 1 12 22 0 Y Y V V Y By Y Y I I V V v u q p ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ v u ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − = = ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − =⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − − 22 21 1 11 12 22 2 2 1 1 12 22 0 0 , ( ) ' V V V V V BVar y B v u Var V V v u E µ µ µ 2 2 1 Y1 V12V22 Y 和Y − − Y2 Y1 p 维正态分布。 ( ) | ( | ) ( | ) ( ) 2 2 1 2 2 1 2 1 12 22 1 2 12 22 1 1 2 = 1 − 12 22 + = = µ + − µ − − − E Y Y E Y V V Y V V Y Y E Y Y V V Y ( ) ( ) ( ) 21 1 2 2 11 12 22 1 2 2 1 12 22 1 2 12 22 1 1 2 1 12 22 Var Y | Y Var Y V V Y V V Y | Y Var Y V V Y | Y V V V V − − − − = − + = − = − 8.2 Gauss 过程 定义 8.2.1:随机过程 X (t),t ∈T 称为 Gauss 过程,如果对于任意的 , ,随机向量 n t1 ,t2 ,Ltn ∈T ( ( ) ( ) ( ) T n X t , X t ,LX t 1 2 ) 为 维正态分布,即其特征函数 为 n ( ) ( ) l k n l n k l k n l t t t n l l u u j u t u u t t n , 2 1 , , exp ( ) 1 1 1 1 , 2 , 1 = ∑ − ∑∑ Γ = = = ϕ L L µ 其中µ( )t = EX (t) ,Γ( ) s,t = cov(X (s), X (t))。 引理 8.2.1: ( ) 是 维实正态随机序列,如果 T n k n n X X X ( ) ( ) 1 ( ) = ,L, k (n) X 均方收敛于 X ,即 Xi (n) 均方收敛于 X i k i ,1 ≤ ≤ ,那么 X 是k 维正态随机向量。 证明:令 ( ) , , , ( ) , , , 1 ( ) ( ) 1 ( ) ( ) T k T n k n n n µ = EX = µ L µ µ = EX = µ L µ ( ) n n n n T k k n ij n E(X )(X ) ( ) ( ) ( ) ( ) ( ) ( ) Σ = σ = − µ − µ × , ( ) T k k Σ = σ ij = E(X − µ)(X − µ) × , 由 (n) X 均方收敛于 X 知 l n l il n il i l k ,从而 n = = ≤ ≤ →∞ lim ,lim ,1 , ( ) µ µ σ σ ( ) } 2 1 } exp{ 2 1 exp{ ( ) ( ) u ju u u ju u u T n T n T T ϕ n = µ − Σ → µ − Σ ,因为 X X u u (u) n n n n lim ,limϕ ( ) ϕ( ),ϕ ( ) = = →∞ →∞ 是 X 的特征函数,故 ( ) } 2 1 u exp{ ju u u T T ϕ = µ − Σ ,从而 X 是k 维正态随机向量。 3
定理821:设{X(1),【∈T}是均方可微的Guss过程,那么{X(,t∈T}也是 Gauss 过程 推论8.,1:设Y()是 Gauss过程,EX()=p(t),cov(X(s),X()=I(s,1),X(1)均 方可微。那么对于任意的n,l,12…,n∈T,(X'(1)…,X(tn)的特征函数是 ar(s, la n)=eXp{∑n1(4)-∑∑1 定理822:如果{X()∈T}是均方可积Gas过程,那么Y()=「X(s)d也是 auss过程。 证明:只需证明对于任意的n,1,l2…,tn∈T,(Y(1)…,Y(n)是n维正态随机向 由于 Y()=∫X(s)ds=m∑X(s)A△s,其中a=S0<…<Sn=1 =max(S,-S,1), 再由引理82.1和上一节定理8.1.2立得。 推论82.2:{X()t∈}是均方可积 Gauss过程,EX(l)=p(t), r(s)=co(x(s)xX(),YO)=X(s)。那么对于任意的n12…;tn∈T, (Y(1),…Y(tn)的特征函数是 u(s)ds 83 Brown运动
定理8.2.1:设 是均方可微的Gauss过程,那么 也是Gauss 过程。 {X (t),t ∈T} {X '(t),t ∈T} 推论 8.2.1:设 X (t)是 Gauss 过程,EX (t) = µ(t),cov(X (s), X (t)) = Γ(s,t), 均 方可微。那么对于任意的 X (t) n,t1 ,t2 ,L,tn ∈T , 的特征函数是 T n (X '(t ), , X '(t )) 1 L ( )} , 2 1 ( , , , ) exp{ '( ) 2 1 1 1 1 , 2 , , 1 2 s t s t u u u j u t u u n l n k l k n l n l l X t t tn ∂ ∂ ∂ Γ = ∑ − ∑∑ = = = ′ ϕ L L µ 。 定理 8.2.2:如果 是均方可积 Gauss 过程,那么 也是 Gauss 过程。 {X (t),t ∈T} ∫ = t a Y(t) X (s)ds 证明:只需证明对于任意的 是 维正态随机向 量。由于 T n n n,t ,t , ,t T,(Y(t ), ,Y(t )) 1 2 L ∈ 1 L n ∫ ∑= − ∆ → ≤ ≤ = = ∆ = < < = ∆ = − k k k k k kn t a n l i i i n k l l n k n Y t X s ds X s s a s s t s s 1 1 1 0 0 ( ) ( ) lim ( ) ,其中 L , max( ), 再由引理 8.2.1 和上一节定理 8.1.2 立得。 推 论 8.2.2 : {X (t),t ∈T} 是均方可积 Gauss 过 程 , EX (t) = µ(t) , , 。那么对于任意的 , 的特征函数是 Γ( ) s,t = cov(X (s), X (t)) ∫ = t a Y(t) X (s)ds n,t1 ,t2 ,L,tn ∈T T n (Y(t ), ,Y(t )) 1 L ( ) ( ) ⎪⎭ ⎪ ⎬ ⎫ ⎪⎩ ⎪ ⎨ ⎧ = − Γ ∑ ∫ ∑∑ ∫∫ = = = l l k n t a t a n l t a n l n k n l l k Y t t t u u j u s ds u u s,t dsdt 2 1 , , exp ( ) 1 1 1 ϕ 1 , 2 ,L, 1 L µ 。 8.3 Brown 运动 4
先讨论对称随机游动,此游动每个单位时间等可能地向左或向右走一个单位 格子。现假设每隔Mt的时间等概率地向左或向右走Ax大小的位移。若以X()记 t时刻的位置,则 X(O)=△(x1+x2+…Xba) 1若长Δx的第步向右 其中X -1若长Ax的第涉步向左且诸x,相互独立 P(X2=1)==P(X1=-1)。因此EX()=0,ar(X()=(△x)1 现令 △,△x→>0,但要使得极限过程是非平凡的(例如若令Ax=M,再令M→0,则 EX(O)a(X(1)→>0,从而X()→0,as),故令Ax=a√M,再令M→0时 EX()=0,ar(Y()→σ2t。再由中心极限定理可见此极限过程 X()~N(0,a2t)。此外随机游动的值在不相重叠的时间区间中的变化是独立的, 位置的变化分布只依赖于区间长度,因此X(),t≥0应有平稳的独立增量。 定义831:随机过程{W(),t≥0}称为 Brown运动或 Wiener过程,若: 1)W(0)=0 2)W()为平稳的独立增量过程 3)对任意t>0,W()~N(0,a2t)。 如果W()~N(μ·t,σ21),则称W()为有漂移的 Brown运动,μ为漂移参数。对 于BoW运动W(),如果a2=1,称为标准 Brown运动。如果a2≠1f() 标准 Brown运动。 Brown运动是 Gauss过程,均值函数为0,协方差函数 I(s,)=EW(s()=a2min(s,l),而且它还是 Markov过程 定义8.3.2:随机过程{()t≥0}称为 Brown运动,若 1)W(0)=0
先讨论对称随机游动,此游动每个单位时间等可能地向左或向右走一个单位 格子。现假设每隔∆t 的时间等概率地向左或向右走∆x 大小的位移。若以 记 时刻的位置,则 X (t) t [ ] ( ) X X X t t X t x = ∆ 1 + 2 + / ∆ ( ) L , 其 中 且 诸 相互独立, ⎩ ⎨ ⎧ − ∆ ∆ = 1 若长 的第 步向左 1 若长 的第 步向右 x i x i Xi Xi ( 1) 2 1 P(Xi = 1) = = P Xi = − 。因此 EX (t) = 0 , ( ) ( ) ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∆ = ∆ t t Var X t x 2 ( ) 。现令 ∆t,∆x → 0 ,但要使得极限过程是非平凡的(例如若令∆x = ∆t ,再令 ,则 ,从而 ),故令 ∆t → 0 EX (t),Var( ) X (t) → 0 X (t) → 0,a.s. ∆x = σ ∆t ,再令 时 , ∆t → 0 EX (t) = 0 Var(X t ) t 2 ( ) → σ 。再由中心极限定理可见此极限过程 。此外随机游动的值在不相重叠的时间区间中的变化是独立的, 位置的变化分布只依赖于区间长度,因此 应有平稳的独立增量。 ( ) ~ (0, ) 2 X t N σ t X (t),t ≥ 0 定义 8.3.1:随机过程{W (t),t ≥ 0}称为 Brown 运动或 Wiener 过程,若: 1) W (0) = 0; 2) W (t) 为平稳的独立增量过程; 3) 对任意t > 0,W (t) ~ N(0,σ2 t)。 如果W (t) ~ N(µ ⋅t,σ2 t) ,则称W (t) 为有漂移的 Brown 运动,µ 为漂移参数。对 于 Brown 运动W (t) ,如果σ2 = 1,称为标准 Brown 运动。如果 σ σ ( ) 1, 2 W t ≠ 是 标准 Brown 运动。Brown 运动是 Gauss 过程,均值函数为 0,协方差函数 Γ( ) s,t = EW (s)W (t) = σ2 min(s,t),而且它还是 Markov 过程。 定义 8.3.2:随机过程{W (t),t ≥ 0}称为 Brown 运动,若: 1) W (0) = 0; 5