在Octave中创建Halton序列

pkg load statistics

function sequence = halton_sequence(n, base)
    sequence = zeros(n, 1);
    for i = 1:n
        x = i;
        f = 1 / base;
        while x > 0
            sequence(i) = sequence(i) + f * mod(x, base);
            x = floor(x / base);
            f = f / base;
        end
    end
end

n = 100;  % 序列长度
base = 2; % Halton序列的基数
sequence = halton_sequence(n, base);