请教随机数？

seawinds

<>请教如何编写一个产生正态分布的随机数的生成器呢？？？？？？</P>
<>本人有急用</P>

dormouse_buaa

<>经典的数值算法书&lt;&lt;numerical recipe in C&gt;&gt; 中给出的代码如下，我就不解释了。如果不明白请去图书馆查一查就行了，这么经典的书图书馆都会有的，而且还有中文版，不过比较的旧了，不知有没有新版本的。</P><>float gasdev(int *idum)
// Returns a normally distributed deviate with zero mean and unit variance,
// using ran1(idum) as the source of uniform deviates.
{
static int iset=0;
static float gset;
float fac,r,v1,v2;</P><> if  (iset == 0) {
  do {
   v1=2.0*ran1(idum)-1.0;
   v2=2.0*ran1(idum)-1.0;
   r=v1*v1+v2*v2;
  } while (r &gt;= 1.0);
  fac=sqrt(-2.0*log(r)/r);
  gset=v1*fac;
  iset=1;
  return v2*fac;
} else {
  iset=0;
  return gset;
}
}
</P>

dormouse_buaa

<>刚才忘记了还需要ran1().现在补上
#define M1 259200
#define IA1 7141
#define IC1 54773
#define RM1 (1.0/M1)
#define M2 134456
#define IA2 8121
#define IC2 28411
#define RM2 (1.0/M2)
#define M3 243000
#define IA3 4561
#define IC3 51349
float ran1(int *idum)
{
static long ix1,ix2,ix3;
static float r[98];
float temp;
static int iff=0;
int j;
void nrerror();</P><> if (*idum &lt; 0 || iff == 0) {
  iff=1;
  ix1=(IC1-(*idum)) % M1;
  ix1=(IA1*ix1+IC1) % M1;
  ix2=ix1 % M2;
  ix1=(IA1*ix1+IC1) % M1;
  ix3=ix1 % M3;
  for (j=1;j&lt;=97;j++) {
   ix1=(IA1*ix1+IC1) % M1;
   ix2=(IA2*ix2+IC2) % M2;
   r[j]=(ix1+ix2*RM2)*RM1;
  }
  *idum=1;
}
ix1=(IA1*ix1+IC1) % M1;
ix2=(IA2*ix2+IC2) % M2;
ix3=(IA3*ix3+IC3) % M3;
j=1 + ((97*ix3)/M3);
assert(j &gt; 97 || j &lt; 1) ;
temp=r[j];
r[j]=(ix1+ix2*RM2)*RM1;
return temp;
}
#undef M1
#undef IA1
#undef IC1
#undef RM1
#undef M2
#undef IA2
#undef IC2
#undef RM2
#undef M3
#undef IA3
#undef IC3</P>

vckid

<>b4楼上的算法</P><>用M函数中的normrnd(mu,sigma,m,n)</P>