1 次元配列を作成し、乱数ジェネレーター (70 の手段と 10 の標準偏差で乱数を生成するガウス ジェネレーター) を使用して、0 から 100 までの範囲の少なくとも 100 個の数値を配列に入力しようとしています。
C++でこれを行うにはどうすればよいですか?
1 次元配列を作成し、乱数ジェネレーター (70 の手段と 10 の標準偏差で乱数を生成するガウス ジェネレーター) を使用して、0 から 100 までの範囲の少なくとも 100 個の数値を配列に入力しようとしています。
C++でこれを行うにはどうすればよいですか?
C++11 では、ランダム ヘッダーとstd::normal_distributionを使用すると、これは比較的簡単です( live example )。
#include <iostream>
#include <iomanip>
#include <string>
#include <map>
#include <random>
int main()
{
std::random_device rd;
std::mt19937 e2(rd());
std::normal_distribution<> dist(70, 10);
std::map<int, int> hist;
for (int n = 0; n < 100000; ++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';
}
}
C++ 11がオプションでない場合は、boostもライブラリを提供します ( live example ):
#include <iostream>
#include <iomanip>
#include <string>
#include <map>
#include <random>
#include <boost/random.hpp>
#include <boost/random/normal_distribution.hpp>
int main()
{
boost::mt19937 *rng = new boost::mt19937();
rng->seed(time(NULL));
boost::normal_distribution<> distribution(70, 10);
boost::variate_generator< boost::mt19937, boost::normal_distribution<> > dist(*rng, distribution);
std::map<int, int> hist;
for (int n = 0; n < 100000; ++n) {
++hist[std::round(dist())];
}
for (auto p : hist) {
std::cout << std::fixed << std::setprecision(1) << std::setw(2)
<< p.first << ' ' << std::string(p.second/200, '*') << '\n';
}
}
何らかの理由でこれらのオプションのどちらも使用できない場合は、独自のBox-Muller transformをロールすることができます。リンクで提供されているコードは合理的に見えます。
Box Muller ディストリビューションを使用します (ここから):
double rand_normal(double mean, double stddev)
{//Box muller method
static double n2 = 0.0;
static int n2_cached = 0;
if (!n2_cached)
{
double x, y, r;
do
{
x = 2.0*rand()/RAND_MAX - 1;
y = 2.0*rand()/RAND_MAX - 1;
r = x*x + y*y;
}
while (r == 0.0 || r > 1.0);
{
double d = sqrt(-2.0*log(r)/r);
double n1 = x*d;
n2 = y*d;
double result = n1*stddev + mean;
n2_cached = 1;
return result;
}
}
else
{
n2_cached = 0;
return n2*stddev + mean;
}
}
詳細については、wolframe math worldを参照してください。
<random>
C++11 では、ヘッダーによって提供される機能を使用します。ランダム エンジン (必要に応じてstd::default_random_engine
またはstd::mt19937
で初期化) と、パラメータで初期化されたオブジェクトを作成します。次に、それらを一緒に使用して数値を生成できます。ここで完全な例を見つけることができます。std::random_device
std::normal_distribution
代わりに、以前のバージョンの C++ では、 [0, MAX_RAND] の範囲で単純な整数分布を生成するだけの「古典的な」C LCG ( srand
/ ) しかありません。Box-Muller transformrand
を使用してガウス乱数を生成できます。(ここに示すように、C++11 GNU GCC libstdc++はMarsaglia polar メソッドを使用することに注意してください。)std::normal_distribution
と#include <random>
std::default_random_engine de(time(0)); //seed
std::normal_distribution<int> nd(70, 10); //mean followed by stdiv
int rarrary [101]; // [0, 100]
for(int i = 0; i < 101; ++i){
rarray[i] = nd(de); //Generate numbers;
}