主要内容

simByTransition

用过渡密度模拟贝茨样例路径

描述

例子

路径) = simByTransition (MDLNPeriods模拟NTrials贝茨二元模型由两个布朗运动风险源和一个复合泊松过程驱动,表示重要事件的到来NPeriods连续观察时间。simByTransition用过渡密度近似连续时间随机过程。

例子

路径) = simByTransition (___名称,值除了前面语法中的输入参数外,还使用一个或多个名称-值对参数指定选项。

例子

全部折叠

用过渡密度模拟贝茨样例路径。

对象的参数贝茨对象。

AssetPrice = 80;返回= 0.03;JumpMean = 0.02;JumpVol = 0.08;JumpFreq = 0.1;= 0.04;水平= 0.05;速度= 1.0;波动率= 0.2;ρ= -0.7; StartState = [AssetPrice;V0]; Correlation = [1 Rho;Rho 1];

创建一个贝茨对象。

batesObj = bates(返回,速度,级别,波动性,...JumpFreq、JumpMean JumpVol,“startstate”StartState,...“相关”、相关)
贝茨batesObj =类贝茨:二元随机波动  -------------------------------------------------- 维度:状态= 2,布朗= 2  -------------------------------------------------- 开始时间:0 StartState: 2 x1双重数组关联:2 x2双重数组漂移:漂移率函数F (t) X (t))扩散:扩散速率函数G(t,X(t))模拟:模拟方法/函数simByEuler Return: 0.03 Speed: 1 Level: 0.05 Volatility: 0.2 JumpFreq: 0.1 JumpMean: 0.02 JumpVol: 0.08

定义仿真参数。

nPeriods = 5;%模拟未来5年的样本路径路径= simByTransition (batesObj nPeriods);路径
路径=6×280.0000 0.0400 66.0422 0.1012 73.8079 0.1243 48.9742 0.0571 49.9649 0.0638 58.9553 0.0467

输入参数

全部折叠

随机微分方程模型,指定为贝茨对象。有关创建贝茨对象,看到贝茨

数据类型:对象

模拟周期的数目,指定为一个正标量整数。的价值NPeriods确定模拟输出系列的行数。

数据类型:

名称-值对的观点

指定可选的逗号分隔的对名称,值参数。的名字参数名和价值为对应值。的名字必须出现在引号内。可以以任意顺序指定多个名称和值对参数Name1, Value1,…,的家

例子:(路径,次)= simByTransition(贝茨,NPeriods DeltaTimes, dt)

的模拟试验(样本路径)NPeriods观察的每个,指定为逗号分隔对组成“NTrials”一个正的标量整数。

数据类型:

观察之间的正时间增量,指定为逗号分隔的对,包括“DeltaTimes”标量或NPeriods——- - - - - -1列向量。

DeltaTime代表了熟悉的dt在随机微分方程中找到,并确定输出状态变量的模拟路径被报告的时间。

数据类型:

在每个时间增量内的中间时间步长的数目dt(定义为DeltaTimes),指定为逗号分隔的对,由“NSteps”一个正的标量整数。

simByTransition函数对每个时间增量进行分区dtNSteps小区间的长度dt/NSteps,并通过计算仿真状态向量来细化仿真NSteps−1中间点。虽然simByTransition不报告这些中间点的输出状态向量,通过使模拟更接近底层连续时间过程,改进提高了精度。

数据类型:

标志存储和返回方法,指示输出数组如何路径存储并返回,指定为逗号分隔的对,由“路径”和值为的标量逻辑标志真正的

  • 如果路径真正的(默认值)或未指定,则simByTransition返回路径作为一个三维时间序列阵列。

  • 如果路径(逻辑0),然后simByTransition返回路径输出数组作为一个空矩阵。

数据类型:逻辑

句号结束过程或状态矢量调整序列,指定为逗号分隔对,由“过程”和函数或单元格数组的形式

X t P t X t

simByTransition在每个观测周期结束时应用处理函数。处理函数接受当前的观测时间t和当前状态向量Xt,并返回可能调整输入状态的状态向量。

如果指定多个处理函数,simByTransition按函数在单元格数组中出现的顺序调用函数。

数据类型:细胞|函数

输出参数

全部折叠

相关状态变量的模拟路径,返回为(NPeriods + 1)——- - - - - -据nvar——- - - - - -NTrials三维时间序列阵列。

对于给定的试验,每一行路径是状态向量的转置吗Xt在时间t.当输入标志路径simByTransition返回路径作为一个空矩阵。

与模拟路径相关的观测时间,返回为(NPeriods + 1)——- - - - - -1列向量。的每个元素与对应的行路径

更多关于

全部折叠

过渡密度模拟

CIR SDE没有这样的解决方案rt) =fr(0)⋯)。

换句话说,这个方程不是显式可解的。然而,这个过程的跃迁密度是已知的。

分布的精确模拟rt_1)⋯,rt_n)是过程的一部分t_1,⋯,t_n对于相同的值r(0).这个过程的跃迁密度已知,表示为

r t σ 2 1 e α t u 4 α x d 2 4 α e α t u σ 2 1 e α t u r u t > u 在哪里 d 4 b α σ 2

贝茨模型

贝茨模型是二元复合模型。

每个Bates模型由两个耦合的单变量模型组成:

  • 几何布朗运动(“绿带运动”)模型具有随机波动函数和跳跃。

    d X 1 t B t X 1 t d t + X 2 t X 1 t d W 1 t + Y t X 1 t d N t

    这个模型通常对应一个价格过程,其波动率(方差率)由第二个单变量模型控制。

  • Cox-Ingersoll-Ross (圆形的)平方根扩散模型。

    d X 2 t 年代 t l t X 2 t d t + V t X 2 t d W 2 t

    该模型描述了耦合贝茨价格过程的变异率的演化。

参考文献

[1] Glasserman,保罗金融工程中的蒙特卡罗方法。纽约:Springer-Verlag, 2004。

Van Haastrecht, Alexander和Antoon Pelsser。Heston随机波动率模型的高效、几乎精确模拟国际理论与应用金融杂志。13日,没有。01(2010): 1。

介绍了R2020b