【文章內(nèi)容簡介】 )。 if (maxlen = sizeof(DWORD)) { memcpy(amp。buf[ni], amp。dwRes, sizeof(dwRes))。 ni += sizeof(dwRes)。 maxlen = sizeof(dwRes)。 } dwRes = GetTickCount()。 if (maxlen = sizeof(DWORD)) { memcpy(amp。buf[ni], amp。dwRes, sizeof(dwRes))。 ni += sizeof(dwRes)。 maxlen = sizeof(dwRes)。 } = sizeof(MEMORYSTATUS)。 GlobalMemoryStatus(amp。ms)。 if (maxlen = sizeof(MEMORYSTATUS)) { memcpy(amp。buf[ni], amp。ms, sizeof(MEMORYSTATUS))。 ni += sizeof(MEMORYSTATUS)。 maxlen = sizeof(MEMORYSTATUS)。 } return ni。 } void CCreRndNum::GetKey() { int sysInfoLen = get_SystemStatus_W32(300)。 CString sKey32。 //字符串型 int iKey32[32]。 //轉換后的整型 Confidential PagePage 10 of 50 10 1/7/2021 sKey32 = CMD5Checksum::GetMD5(buf,sysInfoLen)。//字符串 for(int i = 0。 i 32。 i++)//字符串轉化為整型 { iKey32[i] = CharToInt(sKey32[i])。 } for(i = 0。 i 64。 i++) //取出每一位 { KK1[i] = iKey32[i / 4] (i % 4) amp。 0x01。 KK2[i] = iKey32[16 + i / 4] (i % 4) amp。 0x01。 } } 函數(shù)名 : GetTime 函數(shù)功能 : 得到 64bit 的當前日期和時間的二進制 輸入 : 系統(tǒng)時間 輸出 : btime[64] 算法示意圖： 圖 3 分別得到產(chǎn)生大隨機數(shù)的輸入隨機數(shù) void CCreRndNum::GetTime() { int i。 FILETIME ft。 GetSystemTimeAsFileTime(amp。ft)。 for (i = 0。 i 32。 i++) { DTi[31 i] = (char)(( i) amp。 0x01)。 // 取得各通過移位異或把系統(tǒng)時間轉換成二進制數(shù)并放到 DTi 的高 32 和低 32 位中去 通過摸除 2得到 64位的二進制隨機數(shù) Vi DTi 和 Vi 都將作為主程序的輸入隨機數(shù) 輸入系統(tǒng)時間 Rand()產(chǎn)生的隨機數(shù) Confidential PagePage 11 of 50 11 1/7/2021 位的值 DTi[63 i] = (char)(( i) amp。 0x01)。 } } void CCreRndNum::GetVi() //得到 64 位的二進制隨機數(shù) {int i。 srand((unsigned)time(NULL))。 for (i = 0。 i 64。 i++) { Vi[i] = rand()%2。 } } 函數(shù)名 : Step_1 函數(shù)功能 : 產(chǎn)生隨機數(shù)第 1步 (見 圖 1) 輸入 : DTi(時間 ),K1,K2(密鑰 ) 輸出 : Result_1[64] 調(diào)用函數(shù) : MyDes::EDE() void CCreRndNum::Step_1() { GetTime()。 //得到時間二進制并保存在 DTi 中 GetKey()。 //產(chǎn)生密鑰并保存在 KK1,KK2 中 CDES MyDes。 (DTi, KK1, KK2, Result_1)。 //三重 DES 加密 (EDE) } 函數(shù)名 : Step_2 函數(shù)功能 : 產(chǎn)生隨機數(shù)第 2步 (見 圖 1) 輸入 : Result_1,Vi(第一次輸入的隨機數(shù) ) 輸出 : Result_2[64] void CCreRndNum::Step_2() { int i。 //之前已得到 V0 for (i = 0。 i 64。 i++) { Result_2[i] = Vi[i] ^ Result_1[i]。 } } 函數(shù)名 : Step_3 函數(shù)功能 : 產(chǎn)生隨機數(shù)第 3步 (見圖 1) 輸入 : Result_2,K1,K2(密鑰 ) 輸出 : Result_3[64] Confidential PagePage 12 of 50 12 1/7/2021 調(diào)用函數(shù) : MyDes::EDE() void CCreRndNum::Step_3() {CDES MyDes。 (Result_2, KK1, KK2, Result_3)。 } 函數(shù)名 : Step_4 函數(shù)功能 : 產(chǎn)生隨機數(shù)第 4步 異或 (見 圖 1) 輸入 : Result_1,Result_3 輸出 : Result_4[64] void CCreRndNum::Step_4() { int i。 for (i = 0。 i 64。 i++) { Result_4[i] = Result_1[i] ^ Result_3[i]。 } } 函數(shù)名 : Step_5 函數(shù)功能 : 產(chǎn)生隨機數(shù)第 5步 (見圖 1) 輸入 : Result_4,K1,K2(密鑰 ) 輸出 : Result_5[64] 調(diào)用函數(shù) : MyDes::EDE void CCreRndNum::Step_5() {CDES MyDes。 (Result_4, KK1, KK2, Result_5)。} 函數(shù)名 : GetRandNumber 函數(shù)功能 : 產(chǎn)生隨機數(shù) ,保存在 dwRndNum[33]中 . 輸入 : 無 輸出 : dwRndNum[33] 備注 : 輸出 1025 位二進制 ,故用 33 個 DWORD 型的數(shù)存放隨機數(shù) . BOOL CCreRndNum::GetRandNumber(DWORD dwRndNum[]) { GetVi()。 int i,m。 for (m = 0。 m 17。 m++) //17 輪產(chǎn)生 1088(取 1025)位的 Ri { Step_1()。 Step_2()。 Step_3()。 Confidential PagePage 13 of 50 13 1/7/2021 for (i = 0。 i 64。 i++) { Ri[(m 6) + i] = Result_3[i]。 //保存到 1024bit 的 Ri 中 } Step_4()。 Step_5()。 for (i = 0。 i 64。 i++) //Vi+1 自動賦給 Vi { Vi[i] = Result_5[i]。 } } // end for UINT temp。 for ( i = 0。 i 1024。 i++) //取第 1~1024 位 { temp = (UINT)(i / 32)。 if(i % 32 == 0) dwRndNum[temp] = dwRndNum[temp] | (DWORD)Ri[i]。 else dwRndNum[temp] = dwRndNum[temp] 1 | (DWORD)Ri[i]。 } for(i=0。i32。i++) { while(dwRndNum[i]2147483648) {dwRndNum[i]+=dwRndNum[i]。} } dwRndNum[32] amp。= (DWORD)Ri[1024]。 //取第 1025 位 ,其余位為 0 return TRUE。 } 6 檢驗隨機數(shù) 隨機數(shù)檢驗方法有如下幾種： 參數(shù)檢驗，檢驗其分布 參數(shù)的觀察值與理論值的差異顯著性。 均勻性檢驗，又稱頻率檢驗，意在檢驗偽隨機數(shù)的經(jīng)驗頻率與理論頻率的差異是否顯著。 獨立性檢驗，即檢驗所產(chǎn)生的偽隨機數(shù)的獨立性和統(tǒng)計相關是否異常，包括相關關系檢驗和聯(lián)列表檢驗等。 組合規(guī)律檢測，按隨機數(shù)出現(xiàn)的先后次序，根據(jù)一定的規(guī)律組合，檢驗其組合的觀察值與理值是不否有顯著差異，包括距離檢驗和配套檢驗等。 游程檢驗，把隨機數(shù)序列按一定的規(guī)則進行分類，分為正負游程檢驗和升降Confidential PagePage 14 of 50 14 1/7/2021 游程檢驗等。 本程序采用了均勻性檢測，即頻率檢測。 由于隨機數(shù)是由 16 次循環(huán)得來。并且存放在 32 個數(shù)組中， 所以取其中高位數(shù)組的值就可以判斷是否滿足均勻性檢測。 算法示意圖： 圖 4 檢驗隨機數(shù)足夠隨機 Int bb=0。 bb++。 char s[40]。 for (int i =0。i33。i++) { sprintf(s,(%2d) %10lX\r\n,i,dwn[i])。 (s)。 } cc=cc/bb*(bb1)+dwn[31]/bb。 sprintf(s,第 %2d 次前位平均值為： \n%10lX\n,bb,cc)。 (s)。 7 系統(tǒng)測試 第一次 多次運行程序，如圖： 隨機數(shù)除于當前隨機數(shù)產(chǎn)生的次數(shù) 如果第一次產(chǎn)生隨機數(shù)，平均值為隨機數(shù)本身 平均值除于隨機數(shù)產(chǎn)生次數(shù)與次數(shù) 1 的積 兩者相加得到隨機數(shù)的平均值 得到生成器產(chǎn)生的隨機數(shù) Confidential PagePage 15 of 50 15 1/7/2021 圖 5 第一次測試結果 第二次多次運行程序，如圖： 圖 6 第 2 次測試結果 第三次多次運行程序，如圖： Confidential PagePage 16 of 50 16 1/7/2021 圖 7 第 3 次測試結果 結 論 因為主導程序對本身就是隨機的數(shù)經(jīng)過多次的三重 DES 和異或等運算，所以能保證隨機數(shù)的足夠隨機性，通過 16 次的循環(huán)也能得到一個組合的大隨機數(shù)，所 以能滿足設計的需求 界面是用 VC++實現(xiàn)與設計的一個隨機數(shù)發(fā)生器，該隨機數(shù)發(fā)生器能夠通過點擊產(chǎn)生隨機數(shù)按鈕而相映產(chǎn)生 1024 位的隨機數(shù)，并且該產(chǎn)生器還記錄每次產(chǎn)生的隨機數(shù)而求其平均值，能夠很直觀的讓我們檢測到隨機數(shù)是否隨機。 隨機數(shù)產(chǎn)生器界面簡單，隨機數(shù)足夠隨機并且保證隨機數(shù)范圍滿足 1024 位2進制數(shù)。由于時間關系，程序也有一定瑕疵。 對程序改進的方法：對檢測的方法增加其他的種類，例如 參數(shù)檢驗 ， 游程檢驗 等方法，畢竟，通過的檢測越多隨機數(shù)的可靠性就越高。 Confidential PagePage 17 of 50 17 1/7/2021 參考文獻 [1]王銳 .網(wǎng)絡最高安全技術指南 [M].北京 :機械工業(yè)出版社 ,1998。 [2] (美 )克魯格林斯基 (Kruglinski,.).Visual C++ 技術內(nèi)幕 [M].北京 :清華大學出版社 ，1999。 [3]云巔工作室 .Visual C 中文版全面剖析 [M].北京 ： 中國水利水電出版社 ， 2021。 [4] Bruce ，算法與 C源程序 [M]. 北京 ： 機械工業(yè)出版社 ， 2021。 [5]馮登國著 .密碼分析學 [M].北京 ： 清華大學出版 ， 2021。 [6]（ 美 ) Steven John 設計模式 [M]. 北京 ： 中國電力出版社 ， 2021。 [7] 林元烈，梁宗霞 .隨機數(shù)學引論 [M].北京 ： 清華大學出版社 ， 2021。 [8] . 貝達特， . 皮爾索 .隨機數(shù)據(jù)分析方法 [M].北京： 北京國防工業(yè)出版社 ， 1976。 Confidential PagePage 18 of 50 18 1/7/