python,C#等語言計算正態分布積分的非常的方便。
提到C++,很多人就會覺得這個語言非常的笨, 如果不是為了追求效率,很少人會用C++寫程序。 其實最近發現,C++有很多封裝好的庫, 比如矩陣計算有eigen庫等。
今天想用C++計算一個標准正態分布的積分。 查了很多東西,包括boost等。
最后發現C++11可以直接算erf(x),稍微變換一下就可以計算高斯函數的積分了。erf的介紹在 https://en.wikipedia.org/wiki/Error_function
具體代碼如下:
const double M_SQRT1_2 = sqrt(0.5);
double normalCFD(double value)
{
return 0.5 * erfc(-value * M_SQRT1_2);
}
主要參考資料: https://stackoverflow.com/questions/2328258/cumulative-normal-distribution-function-in-c-c