c++ - 实现 4 维 Halton 序列

Question

可以找到 Halton 序列的伪代码 here .我编写了一个执行此操作的函数，但出于某种原因检查 4 维 Halton 序列的 Matlab 结果，我的数字不匹配，我不确定为什么。这是我的代码:

double Halton_Seq(int index, double base){
    double f = 1, r;
    while(index > 0){
        f = f/base;
        r = r + f*(fmod(index,base));
        index = index/base;
    }
    return r;
}

这是我得到的前 10 个结果:

1
0.25
0.5
0.75
0.0625
0.3125
0.5625
0.8125
0.125
0.375

这是 MATLAB 获得的前 10 个结果:

  Columns 1 through 2

         0    0.5000

  Columns 3 through 4

    0.2500    0.7500

  Columns 5 through 6

    0.1250    0.6250

  Columns 7 through 8

    0.3750    0.8750

  Columns 9 through 10

    0.0625    0.5625

Answer 1

你忘记在第 2 行初始化 r。

r = 0;

double Halton_Seq(int index, int base){
  double f = 1, r = 0;
  while(index > 0){
    f = f/base;
    r = r + f* (index% base);
    index = index/base;
  }
  return r;
}
// Output for 10 (base 2)
0.000000
0.500000
0.250000
0.750000
0.125000
0.625000
0.375000
0.875000
0.062500
0.562500