隨機過程的通常概念
隨機過程:
T是
(−∞,∞)的子集,對於每一個
t∈T,
Xt是隨機變量,則稱隨機變量的集合
{Xt∣t∈T}是隨機過程,
T是隨機過程的
指標集
隨機過程的觀測或實現:對每一個
t∈T,都獲得了
Xt的觀測值
xt,稱
{xt∣t∈T}是
{Xt∣t∈T}是隨機過程的一次實現或一次觀測
軌跡或軌道:當
T是
[0,∞)或
(−∞,∞)時,
{Xt∣t∈T}的一次觀測又稱爲是一條軌道或一條軌跡
時間序列:
T=N或
T=Z時,隨機過程又稱爲時間序列html
計數過程的通常概念
計數過程:用來描述一段時間內發生的某類事件的次數的隨機過程,計數過程
{N(t)}定義爲:
N(t)是
(0,t]時間內某事件發生的次數,知足:
(1)對
t≥0,
N(t)是取非負整數值的隨機變量
(2)
t≥s時,
N(t)≥N(s)≥0
(3)
N(t)−N(s)表示
(s,t]時間段內該事件發生的次數
(4)
N(t)的軌跡是單調不減右連續的階梯函數
舉例:計數過程能夠表示某段時間內收到的手機短信數,某段時間內到達某個車站的乘客數目,某段時間內商場經理收到商品質量的投訴數,某段時間內股票市場股票成交的次數web
獨立增量性:對於計數過程
{N(t)},若是任取
0<t1<t2<⋯<tn,隨機變量
N(t1),N(t2)−N(t1),⋯,N(tn)−N(tn−1)是相互獨立的,則稱計數過程
{N(t)}具備獨立增量性,該計數過程稱爲獨立增量過程app
平穩增量性:若是對計數過程
{N(t)},對任何的非負實數
s,隨機變量
N(t2)−N(t1),N(t2+s)−N(t1+s)同分布,則稱該計數過程具備平穩增量性,該計數過程稱爲平穩增量過程svg
泊松過程
泊松過程的定義
泊松過程 稱知足如下條件的計數過程
{N(t)}是強度爲
λ的泊松過程:
(1)
{N(t)}是知足
N(0)=0的獨立增量過程和平穩增量過程
(2)
N(t)知足參數爲
λt的泊松分佈函數
強度的含義:
EN(t)=λt,故單位時間內發生的平均事件數爲
tEN(t)=λ所以,強度的實際含義是單位時間內發生時間的平均次數spa
等價定義:若是計數過程
{N(t)}知足:
(1)
N(0)=0
(2)
{N(t)}是獨立增量過程和平穩增量過程
(3)
{N(t)}知足普通性,即:
{P(N(h)=1)=λh+o(h)P(N(h)≥2)=o(h)則
{N(t)}是強度爲
λ的泊松過程(通俗地講是,在極短的時間內,發生一次時間的機率同時間幾乎成正比,幾乎不可能發生兩次或以上事件)orm
證:
令
f0(h)=P{N(h)=0}=1−P{N(h)=1}−P{N(h)≥2}=1−λh+o(h),則當
Δh>0時
===f0(h+Δh)=P{N(h+Δh)=0}P{N(h)=0,N(h+Δh)−N(h)=0}P{N(h)=0}P{N(Δh)=0}f0(h)f0(Δh)上式的第二個等號,能夠從兩個角度理解
第一個角度,計數過程的取值都是非負整數,且知足單調性,
Δh>0,而
N(h+Δh)=N(h)+N(h+Δh)−N(h),因爲
N(h),N(h+Δh)−N(h)都取非負整數值,二者之和爲零的充要條件是二者均爲0,若是
Δh<0
第二個角度,從計數過程的定義來看,
(0,h+Δh]這段時間內沒有任何事件發生的充要條件固然是
(0,h],(h,h+Δh]時間段內都沒有任何事件發生,這兩段時間內只要有任何一段發生了一次事件,那麼
(0,h+Δh]這段時間內都至少會發生1次事件
因而,當
Δh>0時,就有
f0(h+Δh)−f0(h)=[f0(Δh)−1]f(h)=[−λΔh+o(Δh)]f0(h)兩邊除以
Δh,再令
Δh→0+,便可證得
f(h)的右導數存在,而且
f0+(h)=−λf0(h)若是
Δh<0,那麼
f0(h)=f0(h+Δh)f0(−Δh)所以
f0(h+Δh)−f0(h)=f0(h+Δh)[1−f0(−Δh)]=f0(h+Δh)[−λΔh+o(Δh)]由上式不可貴出
Δ→0−limf0(h+Δh)=f(h),兩邊除以
Δh,再令
Δh→0−,便可證得
f(h)的左導數也存在,而且
f0−(h)=−λf0(h)故
f(h)在
h=0處可導,且知足微分方程
f0′(h)=−λf0(h)方程的通解爲
f0(h)=Ce−λh令
h=0,
f0(h)=1,故
C=1,從而獲得
P{N(h)=0}=e−λh。接下來就能夠求解
P{N(h)=1},一樣地,咱們令
f1(h)=P{N(h)=1},當
Δh>0時
f1(h+Δh)=f1(h)f0(Δh)+f1(Δh)f0(h)所以
=f1(h+Δh)−f1(h)=f1(h)[f0(Δh)−1]+f1(Δh)f0(h)f1(h)[−λΔh+o(Δh)]+[λΔh+o(Δh)]e−λh兩邊除以
Δh,令
Δ→0+,便可證得
f1(h)在
h處存在右導數,而且
f1+(h)=−λf1(h)+λe−λh若是
Δh<0,有
==f1(h+Δh)−f(h)f1(h+Δh)[1−f0(−Δh)]−f1(−Δh)f0(h+Δh)[−λΔh+o(Δh)]f1(h+Δh)−[−λΔh+o(Δh)]f0(h+Δh)同理可得
f1−(h)存在,而且
f1−(h)=−λf1(h)+λe−λh因而
f1(h)在
h處可導,而且知足微分方程
f1′(h)=−λf1(h)+λe−λh同時
f1(0)=0,方程組的通解爲
f1(h)=λhe−λh+Ce−λh由
f1(0)=0,獲得
C=0,
f1(h)=λhe−λh,接下來咱們用數學概括法完成接下來的證實,假設對
k≤n,都有
P{N(h)=k}=k!(λh)ke−λh咱們來求
fn+1(h)=P{N(h)=n+1},咱們這裏設
fk(h)=k!(λh)ke−λh,k=1,⋯,n,一樣地,當
Δh>0時,有
fn+1(h+Δh)=k=0∑n+1fk(h)fn+1−k(Δh)故
fn+1(h+Δh)−fn+1(h)=fn+1(h)(1−f0(Δh))+k=0∑nfk(h)fn+1−k(Δh)兩邊除以
Δh,再令
Δh→0+,可證得
fn+1(h)在
h處右導數存在,而且
fn+1+(h)=λfn+1(h)+λfn(h)模仿上面的證實,能夠證得
fn+1在
h處左導數也存在,而且
fn+1−(h)=λfn+1(h)+λfn(h)故
fn+1′(h)=λfn+1(h)+n!λn+1hne−λh該微分方程的通解爲
fn+1(h)=(n+1)!(λh)n+1e−λh+Ce−λh再由
fn+1(0)=0,得
C=0,故
fn+1(h)=(n+1)!(λh)n+1e−λh由概括法,對任意的
n=1,2,⋯,就有
P{N(h)=n}=n!(λt)ne−λt故
{N(t)}是泊松過程xml
泊松呼叫流、等待時間、年齡、壽命
先給出這幾個概念,再給出他們的分佈。假設對
{N(t)}對應的事件流,第一次事件發生的時間點爲
S1,
S2,
S3,
⋯,咱們將其畫在時間軸上,即

泊松呼叫流:計數過程
{N(t)}對應的各事件發生的時刻
{Sn},
{Sn}是單調遞增的非負時間序列,由泊松呼叫流的一次實現能夠獲得
{N(t)}的一次實現,因爲
N(t)的含義是在
(0,t]內發生的事件的次數,故咱們只須要數
(0,t]時間內有多少個
Si,若是
N(t)=n說明:
Sn≤t,
Sn+1>t,咱們獲得如下的基本關係式
{N(t)<n}={Sn>t}{N(t)=n}={Sn≤t<Sn+1}
等待時間:即相鄰兩次事件發生的時間間隔,即
Xi=Si−Si−1,i=1,2,⋯,這裏定義
S0=0
年齡和壽命:咱們把
S1,S2,⋯想象成更換某種零件的時間點,具體地講,在0時刻第一個零件開始運做,在
S1時刻第一個零件損壞,當即更換上第二個零件,在
S2時刻,第二個零件損壞,當即更換上第三個零件,以此類推。假設
N(t)=n,此時已經更換了
n次零件,即此時運做的是機器的第
n+1個零件,那麼
Sn是這個零件開始運做的時刻,
Sn+1是這個零件損壞的時刻,那麼,這個零件的年齡是
t−Sn,剩餘壽命是
Sn+1−t,將
n換成
N(t),年齡就是
t−SN(t),壽命就是
SN(t)+1−t,前者咱們記爲
A(t),後者咱們記爲
R(t)htm
呼叫流的聯合分佈
如今咱們來求呼叫流
(S1,⋯,Sn)的分佈,首先求單個
Sn的分佈:實際上
{Sn≤t}={N(t)≥n},因而就獲得
Sn的分佈函數
Fn(t)=k=n∑∞k!(λt)ke−λt求導便可獲得其密度函數
pn(t)===k=n∑∞kk!λktk−1e−λt−λk=n∑∞k!(λt)ke−λtλk=n∑∞(k−1)!(λt)k−1e−λt−λk=n∑∞k!(λt)ke−λtΓ(n)λntn−1e−λt可見
Sn服從伽馬分佈,參數爲
n,λ,其均值爲
λn,接下來咱們來求
(S1,S2,⋯,Sn)的聯合分佈:咱們假設其分佈函數爲
F(s1,s2,⋯,sn)=P{S1≤s1,S2≤s2,⋯,Sn≤sn}若是
i<j,但
si>sj,那麼
=P{S1≤s1,⋯,Si>si,⋯,Sj≤sj,⋯,Sn≤sn}P{S1≤s1,⋯,Si−1≤si−1,Si+1≤si+1,⋯,Sn≤sn}−F(x1,x2,⋅,xn)=0兩邊對
x1,⋯,xn依次求偏導,則
p(x1,⋯,xn)=0所以,咱們假設
s1≤s2≤s3⋯≤sn,咱們求分佈函數
F(s1,⋯,sn),以
n=2爲例,
s1<s2,則
F(s1,s2)==P{S1≤s1,S2≤s2}P{S2≤s2}−P{S1>s1,S2≤s2}而事件
{S1>s1,S2≤s2}至關於在
(0,s1]段事件發生了0次,
(s1,s2]段事件至少發生了兩次,有
P{S1>s1,S2≤s2}====P{N(s1)=0,N(s2)−N(s1)≥2}P{N(s1)=0}P{N(s2−s1)≥2}e−λs1(k=2∑∞k!λk(s2−s1)k)e−λ(s2−s1)(k=2∑∞k!λk(s2−s1)k)e−λs2則密度函數
p(s1,s2)=∂s1∂s2∂2F(s1,s2)=λ2e−λs2所以
(S1,S2)的密度函數爲
p(s1,s2)=λ2e−λs2Is1<s2,再求
n=3的情形,設
s1<s2<s3,那麼
F(s1,s2,s3)===P{S1≤s1,S2≤s2,S3≤s3}P{S2≤s2,S3≤s3}−P{S1>s1,S2≤s2,S3≤s3}P{S2≤s2,S3≤s3}−P{S1>s1,S2≤s2}+P{S1>s1,S2≤s2,S3>s3}從下圖能夠看出
{S1>s1,S2≤s2,S3>s3}={N(s1)=0,N(s2)−N(s1)=2,N(s3)−N(s2)=0}
故
==P{S1>s1,S2≤s2,S3>s3}P{N(s1)=0}P{N(s2)−N(s1)=2}P{N(s3)−N(s2)=0}e−λs12λ2(s2−s1)2e−λ(s2−s1)e−λ(s3−s2)=2λ2(s2−s1)2e−λs3一樣地方法能夠獲得
(S1,S2,S3)的密度函數爲
p(s1,s2,s3)=λ3e−λs3Is1<s2<s3
n=2及
n=3情形已經給出求
(S1,⋯,Sn)的通常方法,當
n=2k+1,
s1<s2<⋯<s2k時,首先,比較容易求出的是如下機率
==×=Gk(s1,s2,⋯,s2k)P{S1>s1,S2≤s2,S3>s3,S4≤s4,⋯,S2k−1>s2k−1,S2k≤s2k}i=1∏nP{N(s2i−1)−N(s2i−2)=0}.i=1∏n−1P{N(s2i)−N(s2i−1)=2}P{N(s2k)−N(s2k−1)≥2}i=1∏n−12λ2(s2i−s2i−1)2(j=2∑∞j!λj(s2k−s2k−1)j)e−λs2k其中規定
s0=0,先對
s1,⋯,s2k−2依次求偏導
∂s1∂s2⋯∂s2k−2∂2k−2Gk=λ2k−2(j=2∑∞j!λj(s2k−s2k−1)j)e−λs2k這樣對其他兩個變元求偏導就輕鬆不少,具體過程再也不詳述,最終獲得
∂s1∂s2⋯∂s2k∂2kGk=(−1)kλ2ke−λs2k如今,咱們令
G1(s1,s2,⋯,sn)=G2(s1,s2,⋯,sn)=⋯Gi(s1,s2,⋯,sn)=⋯Gk(s1,s2,⋯,sn)=P{S1>s1,S2≤s2,S3≤s3,⋯,S2k≤s2k}P{S1>s1,S2≤s2,S3>s3,S4≤s4,S5≤s5,⋯,S2k≤s2k}P{S1>s1,S2≤s2,⋯,S2i−1>s2i−1,Sj≤sj,n≥j>2i−1}P{S1>s1,S2≤s2,S3>s3,S4≤s4,⋯,S2k−1>s2k−1,S2k≤s2k}因而
F(s1,s2,⋯,sn)=P{S2≤s2,S3≤s3,⋯,S2k≤s2k}−G1(s1,⋯,sn)G1(s1,⋯,sn)=P{S1>s1,S2≤s2,S4≤s4,⋯,S2k≤s2k}−G2(s1,⋯,sn)⋯這樣看
F和
Gk的偏導數的關係應該是
∂s1∂s2⋯∂s2k∂2kF=(−1)k∂s1∂s2⋯∂s2k∂2kGk從而獲得
(S1,⋯,S2k)的聯合密度應該是
p(s1,⋯,sn)=λ2ke−λs2kI0<s1<⋯<s2k當
n爲奇數的時候也進行相似的討論,就獲得
(S1,⋯,Sn)的聯合密度爲
p(s1,⋯,sn)=λne−λsnI0<s1<s2<⋯<sn
總結:
{Sn}是強度爲
λ的泊松過程
{N(t)}的呼叫流,則
(1)
Sn服從
Γ(n,λ)
(2)
(S1,⋯,Sn)的聯合密度爲
p(s1,⋯,sn)=λne−λsnI0<s1<s2<⋯<snblog
呼叫流的條件分佈
下面給定
N(t)=n的條件,
(S1,⋯,Sn)的分佈:
按照條件機率的定義
P(S1≤s1,⋯,Sn≤sn∣N(t)=n)=P{N(t)=n}P{Si≤si,i=1,2,⋯,n,N(t)=n}同非條件的聯合分佈的討論同樣,若是
(s1,⋯,sn)不知足
0<s1<⋯<sn<t,則對
s1,⋯,sn依次求偏導等於0,所以咱們假設
0<s1<⋯<sn<t。若是
n=2k,比較容易求解的一個機率是
P{S1>s1,S2≤s2,S3>s3,S4≤s4,⋯,S2k−1>s2k−1,S2k≤s2k,N(t)=n}咱們把這些時刻畫在時間軸上

這樣,如何求解這個機率值就一目瞭然了,這個事件等價於在
(0,s1]段事件發生了0次,在
(s1,s2]段事件發生了2次,
(s2,s3]段事件發生了0次,
(s3,s4]段事件發生了2次,
⋯,在
(s2k−1,s2k]段事件發生了兩次,在
(s2k,t]事件發生了0次。所以
===G(s1,⋯,sn)P{S1>s1,S2≤s2,S3>s3,S4≤s4,⋯,S2k−1>s2k−1,S2k≤s2k,N(t)=n}i=1∏k(2λ2(s2i−s2i−1)2e−λ(s2i−s2i−1))i=1∏ke−λ(s2i−1−s2i−2).e−λ(t−s2k)i=1∏k[2λ2(s2i−s2i−1)2]e−λt其對
s1,⋯,s2k依次求偏導,獲得
∂s1⋯∂s2k∂2kG=(−1)kλ2ke−λt相似於非條件分佈的討論,令
F(s1,⋯,sn)=P{S1≤s1,⋯,Sn≤sn,N(t)=n},就有
∂s1⋯∂s2k∂2kF=λ2ke−λt因而,在給定
N(t)=n條件下,條件密度爲
p(s1,⋯,sn)=(2k)!(λt)2ke−λtλ2ke−λt=t2k(2k)!=tnn!I0<s1<⋯<sn<t當
n爲奇數時討論也相似,獲得的條件密度函數也是上式。實際上,這也是來自均勻分佈
U(0,t)的樣本
U1,⋯,Un的順序統計量
U(1),⋯,U(n)的聯合分佈。也就是說,在給定
N(t)=n的條件下,這
n個事件發生的時間點在
(0,t]內均勻分佈,再進行一個排列就獲得
S1,⋯,Sn。
結論: