#include<iostream>
#include <iostream>
#include<cmath>
#include <random>
#include<ctime>
using namespace std;
double normalf(double E, double sigema2)
{
double y;
std::random_device rd{};
std::mt19937 gen(rd());
std::normal_distribution<> d{ E,sqrt(sigema2) }; //均值, 标准差
//std::normal_distribution<> d{ E,sigema2};
y = d(gen);
return y;
}
//const double M_SQRT1_2 = sqrt(0.5);
//double normalCFD(double value)
//{
// return 0.5 * erfc(-value * M_SQRT1_2); // fai(x)
//}
const double M_SQRT1_2 = sqrt(0.5);
double normpdpdf(double value)
{
return 0.5 * erfc(-value * M_SQRT1_2);
}
double qfunc(double value)
{
return 0.5 * erfc(value * M_SQRT1_2); // Q(x) only "-
#include <iostream>
#include<cmath>
#include <random>
#include<ctime>
using namespace std;
double normalf(double E, double sigema2)
{
double y;
std::random_device rd{};
std::mt19937 gen(rd());
std::normal_distribution<> d{ E,sqrt(sigema2) }; //均值, 标准差
//std::normal_distribution<> d{ E,sigema2};
y = d(gen);
return y;
}
//const double M_SQRT1_2 = sqrt(0.5);
//double normalCFD(double value)
//{
// return 0.5 * erfc(-value * M_SQRT1_2); // fai(x)
//}
const double M_SQRT1_2 = sqrt(0.5);
double normpdpdf(double value)
{
return 0.5 * erfc(-value * M_SQRT1_2);
}
double qfunc(double value)
{
return 0.5 * erfc(value * M_SQRT1_2); // Q(x) only "-