需要进一步研究真解的正则性及其近似。6.2利普希茨BSDEs6。2.1关于Lipschitz-BSDE的收敛估计在提醒中，让我们测试Lipschitz-BSDE的拟议方案的完整性。前几节中给出的模式同样适用于标准Lipschitz-BSDES，并得出相同的收敛估计。在提供数值结果之前，让我们首先解释允许的设置以及相关结果中的相关变化。设置中的主要变化可总结如下：oξ（x）和f（t，x，0，0）在x中最多有多项式增长。oξ∈ C（Rd）和f对于空间变量（x，y，z）是一次性连续可微的，其中|yf ||zf |有界且|xf |在x中最多有多项式增长。o假设3.1（ii）被删除，（i）修改为max1≤我≤n个|mxbui+1（x）|≤ K′（1+| x）|(-m）∨0），十、∈ Rd，m=0，···，3，带某个正整数 ≥ 1、现在允许扰动（δi）具有多项式增长相应修改假设5.1，以检查上述多项式增长条件。注意，用于qg BSDE的终端函数的有界性是保证导数过程(xΘi，t，x）由结果适定[13]。因为衍生过程xΘi，t，xnow遵循Lips chitz BS DE，有界条件是不必要的，命题3.1中的类似分析会产生| Zit |≤ C（1+| Xt）|). 定理3.1同样的分析仍然成立。请注意，|γit |现在以Lipschitz常数为界，因此命题3.2不再需要。

自Xt起，x∈ 任何p的SPP≥ 1、hi的缩放（如Ch3p/2等）附录B和附录C中的短期扩展保持不变。从这些观察结果可以看出，定理4.1中的估计仍然成立。自(|mxbui |）i，mhave多项式增长，第5节中出现的常数C（特别是引理5.1和定义5.1中的常数C）通常依赖于x的多项式。这导致了C→ C（1+| x|). 然而，可以通过在每个地方适当调整第5节中使用的“ρ”来吸收影响。（5.3）也是如此，因为短期扩展的初始值依赖性最多是多项式。命题5.1变得不必要（仅适用于qg BSDE），因此网格大小M的比例因子δ>0不限于δ<1/ρ。命题5.2与修正ρ的证明方式完全相同，通过简单应用切比雪夫不等式，证明了推论5.1的主要结果。它也被用来证明解的唯一性。6.2.2数值示例：不同利率的期权定价使用相同的比例规则 = ζ|π| 1/2但当然，不截断驱动程序。我们考虑两种不同利率下的欧式期权的一个非常流行的估值问题，r代表投资，r（6=r）代表借贷。由于此问题经常用于测试Lipschitz BSDE的数值格式，因此将当前格式与基于蒙特卡罗模拟的现有数值示例进行比较将具有一定的指导意义。让我们假设证券价格的动态asXt=x+ZtuXsds+ZtσXsdWs，其中d=1和u，σ是正常数。对于到期日T时的期权支付Φ（XT），自融资复制的期权价格Y由YT=Φ（XT）给出-ZTtnrYs+u- rσZs-Ys公司-Zsσ-（R）- r） ods-ZTtZsdWs。

（6.10）尽管终端和驱动器功能都不顺畅，但考虑到qg BSDE的图3中的结果，我们可以预期相当准确的结果。图4：认购期权拟议方案（6.10）的经验收敛性。首先，我们研究了Payoff函数等于看涨期权函数的情况：Φ（x）=（x- K） +，其中K>0是履约价格。正如[32]所指出的，这个例子提供了一个非常有趣的测试，因为价格必须与利率为R的Black Scholes模型的价格完全相等。这是因为复制投资组合由多头唯一头寸组成，因此投资者必须始终借钱为其头寸融资。我们选择了{T=1，r=0.01，r=0.06，u=0.06，X=100}等常见参数，并在图4中测试了以下五组（K，σ），其中n=10到n=3000：set={K=106，σ=0.3}，set={K=166，σ=0.3}，set={K=106，σ=1.0}，set={K=306，σ=1.0}，set={K=106，σ=2.0}。虽然我们用molli-fied函数测试了相同的模型，但我们发现在经验收敛速度上没有显著差异。对于T=1，利率为6%时，K=106接近于货币远期利率。较大的罢工对应于货币的2σ。每套Black-Scholes价格分别为BS={12.000、1.117、38.346、11.662、68.296}。尽管OTM期权的相对误差略高，但在每种情况下，精确BS价格的收敛速度都接近1。令人惊讶的是，尽管存在非光滑函数和较大的波动性，但我们并没有看到收敛速度的任何恶化。观测到的误差大小的不规则性可能是由于靠近终端时间的网格的配置相对于终端功能的不可区分点发生了变化。接下来，让我们考虑一个调用扩展情况：Φ（x）=（x-K）+-2（x-K） +。

这与[32]中研究的设置完全相同，因此我们可以测试与基于标准回归的蒙特卡罗模拟相关的ou r sch eme性能。让我们选择[32]中的相同参数集：{r=0.01，r=0.06，u=0.05，X=100，T=0.25，σ=0.2，K=95，K=105}（6.11）。[32]的结果表明，Y=2.96±0.01或Y=2.95±0.01，一个标准偏差取决于r回归的基函数选择。在图e 5中，我们将我们的方案的估计值与[32]中的估计值进行了比较。虚线代表2。96±0.01，便于比较。在我们的s模式中，Yconverge接近2.96。事实上，Bender&Steiner（2012）[5]提出的使用鞅基函数的回归方法[32]的改进建议为2.96，这与我们的结果完全一致。图5：建议的（6.10）通话价差方案的经验收敛性。图6:T=1.0且波动率较高的看涨期权价差的（6.10）预测方案的经验收敛性。我们还测试了最终支付Φ（x）=（x）的更长期限和更高波动性的收敛性- K）+- （十）- K） +。我们像以前一样使用了{r=0.01，r=0.06，u=0.05，X=100，K=95，K=105}，但成熟度更长的是T=1.0和set：={σ=0.3}，set：={σ=0.5}和set：={σ=1.0}。从图6可以看出，所有情况下都有一个ob-servesmooth收敛。波动率较高时价格下降是自然的，因为Kis更接近货币远期点，因此短期头寸对波动率具有更高的敏感性。一个大型Lipschitz-constantBender和Steiner的例子[5]测试了一个极端情况，p参数集（6.11）替换为byR=3.01。

在这种情况下，d河的非线性具有Lipschitz常数（R-r） /σ=15。他们的实验表明，在他们尝试的模拟设置下，[32]的标准方法在这个例子中无法收敛。他们用鞅基函数改进的方法（见[5]中的表3）给出了Y 6.47，n=128，Y 6.44，最小分区n=181。图7:R=3.01时（6.10）建议方案的经验收敛性。在图7中，我们绘制了n=10到n=3000的方案的估计yf。T虚线对应于【5】中给出的值6.44。在我们的方案中，我希望收敛到6.38。特别是，使用相同的离散化n=181，我们的方案产生Y 6.43表现出良好的一致性。注意，方法【5】要求根据X定律更改基函数。备注6.2。数值计算在Microsoft Excel VBA中实现，XeonX5570 cpu@2.93GHz。对于n=100个时间步的qg-bsdew二维示例大约需要8秒，对于n=1000个时间步的Lipschitz-bsdew一维示例大约需要0.8秒。7结论在本文中，我们通过连接由渐近展开技术以显式形式给出的短期表达式，开发了马尔可夫qgBSDEs的半解析计算方案。该方法也可以以几乎相同的方式应用于标准Lipschitz BSDEsin。至少对于低维设置，该方案作为一种略微精细的TREE方法非常容易实现，并且即使对于非常大的qu系数和Lipschitz常数也具有较高的精度。对于高维设置，我们建议使用sparsegrid方案作为有希望的候选方案，以至少在一定程度上克服维度诅咒。

测试这个有趣的想法是留给未来的工作。从理论上看，主要困难在于控制导数的界限|mxbui |统一，这仍然是一个悬而未决的问题。这迫使我们需要对假设5.1中的这些边界进行后验检查，以保证agiven范围的收敛性。在适当的假设下，相应的半线性模型的性质可以提供必要的规律性。这个问题的解决留待将来的工作。BMO鞅及其性质在本节中，让我们总结BMO鞅的性质、相关的HBMO空间及其在讨论中起重要作用的性质。定义A.1。BMO鞅M是满足M=0且| | M | | BMO：=supτ的平方可积鞅∈TTEhhMiT公司- hMiτ| Fτi∞< ∞ ,其中最高值接管所有停车时间τ∈ TT。定义A.2。HBMO（Rk）是Rk值渐进可测过程Z的集合，满足| | Z | | HBM O：=supτ∈TTEhZTτ| Zs | dsFτi∞< ∞ .注意，如果Z∈ HBMO（R1×d），我们有Z·ZsdWsBMO=supτ∈TTEhZTτ| Zs | dsFτi∞= ||Z | | HBM O<∞ ,因此Z*W：=R·ZSDW是BMO鞅。下一个结果就是众所周知的能源不平等。引理A.1。让Z在HBMO中。那么，对于任何n∈ N、 呃ZT | Zs | dsni公司≤ N||Z | | HBM On、 证明。参见[20]中引理9.6.5的证明。设E（M）是M.引理a.2的Dol\'eans-Dade指数。（逆H¨older不等式）设M是BMO鞅。然后Et（M），t∈ [0，T]是一致可积鞅，对于每个停止时间τ∈ TT，存在一些正常量r*> 1使得不等式het（M）r | Fτi≤ Cr，MEτ（M）r，每1<r保持一次≤ R*对于某些正常数Cr，Mdepending仅在r和| | M | | BMO上。证据见Kazamaki的Theo rem 3.1（1994）[34]。引理A.3。

设M是平方可积鞅，cm：=hM i-M那么，M∈ BM O（P）ifand only ifcM∈ BM O（Q），dQ/dP=ET（M）。此外，| | cM | | BMO（Q）由| | M | | BMO（P）的某些函数确定。证据参见【34】中的第2.4条和第3.3条。备注A.1。定理3.1[34]还指出，存在一些递减函数Φ（r），Φ（1+）=∞ 和Φ(∞) = 0这样，如果| | M | | BMO（P）满足| M | | BMO（P）<Φ（r），则E（M）允许具有幂r的反向H¨older不等式。这意味着与引理A.3一起，可以取一个满足1<r的公共正常数≤ R*这样，E（M）和E（^M）在各自的概率测度P和Q下都满足具有幂的反H¨older不等式。此外，上界r*仅由| | M | | BMO（P）（或等效于| | M | | BMO（Q））确定）。B短期扩展：步骤1在以下两部分中，我们近似BSDE（3.1）半解析解（Yi，Zi），并获得其误差估计。我们需要两个步骤来实现这一目标，分别涉及Fujii和Takahashi（2012）[27]和（2015）[31]中提出的BSDE线性化方法和小方差扩展方法。在没有混淆的情况下，我们采用所谓的爱因斯坦约定，假设重复索引（如∈ xi的{1，···，d}），而不显式使用summationsymbolP。例如xi，xjξ（XT）xXiT公司XXJT假设指数i和j的总和表示SPDI，j=1xi，xjξ（XT）xXiT公司xXjT。（B节的持续假设）我们做出假设2.1、2.2和假设3.1（i）本节的持续假设。让我们为每个间隔t介绍BSDE（3.1）的下一个分解∈ 二、一∈ {1，···，n}：Yi，[0]t=bui+1（Xti）-Ztitz，[0]rdWr，（B.1）Yi，[1]t=Ztitfr、 Xr，Yi，[0]r，Zi，[0]r博士-Ztitz，[1]rdWr。

（B.2）它们是线性化方法的主要贡献【27，48】。引理B.1。对于每个间隔Ii，i∈ {1，···，n}，存在一个唯一的解（Yi，[0]，Zi，[0]），以满足BSDE（B.1），具有一些（i，n）独立的正常数C和Cp，即| | Yi，[0]| | S∞[技术信息-1，ti）+| | Zi，[0]| | HBM O[ti-1，ti]≤ C、 还有| | Zi，[0]| | Sp[ti-1，ti]≤ 任何p的CPP≥ 2.证明。有界性| | Yi，[0]| | S∞≤ C很容易遵循假设3.1（i），这意味着| | Zi，[0]| | HBM O≤ C、 第二项主张源自命题3.1中使用的类似论据。引理B.2。对于每个间隔Ii，i∈ {1，···，n}，存在一个唯一的解（Yi，[1]，Zi，[1]），以满足BSDE（B.2），且一些（i，n）独立的正常数Cp，即| | Yi，[1]| | Sp[ti-1，ti]+| | Zi，[1]| | Hp[ti-1，ti]≤ 任何p的CPP≥ 2.证明。由于它是一个Lipschitz BSDE（Lipschitz常数为零），唯一解的存在性注意，自Takahashi（1999）[47]和Kunitomo&Takahashi（2003）[37]首次尝试以来，小方差渐近展开已广泛应用于Europeancontin-gent索赔的定价。易于遵循。标准估计值（例如参见[12]）和假设2.2（i）暗示（一，字，一）pKp【ti】-1，ti]≤ CpEh公司Ztiti公司-1 | f（r，Xr，Yi，[0]r，Zi，[0]r）| dr圆周率≤ CpEh公司Ztiti公司-1.lr+β| Yi，[0]r |+γ| Zi，[0]r|博士圆周率≤ Cp公司||l | | pT+| | Yi，[0]| | pSp[ti-1，ti]+| | Zi，[0]| | 2pS2p[ti-1，ti].因此，我们可以通过引理B.1得到期望的结果。现在，我们定义了每个时间间隔t的过程（eYi，eZi）∈ 二、一∈ {1，···，n}byeYit：=Yi，[0]t+Yi，[1]t，eZit：=Zi，[0]t+Zi，[1]t.命题B.1。存在一些（i，n）独立的正常数Cpsuch，不等式易-eYi公司p【ti】-1，ti]+Ztiti公司-1.齐尔-埃齐尔博士圆周率≤ Cph3p/2适用于每个间隔Ii、i∈ {1，···，n}与任意p≥ 2.证明。

为了符号的简单性，让我们把δYi，[0]t：=Yit- Yi，[0]t，δZi，[0]t：=Zit- Zi，[0]tδYi，[1]t：=Yit-eYit，δZi，[1]t：=Zit-eZitfor每个间隔t∈ 二、一∈ {1，···，n}。然后，它们分别由以下BSD的解给出：δYi，[0]t=zitf（r，Xr，Yir，Zir）dr-ZtitδZi，[0]rdWr，δYi，[1]t=Ztitf（r、Xr、Yir、Zir）- f（r，Xr，Yi，[0]r，Zi，[0]r）博士-ZtitδZi，[1]rdWr。根据Lipschitz BSDE的稳定性结果（例如，参见[12]），假设2.2（i），（3.3）a ndProposition 3.1，可以得到（δYi，[0]，δZi，[0]）pKp【ti】-1，ti]≤ CpEh公司Ztiti公司-1.lr+β| Yir |+γ| Zir|博士圆周率≤ Cphpi公司||l | | pT+| | Yi | pS∞[技术信息-1，ti]+| | Zi | | 2pS2p[ti-1，ti]≤ Cphpi，（B.3）和一些（i，n）-独立正常数CphpiP≥ 使用假设2.2（ii）yie lds对（δYi，[1]，δZi，[1]）进行类似分析（δYi，[1]，δZi，[1]）pKp【ti】-1，ti]≤ CpEh公司Ztiti公司-1h |δYi，[0]r |+（1+| Zir |+| Zi，[0]r |）|δZi，[0]r | idr圆周率≤ Cp公司hpi | |δYi，[0]| | pSp[Ii]+Eh1+| | Zi | 2pIi+| | Zi，[0]| | 2piieh希兹提提-1 |δZi，[0]r | dr圆周率.通过应用属性3.1、emma B.1和之前的估计（B.3），可以得到期望的结果。C短期展开：步骤2在第二步中，我们获得了BSDE（B.1）和（B.2）的简单解析近似，同时保持命题B.1中相同的精度顺序。我们使用[31]中提出的BSDE小方差展开方法，将所有问题转化为一组简单的ODE。此外，我们将看到，这些常微分方程可以用单步Euler方法近似每个中心值Ii。（C节的长期假设）与上一节类似，我们制定了假设2.1、2.2和假设3.1（i）本节的长期假设。C、 （Yi，[0]，Zi，[0]）的1近似对于每个区间，我们引入一个新的参数successing∈ (-c、 c）对于一些常数c>1 toperturb（2.1）和（B.1）：Xt=Xti-1+Ztti-1b（r，Xr）dr+Ztti-1σ（r，Xr）dWr。

（C.1）Yi，[0]，t=bui+1（Xti）-Ztizi，[0]，rdWr。（C.2）对于t∈ Ii=【ti】-1，ti]，i∈ {1，···，n}。请注意，被引入X的方式，通过这种方式，我们对每个间隔Ii有不同的过程。在下文中，为了避免混淆指定区间的指数和x分量的指数∈ 对于后者，我们使用粗体哥特式符号s，如{i，j，····}，每个符号都贯穿1到d。引理C.1。（X，Yi，[0]，，Zi，[0]，）关于的经典导数kXt：=KkXt，kYi，[0]，：=KkYi，[0]，t，kZi[0]，：=KkZi，[0]，tfor k={1，2，3}由以下正向和反向SDE的解给出：X，它=Ztti-1.xjbi（r，Xr）X，jrdr+Ztti-1.σi（r，Xr）+(X，jr）xjσi（r，Xr）dWr，X，它=Ztti-1.xj bi（r，Xr）X，jr+xj，xkbi（r，Xr）X，jrX，krdr+Ztti-1.2(X，jr）xjσi（r，Xr）+(X，jr）xjσi（r，Xr）+(X，jr）(X，kr）xj，xkσi（r，Xr）dWr，X，它=Ztti-1.xj bi（r，Xr）X，jr+3xj，xkbi（r，Xr）X，jrX，kr+xj、xk、xmbi（r、Xr）X，jrX，krX，先生dr+Ztti-1.3(X，jr）xjσi（r，Xr）+3(X，jr）(X，kr）xj，xkσi（r，Xr）+(X，jr）xjσi（r，Xr）+3(X，jr）(X，kr）xj，xkσi（r，Xr）+(X，jr）(X，kr）(X，先生）xj，xk，xmσi（r，Xr）dWr，写Xi更合适，强调对区间t的依赖性∈ Ii，但我们省略了“i”，以减轻符号的重量。Yi，[0]，t=xj bui+1（Xti）X，jti-中兴通讯Zi，[0]，rdWr，Yi，[0]，t=xj bui+1（Xti）X，jti+xj、xkbui+1（Xti）X，jtiX，kti-中兴通讯Zi，[0]，rdWr，Yi，[0]，t=xj bui+1（Xti）X，jti+3xj、xkbui+1（Xti）(X，jti）(X，kti）+xj、xk、xmbui+1（Xti）(X，jti）(X，kti）(X，mti）-中兴通讯Zi，[0]，rdWr，用于t∈ Ii=【ti】-1，ti]。爱因斯坦约定用于{i，j，···}遍历1到d证明。

经典的可微性可以通过遵循[40]中定理3.1的公式来表示。有关更多详细信息，请参见[3 1]第6节。引理C.2。对于k={1，2，3}，存在一些（i，n）-独立的正常数Cp，ksuchkXp【ti】-1，ti]i≤ Cp、khkp/2适用于每个间隔Ii、i∈ {1，···，n}与任意p≥ 2.证明。这可以通过应用[31]附录A中给出的Lipschitz SDE的标准估计来说明。对于k=1，Eh||X| | p[ti-1，ti]i≤ CpEh公司Ztiti公司-1 |σ（r，Xr）| drp/2i≤ Cphp/2i。k=2时，一个为SEH||X| | p[ti-1，ti]i≤ CpEh公司Ztiti公司-1个|Xr | drp+Ztiti公司-1.|Xr |+|Xr|博士圆周率≤ Cp公司hpiEh公司||X| | 2 III+hp/2 EH||X| | pIi+||X| | 2pIii≤ Cphpi，根据需要。我们可以用类似的方式确定最后一种情况k=3。让我们用k介绍以下过程∈ {0，1，2}，X[k]t：=KkXt=0，Yi，[0]，[k]t：=KkYi[0]，t=0，Zi，[0]，[k]t：=KkZi，[0]，t=0和a lsoeYi，[0]t:=Xk=0k！Yi，[0]，[k]t，eZi，[0]t:=Xk=0k！Zi，[0]，[k]t（C.3），对于每个间隔t∈ 二、一∈ {1，···，n}。引理C.3。存在一些（i，n）独立的正常数Cpsuch，不等式彝文[0]-eYi，[0]p【ti】-1，ti]+Ztiti公司-1 | Zi，[0]r-eZi，[0]r | drp/2i≤ Cph3p/2适用于每个间隔Ii、i∈ {1，···，n}与任意p≥ 2.证明。

由于Θi，[0]，相对于的经典微分性，我们可以使用泰勒展开的残差公式；呃彝文[0]-eYi，[0]p【ti】-1，ti]+Ztiti公司-1 | Zi，[0]r-eZi，[0]r | drp/2i≤ CpEhsupr公司∈二、Z（1- ）Yi，[0]，rdp+Ztiti公司-1.Z（1- ）Zi，[0]，rd博士p/2i≤ CpZ公司呃Yi，[0]，p【ti】-1，ti]+Ztiti公司-1.Zi，[0]，r博士p/2id。应用Lipschitz BSDE的标准估计值（例如，参见[12]），以下各项的边界kxbui+1以及引理C.2，一个获得SEH彝文[0]-eYi，[0]p【ti】-1，ti]+Ztiti公司-1.Zi，[0]r-eZi，[0]r博士p/2i≤ CpZ公司呃||X| | pIi+||X| | pIi||X| | pIi+||X| | 3pIiid≤ Cph3p/2根据需要。最后一个引理意味着我们需要得到（eYi，[0]，eZi，[0]），这是（Yi，[0]，Zi，[0]）的二阶近似值。此外，正如我们接下来将看到的，由于共形展开引入的分级结构，这些BSDES的解可以通过简单的常微分方程显式获得。FBSDEs的相关系统总结如下：X【0】t=Xti-1+Ztti-1b（r，X[0]r）dr，X[1]，it=Ztti-1.xj bi（r，X[0]r）X[1]，jrdr+Ztti-1σi（r，X[0]r）dWr，X[2]，它=Ztti-1.xj bi（r，X[0]r）X[2]，jr+xj，xkbi（r，X[0]r）X[1]，jrX[1]，krdr+Ztti-12倍[1]，jrxjσi（r，X[0]r）dWr，Yi，[0]，[0]t=bui+1（X[0]ti）-Ztizi，[0]，[0]rdWr，（C.4）Yi，[0]，[1]t=xjbui+1（X[0]ti）X[1]，jti-Ztizi，[0]，[1]rdWr，（C.5）Yi，[0]，[2]t=xjbui+1（X[0]ti）X[2]，jti+xj、xkbui+1（X[0]ti）X[1]，jtiX[1]，kti-Ztizi，[0]，[2]rdWr，（C.6）用于t∈ 二、一∈ {1，···，n}与{i，j，···}的爱因斯坦约定。定义C.1。

（系数函数）我们定义了函数集χ：Ii×Rd→ Rd，y:Ii×Rd→ R、 y【1】：Ii×Rd→ Rd，y【2】：Ⅱ×Rd→ Rd，G【2】：Ⅱ×Rd→ Rd×d，y[2]：Ii×Rd→ R乘以χ（t，x）：=x+Ztti-1b（r，χ（r，x））dr，y（t，x）：=bui+1（χ（ti，x）），见备注3.3。y[1]j（t，x）：=xj bui+1（χ（ti，x））+Ztitxj bk（r，χ（r，x））y[1]k（r，x）dr，G[2]j，k（t，x）：=xj，xkbui+1（χ（ti，x））+Ztitnxb（r，χ（r，x））G[2]（r，x）<->j、 k级+xj，xkbm（r，χ（r，x））y[2]m（r，x）odr，y[2]（t，x）：=ZtitTrG[2]（r，x）[σσ]（r，χ（r，x））dr，y[2]=y[1]，表示（t，x）∈ Ii×Rd，i∈ {1，···，n}。我们使用了爱因斯坦约定和符号[xb（r，x）]i，j=xi bj（r，x），i，j∈ {1，···，d}. 我们用A表示s ymm量化<->:= A+A对于ad×d矩阵A。请注意，上述系数函数由给定x的ODE给出∈ Rdin每个时段。这些函数以以下方式表示BSDE的解：引理C.4。对于每个周期t∈ 二、一∈ {1，····，n}，BSDE（C.4），（C.5）和（C.6）的溶液由，根据爱因斯坦约定，Yi，，[0]t=y（t，Xti）给出-1） ，Zi，[0]，[0]t≡ 0（0阶）Yi，[0]，[1]t=y[1]j（t，Xti-1） X[1]，jt，Zi，[0]，[1]t=y[1]j（t，Xti-1） σj（t，χ（t，Xti-1） ），（一阶），最后，对于二阶yi，[0，[2]t=y[2]j（t，Xti-1） X[2]，jt+G[2]j，k（t，Xti-1） X[1]，jtX[1]，kt+y[2]（t，Xti-1） ，Zi，[0]，[2]t=2y[2]j（t，Xti-1） X【1】，ktxkσj（t，χ（t，Xti-1） ）+G[2]j，k（t，Xti-1） X[1]，jtσk（t，χ（t，Xti-1） （）.证据这是【31】第8节结果的特例。BSDE（C.4）、（C.5）和（C.6）的唯一解的存在性是非常明显的。

使用定义C.1中给出的ODE，将IT^o公式应用于建议的表格，并将结果与BSDE进行比较，可以直接检查表达式。由于每个区间II的跨度hi都很短，我们希望我们可以通过Euler方法的一个步骤来近似上述ODS，而不影响Le mma C.3中给出的误差顺序。定义C.2。（近似系数函数）我们定义了函数集；χ：Ii×Rd→ Rd，y:Ii×Rd→ R、 y【1】：Ii×Rd→ Rd，y【2】：Ⅱ×Rd→ Rd，G【2】：Ⅱ×Rd→ Rd×d，y[2]：Ii×Rd→ R乘以χ（t，x）：=x+（t） b（ti-1，x），y（t，x）：=bui+1（χ（ti，x）），y[1]j（t，x）：=xj-bui+1（χ（ti，x））+δ（t）xj bk（ti，χ（ti，x））xk bui+1（χ（ti，x）），G[2]j，k（t，x）：=xj，xkbui+1（χ（ti，x））+δ（t）n[xb（ti，χ（ti，x））]x、 xbui+1（χ（ti，x））<->j、 k级+xj，xkbm（ti，χ（ti，x））xmbui+1（χ（ti，x））o，y[2]（t，x）：=δ（t）TrG[2]（ti，x）[σσ]（ti，χ（ti，x））,andy[2]=y[1]，表示（t，x）∈ Ii×Rd，i∈ {1，···，n}。我们有爱因斯坦公约和符号（t） ：=t- ti公司-1，δ（t）：=ti- t、 因此，G[2]是对称矩阵值。定义C.2中的函数在以下意义上为定义C.1中的系数函数提供了良好的近似值：引理C.5。存在一些（i，n）独立的正常数CpsatisfyingE支持∈二、χ（t，Xti-（1）-χ（t，Xti-（1）p+支持∈二、y（t，Xti-（1）-y（t，Xti-（1）p+Xk=1上升∈二、y[k]（t，Xti-（1）-y[k]（t，Xti-（1）p+支持∈二、G[2]（t，Xti-（1）-G[2]（t，Xti-（1）p+支持∈二、y[2]（t，Xti-（1）-y[2]（t，Xti-（1）P≤ Cph3p/2i，对于每个间隔Ii，i∈ {1，···，n}与任意p≥ 2.证明。首先，让我们考虑一下（χ，χ）。

利用t中的1/2-H¨older连续性、x中b的全局Lipschitz和lineargrowth性质，我们得到了|χ（t，x）-χ（t，x）|≤Ztti公司-1 | b（r，χ（r，x））- b（ti-1，x）| dr≤ KZtti公司-1小时（r） 1/2+|χ（r，x）-χ（r，x）|+（r） | b（ti-1，x）|印尼盾≤ C（1+| x | h1/2i）h3/2i+KZtti-1 |χ（r，x）-χ（r，x）|由Gronwall不等式导出，支持∈Ii |χ（t，x）-χ（t，x）|≤ eKhiC（1+| x|√hi）h3/2i。ThusEhsupt公司∈Ii |χ（t，Xti-（1）-χ（t，Xti-1） | pi≤ Cph3p/2i1+马力/2 EH | Xti-1 | pi≤ Cph3p/2i（C.7），具有一些（i，n）独立的正常数Cp.自|xbui+1 |≤ 假设3.1（i）下的K′，| y（t，x）-y（t，x）|=| bui+1（χ（ti，x））- bui+1（χ（ti，x））|≤ K′|χ（ti，x）- χ（ti，x）|。因此，根据（C.7），Ehsupt∈Ii | y（t，Xti-（1）-y（t，Xti-1） | pi≤ Cph3p/2i（C.8），具有一些（i，n）独立的正常数Cp。现在让我们考虑一下[1]（t，x）=xbui+1（χ（ti，x））+δ（t）xb（ti，χ（ti，x））xbui+1（χ（ti，x））。因为两者|xbui+1 |和|xb |有界，很容易看到sup（t，x）∈Ii×Rd | y[1]（t，x）|≤ C（C.9）具有一些正常数C。

对于t∈ IIx给定∈ Rd，我们有[1]（t，x）-y[1]（t，x）=xbui+1（χ（ti，x））- xbui+1（χ（ti，x））+Ztitxb（r，χ（r，x））y[1]（r，x）- xb（ti，χ（ti，x））y[1]（ti，x）博士自（C.9），1/2-H¨older连续性和的全局Lipschitz性质xb，我们得到| y[1]（t，x）- y[1]（t，x）|≤ |xbui+1（χ（ti，x））- xbui+1（χ（ti，x））|+Ztitn|xb（r，χ（r，x））| | y[1]（r，x）-y[1]（r，x）|+|xb（r，χ（r，x））| | y[1]（r，x）- y[1]（ti，x）|+|xb（r，χ（r，x））- xb（ti，χ（ti，x））| | y[1]（ti，x）| odr≤ K′|χ（ti，x）- χ（ti，x）|+KZtit | y[1]（r，x）-y[1]（r，x）| dr+Chi+CZtitδ（r）1/2+|χ（r，x）-χ（r，x）|+|χ（r，x）- χ（ti，x）|博士≤ KZtit | y[1]（r，x）-y[1]（r，x）| dr+Ch3/2i1+| x |φ）。因此，后向Gronwall不等式（例如，参见[44]中的推论6.62）给出了假设∈Ii | y[1]（t，x）-y[1]（t，x）|≤ Ch3/2i（1+| x | phi）eKhi和henceEhsupt∈Ii | y[1]（t，Xti-（1）-y[1]（t，Xti-1） | pi≤ Cph3p/2i，（C.10），需要一些（i，n）独立的常数cpa。根据的有界性|mxbui+1（x）|和|mxb |带m∈ {1，2}，很容易看出| G[2]|是一个lsobonddsup（t，x）∈Ii×Rd | G[2]（t，x）|≤ C（C.11）和一些正常数C。使用（C.11）、1/2-H¨older和Lipschitz连续性对[1]进行了类似分析xb，xb，后向Gronwall不等式yieldssupt∈Ii | G[2]（t，x）-G[2]（t，x）|≤ Ch3/2i（1+| x | phi）和henceEhsupt∈Ii | G[2]（t，Xti-（1）-G[2]（t，Xti-1） | pi≤ Cph3p/2i（C.12），需要一些（i，n）独立的正常数cpa。最后，我们考虑[2]（t，x）=δ（t）TrG[2]（ti，x）[σσ]ti，χ（ti，x）.根据（C.11）和σ的线性增长性质，| y[2]（t，x）|≤ Cδ（t）1+| x|, （C.13）满足每（t，x）∈ Ii×Rd具有一些正常数C。

我们有[2]（t，x）-y[2]（t，x）=ZtitTrG[2]（r，x）[σσ]（r，χ（r，x））-G[2]（ti，x）[σσ]（ti，χ（ti，x））因此，drand | y[2]（t，x）- y[2]（t，x）|≤ZtitTrn公司|G[2]（r，x）-G[2]（r，x）|+| G[2]（r，x）- G[2]（ti，x）|[σσ]（r，χ（r，x））+G[2]（ti，x）[σσ]（r，χ（r，x））- [σσ]（ti，χ（ti，x））odr公司≤ Ch3/2i（1+| x |）+Chi（1+hi | x |）和一些（i，n）-独立常数C。因此，我们得到，对于任何p≥ 2、Ehsupt∈Ii | y[2]（t，Xti-（1）-y[2]（t，Xti-1） | pi≤ Cph3p/2i（C.14），根据需要。根据（C.7），（C.8），（C.10），（C.12），（C.14）和y[2]=y[1]，该权利要求得到了证明。现在我们介绍每个周期t的过程（bYi，[0]t，bZi，[0]t）∈ 二、它们由（C.3）中的（eYi，[0]t，eZi，[0]t）定义，定义C.1中的系数函数被定义C.2中的近似函数取代，即：。；bYi，[0]t：=y（t，Xti-1） +（X[1]t）y[1]（t，Xti-（1）+（X【2】t）y[2]（t，Xti-1） +（X[1]t）G[2]（t，Xti-1） X[1]t+y[2]（t，Xti-（1）, （C.15）bZi，[0]t：=y[1]（t，Xti-1） σ（t，χ（t，Xti-i） （）+（X【1】t）xσ（t，χ（t，Xti-1） ）y[2]（t，Xti-1） +（X[1]t）G[2]（t，Xti-1） σ（t，χ（t，Xti）-1） （）, （C.16）为了简单起见，我们使用了矩阵表示法。索引的细节可以从引理C.4中的thosegiven检查。引理C.6。存在一些（i，n）独立的正常数Cpsuch，不等式eYi，[0]-bYi，[0]p【ti】-1，ti]i+EheZi，[0]-bZi，[0]p【ti】-1，ti]i≤ Cph3p/2适用于每个间隔Ii、i∈ {1，···，n}对于任何p≥ 2.证明。它可以很容易地从引理C.2和C.5中显示出来。推论C.1。存在一些（i，n）独立的正常数Cpsuch thatEh | | Yi，[0]-bYi，[0]| | p[ti-1，ti]+Ztiti公司-1 | Zi，[0]r-bZi，[0]r | drp/2i≤ Cph3p/2适用于每个间隔Ii、i∈ {1，···，n}与任意p≥ 2.证明。

它直接来自引理C.3和C.6。自X[1]ti起-1=X【2】ti-1=0，我们有一个非常简单的表达式-1每个周期II=[ti-1，ti）：bYi，[0]ti-1=y（ti-1，Xti-1） +y[2]（ti-1，Xti-1） ，bZi，[0]ti-1=y[1]（ti-1，Xti-1） σ（ti-1，Xti-1） 。近似解（bYi，[0]，bZi，[0]）具有以下连续性：引理C.7。存在一些（i，n）独立的正常数Cpsuch，不等式支持∈二、bYi，[0]t-bYi，[0]ti-1.pi+Ehsupt∈二、bZi，[0]t-bZi，[0]ti-1.圆周率≤ Cphp/2i，每间隔Ii、i保持∈ {1，···，n}与任意p≥ 2.证明。自y（t，x）=y（ti-1，x）表示（t，x）∈ Ii×Rd，我们有-bYi，[0]ti-1=X【1】ty[1]（t，Xti-（1）+X【2】ty[2]（t，Xti-1） +X[1]tG[2]（t，Xti-1） X【1】t+y[2]（t，Xti-（1）-y[2]（ti-1，Xti-（1）.（C.9）（记住y[1]=y[2]）、（C.11）和（C.1 3）中的界限，以及表C.2中的估计值，意味着∈二、bYi，[0]t-bYi，[0]ti-1.圆周率≤ CpEh | | X[1]| | pIi+| | X[2]| | pIi+| | X[1]| | 2pIi+hpi1+| Xti-1 | 2p我≤ Cphp/2根据需要。同样，我们有| bZi，[0]t-bZi，[0]ti-1 |≤ |y[1]（t，Xti-（1）-y【1】（ti-1，Xti-1） | |σ（t，χ（t，Xti-1） ）|+| y[1]（ti）-1，Xti-1） | |σ（t，χ（t，Xti-1） （）- σ（ti-1，Xti-1） |+| X[1]t|xσ（t，χ（t，Xti-1） ）y[2]（t，Xti-1） +G[2]（t，Xti）-1） σt、 χ（t，Xti-（1）≤ C（t）1+| Xti-1个|+ C（t） 1/2+（t） （1+| Xti-1 |）+ C | X[1]t|1+| Xti-1个|,通过一些正常数C，我们得到了Ehsupt∈二、bZi，[0]t-bZi，[0]ti-1.圆周率≤ Cphp/2根据需要。C、 2（Yi，[1]，Zi，[1]）的近似值我们现在想要近似（Yi，Zi）分解中出现的剩余BSDE（B.2）。我们将在下文中看到，这可以通过非常简单的fa shion实现。我们定义了过程（bYi，[1]，bZi，[1]），bybYi，[1]t：=δ（t）fti公司-1，Xti-1，bYi，[0]ti-1，bZi，[0]ti-1., （C.17）bZi，[1]t：=0，（C.18）对于每个周期t∈ 二、一∈ {1，···，n}。这里，δ（t）=ti- t如前所述。引理C.8。

存在一些（i，n）独立的正常数Cpsuch，不等式易[1]-bYi，[1]p【ti】-1，ti]+Ztiti公司-1 | Zi，[1]r | drp/2i≤ Cph3p/2适用于每个间隔Ii、i∈ {1，···，n}与任意p≥ 2.证明。让我们把δYi，[1]t：=Yi，[1]t-bYi，[1]t，δZi，[1]t：=Zi，[1]t对于t∈ 二、那么，（δYi，[1]，δZi，[1]）是以下Lipschitz BSDE的解：δYi，[1]t=zitδf（r）dr-ZtitδZi，[1]rdWr，其中δf（r）：=f（r，Xr，Yi，[0]r，Zi，[0]r）-f（ti-1，Xti-1，bYi，[0]ti-1，bZi，[0]ti-1） 。根据假设2.2（ii），它满足正常数K |δf（r）|≤ |f（r，Xr，Yi，[0]r，Zi，[0]r）- f（ti-1，Xti-1，Yi，[0]r，Zi，[0]r）|+| f（ti-1，Xti-1、易、[0]r、子、[0]r）- f（ti-1，Xti-1，bYi，[0]r，bZi，[0]r）|+| f（ti-1，Xti-1、bYi、[0]r、bZi、[0]r）- f（ti-1，Xti-1，bYi，[0]ti-1，bZi，[0]ti-（1）|≤ K（r） 1/2+（1+| Yi，[0]r |+| Zi，[0]r |）| Xr- Xti公司-1个|+K | Yi，[0]r-bYi，[0]r |+K（1+| Zi，[0]r |+| bZi，[0]r |）| Zi，[0]r-bZi，[0]r |+K | bYi，[0]r-bYi，[0]ti-1 |+K（1+| bZi，[0]r |+| bZi，[0]ti-1 |）| bZi，[0]r-bZi，[0]ti-1 |。从引理B.1，我们知道| | Yi，[0]| | | S∞[技术信息-1，ti]+| | Zi，[0]| | Sp[ti-1，ti]≤ 任何p的CPP≥ 从（C.15），（C.16），引理C.2，以及引理C.5的证明中所示的（y，y[i]，G[2]）的有界性，一个类似的不等式| | bYi，[0]| | Sp[ti-1，ti）+| | bZi，[0]| | Sp[ti-1，ti]≤ Cpholds。Lipschitz SDE Ehsupt的连续性∈Ii | Xt- Xti公司-1 | pi≤ Cphp/2对于任何人来说都是众所周知的≥ 2.

那么，Eh | |δYi，[1]| | pIi+Ztiti公司-1 |δZi，[1]r | drp/2i≤ CpEh公司Ztiti公司-1 |δf（r）| dr圆周率≤ Cp（h3p/2i+hpiEh1+| | Yi，[0]| | 2II+| | Zi，[0]| | 4PIIHSUPT∈Ii | Xt- Xti公司-1 | 2pi+hpi呃| |易[0]-bYi，[0]| | pIii+Ehsupt∈二、bYi，[0]t-bYi，[0]ti-1.圆周率+Eh1+| | Zi，[0]| | 2pIi+| | bZi，[0]| | 2piieh希兹提提-1 | Zi，[0]r-bZi，[0]r | drpi+hpiEh1+| | bZi，[0]| | 2PI+| bZi，[0]ti-1 | 2piEhsupt∈二、bZi，[0]t-bZi，[0]ti-1.2pi）≤ Cph3p/2来自推论C.1和引理C.7。关于短期近似的主要结果可以总结为下一个定理C.1。根据假设2.1、2.2和假设3.1（i），过程（bYi、bZi）定义为（bYit：=bYi，[0]t+bYi，[1]t，bZit：=bZi，[0]t，t∈ Ii）是qg BSD E（3.1）和满意度的解（Yi，Zi）的短期近似值，一些（i，n）独立的正常数Cp，thatEh | | Yi-bYi | | p[ti-1，ti]+Ztiti公司-1 | Zir-bZir | drp/2i≤ Cph3p/2i，对于每个周期Ii，i∈ {1，···，n}和P≥ 2.证明。它直接来自命题B.1、推论C.1和引理C.8。确认该研究得到了金融高级研究中心（CARF）的部分支持。参考文献[1]Ankirchner，S.、Imkeller，P.和Dos Reis，G.，2007二次增长SDES的经典和变分可微性，概率电子杂志，第12卷，1418-1453。[2] Barrieu，P.和El Karoui，N.，2013，《二次半鞅的单调稳定性及其在无界广义二次BSD Es中的应用》，《概率年鉴》，第41卷，第3B期，1831-1863年。[3] Barthelmann，V.、Novak，E.和Ritter，K.，2000，《解析网格上的高维多项式插值》，计算数学进展12273-288。[4] Bender，C.和Denk，R.，2007，《反向SDE的正向方案，随机过程及其应用》，1171793-1812。[5] Bender，C。

和Steiner，J.，2012，《反向SDE的最小二乘蒙特卡罗》，《金融中的数值方法》（由Carmona e t al.编辑），257-289，柏林斯普林格。[6] Bally，V.，和Pag\'es，G.，2 003，《解决离散时间多维最优停止问题的量化算法》，Bernoulli，6，1003-1049。[7] Bianchetti，M和M Morini（编辑）（201 3），《金融危机后的利率建模》，RiskBooks，伦敦。[8] Bimit，J.M.，1973，《最优随机控制中的共轭凸函数》，数学杂志。肛门。Apl公司。44384-404。[9] Bouchard，B.和Chassagneux，J-F.，2008，《连续和离散反射BSDE的离散时间近似，随机过程及其应用》，1182269-2293。[10] Bouchard，B.和Elie，R.，2008，《带跳跃的解耦前向反向离散时间近似，随机过程及其应用》，118，53-75。[11] Bouchard，B.和Touzi，N.，2004，《反向随机微分方程的离散时间近似和蒙特卡罗模拟，随机过程及其应用》，111175-206。[12] Briand，P.、Delyon，B.、Hu，Y.、Pardoux，E.和Stica，L.，2003，《反向随机微分方程的Lpsolutions，随机过程及其应用》，108，109-129。[13] Briand，P.和Confortola，F.，2008，《希尔伯特空间中具有st-ochastic-Lipschitz条件和二次偏微分方程的BSDE，随机过程及其应用》，118，第8 18-838页。[14] Brigo，D，M Morini和A Pallavicini（2013年），《交易对手信用风险、抵押品和融资》，威利，西南苏联。[15] Chassagneux，J.F.和Richou，A.，2016年，Qu-adratic BSDE的数值模拟，应用概率年鉴，第26卷，第1期，262-304。[16] Chassagneux，J.F。

Richou，A.，2016，《切换问题中产生的反射BSDE的离散时间近似收敛速度》，工作论文，arXiv:1602:00015。[17] Cr'epey，S，T Bielecki与D Brigo的介绍性对话（201 4），交易对手风险与融资，CRC出版社，纽约。[18] Cr'epey，S.和So ng，S.，2016，《交易对手风险与风险：沉浸与超越》，《金融与随机》，20901-930。[19] Crisan，D.和Monolarakis，K.，2014，《反向SDE的二阶离散化和容积法模拟》，《应用概率年鉴》，第24卷，第2期，652-678。[20] Cvitani\'c，J.和Zhang，J.，2013，《连续时间方法中的契约理论》，柏林斯普林格。[21]德隆，L。，2013年，《带跳跃的倒向随机微分方程及其精算和金融应用》，斯普林格·维拉格，LN。[22]Delarue，F.和Menozzi，S.，2006年，《准线性模型的前向-后向随机算法》，《应用概率年鉴》，第16卷，第1期，140-184。【23】Delarue，F.和Menoz zi，S.，2008，《拟线性偏微分方程的插值随机算法》，《计算数学》，第77卷，261页，第125-158页。[24]E l Karoui，N.和Mazliak，l.（编辑），1997年，反向随机微分方程，AddisonWesley Longman有限公司，美国。。[25]E l Karoui，N.，Peng，S.和Quenez，M.C.，1997，《金融中的反向随机微分方程》，数学金融，第7卷，第1期，1-71。[26]E pstein，L.和Zin，S.，1989，《消费和资产回报的替代、风险规避和时间行为：理论框架》，生态计量学，57937-969。[27]藤井，M。

和Takahashi，A.，2012，《带摄动方案的非线性FBSDE解析近似》，国际理论与应用金融杂志，15，5125 0034（24）。[28]Fujii，M.，2014，《s-Tocastic filtering渐近展开的动量空间方法》，统计数学研究所年鉴，第66卷，93-120页。[29]Fujii，M.和Takahashi，A.，2015，《非线性FBSDE的微扰扩展技术与相互作用粒子方法》，亚太金融市场，第22、3、283-304卷。[30]Fujii，M.和Takahashi，A.，2017年，《带跳跃的二次指数增长BSDE及其马尔可夫差异性、随机过程及其应用》，2017年9月21日在线发布，https://doi.org/10.1016/j.spa2017.09.002。[31]Fujii，M.和Takahashi，A.，2015，《带跳跃的前后向SDE的渐近展开》，工作论文，CARF-F-372。arXiv中提供。【32】Gobet，E.，Lemor，J-P.a and Warin，X.，2005，《基于回归的蒙特卡罗方法求解反向随机微分方程》，《应用概率年鉴》，第15卷，第3期，21722202。[33]Imkeller，P.和Dos Reis，G.，2010，《截断二次增长BSDE的路径正则性和显式收敛速度，随机过程及其应用》，120348-379。定理5.520101202286-2288勘误表。[34]K azamaki，N.，1994，《连续指数鞅与BMO》，《数学课堂讲稿》，第1579卷，柏林斯普林格·维拉格出版社。[35]K loeden，P.和Platen，E.，1992，《随机微分方程的数值解》，数学应用（纽约）23。柏林斯普林格。【36】K obylanski，M.，2000，《具有二次增长的倒向随机微分方程和偏微分方程》，《生育年鉴》，第28卷，第2期，558-602。[37]K unito mo，N。

和Takahashi，A.，2003，《关于不一致索赔分析中渐近展开法的有效性》，《应用概率年鉴》，第13期，第3期，914-952页。[38]Ma，J.，Protter，P.和Yong，J.，1994，求解正反向随机微分方程多重-四步方案，概率论和相关领域，98339-359。[39]Ma，J.和Yong，J.，2000，《向前向后随机微分方程及其应用》，柏林斯普林格。[40]马，J.和张，J.，2002，倒向随机微分方程的表示定理，应用概率年鉴，12，41390-1418。【41】Ma，X.和Zabaras，N.，2009，《用于求解随机微分方程的自适应分层s解析网格配置算法》，计算物理杂志，2283084-3113。【42】Pag\'es，G.和Sagna，A.，2017年，《BSDE和非线性滤波的基于量化的数值模式、随机过程及其应用的改进误差界》，2017年7月5日发布，https://doi.org/10.1016/j.spa.2017.05.009.【43】Pardoux，E.和Peng，S.，1990年，《倒向随机微分方程的自适应解》，系统控制Lett。，14、55-61。【44】Pardoux，E.和Ra scanu，A.，2014，《随机微分方程，反向SDE，偏微分方程》，瑞士斯普林格国际出版社。【45】Sauer，T.，1995，《最小次数的多项式插值》，Numerische Mathematik，78，59-85。[46]Yong，J.a and Zhou，X.Y.，1999，随机控制：哈密顿系统和HJB方程，纽约州斯普林格。【47】Takahashi，A.，1999，《未定权益定价的渐进扩展方法》，亚太金融市场，6115-151。[48]Takahashi，A.和Yamada，T.，2015年，《具有扰动驱动的前向-后向SDE的渐近扩展》，即将发表在《国际金融工程杂志》上。【49】高桥，A。

和Yamada，T.，2016，《带渐近展开和多维Malliavin权重的弱近似》，《应用可能性年鉴》，第26卷，第2期，818-856。[50]Zhang，G.，Gunzburger，M.和Zhao，W.，2013，多维后向随机微分方程的稀疏网格方法，计算数学杂志，31，221248页。[51]Zhang，J.，2004，BSD Es的数值格式，应用概率年鉴，第14卷，第1459-488号。[52]Zhang，J.，2001，倒向随机微分方程的一些性质，普渡大学博士论文。