实现4维Halton序列

Question

可以找到Halton序列的伪代码here。我写了一个函数来做这个，但由于某种原因检查Matlab结果的第四维Halton序列我的数字不匹配，我不知道为什么。这是我的代码：

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