C++产生正态分布

1.求均值

1. double sum = std::accumulate(std::begin(resultSet), std::end(resultSet), 0.0);  
2. double mean =  sum / resultSet.size(); //均值

 2.求方差

1. double accum  = 0.0;  
2. std::for_each (std::begin(resultSet), std::end(resultSet), [&](const double d) {  
3.     accum  += (d-mean)*(d-mean);  
4. });  
5.   
6. double stdev = sqrt(accum/(resultSet.size()-1)); //方差 

3.求正太分布值密度函数

double  LogNormal( double  x, double  miu, double  sigma)  //对数正态分布概率密度函数

{

         return  1.0/(x* sqrt (2*PI)*sigma) *  exp (-1*( log (x)-miu)*( log (x)-miu)/(2*sigma*sigma));

}




4.求正态分布随机数




1.     std::default_random_engine e; //引擎 
2.     std::normal_distribution<double> n(mean, stdev); //均值, 方差  
3.     std::vector<unsigned> vals(9);  
4.     for(std::size_t i=0; i != 200; ++i) {  
5.         unsigned v = std::lround(n(e)); //取整-最近的整数  
6.         if (v < vals.size())  
7.             ++vals[v];  
8.     }  
9.     for (std::size_t j=0; j != vals.size(); ++j)  
10.         std::cout << j << " : " << vals[j] << std::string(vals[j], '*') << std::endl;