MATLAB循环扩展迭代范围后部分结果异常，求原因分析

问题根源 #

你的代码中，cell数组yv在外层kx循环中未被重置。首次迭代后yv被赋值，后续循环若不主动清空，会保留之前迭代的内容。当迭代到p=0.10时，yv{1}的长度可能与之前迭代不同，导致后续ndgrid、combinations计算混入旧数据，进而引发概率计算错误。而单独运行某一步时，yv是全新创建的，无旧数据干扰，因此结果正确。

修复方案 #

在每次外层kx循环的起始位置，添加yv = cell(size(v));重置yv，确保每次迭代都用全新的cell数组存储当前p对应的结果。

修改后的完整代码 #

v = [0.5];
pv = 0.01:0.01:0.1;  % 原问题版本
n = 1000;
h = 3.841;
fp = zeros(size(pv)); % 预分配数组避免尺寸问题

for kx = 1:length(pv)
    p = pv(kx); 
    yv = cell(size(v)); % 新增：重置yv，清除旧数据干扰

    % 重新计算依赖p的函数
    a = @(y) y;
    b = @(x) x;
    c = @(x, y) p*n - y;
    d = @(x) n*(1 - p) - x;

    new = @(x, y) (n.*(a(y).*d(x) - b(x).*c(x, y)).^2) ./ ...
        ((a(y) + b(x)) .* (c(x, y) + d(x)) .* (a(y) + c(x, y)) .* (b(x) + d(x))) - h;

    for j = 1:length(v)    
        yp = 0:p*n;
        xp = v(j).*n - yp;    
        f = new(xp,yp);
        fin = find(f<0);
        yv{j} = yp(fin);    
    end

    % 动态为所有cell元素生成网格
    grid_cells = cell(1, numel(yv));
    [grid_cells{:}] = ndgrid(yv{:});

    % 将网格转换为组合矩阵
    combinations = cell2mat(cellfun(@(x) x(:), grid_cells, 'UniformOutput', false));

    % 检查每行元素是否非递减
    valid_rows = all(diff(combinations, 1, 2) >= 0, 2);

    % 保留有效组合
    filtered_combinations = combinations(valid_rows, :);

    i = 1;
    probs = zeros(size(filtered_combinations)); 
    probs(:,1) = hygepdf(filtered_combinations(:,1),n,round(n.*p),round(v(1).*n));

    while i < (length(v))
        remy = round(n.*p - filtered_combinations(:,i));
        yhits = filtered_combinations(:,i+1) - filtered_combinations(:,i);
        npick = round(n.*(v(i+1) - v(i)));
        nrem = round(n.*(1-v(i)));    
        probs(:,i+1) = hygepdf(yhits,nrem,remy,npick);       
        i = i +1;
    end

    full = prod(probs,2); 
    falsepositive = 1 - sum(full);

    % 将结果存入fp数组
    fp(kx) = falsepositive;
end

补充说明 #

• 即便你避免使用clearvars，对于yv这类需要每次迭代重新初始化的cell数组，显式重置是必要的。
• 验证时可在每次循环结束后查看yv内容，确认其仅包含当前p对应的结果，无旧数据残留。

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