美式期权不能用蒙特卡洛模拟（蒙特卡洛模拟期权定价）

内容导航：
求程序对美式亚式期权的蒙特卡罗定价的代码！谢谢啦！！
为什么b-s定价模型不适用于美式期权？
怎么用 Excel 做蒙特卡洛模拟
如何在matlab中用蒙特卡洛模拟计算欧式期权价格
蒙特卡洛期权定价过程是二叉树吗
求程序对美式亚式期权的蒙特卡罗定价的代码！谢谢啦！！
看涨期权的定价公式

Q1：求程序对美式亚式期权的蒙特卡罗定价的代码！谢谢啦！！

用matlab的包吧，这里是官方介绍Matlab__asiansensbyls

asiansensbyls：“Calculate price and sensitivities for European or American Asian options using Monte Carlo simulation”

PriceSens = asiansensbyls(RateSpec,StockSpec,OptSpec,StrikeSettle,ExerciseDates)

代码那里没得选matlab，下面是我的例子：

Rates = 0.06;

StartDate = 'nov-1-2018';

EndDate = 'nov-1-2020';

RateSpec = intenvset('ValuationDate', StartDate, 'StartDates', StartDate, 'EndDates', EndDate,'Compounding', -1, 'Rates', Rates)

AssetPrice = 100;

Sigma = 0.25;

StockSpec = stockspec(Sigma, AssetPrice)

Settle = 'nov-1-2018';

ExerciseDates = 'nov-1-2020';

Strike = 100;

OptSpec = 'call';

NumTrials = 10000;

NumPeriods = 60;

AvgType = 'arithmetic';

Antithetic = true;

AmericanOpt = 1;

OutSpec = 'Price';

Price = asiansensbyls(RateSpec, StockSpec, OptSpec, Strike, Settle, ExerciseDates, 'NumTrials', NumTrials, 'NumPeriods', NumPeriods,'Antithetic', Antithetic, 'AvgType', AvgType, 'AmericanOpt', AmericanOpt, 'OutSpec',OutSpec)

另外，美亚式可以自己写二叉树来定价，参考Hull的“ Efficient procedures for valuing european adn american path-dependent options”。

Q2：为什么b-s定价模型不适用于美式期权？

因为BS模型是以行权日价格作为一个重要参数的，而与欧式期权确定在某一天行权不同，美式期权是可以在一段时间以内行权的，其行权期限较为灵活。

Q3：怎么用 Excel 做蒙特卡洛模拟

下面是在Excel中模拟一只股票价格的例子。假设股票价格
的对数收益率服从正态分布，均值为0，每日变动标准差为0.1，
模拟股票价格1年的路径，过程如下：
用到两个内置函数，即用rand()来产生0到1之间的随机数，然后用norminv()来获得服从既定分布的随机数，即收益率样本=norminv(rand(), 0, 0.1)。假定股票价格的初始值是100元，那么模拟的价格就是 S=100 * exp(cumsum(收益率样本))。
其中的cumsum（）不是Excel的内置函数，其意思就是收益率样本的累积，每个时刻的值都是当前样本及此前所有样本的和，如，收益率样本从单元格C3开始，当前计算C15对应的模拟价格，则模拟价格计算公式是：100 * exp(sum($C$3:C15))。
由此可以得到股票价格的一条模拟路径。

其他非正态分布也可以通过类似方式得到分布的抽样，即分布函数的逆函数，这些函数Excel都内置了。所以，做蒙特卡洛模拟的时候，关键是先确定所需模拟的分布，然后进行抽样，然后应用层面的各种公式就可以在抽样的基础上进行计算了。

--------以下是补充的--------

根据上面提到的思路，其实可以很便捷地为期权做定价。下面就用蒙特卡洛方法为一个普通的欧式看涨期权定价（蒙特卡洛在为普通期权plain vanilla option定价时不占优势，因为相对于解析法而言计算量很大。但是，如果要给结构比较复杂的奇异期权定价时，可能蒙特卡洛法就比较实用，有时可能成为唯一的方法）。

1）假设这个期权是欧式看涨期权，行权价格为50元，标的股票当前的价格也是50元，期权剩余时间是1天。

2）假设标的股票的价格服从对数正态分布，即股票的每日收益率服从正态分布，均值为0，每日标准差为1%。

根据分布假设，首先用rand()函数产生在0到1之间的均匀分布样本。为了提高精确度，这里抽样的数量为1000个（其实1000个是很少的了，通常需要10万个甚至50万个，但是在Excel表格中操作这么多数字，不方便，这是Excel的不足之处）。

下一步，用norminv(probability, mean, std)函数来获得股票收益率分布的1000个抽样，其中的probability参数由rand()产生的抽样逐个代入，mean=0.0, std = 0.01。注意这里抽样得到的日度收益率。也就是说，这个样本对应的下一个交易日股票价格的收益率分布。

下一步，股票价格=50×exp(收益率样本)，得到股票价格分布的抽样，有1000个样本。

根据我做的实验，这1000个样本的分布图形（histogram）跟对数正态分布是比较接近的，如下图所示：

图的横轴是股票价格，纵轴是样本中出现的频率。

得到了股票价格未来一天分布的样本之后，就可以以此样本来计算期权的价格了。

欧式看涨期权的定义为：

C=max(S-K,0)

所以，根据这个计算公式可以计算出在到期那天在特定的价格下期权的价值。在Excel中，相当于期权价值=max(股票价格样本 - 50，0)。由此就可以得到了该期权未来1天价值的样本。

然后，将未来价值贴现回来（用无风险利率贴现，假设无风险利率为0.05，则贴现公式是=exp(-0.05/360）×期权价值，得到期权价格的1000个样本。

最后，对期权价格的1000个样本求平均，Excel函数average(期权价格样本)，就可以得到期权的价格了。

我这里算出来的是：0.2015元。

而根据Black-Scholes期权定价公式算出来的理论价格则是0.2103元。二者比较接近，但是还是有差距。

而且，每次刷新Excel表格，就重新做一次模拟，得到的模拟价格变动比较大，有时是0.2043元，有时是0.1989元。由于这个抽样的数量比较小（1000个样本），所以估算的结果受到样本的影响会比较大。如果把抽样数量提高100倍甚至500倍，那么样本变动的影响可能会小一个或者两个数量级。但是计算量就大了，如果计算机性能不够高，那么利用Excel来做的话，比较困难。

这就是我的工作台：

------ 再来一个 --------

看到有人提到利用蒙特卡洛方法来估计圆周率Pi，挺有意思，也简单，所以就在Excel中做了一个实验。

基本原理在于在直角坐标系中的第一个象限中的一个单位圆，如下图所示：

在这个面积为1的正方形中，有四分之一的圆，圆的半径与正方向的边长都是1。那么根据圆的面积公式，这个图形中阴影部分的面积应该是 Pi/4。

下面开始进入蒙特卡洛的解法。

即，如果我们对这个正方形平面中的点进行均匀地抽样，随着抽样点的增多，那么落入阴影内的点的数量与总抽样数量的比，应该基本上等于阴影的面积Pi/4与整个正方形面积1的比，即Pi/4。用数学表示，就是

阴影内的样本点数量 ÷ 总数量 = Pi/4

所以，Pi = 4 × 阴影内的样本点数量 ÷ 总数量。

下面就在Excel中进行实验。

用rand()函数生成2000个随机数，作为随机样本点的X轴坐标，

再用rand()函数生成2000个随机数，作为随机样本点的Y轴坐标。

如此就得到了2000个随机样本点，这些点的X轴坐标和Y轴坐标都大于零且小于1，所以是在前面所说的正方形之中的点。

下一步，判断样本点是否处于阴影之内，由于这个阴影就是单位圆在直角坐标系第一想象的四分之一，所以圆阴影内的点都符合如下不等式：

翻译到Excel中，就是用IF函数来判断，例如：

IF(A2^2 + B2^2 <=1, 1, 0)

即，如果样本点在阴影中，得到1，否则得到0。这样就把样本点区分开来了。

最后，把所有得到的1和0加总，就知道所有样本点中处于阴影中样本点的数量了。

最后根据

Pi = 4 × 阴影内的样本点数量 ÷ 总数量

就可以算出Pi来了。

我这个试验中算出来的 Pi=3.142。

以下是样本点的散点图：

由于样本数量有限，所以计算出来的Pi的精度并不高。

以下是工作界面，挺简单的。

来源：知乎

Q4：如何在matlab中用蒙特卡洛模拟计算欧式期权价格

function [c,p]=ucoption(S,X,sigma,r,T,M)
sig2=sigma^2;
srT=sqrt(T);
srTa=sigma*srT;
c=0;
p=0;
for i=1:M
ST=S*exp((r-0.5*sig2)*T+srTa*randn);
c=c+max(ST-X,0);
p=p+max(X-ST,0);
end
c=c/M;
p=p/M;
[Call,Put] = blsprice(S, X, r, T, sigma);
error=[c,p]-[Call,Put]
%可以试试 [c,p]=ucoption(10,10,0.3,0.05,0.5,10^4*100);

Q5：蒙特卡洛期权定价过程是二叉树吗

蒙特卡罗方法的解题过程可以归结为三个主要步骤:构造或描述概率过程;实现从已知概率分布抽样;建立各种估计量。
蒙特卡罗方法解题过程的三个主要步骤:
(1)构造或描述概率过程
(2)实现从已知概率分布抽样
(3)建立各种估计量
应用到期权上一定程度上你可以这么理解，但不完全相同，因为有的时候会过于简单，蒙特卡罗过程如果本身的设定是偏离实际的，会没有意义，所以二叉树是一种比较理想的状态而已。

Q6：求程序对美式亚式期权的蒙特卡罗定价的代码！谢谢啦！！