高斯随机数产生原理及代码_笔记.docxVIP

一:随机数的产生 C++中不提供random函数,但是提供了rand函数,产生0~RAND_MAX之间的整数,但严格意义上来讲生成的只是伪随机数(pseudo-random integral number).生成随机数时需要指定一个种子，如果在程序内循环，那么下一次生成随机数时调用上一次的结果作为种子。但如果分两次执行程序，那么由于种子相同，生成的“随机数”也是相同的。rand()函数不接受参数，默认以1为种子（即起始值）。若随机数生成器总是以相同的种子开始，所以形成的伪随机数列也相同，失去了随机意义。（但这样便于程序调试）. C++中另一函数srand（），可以指定不同的数（无符号整数变元）为种子。但是如果种子相同，伪随机数列也相同。一个办法是让用户输入种子，但是仍然不理想。 比较理想的是用变化的数，比如时间来作为随机数生成器的种子。 time的值每时每刻都不同。所以种子不同，所以，产生的随机数也不同。 rand函数产生随机数的方法: 1如果要产生0~10的10个整数，可以表达为： 　　int N = rand() % 11; 这样，N的值就是一个0~10的随机数，如果要产生1~10，则是这样： 　　int N = 1 + rand() % 11; 总结来说，可以表示为： 　　a + rand()%(b-a+1) 其中的a是起始值，(b-a+1)是整数的范围。 若要0~1的小数，则可以先取得0~10的整数，然后均除以10即可得到随机到十分位的10个随机小数，若要得到随机到百分位的随机小数，则需要先得到0~100的10个整数，然后均除以100，其它情况依 此类推。当要求的精度较高的时候，使用RAND_MAX作为分母。如果要求左闭右开的话，分母设置为RAND_MAX+1即可。精度要求高时形式如下： x=a+((rand()%RAND_MAX)/(double)RAND_MAX)*(b-a); /*x belong [a,b] */ 通常rand()产生的随机数在每次运行的时候都是与上一次相同的，这是有意这样设计的，是为了便于程序的调试。若要产生每次不同的随机数，可以使用srand(seed)函数进行随机化，随着seed的不同，就能够产生不同的随机数。还可以包含time.h头文件，然后使用srand(time(0))来使用当前时间使随机数发生器随机化，这样就可以保证每两次运行时可以得到不同的随机数序列(只要两次运行的间隔超过1秒)。 后一种方法性能更优！ 在工程应用时，我们一般将系统当前时间(Unix时间)作为种子，这样生成的随机数更接近于实际意义上的随机数。给一下例程如下： #include iostream #include time.h using namespace std; void main(){ int i=0; int count[6]={0}; srand(time(0)); /* 根据系统时间设置随机数种子,只需要语句一次 */ for(int k=0;k100000;k++){ i=rand(); i=i%5+1; /* 产生1~5之间的整数,均匀产生 */ switch(i){ case 1:count[1]++; break; case 2:count[2]++; break; case 3:count[3]++; break; case 4:count[4]++; break; case 5:count[5]++; break; } } for (i=1;i6;i++) { coutcount[i]endl; } } 输出结果如下所示： 可以看出，随机产生的数，是符合均匀分布的。 Pay attention: 如果想在一个程序中生成随机数序列，需要至多在生成随机数之前设置一次随机种子。同时，如果程序在执行很短的时间内，连续执行两次随机数产生器(函数)，即使使用srand函数把系统时间做为种子(time(NULL)所取得的时间只能精确到秒)，会得到相同的时间，因此产生的第二组随机数将会和第一组随机数相同，或者产生的后一组随机数比较有规律。这是因为计算机计算能力强，此时系统的时间并没有发生变化。另一种情况是系统的时间呈线性变化的，使