MATLABでは、関数normrnd(mu、sigma)を使用して通常の乱数を生成します。C ++およびそのライブラリでの同等の関数は何ですか?
1214 次
3 に答える
2
randomライブラリを参照してください。例を次に示します。
#include <iostream>
#include <iomanip>
#include <string>
#include <map>
#include <random>
int main()
{
std::random_device rd;
//
// Engines
//
std::mt19937 e2(rd());
//std::knuth_b e2(rd());
//std::default_random_engine e2(rd()) ;
//
// Distribtuions
//
std::normal_distribution<> dist(2, 2);
//std::student_t_distribution<> dist(5);
//std::poisson_distribution<> dist(2);
//std::extreme_value_distribution<> dist(0,2);
std::map<int, int> hist;
for (int n = 0; n < 10000; ++n) {
++hist[std::round(dist(e2))];
}
for (auto p : hist) {
std::cout << std::fixed << std::setprecision(1) << std::setw(2)
<< p.first << ' ' << std::string(p.second/200, '*') << '\n';
}
}
これは、 normal_distributionの特定のページです。
于 2013-03-19T19:51:33.923 に答える
0
自分で書きたいだけなら…
double normrnd(double mean, double stdDev) {
double u, v, s;
do {
u = ((double)rand()/(double)RAND_MAX) * 2.0 - 1.0;
v = ((double)rand()/(double)RAND_MAX) * 2.0 - 1.0;
s = u * u + v * v;
} while (s >= 1 || s == 0);
double mul = sqrt(-2.0 * log(s) / s);
return mean + stdDev * u * mul;
}
http://en.wikipedia.org/wiki/Marsaglia_polar_method#Implementationのアルゴリズムから簡略化
于 2013-03-20T11:34:46.167 に答える
0
正規分布で Boost.Random を使用できます。
于 2013-03-19T19:50:40.197 に答える