带奇异摄动马氏链的倒向随机微分方程及其应用

【摘要】本篇论文主要讨论带奇异摄动马氏链的倒向随机微分方程(BSDEs)及相应偏微分方程(PDEs)的渐进性质和在随机控制及金融数学中的应用.论文包括以下三个部分：第一部分研究带奇异摄动马氏链的倒向随机微分方程在Meyer-Zheng拓扑下的弱收敛问题；第二部分研究随机体制转换系统的最优转换问题,通过带马氏链的斜反射倒向随机微分方程得到转换问题的最优解,并在奇异摄动马氏链情形下,得到相应变分不等式的渐进性质；第三部分给出带马氏链的倒向随机微分方程在随机控制中的一个应用,即正倒向体制转换系统的随机最大值原理.本文中奇异摄动马氏链指转移速率为多时间尺度的马氏链.在许多物理模型中,系统中的组成元素以不同的速率发生变化,为了定量的研究这种变化,并且降低模型的复杂度,我们运用奇异摄动的方法,即假设马氏链的转移速率为两个时间尺度,进而把与快速变换相关的变量平均化,得到与平稳分布相关的极限问题.具体来说,我们假设马氏链at=αtε的生成矩阵为Qε(t)＝Q(t)/ε＋Q(t).这里,ε>0是一个时间参数,并且Q(t)和Q(t)都是一个马氏链的转移矩阵,其中Q(t)表示快速转换部分,Q(t)表示慢速部分.我们将研究在ε→0时,相应方程及优化控制问题的渐进性质.通过得到简化的极限问题,可以降低问题的计算复杂度.下面将进一步介绍论文的内容及结构.第一章,介绍本文中所研究问题的背景及基础知识.第二章,我们研究了奇异摄动马氏链的倒向随机微分方程在Meyer-Zheng拓扑下的弱收敛问题.我们首先通过估计得到带参数ε的一族倒向随机微分方程的紧性,然后通过鞅问题刻画了极限过程.通过倒向随机微分方程给出的概率表示,我们得到相应偏微分方程在ε→0时的收敛性,并给出了数值计算例子.本章主要来自于：R.Tao,Z.WuandQ.Zhang,BSDEswithregimeswitching:Weakconvergenceandapplications,JournalofMathematicalAnalysisandApplications,407(1),97-111,2013.第三章,我们研究了体制转换系统的最优转换问题,其中体制转换系统由带马氏链αt的随机微分方程描述,决策者可以在有限的控制集I={1,2，…,N}中转换状态以最大化效用函数.此问题存在两个“状态转换”,at的状态转换由市场确定,另一个由决策者根据市场体制及市场价格决定.我们通过带马氏链的斜反射倒向随机微分方程刻画了值函数,并构造出最优的转换策略.在马氏链的转移速率为两尺度的情形下,我们通过倒向随机微分方程方法证明了相应变分不等式的收敛性,从而证明了值函数的收敛性.我们也给出了一个数值计算的例子.本章主要来自于：R.Tao,Z.WuandQ.Zhang,Optimalswitchingunderaregimeswitchingmodelwithtwo-time-scaleMarkovchains,已投稿.第四章,我们研究了正倒向随机控制系统的最大值原理.假设系统由带马氏链的耦合的正倒向随机控制方程描述,控制域为凸集,通过凸变分方法,我们给出了最优控制存在的必要和充分条件,并给出一个在投资消费问题中的应用.本章主要来自于：R.TaoandZ.Wu,Maximumprincipleforoptimalcontrolproblemsofforward-backwardregime-switchingsystemandapplications,SystemandControlLetters,61(9),911-917,2012.下面我们给出本论文的主要结论.1.带奇异摄动马氏链的倒向随机微分方程弱收敛问题我们的目的是研究如下倒向随机微分方程的弱收敛：其中,马氏链αε(t)的生成矩阵为马氏链的状态空间为M=M1∪…∪Ml,其中Mk={Skl,…,skmk},这里K∈{1,…,L},M=m1+…+ml.另外,Q(t)有对角线结构其中,(?)k∈{1,…l},Qk(t)是一个状态空间是Mk的马氏链的生成矩阵.Xεt是带奇异摄动马氏链αεt的随机微分方程的解马氏链αt或者它的生成矩阵Q(t)称为弱不可约,如果有唯一的非负解v(t)=(v1(t),…,vm0(t)).此解称为马氏链αt或者它的生成矩阵Q(t)的拟平稳分布.这一部分的主要结果是如下定理：定理0.1.Yεt是BSDEs(0.0.1)的解.在假设条件似2.3-A2.7)下,随机过程Yεt弱收敛于过程Yt,其中Yt是BSDE的解.这里,Bt是一个布朗运动,Vt(j)是马氏链αt的补偿鞅测度.αt的生成矩阵为Q(t)=diag(v1(t),…,v1(t))Q(t)diag(1m1,…,1ml).这里,1mk.=(1,…,1)’∈Rmk是一个mk维列向量,vk(t)=(v1k(t),…,vmkk(t))∈R1×mk是QK(t)的拟平稳分布..f的定义为F(t,i,x,y)=∑jmi=1v3i(t)f(t,sij,x,y).为证明此定理,首先根据Meyer-Zheng判定准则得到Yεt的紧性.为了刻画极限过程,考虑如下算子其中这里Q(t)ij=(λij(t)),b,a的定义为通过算子g对应鞅问题解的唯一性,我们可以刻画Yεt的极限过程.然后,我们分别在粘性解和经典解意义下,通过BSDE(0.0.1)给出对应耦合PDE方程组的概率表示,进而得到PDE方程组的收敛.这部分的主要结果是：定理0.2.uε是如下反应扩散方程的唯一粘性解则Vt∈[0,T]和x∈Rn,当ε→0,uε(t,x)收敛于u(t,x),其中u是下列方程唯一的粘性解注0.1.uε(t,x)收敛于u(t,x)表示对于任给的(t,x)∈[0,T]×Rn和i∈Mk,uε(t,i,x)→u(t,k,x).定理0.3.在假设条件(A2.3-A2.4)和(A2.6-A2.11)下,方程(0.0.3)有唯一的C1,2b解uε(>)t[0,T]和x∈Rn,当ε→0uε(t,x)收敛于u(t,x),这里u是极限方程(0.0.4)的唯一C1,2b解.2.带奇异摄动马氏链的随机系统的最优转换问题及其渐进性质在这一部分,我们主要研究如下体制转换系统的最优转换问题,其中α(s)是连续时间有限状态时齐的马氏链.决策者可以在有限的转换控制集N={1,…,N}中选择,即一个转换控制过程是一列(τn,ξn)n≥1,其中τn是一列停时,表示转换的时间,ξn是一列取值于N的随机变量,表示转换之后的状态.给定初始时间t和初始状态i,一个转换控制过程可以表示为这里1是一个示性函数.我们的目标是选择一个转换控制过程Iti*来最大化效用函数其中gij是从状态i变换到j的固定成本.Vi,p(t,x):=J(i,t,p,x,It,i,*)称为最优转换问题的值函数.为解决此问题,我们运用BSDE方法.考虑如下带马氏链的斜反射倒向随机微分方程：通过惩罚函数方法,BSDE(0.0.6)解的存在性结果由如下定理得到：定理0.4.假设(A3.1-A3.3)成立,那么,BSDE(0.0.6)存在一个解(Yt,p,x,Zt,p,x,Wt,p,x,Kt,p,x)∈S2×M2×H2×N2.下面,我们由验证定理来说明BSDE(0.0.6)解的唯一性.对任意转换控制过程我们定义如下增过程4I：引入下列带转换控制的BSDE:上述BSDE的唯一解记为(YiI,ZiI,WiI).BSDE(0.0.6)解的唯一性及最优策略的选择可由下列验证定理得到：定理0.5.假设(A3.1-A3.3)成立.设(Yt,p,xZt,p,x,Wt,p,x,Kt,p,x)是BSDE(0.0.6)在S2×M2×H2×N2中的一个解.那么,(1)任给I∈Ait,我们有Yit,p,x(s)≥YiI(s),s∈[t,T].(2)令τ0*=t,ξ0=i.定义下列序列{τj*,ξj*}:其中ξ*j是一个随机变量使得那么,是最优转换问题的最优策略,并且另外,我们有Yit,p,x(s)=YiI*(s),s∈[t,T],从而说明了BSDE(0.0.6)解的唯一性.进一步,通过以下定理,我们给出最优转换问题的值函数对相应变分不等式的概率表示：定理0.6.假设(A3.1-A3.3)成立.最优转换问题的值函数V(t，x)是如下变分不等式系统的唯一粘性解终端条件为Vi,p(T,x)=Φ(x),下一部分,我们假设马氏链αε的结构为双时间尺度,即生成矩阵为αε的状态空间为M=M1∪…∪ML,其中Mk={sk1…,skmk},M=m1+…mL.生成矩阵有与第一部分相同的对角线结构.我们将研究当ε→0时最优转换问题对应变分不等式的收敛.首先,我们给出关于增过程Kit,p,x的一个估计：引理0.1.Kit,p,x对Lebesgue测度绝对连续,并且考虑如下带参数ε的变分不等式系统：下面,对于k=1,…,L,令Vi.skl(t,x)=Vi,k(t,x).定义一个平均系数的极限最优转换问题.vk=(v1k,…,vmkk)记为Qk的平稳分布.令记α为一个马氏链,生成矩阵为Q=diag(v1,…,vL)Qdiag(1m1,…,1mL),其中vk为Qk的平稳分布,1n=(1,…,1)’∈Rn.令Q=[λpq](p,q∈{1,…,L}).考虑系数为b,σσ’,f,马氏链为α的最优转换问题.对应变分不等式系统为终端条件为Vi,k(T,x)=Φ(x)其中本部分的主要结果是：定理0.7.对于k=1,…,L和l=1,…,mk,我们有这里,Vik(t,x)是极限变分不等式系统(0.0.9)的唯一粘性解.3.正倒向体制转换系统的随机最大值原理这一部分,我们主要研究如下正倒向随机控制系统的最优控制问题：其中,马氏链的状态空间为M={1,…,k}.Wt=(Wt(1),…,Wt(k)),nt=(nt(l),…,nt(k)),其中,nt(j)=lαt-≠j}λ(αt-,j).记U为取值与凸集U的可行控制集并且满足E∫0t|ut|2dt定义如下效用函数这里,l,h,r是确定性函数.我们的目标是在U中寻找最优控制来最大化效用函数.首先考虑最优控制存在的必要条件.设u(·)是最优控制问题的一个最优控制,对应的系统记为(X(·),Y(·),Z(·),W(·)).设υ(·)是另一个控制过程(不一定取值与U)并且满足u(·)+v(·)∈u.因为控制域U为凸,我们有(?)0≤p≤1,uρ(·):=u(·)+pv(·)∈u.引入如下变分方程：首先,我们可以得到关于变分不等式的如下引理：引理0.2.假设(A4.1-A4.3)成立,如下变分不等式成立定义如下汉密尔顿函数H：[0,T]×M×R×R×Rl×dMρ×U×R×R1×d×R:H(t,i,x,y,z,w,u,p,k,q)=(p,b(t,i,x,u))+(k,σ(t,i,x,u))-(g,g(t,i,x,y,z,wn,u))+l(t,i,x,y,z,wn,u),(0.0.14)这里,wn=(w(l)n(1),…,w(k)n(k)),n(j)=1{i≠j}λij.引入对偶方程：运用Ito公式,可以得到本节的主要结果定理0.8(最大值原理).设u(·)是一个最优控制,(X(·),Y(·),Z(·),W(·))是相应系统方程的解.(P(·),K(·),Q(·))是对偶方程的唯一解.那么,Vv∈U,我们有Hu·(v-vt)≤0,a.e.,a.s..在一定的凸性条件下,我们还可以得到最优控制存在的充分条件：定理0.9.设(A4.1-44.3)成立.此外,我们假设h,r,H对于变量(X,Y,Z,u,W)是凹的(Concave),并且YT=Φ(XT)是YT=K(αT)XT这种特殊形式,这里K(i)是一个确定性函数.(P,Q,K,M)是对偶方程控制取u(·)时的唯一解.那么,u(·)是最优控制如果满足(0.0.16).在论文的最后部分,我们给出此最大值原理在一个投资消费问题中的应用.
【作者】陶然；
【导师】吴臻；
【作者基本信息】山东大学，金融数学与金融工程，2014，博士
【关键词】倒向随机微分方程；双时间尺度马氏链；奇异摄动；弱收敛；Meyer-Zheng拓扑；最优转换；斜反射；变分不等式；粘性解；随机最大值原理；凸变分；

【参考文献】
[1]侯海军.中外房地产税收制度对比研究[D].西安科技大学,企业管理,2013,硕士.
[2]刘晓茹.我国公用企业滥用垄断力行为的法律规制[D].山西财经大学,经济法学,2014,硕士.
[3]解建芸.清政府在抵制美货运动期间的对美外交[D].安徽大学,专门史,2014,硕士.
[4]高精莹.合作学习法在小学英语教学中的应用[D].吉林大学,外国语言学及应用语言学,2013,硕士.
[5]李清杰,钟鹰,康应业.504所天线近场测量实验室(二)[J].空间电子技术.2000(03)
[6]杨晶.人工流产术后三种避孕方法的临床效果比较[D].山西医科大学,妇产科学（专业学位）,2013,硕士.
[7]李荣.共情能力、情感风格对情绪调节的影响[D].湖南师范大学,应用心理学,2013,硕士.
[8]陆邦干.百尺竿头更上—层楼[J].石油地球物理勘探,1990,06:649.
[9]李意德,吴仲民,曾庆波,周光益,陈步峰,方精云.尖峰岭热带山地雨林群落生产和二氧化碳同化净增量的初步研究[J].植物生态学报,1998,02:32-39.
[10]杨远舟.快速公交时空优先策略及其建模研究[D].北京交通大学,2011.
[11]赵利明.机电产品供销存系统的设计与实现[D].大连理工大学,软件工程（专业学位）,2012,硕士.
[12]施慧洪.农户经济行为与农户金融体制研究[D].首都经济贸易大学,金融学,2004,硕士.
[13]钟兴龙.《唐六典》撰修研究[D].吉林大学,历史文献学,2004,硕士.
[14]周玉叶.RAGE在SAP大鼠胰腺和肾脏中的表达及其作用的研究[D].广东药学院,内科学,2013,硕士.
[15]芮洲.ZnO纳米线阵列的制备及其光电化学性能研究[D].哈尔滨工业大学,材料科学与工程,2014,硕士.
[16]罗晓莉.论诗人沈浩波的创作演变[D].兰州大学,中国现当代文学,2013,硕士.
[17]孙小凯.基于RGB-D信息的物体定位与识别[D].浙江大学,2014.
[18]阎婕.脂多糖对小鼠心肌TMEM16A表达的影响[D].河北医科大学,内科学（专业学位）,2014,硕士.
[19]窦仁庆.安徽省教育技术学本科专业发展的问题及对策研究[D].安徽大学,高等教育学,2013,硕士.
[20]蔡晓炜.智能在线PH检测记录仪的研究和设计[D].浙江工业大学,2009.
[21]王冬梅.理性选择视角下信息安全违背行为影响因素实证研究[D].江苏科技大学,管理科学与工程,2014,硕士.
[22]董雪.血链球菌细菌素对牙龈卟啉单胞菌生长曲线及细胞内钙离子浓度的影响[D].佳木斯大学,口腔临床医学（专业学位）,2014,硕士.
[23]冯伟宇.应用巴马香猪建立腹腔热灌注模型安全性评估[D].郑州大学,临床医学（专业学位）,2013,硕士.
[24]王骏.武维屏教授治疗肺纤维化临床资料回顾及经验总结[D].北京中医药大学,中医内科学,2014,硕士.
[25]徐硕.中国影视植入式广告的规制研究[D].中南民族大学,传播学,2013,硕士.
[26]陆广欣.彩色马铃薯试管薯诱导及遗传转化研究[D].浙江大学,植物保护,2012,硕士.
[27]王礼信.阶跃输入作用下单输入单输出MRAC的简化设计方案[J].自动化学报,1987,06:473-474.
[28]肖丽丽.新型高性能Li_4Ti_5O_(12)负极材料的合成及电化学性能研究[D].哈尔滨工业大学,2013.
[29]顾玲玲.我国虚拟经济和实体经济的关系研究[D].西北师范大学,政治经济学,2013,硕士.
[30]李顺光.镍基微结构的飞秒激光烧蚀机理与工艺研究[D].吉林大学,机械制造及其自动化,2014,硕士.
[31]王志伟,孟玲琴,刘丹,薄林海.振动筛选机的优化设计[J].农业机械学报.2006(05)
[32]李松.PS/PtOEP复合微球的制备及其氧传感性能的研究[D].哈尔滨工业大学,无机化学,2014,硕士.
[33]梁婷婷.单层球面网壳在外爆荷载下的可靠度分析及风险评估研究[D].哈尔滨工业大学,结构工程,2013,硕士.
[34]姜岳庆.水平井测井及仪器输送方法研究[D].东北石油大学,地质工程（专业学位）,2013,硕士.
[35]卢鸯.新型样品前处理技术在纺织品有害物质检测中的应用研究[D].浙江大学,2014.
[36]孙兆勇.石墨烯与复合物的制备以及铀吸附性能研究[D].哈尔滨工程大学,2013.
[37]黄晓阳.从铝土矿直接制备氧化铝载体的研究[D].华东理工大学,化学工程,2013,硕士.
[38]罗菲菲,刘贵全,安景琪,张婷慧.一种分层聚类方法及其应用研究[J].成都理工大学学报(自然科学版),2005,06:649-652.
[39]杨亚萍,谢颖.开征碳税对建筑业发展的预期影响[J].经济师,2014,11:133-134.
[40]于倩倩.沟眶象生物学特性及人工饲养初探[D].北京林业大学,森林保护学,2013,硕士.
[41]刘永伟.基于空间统计的中国碳排放时空动态研究[D].中国矿业大学,自然地理学,2014,硕士.
[42]刘杨.简青霉降解木质素的机理及其对碱木质素生物吸附的研究[D].湖南大学,环境工程,2014,硕士.
[43]梅静梁.淮河支流微量元素分布及环境生物效应研究[D].安徽理工大学,环境工程,2013,博士.
[44]晏波.近代甘肃东南地区政治地理研究[D].复旦大学,历史地理,2012,博士.
[45]张勇.松辽盆地南部小城子地区登娄库组—泉一段储层特征研究[D].吉林大学,矿物学、岩石学、矿床学,2013,硕士.
[46]鲍晓华.论英国侵权法上的私人侵扰[D].对外经济贸易大学,法律,2004,硕士.
[47]王晓璐.企业完善内控体系的案例研究[D].山东大学,国际商务（专业学位）,2013,硕士.
[48]顾云.“两马三车”——谈导控设备的共享[J].广播与电视技术.2003(02)
[49]罗日萍.碳纳米管/石墨烯/环氧树脂复合材料的制备及其性能研究[D].华东理工大学,2013.
[50]郝育晨.政治排斥理论视角下的中国农民工政治参与问题探析[D].吉林大学,政治学理论,2014,硕士.