金融工程-维纳过程与伊藤引理.ppt

* * * * 经济和金融中的模拟方法 Monte Carlo 方法 在计量经济学里，如果我们对某种估计方法的统计性质不是很了解，而又要用到该种方法时，可以用Monte Carlo 方法来解决. 在计量经济学中的例子: 对联立方程偏误的定量研究. 确定Dickey-Fuller 检验的临界值. 确定在自相关检验中样本大小对检验功效的影响. 经济和金融中的模拟方法 Monte Carlo 方法 在金融中的例子: 奇异期权的定价. 确定宏观环境对金融市场的影响. 风险管理建模: 压力测试，例如，确定最小资本要求. 模拟中的“随机数” 进行蒙特卡罗模拟首先要设定数据生成系统。而设定数据生成系统的关键是要产生大量的随机数。例如模拟样本为100的随机趋势过程的DF统计量的分布，若试验1万次，则需要生成200万个随机数。 计量经济学中蒙特卡罗模拟和自举模拟所用到的随机数一般是服从N(0,1)分布的随机数。 计算机所生成的随机数并不是“纯随机数”，而是具有某种相同统计性质的随机数，即某种“伪随机数”（pseudo-random number）。生成随机数的程序称作“伪随机数生成系统”。实际上计算机不可能生成纯随机数。 模拟的计算机实现 蒙特卡罗模拟和自举模拟的实现要通过计算机编程来实现。 常用的软件有Mathematica，Gauss，Ox，EViews，Stata等。其原理基本一样。 若干例子见图。 图1 随机游走序列 图2 带趋势项的随机游走序列 图3 三维图圆环 图4 空间曲面 图5 投币1000次的概率值模拟 图6 生长曲线 图7 二元正态分布 图8 蒲丰问题 Dr. Fan 12.7 蒙特卡罗模拟的实现 我们从几个例子来看 例1：两个I(1)变量相关系数分布的蒙特卡罗模拟 未达到N 图11 蒙特卡罗模拟过程示意图 生成 xt, yt?I(1) 估计相关系数r 分析r的 分布 设定循环次数N 设定 xt,yt ?I(1) EViews程序如下： 定义一个非时间序列（u）工作文件，corr，容量为500。 定义一个空序列result，用来存储相关系数的计算结果。 !i为控制变量，通过一个for循环语句使计算进行500次。 把样本范围设置成100。 生成两个互不相关的白噪声序列x、y，样本容量100。 定义两个空的序列xx和yy，样本容量也是100。 定义两个标量sum1和sum2，初始值为0。 !counter为控制变量，在这个for循环中，分别对序列x和y进行一次累加生成两个一阶单整的序列，将结果分别放到序列xx和yy中。 累加一次。 计算序列xx和yy的相关系数，并将结果放到标量r中。 将相关系数计算结果放到序列result中，在这个for循环中，这个操作要进行500次。 显示序列result的直方图以及有关统计量。 图13 两个非相关I(1) 序列的相关系数的分布 例2：DW统计量分布的蒙特卡罗模拟 生成T=50的相互独立的IN(0,1)序列ut 和vt 用ut 和vt分别生成两个相互独立的I(1)序列 yt = yt-1 + ut , y0 = 0, xt = xt-1 + vt , x0 = 0, 估计模型yt = ?0 + ?1xt + wt 并计算残差 用残差计算DW统计量的值 存储2000个DW值 画DW频数分布直方图。记录T=50条件下DW分布的均值、标准差和第90、95、99百分位数。 分别估计DW均值、标准差和第90、95、99百分位数值对(1/T )的响应面函数 例3（利用模拟方法对欧式期权进行定价） 设股票价格St服从风险中性测度下的几何Brown运动： 其离散化形式为 根据金融工程理论，设现在股票价格为S0，T时刻到期（单位天），敲定价为K的欧式看涨期权的价格为 MC方案：按照（1）递推产生n条风险中性测度下的轨道，提取出ST (n)；（2） 对一个大众型欧式看涨期权的定价. 具体步骤如下: 确定标的资产的数据产生过程. 通常假设该过程为具有漂移的随机游走，即要确定漂移和波动参数. 同时要确定行权价格K 及到期日T. 产生T 个标准正态分布的数据，作为误