生成符合正太分布随机数

1. const int N = 100;  
2.     //const double MAX = 50.0;  
3.     //const double MIN = 40.0;  
4.     //const double MIU = 40.0;  
5.     //const double SIGMA = 1.0;  
6.     static Random aa = new Random((int)(DateTime.Now.Ticks / 10000));  
7.   
8.     public double AverageRandom(double min, double max)//产生(min,max)之间均匀分布的随机数  
9.     {  
10.        int MINnteger = (int)(min * 10000);  
11.        int MAXnteger = (int)(max * 10000);  
12.        int resultInteger = aa.Next(MINnteger, MAXnteger);  
13.        return resultInteger / 10000.0;  
14.      }  
15.   
16.     public double Normal(double x, double miu, double sigma) //正态分布概率密度函数  
17.     {  
18.         return 1.0 / (x * Math.Sqrt(2 * Math.PI) * sigma) * Math.Exp(-1 * (Math.Log(x) - miu) * (Math.Log(x) - miu) / (2 * sigma * sigma));  
19.     }  
20.   
21.     public double Random_Normal(double miu, double sigma, double min, double max)//产生正态分布随机数  
22.     {  
23.         double x;  
24.         double dScope;  
25.         double y;  
26.         do  
27.         {  
28.            x = AverageRandom(min, max);  
29.            y = Normal(x, miu, sigma);  
30.            dScope = AverageRandom(0, Normal(miu, miu, sigma));  
31.          } while (dScope > y);  
32.          return x;  
33.      }  

[csharp]  view plain copy

1. protected void btnGenerate_Click(object sender, EventArgs e)  
2.    {  
3.        //定义这个正态分布随机数的范围、均值和方差  
4.        double min = double.Parse(txtBegin.Text);  
5.        double max = double.Parse(txtEnd.Text);  
6.        double miu = double.Parse(txtMIU.Text);  
7.        double sigma = double.Parse(txtSIGMA.Text);  
8.        //随机差生四个随机数  
9.        TextBox1.Text = Random_Normal(miu,sigma,min,max).ToString();  
10.        TextBox2.Text = Random_Normal(miu,sigma,min,max).ToString();  
11.        TextBox3.Text = Random_Normal(miu,sigma,min,max).ToString();  
12.        TextBox4.Text = Random_Normal(miu,sigma,min,max).ToString();  
13.    }