如何产生高斯随机数

caoyi132UID：2350218

/* boxmuller.c           Implements the Polar form of the Box-Muller
                         Transformation

                      (c) Copyright 1994, Everett F. Carter Jr.
                          Permission is granted by the author to use
                          this software for any application provided this
                          copyright notice is preserved.

*/

#include <math.h>


extern float ranf();         /* ranf() is uniform in 0..1 */


float box_muller(float m, float s)        /* normal random variate generator */
{                                        /* mean m, standard deviation s */
        float x1, x2, w, y1;
        static float y2;
        static int use_last = 0;

        if (use_last)                        /* use value from previous call */
        {
                y1 = y2;
                use_last = 0;
        }
        else
        {
                do {
                        x1 = 2.0 * ranf() - 1.0;
                        x2 = 2.0 * ranf() - 1.0;
                        w = x1 * x1 + x2 * x2;
                } while ( w >= 1.0 );

                w = sqrt( (-2.0 * log( w ) ) / w );
                y1 = x1 * w;
                y2 = x2 * w;
                use_last = 1;
        }

        return( m + y1 * s );
}

搞不清楚这个函数中 static int use_last=0；干什么用呀？ 这样的话不是每次调用是都只执行 if 嘛？

DriverEntryUID：143979

那是 STATIC INT 。 每次进入这个函数的时候，上次的值会保存下来的。自己搜搜 STATIC.

caoyi132UID：2350218

2# DriverEntry
嗯 是的，搞清楚了。谢啦。

DriverEntryUID：143979

不客气