摘 要 嚴格定義了Markov相依風險模型.證明了該模型的一個等價定理,使得Markov相依風險模型中的諸過程之間的關系更清晰. 獲得了Markov相依風險模型的概率結構,構造性地證明了該模型的存在定理.
關鍵詞 Markov相依風險模型;等價定理;Markov鏈;概率結構;組裝法
中圖分類號 F224.7 文獻標識碼 A
Equivalence Theorem and Probabilistic Structure for MarkovDependence Risk Model
MO Xiaoyun1,2,OU Hui1,ZHOU Jieming1
(1. College of Mathematics and Computer Sclence,Hunan Normal University, Changsha,Hunan 410081;
2. Hunan University of Finance and Economics, Changsha,Hunan 410205)
Abstract The Markovdependent risk model was defined rigorously. The equivalence theorem of this model was proved,which makes the relationships of the stochastic processes among Markovdependent risk model more clear. Probabilistic structure of the Markovdependent risk model was obtained,and the existence theorem of this model was proved constructively.
Key words Markovdependence risk model; Equivalence theorem ; Markov chain; Probabilistic Structure
1 引 言
許多文獻將經典風險模型作了多種推廣并研究.特別是將Markov鏈引入模型中,如文獻[1-4]中的Markov相依風險模型,[5,6]中的Markov調制風險模型.文獻[1]將經典風險模型推廣,引進了Markov相依風險模型如下:將考慮下面的半Markov相依結構的模型:“設Wi表示第i-1次到第i次索賠的等待時間,W0=X0a.s.則
P(Wn+1≤x,Xn+1≤y,Zn+1=j|Zn=i,
(Wr,Xr,Zr),0≤r≤n)
=PW1≤x,X1≤y,Z1=j|Z0=i
=(1-e-λix)PijBj(y).(1)
其中,{Zn,n≥0}是不可約的離散時間Markov鏈,具有狀態空間{1,2,…M}和轉移概率P=(Pij,1≤i,j≤M).這樣,在每一個索賠時刻,如Markov鏈跳躍到狀態j,索賠額的分布Bj將依賴于新的狀態j,而下一個索賠到來的等待時間有參數為λj的指數分布.注意在給定Zn-1和Zn的條件下,隨機變量Wn+1和Xn+1是獨立的,但在相鄰的索賠額之間和在相鄰的索賠等待時間之間存在自相關性,如同Wn和Xn之間存在交互相關那樣.”文獻[1]還指出,Markov相依風險模型包含復合Poisson模型(M=1),Sparre Andersen模型(帶廣義Erlang(n)索賠額分布)[7,8]和Albrecher、Boxma研究的模型[9]等.但文獻[1]中模型的上述描述存在下面的問題.第一,模型的描述不很明確.第二,模型涉及3個相依的隨機過程,3個相依隨機過程是否存在,文獻[1]沒有給出明確回答.第三,式(1)后面的幾個論斷不知是對模型的要求還是式(1)的推論.
本文對Markov相依風險模型給出了上面所述的3個問題的明確回答,獲得了Markov相依風險模型的一個等價定理,給出了Markov相依風險模型的概率結構,構造性地證明了這種模型的存在定理,為研究Markov相依風險模型的隨機過程性質和破產問題打下了嚴實的基礎.2 Markov相依風險模型
設M是給定的正整數,E=(1,2,…,M),R+=(0,∞).
定義1 設Ρ=(Pij,i,j∈E)是一個隨機矩陣;正數λi>0(i∈E),Bi(i∈E)是R+上的概率分布函數,稱
A=Ρ,λi,Bi,i∈E (2)
為Markov特征組.
定義2 給定Markov特征組A=Ρ,λi,Bi,i∈E和定義在同一概率空間(Ω,F,P)上的3個隨機過程:R+值過程W={Wn,n≥1}和X={Xn,n≥1},E值過程Z={Zn,n≥0},W0=X0=0 a.s.,如果
P(Wn+1≤x,Xn+1≤y,Zn+1=j|Zn=i,
(Wr,Xr,Zr),0≤r≤n)
=(1-e-λix)·PijBj(y),i,j∈E,x,y∈R+(3)
成立,稱4元組U=(A;W,X,Z)為Markov相依風險模型,并稱A為模型U的特征組.
定義3 設U=(A;W,X,Z)是Markov相依風險模型.對常數x≥0,c>0,令
N(t)=sup {n≥1:∑ni=1Wi≤t},t≥0. (4)
R(t)=x+c t-∑N(t)j=1Zj. (5)
稱R={R(t),t≥0}為由模型U導出的Markov相依風險過程.
3 等價定理
證明Markov相依風險模型的一個等價定理.從該定理看出Markov相依風險模型各要素之間清晰、明確的聯系.
定理1 給定Markov特征組A=(Ρ,λi,Bi,i∈E)和定義在同一概率空間(Ω,F,P)上的3個隨機過程:R+值過程W={Wn,n≥1}和X={Xn,n≥1},及E值過程Z={Zn,n≥0},W0=X0=0 a.s..則式(3)成立的充要條件是下面的(ⅰ) ~(ⅴ)同時成立.下面的n≥0,i,j∈E,x,y∈R+.
(ⅰ)P(Zn+1=j|Zn=i,(Wr,Xr,Zr),0≤r≤n)=Pij.
(ⅱ) Z是一步轉移矩陣為Ρ的時間齊次Markov鏈:P(Zn+1=j|Zr=ir,0≤r≤n-1,Zn=i)=Pij.
(ⅲ)當Pij>0時,
P(Wn+1≤x|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n)
=P(Wn+1≤x,|Zn=i,(Wr,Xr,Zr),0≤r≤n)
=(1-e-λix).
(ⅳ)當Pij>0時,
P(Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n)
=P(Xn+1≤y|Zn+1=j,(Wr,Xr,Zr),0≤r≤n)
=Bj(y).
(ⅴ) 當Pij>0時,在條件Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n下,Wn+1與Xn+1獨立.即當Pij>0時,有
P(Wn+1≤x,Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n)
=P(Wn+1≤x|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n)
P(Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n).
證明 必要性.設式(3)成立.
在式(3)中令x→+∞,y→+∞得(ⅰ) .由于σ代數σ(Zr,0≤r≤n-1)σ((Wr,Xr,Zr),0≤r≤n),故由條件期望的性質,從(ⅰ)推出(ⅱ),即Z是時間齊次Markov鏈,一步轉移矩陣為Ρ.
在式(3)中令y→+∞得
(1-e-λix)Pij=P(Wn+1≤x,Zn+1=j|Zn=i,
(Wr,Xr,Zr),0≤r≤n).(6)
而上式右方又等于
P(Zn+1=j|Zn=i,(Wr,Xr,Zr),0≤r≤n)·
P(Wn+1≤x|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n).
利用(ⅰ),上式第1因子等于Pij,從而再利用式(6)得出:當Pij>0時,(ⅲ)中式子的最左方等于最右方,從而中間的等號也成立.得證式(ⅲ).
在式(3)中令x→+∞得
PijBj(y)=P(Xn+1≤y,Zn+1=j|Zn=i,
(Wr,Xr,Zr),0≤r≤n)
=P(Zn+1=j|Zn=i,(Wr,Xr,Zr),0≤r≤n)·
P(Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n).
利用(ⅰ),上式右方第1因子等于Pij,故當Pij>0時,從上式得(ⅳ)的最左方等于最右方,從而中間的等號也成立.得證式(ⅳ).
由于式(3)成立,式(3)的右方等于左方,利用(ⅰ),得
(1-e-λix)PijBj(y)
=Pij·P(Wn+1≤x,Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n).
故當Pij>0時,有
P(Wn+1≤x,Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n)
=(1-e-λix)Bj(y).
利用已經證明的結論(ⅲ) (ⅳ)得(ⅴ).
充分性.設 (ⅰ)~(ⅴ)成立,來證式(3). 為證明式(3)成立,只需考慮Pij>0的情形.
式(3)的左方等于
P(Zn+1=j|Zn=i,(Wr,Xr,Zr),0≤r≤n)·
P(Wn+1≤x,Xn+1≤y|Zn=i,Zn+1=j,
(Wr,Xr,Zr),0≤r≤n).
利用(ⅰ),上式第1因子等于Pij,利用(ⅴ)和 (ⅲ) (ⅳ),上式第2因子等于(1-e-λix)Bj(y),從而上式等于式(3)的右方,故式(3)成立.
注 從(ⅱ) (ⅲ) (ⅳ)分別地可推得
PZn+1=j|Zn=i=Pij.(7)
PWn+1≤x|Zn=i=(1-e-λix). (8)
PXn+1≤y|Zn+1=j=Bj(y). (9) 4 Markov相依風險模型的概率結構
定理2 給定Markov特征組A=Ρ,λi,Bi,i∈E,則存在概率空間(Ω,F,P)及定義在其上的3個過程:R+值過程W={Wn,n≥1}和X={Xn,n≥1},E值過程Z={Zn,n≥0},W0=X0=0 a.s.,使得U=(A;W,X,Z)是Markov相依風險模型,其特征組就是給定的A.
證明 用文獻[10]中的獨立乘積空間構造相依隨機變量的組裝法來構造Markov相依風險模型U=(A;W,X,Z).
第1步 在不同的概率空間上分別構造隨機變量族.
依照給定的Markov特征組A中的隨機矩陣Ρ,由時間齊次Markov鏈的存在定理,存在概率空間(Ω1,F1,P1)及定義在其上的時間齊次Markov鏈={Zn~,n≥0},Zn~=Zn~(ω1),ω1∈Ω1,其一步轉移矩陣是Ρ.
依照給定的Markov特征組A中的正數λi,i∈E,利用獨立乘積空間的技巧[11],存在概率空間(Ω2,F2,P2)及定義在其上的獨立隨機變量族{Wn~(i,ω2):i∈E,n≥1}(ω2∈Ω2),使得Wn~(i,ω2)
(ω2∈Ω2)有以λi為參數的指數分布.
依照給定的Markov特征組A中的概率分布函數Bj,j∈E,利用獨立乘積空間的技巧,存在概率空間(Ω3,F3,P3)及定義在其上的獨立隨機變量族{Xn~(j,ω3):j∈E,n≥1}(ω3∈Ω3),使得Xn~(j,ω3)
(ω3∈Ω3)的分布是Bj .
第2步 將不同概率空間上的隨機變量族放在同一概率空間上.作乘積概率空間
(Ω,F,P):=(Ω1×Ω2×Ω3,
F1×F2×F3,P1×P2×P3).(10)
對ω=(ω1,ω2,ω3)∈Ω,i,j∈E,令
Zn(ω)=Zn~(ω1)(n≥0);Wn~(i,ω)=Wn~(i,ω2)(n≥1);Xn~(j,ω)=Xn~(j,ω3)(n≥1). (11)
于是有
獨立性質D:定義在(Ω,F,P)上的諸量Z={Zn,n≥0},Wn~(i),i∈E,n≥1,Xn~(j),j∈E,n≥1,是相互獨立的.
第3步 將同一概率空間上的隨機變量族組裝成3個隨機過程.
Z={Zn(ω),n≥0}(ω∈Ω)已在第2步中定義.對ω∈Ω,令
Wn(ω)=Wn~(Zn-1(ω),ω)
(n≥1),W0(ω)=0;
Xn(ω)=Xn~(Zn(ω),ω)
(n≥1),X0(ω)=0.(12)
得到(Ω,F,P)上的3個隨機過程:Z={Zn,n≥0},W={Wn,n≥1},X={Xn,n≥1},且W0=X0=0.
a.s.,于是有
相依性質H:對任意固定的n(n≥1),定義在(Ω,F,P)上的Wn依賴Zn-1和Wn~(i),i∈E,從而Wr,0≤r≤n,依賴于Z0,…,Zn-1,Wr~(i),i∈E,
0≤r≤n.定義在(Ω,F,P)上的Xn(n≥1)依賴于Zn和Xn~(j),j∈E,從而Xr,0≤r≤n依賴于Z0,…,Zn,Xr~(j),j∈E,0≤r≤n.
下面證明:U=(A;W,X,Z)就是Markov相依風險模型,其特征組是給定的A.只需證明式(3).
由于式(12),式(3)的左方等于
P(n+1(Zn)≤x,n+1(Zn+1)≤y,Zn+1=j|
Zn=i,(Wr,Xr,Zr),0≤r≤n)
=P(n+1(i)≤x,n+1(j)≤y,Zn+1=j|
Zn=i,(Wr,Xr,Zr),0≤r≤n)
=P(Zn+1=j|Zn=i,(Wr,Xr,Zr),0≤r≤n)·
P(n+1(i)≤x,n+1(j)≤y,|Zn=i,
Zn+1=j,(Wr,Xr,Zr),0≤r≤n).
由獨立性質D、相依性質H以及Z是Markov鏈,上式最右方的第1因子等于
P(Zn+1=j|Zn=i,Zr,0≤r≤n)=Pij.
由相依性質H,第2因子中的條件只依賴于Z0,…,Zn,Zn+1,Wr~(i),Xr~(j),i,j∈E,0≤r≤n.從而由獨立性質D,第2因子中的條件概率成為絕對概率
P(n+1(i)≤x,n+1(j)≤y)
=P(n+1(i)≤x)P(n+1(j)≤y)
=(1-e-λix)Bj(y).
上式第1個等號是由于n+1(i)和n+1(j)獨立.第2個等號是由于n+1(i,ω)與n+1(i,ω2)同分布,分布為以λi為參數的指數分布.n+1(j,ω)與n+1(j,ω2)同分布,分布為Bj.于是式(3)的左方等于右方,即式(3)成立.
這樣,U=(A;W,X,Z)是Markov相依風險模型,其特征組是給定的A.定理證完.參考文獻
[1] H AIBRECHER,O BOXMA. On the discounted penalty function in a Markovdependent risk model[J]. Insurance: Mathematics and Economics 2005,37:650-672.
[2] W G YU, Y J HUANG.Dividend payments and related problems in a markovdependent insurance risk model under absolute ruin [J]. American Journal of Industrial and Business Management, 2011,1(1):1-9.
[3] I ADAN, V KULKAMI.Singleserver queue with Markov dependent interarrival and service times[J]. Queuing Systems 2003,45:113-114.
[4] J JANSSEN, J REINHARD. Probabilites de ruine pour une classe de modeles de risqué semiMarkoviens[J].ASTIN Bulletin 1985,15(2):123-134.
[5] J Q WEI, H L YANG, R M WANG. On the Markovmodulated insurance risk model with tax. Blatter DGVFM,2010,31(1):65-78-28.
[6] S ASMUSSEN.Risk theory in a Markovian environment[J].Scandinavia Actuarial Journal,1989,2:69-100.
[7] H GERBER, E SHIU.The time value of ruin in sparre andersen model[J].North American Actuarial Journal, 2005,9(2):49-69.
[8] H ALBRECHER. Discussion on “The time value of ruin in a Sparre Andersen Model” by H.Gerber and E Shiu[J]. North American Actuarial Journal,2005 9(2):71-74.
[9] H AIBRECHER, O BOXMA. A ruin model with dependence between claim sizes and cliam intervals[J]. Insurance: Mathematics and Economics 2004,35(2):245-254.
[10]莫曉云. 用獨立乘積空間構造相依隨機變量的組裝法[J].湖南師范大學學報:自然科學版,2010,33(2):3-6.
[11]嚴加安. 測度論講義[M].2版,北京:科學出版社,1998.