金融工程维纳过程与伊藤引理PPT课件

1 教学目的与要求 掌握随机变量的概念，了解马尔科夫过程的特点，掌握维纳过程的特点和性质，掌握一般维纳过程的特征以及其漂移率和方差率，维纳过程的均值和标准差。掌握Ito过程的特征。 第1页/共47页2 教学重点及难点 一、马尔科夫过程与效率市场的关系。 二、维纳过程、一般维纳过程与此同时Ito过程的特征，漂移率和方差率，变量的均值与方差。以及这几种过程的内在联系和变化。 三、Ito定理及其运用。 第2页/共47页3期权的估值 欧式期权的到期收益Max (STX, 0)ST不确定，所以期权到期的收益也不确定。 期权当期的价值？ 风险中性估值 期权当期的价值未来收益折现后的期望值cE Max (STX, 0) 问题ST的分布是怎样的？ 只有确定ST的分布才能确定c的价值第3页/共47页412.1弱式效率市场假说与马尔可夫过程 效率市场假说 1965年，法玛(Fama)提出了著名的效率市场假说。该假说认为： 投资者都力图利用可获得的信息获得更高的报酬； 证券价格对新的市场信息的反应是迅速而准确的，证券价格能完全反应全部信息； 市场竞争使证券价格从一个均衡水平过渡到另一个均衡水平，而与新信息相应的价格变动是相互独立的。 效率市场分类 效率市场假说可分为三类：弱式、半强式和强式。 弱式效率市场假说认为，证券价格变动的历史不包含任何对预测证券价格未来变动有用的信息，也就是说不能通过技术分析获得超过平均收益率的收益。第4页/共47页5 半强式效率市场假说认为，证券价格会迅速、准确地根据可获得的所有公开信息调整，因此以往的价格和成交量等技术面信息以及已公布的基本面信息都无助于挑选价格被高估或低估的证券。 强式效率市场假说认为，不仅是已公布的信息，而且是可能获得的有关信息都已反映在股价中，因此任何信息(包括“内幕信息”)对挑选证券都没有用处。 效率市场假说提出后，许多学者运用各种数据对此进行了实证分析。结果发现，发达国家的证券市场大体符合弱式效率市场假说。第5页/共47页6 马尔可夫过程 弱式效率市场假说可用马尔可夫随机过程(Markov Stochastic Process)来表述。 马尔科夫过程(Markov process)是一种特殊类型的随机过程。 未来的预测只与变量的当前值有关，与变量过去的历史和变量从过去到现在的演变方式不相关。 股价的马尔科夫性质与弱型市场有效性(the weak form of market efficiency)相一致： 一种股票的现价已经包含了所有信息，当然包括了所有过去的价格记录。 如果弱型市场有效性正确的话，技术分析师可通过分析股价的过去历史数据图表获得高于平均收益率的收益是不可能的。 是市场竞争保证了弱型市场有效性成立。 第6页/共47页712.2维纳过程（ Wiener Process ） 布朗运动起源于物理学中对完全浸没于液体或气体中的小粒子运动的描述，以发现这种现象的英国植物学家Robert Brown命名。 描述布朗运动的随机过程的定义是维纳(wiener)给出的，因此布朗运动又称维纳过程 股价行为模型通常用布朗运动来描述。 布朗运动是马尔科夫随机过程的一种特殊形式。第7页/共47页8维纳过程（ Wiener Process ） 维纳过程（Wiener Process） 性质一：股票价格的变动是一个正态变量与时间的乘积 （服从标准正态分布） 性质二：任意两个不重叠时段的股票价格变动相互独立 从性质一，我们知道z服从正态分布，性质2则隐含z遵循马尔科夫过程。 维纳过程/布朗运动的特征 股票价格在任意时段变动的均值都为0。 股票价格在某一时段变动的方差等于时间的长度tz0),(21zzCorr第8页/共47页9程序：维纳过程的模拟 假定股票价格服从普通布朗运动，即dS=dt+dz，其中和均为常数，dz遵循标准布朗运动，也就是说，在短时间t后，S值的变化值S为 假定股票价格服从几何布朗运动，即dS=dt+dz,其中和均为常数，dz遵循标准布朗运动，也就是说，在短时间t后，ln(S)值的变化值ln(S)为Stt 0,TSSNTT0,TSN STT20lnln(/ 2) ,TSSNTT20lnln(/ 2) ,TSNSTT第9页/共47页当股票价格服从普通布朗运动时的走势图104242.54343.5st0.511.52t第10页/共47页当股票价格服从几何布朗运动时的走势图11406080100120140g_st0.511.52t第11页/共47页12股票价格的一般变动 一般化的维纳过程 变量本身随着时间的推移会有定量的增长at 除了时间价值之外的变动为布朗运动zbtax第12页/共47页1312.3 股票价格的一般变动 股票价格的变动 股票价格有随时间推移增长的稳定趋势 股票“实际”价格变动为布朗运动zStSS第13页/共47页布朗运动股票价格tWStt第14页/共47页指数布朗运动股票价格)exp(tWStt第15页/共47页上证指数第16页/共47页1712.4 Itos Lemma Itos Lemma 假设存在一个伊藤过程： 如果G是x和t的函数，即：G=G(x,t) 那么： 期权及其他衍生证券的价格变动 股票价格服从维纳过程： 那么：ztxbttxax),(),(zbxGtbxGtGaxGG22221zStSSzSSftSSftfSSff222221第17页/共47页18证明：如前述，假设标的资产价格变动过程服从：其中其中利用泰勒展开，忽略高阶项，利用泰勒展开，忽略高阶项，G(x,t)可以展开为可以展开为tzztxbttxax),(),(),(212121),(2222222txotxtxGxxGttGxxGttGtxG第18页/共47页1932200limlim0ttx ta tbt 因此，上式可以改写为因此，上式可以改写为22201lim2tfffftxxtxx 保留保留1阶项，忽略阶项，忽略1阶以上的高阶项阶以上的高阶项第19页/共47页202200322222022limlim lim2 tttxa tbtatbtabtbt 其中（忽略高阶项）：其中（忽略高阶项）：第20页/共47页2122222()()()ExE btbtE因此，可得22(0,1),( )(0) ()1NDEE由于则22()Exbt由此得到代入前述公式可得到伊藤引理。第21页/共47页2212.5 股票价格的对数正态特性 对数正态分布 股票价格服从维纳过程 股票价格的分布为对数正态分布 公式ttSS),(ttSStTtTSST),)(2(lnln20第22页/共47页23关于对数正态分布 定义G=lnS，由于： 所以有： 即： 显然G为一个广义维纳过程，其漂移率为常数 ，波动率为常数 。 因此，lnS的变化服从正态分布，不难知道：0,11222tGSSGSSG，dzdtdG)2(2)2(2222lnln(/2) ,0, TsNSTt tT dzdtSd)2(ln2第23页/共47页24对数正态分布第24页/共47页几何布朗运动的深入分析 在很短的时间t后，证券价格比率的变化值 为： 可见，在短时间内， 具有正态分布特征 其均值为 ，标准差为 ，方差为 。 SSttSS (,)SttS tt2t第25页/共47页几何布朗运动的深入分析（2） 但是，在一个较长的时间T后， 不再具有正态分布的性质： 多期收益率的乘积问题，服从正态分布的变量的乘积并不服从正态分布。而由于总的连续复利收益率等于各期收益率的加和，因此仍为正态分布。 因此，尽管是短期内股票价格百分比收益率的标准差，但是在任意时间长度T后，这个收益率的标准差却不再是 。股票价格的年波动率并不是一年内股票价格百分比收益率变化的标准差。SST第26页/共47页几何布朗运动的深入分析（3） 如果股票价格服从几何布朗运动，则可以利用Ito引理来推导证券价格自然对数lnS所遵循的随机过程： 这个随机过程的特征： 普通布朗运动：恒定的漂移率和恒定的方差率。 在任意时间长度T之后，G的变化仍然服从正态分布，均值为 ，方差为 。标准差仍然可以表示为 ，和时间长度平方根成正比。 从自然对数lnS所遵循的这个随机过程可以得到两个结论:2()2dGdtdz2(/2)()Tt2()T ttT22 ()(),GTtTt 第27页/共47页（1）几何布朗运动意味着股票价格服从对数正态分布。 令t时刻G的值为lnS，T时刻G的值为lnST，其中S表示t时刻（当前时刻）的证券价格，ST表示T时刻（将来时刻）的证券价格，则在Tt期间G的变化为： 这意味着： 进一步从正态分布的性质可以得到 也就是说，证券价格对数服从正态分布。如果一个变量的自然对数服从正态分布，则称这个变量服从对数正态分布。这表明ST服从对数正态分布。 这正好与作为预期收益率的定义相符。lnlnTSS22lnln ()(),TSSTtTt 22ln ln()(),TSSTtTt()()T tTE SSe222()()var()1T tT tTSS ee第28页/共47页（2）股票价格对数收益率服从正态分布 由于dG实际上就是连续复利的对数收益率。因此几何布朗运动实际上意味着对数收益率遵循普通布朗运动，对数收益率的变化服从正态分布，对数收益率的标准差与时间的平方根成比例。 将t与T之间的连续复利年收益率定义为，则22t22lnln ()(1e), (n,lt)TTTSSSSSTtTStTt （T- ）， 可得T-由第29页/共47页结论 几何布朗运动较好地描绘了股票价格的运动过程。第30页/共47页3112.6 随机过程的蒙特卡罗模拟 有关蒙特卡罗方法的由来 取名于摩纳哥的著名赌城 掷色子是一个随机事件 蒙特卡罗方法 任何涉及随机采样的数值方法 不仅仅用于有关随机的问题 估计 圆周率 优化问题 40年代美国Los Alamos 实验室的科学家用于核武器的研究 代表人物：冯诺依曼第31页/共47页经济和金融中的模拟方法 Monte Carlo 方法在计量经济学里，如果我们对某种估计方法的统计性质不是很了解，而又要用到该种方法时，可以用Monte Carlo 方法来解决.在计量经济学中的例子:1.对联立方程偏误的定量研究.2.确定Dickey-Fuller 检验的临界值.3.确定在自相关检验中样本大小对检验功效的影响.第32页/共47页经济和金融中的模拟方法 Monte Carlo 方法在金融中的例子:1.奇异期权的定价.2.确定宏观环境对金融市场的影响.3.风险管理建模: 压力测试，例如，确定最小资本要求.第33页/共47页模拟中的“随机数” 进行蒙特卡罗模拟首先要设定数据生成系统。而设定数据生成系统的关键是要产生大量的随机数。例如模拟样本为100的随机趋势过程的DF统计量的分布，若试验1万次，则需要生成200万个随机数。 计量经济学中蒙特卡罗模拟和自举模拟所用到的随机数一般是服从N(0,1)分布的随机数。 计算机所生成的随机数并不是“纯随机数”，而是具有某种相同统计性质的随机数，即某种“伪随机数”（pseudo-random number）。生成随机数的程序称作“伪随机数生成系统”。实际上计算机不可能生成纯随机数。第34页/共47页模拟的计算机实现 蒙特卡罗模拟和自举模拟的实现要通过计算机编程来实现。 常用的软件有Mathematica，Gauss，Ox，EViews，Stata等。其原理基本一样。 若干例子见图。第35页/共47页 图1 随机游走序列 图2 带趋势项的随机游走序列-202-202-1-0.500.51-202 图3 三维图圆环 图4 空间曲面第36页/共47页 图5 投币1000次的概率值模拟 图6 生长曲线 图7 二元正态分布 图8 蒲丰问题第37页/共47页3812.7 蒙特卡罗模拟的实现 我们从几个例子来看 例1：两个I(1)变量相关系数分布的蒙特卡罗模拟 未达到未达到N图图11 蒙特卡罗模拟过程示意图蒙特卡罗模拟过程示意图生成生成 xt, yt I(1) 估计相关估计相关系数系数r 分析分析r的的 分布分布 设定循环设定循环次数次数N 设定设定 xt,yt I(1) 第38页/共47页EViews程序如下：workfile corr u 1 500series resultfor !i=1 to 500smpl 1 100series x=nrndseries y=nrndseries xxseries yyscalar sum1=0scalar sum2=0for !counter=1 to 100sum1=sum1+x(!counter)sum2=sum2+y(!counter) xx(!counter)=sum1yy(!counter)=sum2nextscalar r=cor(xx,yy)result(!i)=rnextresult.hist 定义一个非时间序列（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的直方图以及有关统计量。第39页/共47页图13 两个非相关I(1) 序列的相关系数的分布第40页/共47页例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 )的响应面函数的响应面函数第41页/共47页 例3（利用模拟方法对欧式期权进行定价）设股票价格St服从风险中性测度下的几何Brown运动：ttttdSrS dtS dB其离散化形式为其离散化形式为1(1)(0,1)(1)iiiiiSr SSBBN根据金融工程理论，设现在股票价格为根据金融工程理论，设现在股票价格为S0，T时时刻到期（单位天），敲定价为刻到期（单位天），敲定价为K的欧式看涨期权的欧式看涨期权的价格为的价格为rTTCeESKMC方案：按照（方案：按照（1）递推产生）递推产生n条风险中性测度下的条风险中性测度下的轨道，提取出轨道，提取出ST (n)；（；（2） 11nirTTiCeSKn第42页/共47页43 对一个大众型欧式看涨期权的定价. 具体步骤如下: 确定标的资产的数据产生过程. 通常假设该过程为具有漂移的随机游走，即要确定漂移和波动参数. 同时要确定行权价格K 及到期日T. 产生T 个标准正态分布的数据，作为误差项, ut N(0,1). 构造T 个标的资产的观测值.第43页/共47页44 记录在时刻 T时标的资产的价格ST. 对一个看涨期权如果 ST K, 则价值为ST - K, 然后用无风险利率贴现. 重复 1到 4步 N 次, 取 N 重复的平均值，这个平均值就是该欧式看涨期权的价格. 对于更加复杂期权的定价可以按同样思路进行.第44页/共47页有关随机数发生器 随机数产生器能够产生在区间0,1上均匀分布的随机数 Excel中 在excel中可以采用RAND（）产生一个0和1之间的随机数， 累积正态分布分布的反函数为NORMSINV。 所以，EXCEL中从标准正态分布中随机抽样的方法是：NORMSINV（RAND（） STATA中 生成（0，1）之间伪随机数的命令为uniform 每运行”di uniform()”命令一次，就可以得到一个随机数 . simulate exp_list , reps(#) saving(filename , replace) seed(#) : command 例 simulate max=r(max), reps(10000) nodots:seq3，效果就是重复一万次，取平均数第45页/共47页有关随机数发生器 例：虚假回归。随机生成两个随机游走序列（100个样本），进行回归， 计算估计量和t统计量；如此反复10000次。 （1）利用均匀分布或正态分布随机数生成器生成随机游走过程。 （2）通过如下命令进行模拟。 . simulate beta=(r(b1) se=(r(se1), reps(200) : rdwalk, diff( 1 1) （3）观察其分布，并与标准正态分布相比较。 . gen t=beta/se . histogram t . twoway histogram t | function y=normalden(x,0,1), range(-3 3) 第46页/共47页Dr. Fan47感谢您的观看。第47页/共47页