Octave循环生成v1+v2+v3=1的变量组合时缺失结果，求修正

问题描述 #

有三个变量v1、v2、v3，满足v1 + v2 + v3 = 1，需要通过Octave循环以0.1为步长生成所有可能组合（预期共66种），原代码如下：

v1 = 1;
v2 = 0;
v3 = 0;
        while (v2 < 1)
        v3 = 0;
        v4 = v2;
          for v3 = 0:0.1:v4
          varsum = v1 + v2 + v3;
          data(1,n) = v1;
          data(2,n) = v2;
          data(3,n) = v3;
          data(4,n) = varsum;
          v2 -= 0.1;
          n+=1
          endfor
        v1 -= 0.1;
        v2 = 1 - v1;
        endwhile

实际运行仅生成65种组合，缺失0.9 0 0.1（对应预期格式的9 0 1）这一组合。

问题原因 #

原代码存在两个核心逻辑错误：

• 在for循环内部修改了循环控制变量v2，导致循环逻辑混乱；
• 当v1从1减到0.9后，v2被直接赋值为1-v1=0.1，跳过了v2=0的情况，因此遗漏了对应的v3=0.1组合。

修正代码 #

重新梳理循环逻辑，先遍历v1的所有可能值，再对每个v1遍历v2的合法取值，v3直接通过1-v1-v2计算，避免循环变量被意外修改：

% 初始化存储矩阵和计数变量
data = [];
n = 1;

% 遍历v1从1到0，步长-0.1
for v1 = 1:-0.1:0
    % 遍历v2从(1-v1)到0，步长-0.1
    for v2 = (1 - v1):-0.1:0
        v3 = 1 - v1 - v2;
        varsum = v1 + v2 + v3;
        % 将当前组合存入矩阵
        data(:,n) = [v1; v2; v3; varsum];
        n += 1;
    endfor
endfor

% 可选：将数值放大10倍，匹配预期格式的整数显示
data_scaled = data * 10;
disp(data_scaled);

代码说明 #

• 外层循环控制v1的取值，从1递减到0，步长为-0.1，覆盖所有可能的v1值；
• 内层循环针对每个v1，控制v2从1-v1递减到0，确保v1+v2不超过1；
• v3直接通过等式计算，无需额外循环，保证满足v1+v2+v3=1的约束；
• 最后可选择将数值放大10倍，得到预期的整数格式输出。

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