不使用任何外部库的C++中生成高斯分布的误差值的方法

/**
*@brief 产生高斯噪声
*@param E 噪声期望大小
*@param V 噪声方差大小
*@return 返回期望为E,方差为V的double型高斯值
*/
double gaussrand(double E,double V)
{
	static double V1, V2, S;
	static int phase = 0;
	double X;

	if (phase == 0) {
		do {
			double U1 = (double)rand() / RAND_MAX;
			double U2 = (double)rand() / RAND_MAX;

			V1 = 2 * U1 - 1;
			V2 = 2 * U2 - 1;
			S = V1 * V1 + V2 * V2;
		} while (S >= 1 || S == 0);

		X = V1 * sqrt(-2 * log(S) / S);
	}
	else
		X = V2 * sqrt(-2 * log(S) / S);

	phase = 1 - phase;
	X = X * V + E;
	return X;
}