📄 gaussrand.luben.html
字号:
<html><!-- Mirrored from c-faq.com/lib/gaussrand.luben.html by HTTrack Website Copier/3.x [XR&CO'2008], Sat, 14 Mar 2009 08:03:00 GMT --><head><title></title></head><body>From: Luben<br>Subject: C-FAQ, q. 13.20<br>Message-ID: <Pine.GSO.3.95.980510025115.27152A-100000@chimp><br>Date: Sun, 10 May 1998 03:12:21 -0400<p>The following version uses explictly the inverse of the density functionof the Normal distribution, as well as you can specify mean and standarddeviation. It's no longer than the one at q. 13.20 and simpler.Check it out.<p><pre>/* Generate random numbers with mean mean, standard deviation std_dev, *//* with a Normal (Gaussian) distribution. *//* mean is any real number, std_dev > 0 (or machine epsilon). *//* The formula used is the inverse of the density function of the *//* Normal distribution: *//* x = mean +/- std_dev * sqrt((-2.0) * log(y)), 0 < y <= 1 *//* The client should call srand(int) to initialize the seed, *//* before any calls to gauss_rand(). *//* L.T., May 3, 1998, University of Toronto */#include <math.h>#include <float.h>#include <errno.h>#include <stdio.h>#include <stdlib.h>double gauss_rand(double mean, double std_dev){ double y; unsigned bin; errno = 0; /* std_dev must be greater than 0.0 (or machine epsilon) */ if ( ! (std_dev > DBL_EPSILON)) { errno = EDOM; /* domain error */ perror("gaussrand: std_dev too small or zero"); return mean; } y = (double) rand() / (RAND_MAX + 1.0); /* 0.0 <= y < 1.0 */ bin = (y < 0.5) ? 0 : 1; y = fabs(y - 1.0); /* 0.0 < y <= 1.0 */ y = std_dev * sqrt((-2.0) * log(y)); return bin ? (mean + y) : (mean - y);}</pre></body><!-- Mirrored from c-faq.com/lib/gaussrand.luben.html by HTTrack Website Copier/3.x [XR&CO'2008], Sat, 14 Mar 2009 08:03:00 GMT --></html>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -