732 2

[有偿编程] 想询问关于这段程序的一些语句,有关最小二乘美式卖权的定价 [推广有奖]

  • 0关注
  • 0粉丝

初中生

23%

还不是VIP/贵宾

-

威望
0
论坛币
944 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
182 点
帖子
6
精华
0
在线时间
17 小时
注册时间
2016-12-24
最后登录
2020-10-27

100论坛币
程序如下,是利用最小二乘法和蒙特卡洛计算美式看跌期权的价格

fhandles = {@(x)x, @(x)x.^2};NSteps=5;NRepl=10000;
S0=83.5;K=80;r=0.1;sigma=0.35;T=1;
dt = T/NSteps;
discountVet = exp(-r*dt*(1:NSteps)');
NBasis = length(fhandles);        %基函数个数
% generate sample paths
SPaths=AssetPaths(S0,r,sigma,T,NSteps,NRepl);%价格路径
CashFlows = max(0, K - SPaths(:,NSteps));%
ExerciseTime = NSteps*ones(NRepl,1);%
for step = NSteps-1:-1:1
    InMoney = find(SPaths(:,step) < K);
    XData = SPaths(InMoney,step);
    RegrMat = zeros(length(XData), NBasis);%初始化
    for k=1:NBasis
        RegrMat(:, k) = feval(fhandles{k}, XData);
    end
    YData = CashFlows(InMoney).*discountVet(ExerciseTime(InMoney)-step);
    alpha = RegrMat \ YData;
    IntrinsicValue = K - XData;
    ContinuationValue = RegrMat * alpha;
    Index = find(IntrinsicValue > ContinuationValue);
    ExercisePaths = InMoney(Index);%
    CashFlows(ExercisePaths) = IntrinsicValue(Index);%
    ExerciseTime(ExercisePaths) = step;%
end
price = max(K-S0, mean(CashFlows.*discountVet(ExerciseTime)));

想问的是YData那一句是怎么算,按我的理解是就倒着循环每一列,假如现在step=2,就是找到第二期的inmoney的行,然后将最后一列的值折现到这一期,计算出来的YData就相当于等待的价值。就相当于这个程序只要等待就是等到最后一期,所以等待的价值才是最后一期往回折现,可以这样理解吗?但是美式期权每期都有可能执行或者等待啊?所以想请教一下论坛的大佬这句到底要怎么理解?




因为急着交作业所以悬赏100个币,之前没有学过编程之类的,所以读别人程序有点困难,希望能得到大佬们的帮助,真的很急

使用道具

这个是生成价格路径的代码
function SPaths=AssetPaths(S0,mu,sigma,T,NSteps,NRepl)
%生成价格路径
SPaths = zeros(NRepl, 1+NSteps);
SPaths(:,1) = S0;
dt = T/NSteps;
nudt = (mu-0.5*sigma^2)*dt;
sidt = sigma*sqrt(dt);
for i=1:NRepl
   for j=1:NSteps
      SPaths(i,j+1)=SPaths(i,j)*exp(nudt + sidt*randn);
   end
end
%作图
x1=0:dt:T;
y1=SPaths';
y2=mean(SPaths);
figure
subplot(2,1,1);
plot(x1,y1)
subplot(2,1,2);
plot(x1,y2)
xlabel('期权存续期间')
ylabel('股票价格模拟路径')
SPaths(:,1) = [];

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-25 02:25