摘 要研究了保險公司的最優(yōu)投資和再保險問題.保險公司的盈余通過跳擴(kuò)散風(fēng)險模型來模擬,可以把盈余的一部分投資到金融市場,金融市場由一個無風(fēng)險資產(chǎn)和n個風(fēng)險資產(chǎn)組成.并且保險公司還可以購買比例再保險;在買賣風(fēng)險資產(chǎn)時,考慮了交易費用.通過隨機(jī)控制的理論,獲得了最優(yōu)策略和值函數(shù)的顯示解.
關(guān)鍵詞 隨機(jī)控制; Hamilton-Jacobi-Bellman(HJB);投資;再保險
中圖分類號 F830,O211.3文獻(xiàn)標(biāo)識碼 A
Optimal Investment and Proportional Reinsurance with Transaction Costs
YANGPeng1, LINXiang2
(1.Department of Basic ,Xijing College, Xian,Shannxi 710123,China;
2. School of Mathematics, Central South University, Changsha,Hunan 410075,China)
Abstract This paper studied an optimal investment and reinsurance problem for insurance company, whose surplus was modeled by a jump diffusion risk process.The insurance company can invest part of the surplus in a risk-free asset, and n risky asset and purchase proportional reinsurance for claims. When purchasing risky asset,we assume there exist transaction costs.Thestochastic control theory was appliedto solve this problem,andthe close form expression for vulue function and optimal policy was obtained.
Keywords stochastic control; Hamilton-Jacobi-Bellman(HJB) equation; investment;reinsurance
1 引 言
Browne[1]首先應(yīng)用隨機(jī)控制理論研究了擴(kuò)散風(fēng)險模型的最優(yōu)投資問題.之后,許多學(xué)者在這一領(lǐng)域進(jìn)行了研究.Bai 和Guo[2]對擴(kuò)散風(fēng)險模型,獲得了使終值財富的指數(shù)效用最大的投資和再保險策略.林祥和楊鵬[3]對擴(kuò)散風(fēng)險模型, 考慮投資和再保險對邊界分紅的影響,得到了使得期望邊界紅利最大的最優(yōu)投資和再保險策略,以及最大期望貼現(xiàn)紅利的顯示表達(dá),并通過數(shù)值計算得到投資和再保險對紅利的影響.對跳-擴(kuò)散風(fēng)險模型,Irgen 和Paulsen[4],獲得了使終值財富的指數(shù)效用最大的投資和再保險策略.Yang和Zhang[5],對跳擴(kuò)散風(fēng)險模型獲得了使終值財富的指數(shù)效用最大的投資策略,但沒有考慮再保險.林祥和錢藝平[6]研究了最小化破產(chǎn)概率的最優(yōu)策略,得到了破產(chǎn)概率和最優(yōu)策略的顯示解.
再保險的保費計算通常使用期望值原理.Hald 和 Schimidli[7]研究了最大化調(diào)節(jié)系數(shù)的最優(yōu)投資問題,對再保險的保費他們考慮了方差原理.目前,還沒發(fā)現(xiàn)有文獻(xiàn)應(yīng)用方差原理來研究最優(yōu)投資和再保險問題.
買賣股票是需要傭金、印花稅、過戶費的,即交易股票是有交易費用的,尤其是頻繁的交易時,交易費用是很大的.因此在風(fēng)險投資時,應(yīng)考慮交易費用,考慮交易費用才更符合實際情況.Zhang 等[8]僅對擴(kuò)散風(fēng)險模型,考慮了含交易費用的投資投資和再保險問題.本文對Zhang 等[8]的研究模型進(jìn)行推廣,研究了跳-擴(kuò)散風(fēng)險模型的帶交易費用的最優(yōu)投資和再保險問題.再保險的保費應(yīng)用方差原理.得到了最優(yōu)策略和值函數(shù)的顯示解.
2 模型和Hamilton-Jacobi-Bellman方程
2.1 模型
考慮跳-擴(kuò)散風(fēng)險模型:
dXt=cdt+βdW0t-d∑Nti=1Yi,(1)
X0=x,
其中,
x≥0表示保險公司的初始盈余;
c>0是保險公司單位時間的保費收入;{Yk,k=1,2,…}是一列獨立同分布的(嚴(yán)格)正值隨機(jī)變量,其共同分布為Fy,密度函數(shù)為fy,F(0)=0,Yk表示第k次賠付的大小;{N(t),t≥0}是參數(shù)為λ>0的泊松過程,表示到時刻t為止的總的索賠發(fā)生次數(shù);{W0t,t≥0}是標(biāo)準(zhǔn)的布朗運動,β≥0是常數(shù),表示擴(kuò)散變差參數(shù).此外,假設(shè){Yk,k=1,2,…},{N(t),t≥0}和{W0t,t≥0}之間是相互獨立的.{Xt,t≥0}為保險公司在t時刻的盈余.
下面考慮比例再保險,比例再保險的水平為1-a,即保險公司的自留額為0≤a≤1,分出為1-a,即在每次理賠時保險公司支付100a%,同時再保險公司支付剩余的100(1-a)%.假設(shè)保險公司向再保險公司支付的再保費率為(1-a)λμ1+α(1-a)2λμ2,且c<λ(μ1+αμ2),α為一常數(shù),μ1=EY,μ2=EY2.則考慮比例再保險后,保險公司在t時刻的盈余為
dXat=[c-1-aλμ1-α1-a2λμ2]dt+
βdW0t-d∑Nti=1aYi. (2)
考慮一個金融市場,由n+1個金融資產(chǎn)組成,其中一個是無風(fēng)險資產(chǎn)(債券),時刻t的價格{Bt,t≥0}滿足方程
dBt=r0Btdt,
其中,r0>0為無風(fēng)險利率.n個風(fēng)險資產(chǎn)(股票),在時刻t時的價格Sit滿足下面的隨機(jī)微分方程:
dSit=Sit[ridt+∑nj=1σijdWj(t)],
i=1,2,…,n,
其中,ri≥r0,σij>0為常數(shù),W(t)={W1(t),…,Wn(t)}是n維標(biāo)準(zhǔn)布朗運動,假設(shè)Wjt,j=0,1,2,...,n相互獨立.
買賣風(fēng)險資產(chǎn)都需要交易費用.
θb=[θb1,θb2,…,θbn]′,
θs=[θs1,θs2,…,θsn]′
分別為買賣風(fēng)險資產(chǎn)的交易費用,即買一個單位的風(fēng)險資產(chǎn)i將花費1+θbiSi(t)的資金,賣一個單位的風(fēng)險資產(chǎn)i將得到1+θsiSi(t)的現(xiàn)金.
設(shè)πb,πs分別為買和賣風(fēng)險資產(chǎn)的資金,這里πb=(π1b,π2b,…,πnb), πs=(π1s, π2s,…,πns).因為不能同時買賣風(fēng)險資產(chǎn),所以有πb#8226;π′s=0.
比例再保險水平at和在風(fēng)險資產(chǎn)上的投資πb,πs作為控制變量.在任意時刻t≥0,πb=πb(t),πs=πs(t)和比例再保險水平a=a(t)由保險公司選擇,記π#8226;=(a(#8226;),πb(#8226;),πs(#8226;)).一旦π#8226;被選擇了,則保險公司的財富過程為
dXπt=[c-1-aλμ1-α1-a2λμ2+
r0Xπt+πb(t)(r-r0I-θb)+πs(t)(r+
r0I-θs)]dt+[πb(t)+πs(t)]D′D[πb(t)+
πs(t)]′dW1t+β dW0t-d∑Nti=1aYi.(3)
Xπ0=x,
其中,I是n維單位列向量,
r=(r1,r2,…,rn)',D=σ11 σ21… σn1
σ12 σ22… σn2
σ1n σ2n… σnn.
定義1一個策略π#8226;=(a(#8226;),πb(#8226;),πs(#8226;))稱為可行的,如果π#8226;關(guān)于流Ft是可料的,且對于每個t≥0過程π#8226;滿足下面的條件
1) ∫T0[πibt]2dt<
對所有T<
,
i=1,2,…,n.
2) ∫T0[πist]2dt<
對所有T<
,
i=1,2,…,n
3) 0≤a≤1 .
所有可行的策略記為∏.
2.2 Hamilton-Jacobi-Bellman (HJB) 方程
假設(shè)保險公司的目的是使得時刻T時財富的期望效用最大.設(shè)效用函數(shù)為ux=m-δγe-γx,其中δ>0,γ>0.顯然有u′>0,u″<0.記Vπt,x為時刻t盈余為x時,策略為π時的期望效用,即
Vπt,x=EuXπTXπt=x, 0 目標(biāo)是尋找最優(yōu)的值函數(shù) Vt,x=sup π∈∏Vπt,x(4) 和最優(yōu)的策略π*使得 Vt,x=Vπ*t,x. (5) 采用Flemming和Soner[9]或Schmidli[10]中介紹的標(biāo)準(zhǔn)方法,得到最優(yōu)的期望財富指數(shù)效用Vt,x滿足下面的HJB方程. 定理1 假設(shè)由(4)定義的V關(guān)于t是連續(xù)可微,關(guān)于x是二次連續(xù)可微函數(shù),則V滿足下面的HJB方程: sup π∈∏{Vt+[πb(t)B1+πs(t)B2]Vx+ [c-1-aλμ1-α1-a2λμ2+r0x]Vx+ 12[πb(t)+πs(t)]D′D[πb(t)+πs(t)]′Vxx + 12β2Vxx+λE[V(t,x-aY)-V(t,x)]}=0. (6) 邊界條件 VT,x=ux, (7) 這里Vt,Vx,Vxx分別為V關(guān)于t的一階導(dǎo)數(shù),關(guān)于x的一階導(dǎo)數(shù)和關(guān)于x的二階導(dǎo)數(shù),且 B1=(r1-r0-θb1,r2-r0-θb2,…,rn-r0-θbn)′, B2=(r1+r0-θs1,r2+r0-θs2,...,rn+r0-θsn)′. 由Flemming和Soner[9]或Schmidli[10]有下面的檢驗定理. 定理2 設(shè)W∈C2是一凹函數(shù)是HJB方程(6)的解,滿足邊界條件(7),則式(4)給出的期望財富指數(shù)效用V恰好等于W.進(jìn)一步,若π*使得 Wt+[π*b(t)B1+π*s(t)B2]Wx+12[π*b(t)+ π*s(t)]D′D[π*b(t)+π*s(t)]′Wxx+ [c-1-a*λμ1-α1-a*2λμ2+ r0x]Wx+12β2Wxx+λE[W(t,x-a*Y)- W(t,x)]}=0,(8) 則π*#8226;是最優(yōu)的策略,也就是Wt,x=Vt,x=Vπ*t,x. 3 輔助結(jié)果 引理1 設(shè)s1定義為 s1z=12‖D′-1z+D′-1B1‖2, (9) 其中,z∈[0, )n.則s1有唯一的最小值1∈[0, )n,i.e. ‖D′-11+D′-1B1‖2≤ ‖D′-1z+D′-1B1‖2 ,z∈[0, )n.(10) Kuhn-Tucker條件對于式(9)中最小的1在區(qū)間[0, )導(dǎo)致Lagrange multiplier向量1∈[0, )n使得 1= s11=D′D-11+D′D-1B1, 11=0. 引理2 設(shè)s2定義為 s2z=12‖D′-1z+D′-1B2‖2,(11) 其中,z∈[0, )n,則s2有唯一的最小值: 2∈[0, )n, i.e. ‖D′-12+D′-1B2‖2≤ ‖D′-1z+D′-1B2‖2, z∈[0, )n.(12) Kuhn-Tucker條件對于式(11)中最小的s2在區(qū)間[0, )導(dǎo)致Lagrange multiplier向量2∈[0, )n使得 2= s22=D′D-12+D′D-1B2, 22=0. 引理3 設(shè) l1 定義 l1z=12z′D′Dz-ρ1B′1z,(13) 其中,z∈[0, )n,這里ρ1≥0.則l1有唯一的最小值:ρ1D-11∈[0, )n,這里 1=D′-11+D′-1B1, (14) 其中,1為引理1給出的s1z的唯一最小值.進(jìn)一步′1D-11=0且 l1ρ11=l1ρ1D-1=-12ρ21‖1‖.(15) 引理4設(shè)l2定義為 l2z=12z′D′Dz-ρ2B′2z,(16) 其中,z∈[0, )n, 這里ρ2≥0.則l2有唯一的最小值ρ2D-12∈[0, )n,其中 2=D′-12+D′-1B2,(17) 其中,z2為引理2給出的s2z的唯一最小值.進(jìn)一步′2D-12=0且 l2ρ22=l2ρ2D-1=-12ρ22‖2‖. (18) 3 最優(yōu)投資和比例再保險策略 假設(shè)有一解Wx,b,滿足Wx>0和Wxx<0,則由引理1—引理4的輔助結(jié)果,有 sup πb≥0πbB1Wx+12πbD′Dπ′bWxx =Wxxinf πb≥012πbD′Dπ′b+WxWxxπbB1 =-12W2xWxx‖1‖2, 及 sup πs≥0πsB2Wx+12πsD′Dπ′sWxx =Wxxinf πs≥012πsD′Dπ′s+WxWxxπsB2 =-12W2xWxx‖2‖2, 其中,1和2由式(14)和式(17)給出,所以有 π*b=-D-11WxWxx, π*s=-D-12WxWxx.(19) 把式(19)代入式(6),得到 sup a{[c-1-aλμ1-α1-a2λμ2+r0x]Wx+ 12β2Wxx+λE[W(t,x-aY)-W(t,x))]}+ Wt-12W2xWxx‖1‖2+‖2‖2=0.(20) 由Browne[1]或Yang和Zhang[5],假設(shè)有解: Wt,x=m-δγexp {-xγer0(T-t)+ 12‖1‖2+‖2‖2(T-t)+h(T-t)}.(21) 這里h(#8226;)是一個適合的的函數(shù)使得式(21)是式(20)的一個解,且h(0)=0.因此有 Wt=[W(t,x)-m][xr0γer0(T-t)+12(‖1‖2+ ‖2‖2)-h(huán)′(T-t)], Wx=[W(t,x)-m][-γer0(T-t)], Wxx=[W(t,x)-m][γ2e2r0(T-t)], λE[W(t,x-aY)-W(t,x))] =λ[W(t,x)-m]E[exp {γaYer0(T-t)}-1]. 代入式(20),有 sup a{-h(huán)′(T-t)+[c-1-aλμ1-α1-a2λμ2]#8226; [-γer0(T-t)]+12β2γ2e2r0(T-t)+ λE[exp {γaYer0(T-t)}-1]}=0.(22) 假設(shè) G(a)=-h(huán)′(T-t)+[c-1-aλμ1- α1-a2λμ2]#8226; [-γer0(T-t)]+12β2γ2e2r0(T-t)+ λE[exp {γaYer0(T-t)}-1], 則令Gaa=0,有 μ1+2αμ2(1-a)=E[Yexp {γaYer0(T-t)}].(23) 定理3方程(23)有唯一正根,且0<<1. 證明 設(shè) h(a)=μ1+2αμ2(1-a)-E[Yexp {γaYer0(T-t)}] =μ1+2αμ2(1-a)-∫ 0yeγayer0(T-t)F(dy), 則 h′a=-2αμ2-∫ 0y2γer0T-teγayer0(T-t)F(dy)<0, h″a=-∫ 0y3γ2e2r0T-teγayer0(T-t)F(dy)<0. 因此h(a)關(guān)于a是一單調(diào)遞減、凹函數(shù).又因為 h0=2αμ2>0, h1=μ1-∫ 0yeγyer0(T-t)F(dy)<0. 所以式(23)有唯一正解0<<1. 把代入式(22),得到 h′(T-t)=[c-1-λμ1-α1-2λμ2]#8226; [-γer0(T-t)]+12β2γ2e2r0(T-t)+ λE[exp {γYer0(T-t)}-1].(24) 所以 Wt,x=m-δγexp {-xγer0(T-t)+ 12(‖1‖2+‖2‖2)(T-t)+h(T-t)}. 因此,從定理2有下面的定理. 定理4 對盈余過程(3),最優(yōu)的比例再保險策略a*為 μ1+2αμ2(1-a)=E[Yexp {γaYer0(T-t)}]