編集履歴

回答編集履歴

2018/10/10 09:45

投稿

スコア21956

test CHANGED Viewed

@@ -123,3 +123,147 @@
 正規分布を生成するアルゴリズムは他にもあるようです。
 いろいろなアルゴリズムの実装が提供されている [Github](https://github.com/miloyip/normaldist-benchmark) を見つけたので、見てみると面白いかもしれません。
+## 追記 グラフの生成方法について
+以下のコードでテキストファイルに出力して、Python で読み込んで seaborn というライブラリでプロットしました。
+Cだけでプロットする場合は gunplot とか使うことになるのでしょうけど、それでも結構面倒くさいです。Python だと1行できれいなグラフが出力できます。
+```cpp
+#include <fstream>
+#include <math.h>
+#include <stdio.h>
+#include <stdlib.h>
+/**
+   @brief [-1, 1] の一様分布に従う乱数を生成する。
+   @return [-1, 1] の一様分布に従う乱数
+ */
+double urand()
+{
+    return (double)rand() / RAND_MAX * 2. - 1.;
+}
+double nrand()
+{
+    // 1度に2個の正規分布に従う値 r1, r2 が得られるので
+    // 生成したときは r1、次に呼び出したときは r2 を返す。
+    static int sw = 0;  // switch 用のフラグ
+    static double r1, r2, s;
+    if (sw == 0) {
+        sw = 1;
+        // (x, y) \in [-1, 1]^2 を満たす点をランダムに選ぶ。
+        do {
+            r1 = urand();
+            r2 = urand();
+            s = r1 * r1 + r2 * r2;
+        } while (s > 1.0 || s == 0.0);
+        // 正規分布に従う点に変換する。
+        s = sqrt(-2. * log(s) / s);
+        return r1 * s;
+    }
+    else {
+        sw = 0;
+        return r2 * s;
+    }
+}
+int main()
+{
+    int seed = 1234567;
+    srand(seed);
+    std::ofstream ofs("values.csv");
+    for (int i = 0; i < 1000; i++)
+        ofs << nrand() << std::endl;
+}
+```
+ファイルを読み込んで、出力する。
+```python
+import seaborn as sns
+import matplotlib.pyplot as plt
+import numpy as np
+array = np.loadtxt('values.csv')
+print(array.shape)
+sns.set()
+ax = sns.distplot(array)
+plt.show()
+```