毕业设计（论文）简述大随机数生成程序的开发与测试

1、简述大随机数生成程序的开发与测试摘要大随机数已经在当今社会的各个领域中都频繁使用，特别是在加密技术中已经成了不可缺少的一部分，像rsa，md5中随机数成为加密技术的关键。本设计主要为第3代移动通信系统（3g）提供符合要求的随机数（1024位），首先取得系统时间和rand（）函数所产生的随机数作为最初的随机初值，经过三重des（两密钥通过md5算法得来）和异或的变换，保证其随机数的足够随机，然后通过16次的循环得到一个组合起来的1024位随机数，设计还提供一个检验随机数是否随机的平台，采用了均匀性检测，即频率检测的方法检测随机数的随机性，通过检测发现，所产生的随机数能够达到我们所期望的随机性。设

2、计还对常见的随机数的生成方法进行了检析，提供多种随机数的生成方法，并且也提供了多种随机数的检测方法供大家参考，希望对大家有所帮助。关键字：随机数；rsa；md5；加密技术；均匀性检测confidentialpagepage 2 of 24 210/1/2011big random number generator algorithm research and implementabstractthe big random number is used everywhere in modern society especially in the encryption technology. th

3、e random number is the key technology of the encryption.this design mainly provides the request random number (1024) for 3rd generation of mobile communication system. the way to provide the number is discussed in this article, and the randomness test is discussed too. there are many ways to finish

4、the task which are shown in this paper. we hope these techniques can be useful. key words: random number; rsa; md5; encryption technology; even line of examination目 录论文总页数： 19页1 引言11.1随机数的概念11.2课题背景11.3 国内外研究现状11.4 本课题研究的意义11.5 本课题的研究方法12常见随机数生成方法简析22.1 迭代取中法22.2 乘同余法22.3 混同于法22.4 反变换法32.4.1 平均分布 :3

5、2.4.2 指数分布 :42.4.3 正态分布随机变量的生成 :42.5 离散型随机变量43随机数的检验54 大随机数产生的机理64.1 流程图64.2 des算法简介75 算法实现86 检验随机数137 系统测试14结 论16参考文献17致 谢18声 明191 引言1.1随机数的概念在现今的计算机中所产生的随机数，都是伪随机数。即，可以通过一定手段和方法发现或破译其中的规律。真随机数，也有了一定的研究，比如：通过声音或原子衰变等所产生的随机数。伪随机数可以通过一定的数学算法，近似真随机数但仍然不是真随机数。1.2课题背景随机数已经在当今社会的各个领域中都频繁使用，特别是在加密技术中已经成了不

6、可缺少的一部分，甚至很多加密技术的保密程度就取决于随机数。像rsa，md5需求大量随机数的密码技术正需求一个好的随机数发生器的产生。如今很多随机数产生器已经存在，但那些都存在很多的不足，比如产生的随机数位数不够，不是足够随机等等问题，所以编制一个能够产生我们需要足够大的且足够随机的随机数的随机数产生器就变得很重要。1.3 国内外研究现状通过查阅质料和在网上了解，国外对随机数的研究领先于国人对随机数的研究，但是总体来说对随机数的研究都还不够深入与透彻，都还不能脱离伪随机数的阴影，但的确目前的技术支持与环境配置等方面都还制约着我们大多数只能在研究伪随机数的层面，我们只可能的尽量地做到无限接近真随机

7、数，而不能达到真正的随机。特别在随机数的检测这方面，虽然检测的方法很多，但是都不够完善，没有一个很公用很全面的检测方法诞生，所以在随机数的研究中还有很深的东西需要我们去挖掘。1.4 本课题研究的意义保证我们能够很快速的得到需要的随机数，而且随机数能够足够大足够随机，尽量能够实用在需要用到随机数的任何地方，特别是在科研领域，比如第3代移动通信系统（3g）中需要的1024随机数，就能满足它的要求，我们所要做的就是使产生的随机数尽量的靠近真随机数。1.5 本课题的研究方法工作任务：1.大致了解随机数产生器的发展过程和现阶段的大概情况，认识现阶段随机数产生器的产生方式和所用到的知识，结构体系是怎样的。

8、2.分析他们的优点和缺点，能够保留的优点就要尽量用到，如果有不足应该怎样改正，加上自己的理解和题目的要求做一个满意的随机数产生器。要求： 使用vc+平台，编写一个能产生1024位的随机数发生器，而且随机数还要是足够随机的，并且还要编制一个检验平台，能在该平台上检验该随机数是足够随机的。设计思路：采用vc+6.01使用vc+实现控件的开发与界面的设计，尽量使外观简单容易实用，输出结果方便易看2借鉴其他随机数产生器的产生方法，参阅aes，des中随机数的产生方法，借鉴出其中的精华，补上自己的构思与想法尽量使随机数不出现重复。2常见随机数生成方法简析2.1 迭代取中法这里在迭代取中法中介绍平方取中法

9、 , 其迭代式如下 : xn+1=(xn2/10s)(mod 102s) rn+1=xn+1/102s 其中， xn+1 是迭代算子，而 rn+1 则是每次需要产生的随机数。 第一个式子表示的是将 xn 平方后右移 s 位，并截右端的 2s 位。 而第二个式子则是将截尾后的数字再压缩 2s 倍，显然 :0=rn+1= sizeof(dword)memcpy(&bufni, &dwres, sizeof(dwres);ni += sizeof(dwres);maxlen -= sizeof(dwres);dwres = getcurrentthreadid();if (maxlen = size

10、of(dword)memcpy(&bufni, &dwres, sizeof(dwres);ni += sizeof(dwres);maxlen -= sizeof(dwres);dwres = gettickcount();if (maxlen = sizeof(dword)memcpy(&bufni, &dwres, sizeof(dwres);ni += sizeof(dwres);maxlen -= sizeof(dwres);ms.dwlength = sizeof(memorystatus);globalmemorystatus(&ms);if (maxlen = sizeof(m

11、emorystatus)memcpy(&bufni, &ms, sizeof(memorystatus);ni += sizeof(memorystatus);maxlen -= sizeof(memorystatus);return ni;void ccrerndnum:getkey() int sysinfolen = get_systemstatus_w32(300);cstring skey32; /字符串型int ikey3232; /转换后的整型 skey32 = cmd5checksum:getmd5(buf,sysinfolen);/字符串for(int i = 0; i 32

12、; i+)/字符串转化为整型ikey32i = chartoint(skey32i);for(i = 0; i (i % 4) & 0x01; kk2i = ikey3216 + i / 4 (i % 4) & 0x01; 函数名: gettime函数功能: 得到64bit的当前日期和时间的二进制输入: 系统时间输出: btime64算法示意图：通过移位异或把系统时间转换成二进制数并放到dti的高32和低32位中去通过摸除2得到64位的二进制随机数vidti和vi都将作为主程序的输入随机数输入系统时间rand()产生的随机数图3 分别得到产生大随机数的输入随机数void ccrerndnum:

13、gettime()int i;filetime ft;getsystemtimeasfiletime(&ft);for (i = 0; i i) & 0x01); / 取得各位的值 dti63 - i = (char)(ft.dwlowdatetime i) & 0x01);void ccrerndnum:getvi() /得到64位的二进制随机数int i; srand(unsigned)time(null); for (i = 0; i 64; i+)vii = rand()%2;函数名: step_1函数功能: 产生随机数第1步(见图1)输入: dti(时间),k1,k2(密钥)输出:

14、result_164调用函数: mydes:ede()void ccrerndnum:step_1()gettime(); /得到时间二进制并保存在dti中 getkey(); /产生密钥并保存在kk1,kk2中cdes mydes;mydes.ede (dti, kk1, kk2, result_1); /三重des加密(ede)函数名: step_2函数功能: 产生随机数第2步(见图1)输入: result_1,vi(第一次输入的随机数)输出: result_264void ccrerndnum:step_2()int i;/之前已得到v0for (i = 0; i 64; i+)resu

15、lt_2i = vii result_1i;函数名: step_3函数功能: 产生随机数第3步(见图1)输入: result_2,k1,k2(密钥)输出: result_364调用函数: mydes:ede()void ccrerndnum:step_3()cdes mydes;mydes.ede(result_2, kk1, kk2, result_3);函数名: step_4函数功能: 产生随机数第4步异或(见图1)输入: result_1,result_3输出: result_464void ccrerndnum:step_4() int i; for (i = 0; i 64; i+)

16、 result_4i = result_1i result_3i; 函数名: step_5函数功能: 产生随机数第5步(见图1)输入: result_4,k1,k2(密钥)输出: result_564调用函数: mydes:edevoid ccrerndnum:step_5()cdes mydes;mydes.ede (result_4, kk1, kk2, result_5);函数名: getrandnumber函数功能: 产生随机数,保存在dwrndnum33中.输入: 无输出: dwrndnum33备注: 输出1025位二进制,故用33个dword型的数存放随机数.bool ccrern

17、dnum:getrandnumber(dword dwrndnum) getvi(); int i,m;for (m = 0; m 17; m+) /17轮产生1088(取1025)位的ristep_1(); step_2(); step_3(); for (i = 0; i 64; i+) ri(m 6) + i = result_3i; /保存到1024bit的ri中 step_4(); step_5(); for (i = 0; i 64; i+) /vi+1自动赋给vi vii = result_5i; / end for uint temp;for ( i = 0; i 1024;

18、i+) /取第11024位temp = (uint)(i / 32);if(i % 32 = 0) dwrndnumtemp = dwrndnumtemp | (dword)rii; else dwrndnumtemp = dwrndnumtemp 1 | (dword)rii;for(i=0;i32;i+)while(dwrndnumi2147483648)dwrndnumi+=dwrndnumi; dwrndnum32 &= (dword)ri1024; /取第1025位,其余位为0return true; 6 检验随机数随机数检验方法有如下几种：参数检验，检验其分布参数的观察值与理论值的

19、差异显著性。均匀性检验，又称频率检验，意在检验伪随机数的经验频率与理论频率的差异是否显著。独立性检验，即检验所产生的伪随机数的独立性和统计相关是否异常，包括相关关系检验和联列表检验等。组合规律检测，按随机数出现的先后次序，根据一定的规律组合，检验其组合的观察值与理值是不否有显著差异，包括距离检验和配套检验等。游程检验，把随机数序列按一定的规则进行分类，分为正负游程检验和升降游程检验等。本程序采用了均匀性检测，即频率检测。由于随机数是由16次循环得来。并且存放在32个数组中，所以取其中高位数组的值就可以判断是否满足均匀性检测。算法示意图： 随机数除于当前随机数产生的次数如果第一次产生随机数，平均

20、值为随机数本身平均值除于随机数产生次数与次数-1的积两者相加得到随机数的平均值得到生成器产生的随机数 图4 检验随机数足够随机int bb=0;bb+; char s40;for (int i =0;i33;i+)sprintf(s,(%2d) %10lxrn,i,dwni);m_show.replacesel (s);cc=cc/bb*(bb-1)+dwn31/bb;sprintf(s,第%2d次前位平均值为：n%10lxn,bb,cc);m_show.replacesel (s);7 系统测试 第一次多次运行程序，如图：图5 第一次测试结果第二次多次运行程序，如图：图6 第2次测试结果第三次多次运行程序，如图：图7 第3次