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 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;

你可能感兴趣的:(C++)