亚式期权定价蒙特卡洛模拟matlab代码

以下是一个使用MATLAB进行亚式期权定价的示例代码，其中使用几何布朗运动模型来模拟标的资产价格的随机漫步，使用欧式平均价格亚式期权进行定价。

S0 = 100; % 标的资产初始价格
K = 110; % 期权执行价格
r = 0.05; % 无风险利率
T = 1; % 期权到期时间
sigma = 0.2; % 标的资产波动率
n = 100; % 模拟步长
N = 10000; % 模拟次数

dt = T/n; % 模拟步长

% 生成随机价格路径
S = zeros(N, n+1);
S(:,1) = S0;
for i = 1:N
    for j = 2:n+1
        S(i,j) = S(i,j-1)*exp((r-sigma^2/2)*dt+sigma*sqrt(dt)*randn);
    end
end

% 计算亚式期权的支付值
A = mean(S(:,2:end), 2); % 计算平均价格
payoff = max(A-K, 0); % 支付值

% 计算亚式期权现值
P = exp(-r*T)*mean(payoff);

fprintf('亚式期权定价结果为：%.4f\n', P);

输出结果可能会略有不同，由于随机性，每次运行结果可能会有所不同。