MATLAB循环修改变量异常：European Put估值绘图问题

我来帮你排查这个MATLAB循环的问题——这种遍历M值计算欧式看跌期权价值的场景，我碰到过不少常见坑，咱们一步步来捋：

1. 最常见的坑：未提前存储所有计算结果 #

很多人会犯的错误是：循环里只更新单个变量，最后绘图时只用了循环结束后的最后一组值，导致图像不符合预期。比如错误写法：

for M = 20:5:250
    put_value = calculate_european_put(M); % 你的期权计算函数
end
plot(M, put_value) % 这里只会画出最后一个M对应的点

解决方法：提前初始化数组存储所有M值和对应的期权价值，循环中逐个赋值：

% 先定义M的取值序列
M_list = 20:5:250;
% 初始化数组存储结果，长度和M序列一致
put_values = zeros(size(M_list));

for i = 1:length(M_list)
    current_M = M_list(i);
    put_values(i) = calculate_european_put(current_M); % 传入当前M计算
end

% 绘图时用完整的序列
plot(M_list, put_values, '-o')
xlabel('M')
ylabel('European Put Value')
title('Put Value vs M')

2. 计算函数的隐形依赖问题 #

如果你的calculate_european_put函数里用到了全局变量或者persistent持久化变量，可能会导致后续循环调用时，没有更新与M相关的参数。比如：

function put_val = calculate_european_put(M)
    persistent S K r sigma T; % 持久化了其他参数
    % 如果这些参数只初始化过一次，后续M变化时，内部的离散化步骤（比如dt=T/M）可能没正确更新
end

解决方法：把所有需要的参数（比如标的价格S、执行价K等）都作为函数输入传递，避免依赖全局/持久化变量，确保每次调用都是独立计算：

function put_val = calculate_european_put(M, S, K, r, sigma, T)
    dt = T / M; % 每次用当前M计算步长
    % 后续的二叉树/蒙特卡洛计算逻辑...
end

3. 循环范围或步长的验证 #

先在命令行直接输入disp(20:5:250)，确认生成的序列是否符合预期——比如是否从20开始，步长5，最后一个值是250（20+5*46=250，没问题）。如果你的计算逻辑对M有额外限制（比如必须是偶数、大于某个阈值），要确保序列里的所有值都满足条件。

4. 调试小技巧：打印中间结果 #

在循环里加入打印语句，查看每一组M对应的计算结果，定位是循环逻辑问题还是计算函数的问题：

for i = 1:length(M_list)
    current_M = M_list(i);
    put_values(i) = calculate_european_put(current_M);
    fprintf('M=%d, 期权价值=%.4f\n', current_M, put_values(i));
end

如果某几个M的结果明显异常，就去检查计算函数在该M值下的逻辑。

完整示例框架 #

给你一个可参考的完整代码结构（包含二叉树法计算欧式看跌期权的示例逻辑，你可以替换成自己的实现）：

% 欧式看跌期权计算函数（示例，替换成你基于文献的实现）
function put_val = european_put(M, S, K, r, sigma, T)
    dt = T / M;
    u = exp(sigma * sqrt(dt));
    d = 1 / u;
    p = (exp(r*dt) - d) / (u - d);
    
    % 计算期末期权价值
    ST = K * d.^(M:-1:0) .* u.^(0:M);
    put_T = max(K - ST, 0);
    
    % 倒推当前价值
    for i = M:-1:1
        put_T = exp(-r*dt) * (p * put_T(2:end) + (1-p) * put_T(1:end-1));
    end
    put_val = put_T;
end

% 主程序：遍历M并绘图
clear; clc; close all;

% 固定参数（根据你的需求修改）
S = 100; K = 100; r = 0.05; sigma = 0.2; T = 1;

M_range = 20:5:250;
put_values = zeros(size(M_range));

for idx = 1:length(M_range)
    current_M = M_range(idx);
    put_values(idx) = european_put(current_M, S, K, r, sigma, T);
end

% 绘图美化
plot(M_range, put_values, '-b', 'LineWidth', 1.5);
hold on;
scatter(M_range, put_values, 50, 'r', 'filled');
hold off;

xlabel('Parameter M', 'FontSize', 12);
ylabel('European Put Option Value', 'FontSize', 12);
title('European Put Value vs Parameter M', 'FontSize', 14);
grid on;
set(gca, 'FontSize', 10);

先从检查数组初始化和变量存储开始，再结合打印调试定位问题，应该能快速解决循环的故障。

内容的提问来源于stack exchange，提问作者Axion004