王 麗 燕， 柳 揚(yáng)

( 大連大學(xué) 信息工程學(xué)院， 遼寧 大連 116622 )

?

基于污染二維混沌動(dòng)力系統(tǒng)的加密算法

王 麗 燕*， 柳 揚(yáng)

( 大連大學(xué) 信息工程學(xué)院， 遼寧 大連 116622 )

首先定義了污染動(dòng)力系統(tǒng)，將二維Henon動(dòng)力系統(tǒng)用二維Logistic動(dòng)力系統(tǒng)進(jìn)行污染，用這個(gè)污染的二維混沌動(dòng)力系統(tǒng)構(gòu)造序列密碼體系．這種算法可以產(chǎn)生兩列密鑰，從而有效地解決了輸出結(jié)果對密鑰低bit位變化敏感度較低的問題．計(jì)算機(jī)模擬實(shí)驗(yàn)和游程測試、相關(guān)性分析、靈敏度分析、平衡度檢驗(yàn)等安全實(shí)驗(yàn)分析結(jié)果表明，密文、明文和密鑰之間具有高度的非線性和敏感性，算法的密鑰空間巨大，可以有效防止統(tǒng)計(jì)攻擊、唯密文攻擊和窮舉攻擊．

污染動(dòng)力系統(tǒng)；二維Henon動(dòng)力系統(tǒng)；二維Logistic動(dòng)力系統(tǒng)；序列密碼

0 引 言

為增強(qiáng)混沌動(dòng)力系統(tǒng)的敏感性，本文借用污染分布的思想，給出污染動(dòng)力系統(tǒng)的定義，將二維Henon動(dòng)力系統(tǒng)用二維Logistic動(dòng)力系統(tǒng)進(jìn)行污染，討論污染后動(dòng)力系統(tǒng)的簡單性質(zhì)．基于該污染的二維混沌系統(tǒng)產(chǎn)生兩列密鑰流對明文進(jìn)行加密，通過游程測試、相關(guān)性分析、靈敏度分析、平衡度檢驗(yàn)等計(jì)算機(jī)模擬實(shí)驗(yàn)對安全性進(jìn)行驗(yàn)證．

1 混沌動(dòng)力系統(tǒng)

1.1 污染混沌動(dòng)力系統(tǒng)

定義1 設(shè)f1和f2是兩個(gè)動(dòng)力系統(tǒng)，稱f=αf1+(1-α)f2為污染動(dòng)力系統(tǒng)，其中α(0<α<1)稱為污染系數(shù)．

1.2 污染二維混沌動(dòng)力系統(tǒng)

二維Logistic映射動(dòng)力學(xué)方程為

xn+1=μλ1xn(1-xn)+γyn

yn+1=μλ2yn(1-yn)+γxn

(1)

其中λ1、λ2、γ、μ是參數(shù)．當(dāng)μ=4時(shí)，表1列出了系統(tǒng)會(huì)出現(xiàn)混沌現(xiàn)象的不同條件[12]．

表1 二維Logistic映射混沌條件

圖1分別給出了兩種參數(shù)取值條件下的分叉情況．

二維Henon映射的動(dòng)力學(xué)方程為

(2)

其中p、q為參數(shù)．圖2給出了p=1.4，q=0.3條件下的混沌狀態(tài)圖．

將二維Henon映射用二維Logistic映射進(jìn)行污染，得到污染動(dòng)力系統(tǒng)：

(3)

其中0<α<1．

(a) λ1, λ2=1-λ1, λ=0.3

(b)λ1=λ2,λ=1-λ1

圖1 Logistic映射分叉圖

Fig.1 Branch chart of Logistic mapping

圖2 Henon映射的吸引子

利用Matlab計(jì)算，得到Jacobi矩陣的特征值的絕對值為1.871 9>1，根據(jù)差分方程組計(jì)算Lyapunov指數(shù)定義[13]，可知式(3)的Lyapunov指數(shù)大于零，說明污染的二維動(dòng)力系統(tǒng)(3)為混沌動(dòng)力系統(tǒng)．系統(tǒng)輸出的x(n)和y(n)的值如圖3所示，容易看出輸出值服從均勻分布．

(a) x(n)

(b)y(n)

圖3x(n)和y(n)在其區(qū)間的分布

Fig.3 The distribution ofx(n) andy(n) in their domains

2 序列密碼的加密與解密算法

2.1 加密算法

加密過程如圖4所示．

圖4 加密過程

具體加密算法如下：

(1)通過ASCII碼把明文轉(zhuǎn)化為十六位二進(jìn)制序列{m1m2m3…mn}，其中mi(i=1,2,3,…,n)為0或者1．

(2)確定密鑰．給二維污染混沌動(dòng)力系統(tǒng)(3)中的參數(shù)μ、γ、λ1、λ2、p、q、α取值，選定迭代初始值x0、y0，迭代得到兩個(gè)混沌序列{x(i)}(i=1,2,3,…,n)和{y(j)}(j=1,2,3,…,n)．

(3)由離散化算子Tk(x(i))=[10kx(i)]mod 2和Tk(y(j))=[10ky(j)]mod 2，計(jì)算得到兩個(gè)密鑰序列{k(i)}和{k(j)}，其中k(i)=Tk(x(i))(i=1,2,3,…,n)，k(j)=Tk(y(j))(j=1,2,3,…,n)．

(4)比較兩個(gè)密鑰序列{k(i)}(i=1,2,3,…,n)和{k(j)}(j=1,2,3,…,n)中0的個(gè)數(shù)，個(gè)數(shù)多的取正序列，個(gè)數(shù)少的取逆序列，然后將這兩個(gè)序列異或，得到新的密鑰序列{k(l)}(l=1,2,3,…,n)．

(5)將密鑰序列{k(l)}(l=1,2,3,…,n)與明文序列{m1m2m3…mn}進(jìn)行異或運(yùn)算，得到密文二進(jìn)制序列{c1c2c3…cn}．

(6)由密文序列C=c1c2c3…cn的ASCII值得到最終的密文．

類似地可以給出解密算法．

2.2 算法仿真

“二維污染混沌動(dòng)力系統(tǒng)”明文的二進(jìn)制序列為

01001110100011000111111011110100011011

00011000010110011111010011011011011111

01110110110010001100010100101010100001

01001010011011011111001111101101111110

11011111

不同條件下得到的密文如下：

(1)若取μ=4，γ=0.52，λ1=0.7，λ2=0.3，p=1.4，q=0.3，α=0.02，x0=0.131 4，y0=0.112 3，j=4，得到的密文為

(2)若取α=0.02+10-11，其他條件與(1)相同，密文為

(3)若取y0=0.112 3+10-11，其他條件與(1)相同，密文為

(4)若取x0=0.131 4+10-11，其他條件與(1)相同，密文為

(5)若取x0=0.131 4+10-11,y0=0.112 3+10-11，其他條件與(1)相同，密文為

(6)若取j=7，其他條件與(1)相同，密文為

(7)若取x0=0,y0=0，其他條件與(1)相同，密文為

(8)若取γ=0.52+10-7，x0=0，y0=0，j=4，其他條件與(1)相同，密文為

圖5給出了以上8種條件下密文用0-1序列的圖形化表示, 顯然，密鑰的細(xì)微改變將會(huì)導(dǎo)致密文的顯著改變．

3 安全性分析

3.1 游程測試[13]

游程是指序列中由相同bit所構(gòu)成的不間斷的子序列．該測試可以判斷其是否為隨機(jī)序列．

具體測試方法如下：

步驟4 計(jì)算判斷標(biāo)準(zhǔn)P：

如果P<0.01，斷定測試的序列隨機(jī)性較差；反之，斷定序列具有較好的隨機(jī)性．

若選取μ=4，γ=0.52，λ1=0.7，λ2=0.3，p=1.4，q=0.3，α=0.02，x0=0.131 4，y0=0.112 3，j=4，n=160，計(jì)算得到x(n)序列和y(n) 序列的P都為1.99，遠(yuǎn)大于0.01．因此，可以認(rèn)為混沌序列是隨機(jī)序列．

3.2 相關(guān)性分析

本文選取長度n=10 000的0序列作為明文，以動(dòng)力系統(tǒng)(3)參數(shù)值μ=4，γ=0.52，λ1=0.7，λ2=0.3，p=1.4，q=0.3，α=0.02，x0=0.131 4，y0=0.112 3，j=4為例，對明文進(jìn)行加密，相關(guān)度情況如圖6所示．

顯然隨著n的增大，相關(guān)度逐漸趨近于0，說明密文與明文幾乎不相關(guān)．

圖6 明文與密文的相關(guān)度

3.3 靈敏度分析

如果明文表示為{m1m2m3…mn}，密文表示為{s1s2s3…sn}，其中mi和si只取0或1，i=1,2,3,…,n，稱

為明文與密文間的靈敏度[27]．

仍以長度n=10 000的0序列作為明文，動(dòng)力系統(tǒng)(3)參數(shù)值μ=4，γ=0.52，λ1=0.7，λ2=0.3，p=1.4，q=0.3，α=0.02，x0=0.131 4，y0=0.112 3，j=4為例，對明文進(jìn)行加密，密文與明文間的靈敏度情況如圖7所示．

圖7 靈敏度分析圖

圖7表明，相比明文，大致50%的密文序列將會(huì)改變．

3.4 密文的平衡度檢驗(yàn)

仍以長度n=10 000的0序列作為明文，動(dòng)力系統(tǒng)(3)參數(shù)值μ=4，γ=0.52，λ1=0.7，λ2=0.3，p=1.4，q=0.3，α=0.02，x0=0.131 4，y0=0.112 3，j=4為例，對明文進(jìn)行加密，密文序列中0-1平衡度如圖8所示．

圖8 平衡度分析圖

圖8的結(jié)果表明，序列位數(shù)越多，平衡度的值就越趨近于0，說明序列中1和0的個(gè)數(shù)幾乎相等．

3.5 密鑰空間分析

本文選取的密鑰是污染二維混沌動(dòng)力系統(tǒng)隨機(jī)產(chǎn)生的初值x0、y0和離散化算子以及污染系數(shù)α，假設(shè)計(jì)算的精度為10-5，采用本算法形成的密鑰空間至少為1020．實(shí)際上，動(dòng)力系統(tǒng)本身的參數(shù)μ、γ、λ1、λ2、p、q只要在可以形成混沌的范圍內(nèi)取值，都可以作為密鑰．而且計(jì)算機(jī)的計(jì)算精度遠(yuǎn)遠(yuǎn)超過10-5，這樣密鑰空間將大大超過1020．本算法足夠抵抗由于密鑰空間不大而形成的窮舉攻擊．

4 結(jié) 語

本文給出污染混沌動(dòng)力系統(tǒng)的概念，并用污染系數(shù)α將二維Henon動(dòng)力系統(tǒng)用二維Logistic動(dòng)力系統(tǒng)進(jìn)行污染，并進(jìn)一步用污染后的多參數(shù)動(dòng)力系統(tǒng)構(gòu)造序列密碼的加密算法．所進(jìn)行的各項(xiàng)性能分析，如隨機(jī)性分析、相關(guān)性分析、靈敏度分析、0-1平衡度檢驗(yàn)等都表明污染的混沌映射具有良好統(tǒng)計(jì)特性，而且密鑰空間巨大，可以有效防止統(tǒng)計(jì)攻擊、唯密文攻擊和窮舉攻擊．

[1] Matthews R A J. On the derivation of a ″chaotic″ encryption algorithm [J]. Cryptologia, 1984, 8(1):29-41.

[2] Habutsu T, Nishio Y, Sasase I,etal. A secret key cryptosystem by iterating a chaotic map [C] // Davies D W, ed. Advances in Cryptology - EUROCRYPT ′91, LNCS 547. Berlin:Springer-Verlag, 1991:127-140.

[3] Biham E. Cryptanalysis of the chaotic map cryptosystem suggested at EUROCRYPT′91 [C]// EUROCRYPT′91 Proceedings of the 10th Annual International Conference on Theory and Application of Cryptographic Techniques. Berlin:Springer-Verlag, 1991:532-534.

[4] 周 紅,羅 杰,凌燮亭. 混沌非線性反饋密碼序列的理論設(shè)計(jì)和有限精度實(shí)現(xiàn)[J]. 電子學(xué)報(bào), 1997, 25(10):57-60.

ZHOU Hong, LUO Jie, LING Xie-ting. Generating nonlinear feedback stream ciphers via chaotic systems [J]. Acta Electronica Sinica, 1997, 25(10):57-60. (in Chinese)

[5] 周 紅,俞 軍,凌燮亭. 混沌前饋型流密碼的設(shè)計(jì)[J]. 電子學(xué)報(bào), 1998, 26(1):98-101.

ZHOU Hong, YU Jun, LING Xie-ting. Design of chaotic feed forward stream cipher [J]. Acta Electronica Sinica, 1998, 26(1):98-101. (in Chinese)

[6] 桑 濤,王汝笠,嚴(yán)義塤. 一類新型混沌反饋密碼序列的理論設(shè)計(jì)[J]. 電子學(xué)報(bào), 1999, 27(7):47-50.

SANG Tao, WANG Ru-li, YAN Yi-xun. The theoretical design for a class of new chaotic feedback stream ciphers [J]. Acta Electronica Sinica, 1999, 27(7):47-50. (in Chinese)

[7] 孫 楓,秦紅磊,徐耀群,等. 基于混沌的分組密碼置換網(wǎng)絡(luò)的設(shè)計(jì)[J]. 中國工程科學(xué), 2000, 2(9):47-49.

SUN Feng, QIN Hong-lei, XU Yao-qun,etal. Design of block cipher substitution network on chaos [J]. Engineering Science, 2000, 2(9):47-49. (in Chinese)

[8] 王相生,甘駿人. 一種基于混沌的序列密碼生成方法[J]. 計(jì)算機(jī)學(xué)報(bào), 2002, 25(4):351-356.

WANG Xiang-sheng, GAN Jun-ren. A chaotic sequence encryption method [J]. Chinese Journal of Computers, 2002, 25(4):351-356. (in Chinese)

[9] 翁貽方,鞠 磊. 基于混沌的序列密碼加密算法[J]. 計(jì)算機(jī)工程, 2002, 28(11):79-80, 83.

WENG Yi-fang, JU Lei. Chaotic stream cipher encryption algorithms [J]. Computer Engineering, 2002, 28(11):79-80, 83. (in Chinese)

[10] 李紅達(dá),馮登國. 基于復(fù)合離散混沌動(dòng)力系統(tǒng)的序列密碼算法[J]. 軟件學(xué)報(bào), 2003, 14(5):991-998.

LI Hong-da, FENG Deng-guo. Stream cipher algorithms based on composite nonlinear discrete chaotic dynamical systems [J]. Journal of Software, 2003, 14(5):991-998. (in Chinese)

[11] 李紅達(dá),馮登國. 復(fù)合離散混沌動(dòng)力系統(tǒng)與序列密碼體系[J]. 電子學(xué)報(bào), 2003, 31(8):1209-1212.

LI Hong-da, FENG Deng-guo. Composite nonlinear descrete chaotic dynamical systems and stream cipher systems [J]. Acta Electronica Sinica, 2003, 31(8):1209-1212. (in Chinese)

[12] Gotz M, Kelber K, Schwarz W. Discrete-time chaotic encryption systems. I. Statistical design approach[J]. IEEE Transactions on Circuits and Systems. I. Fundamental Theory and Applications, 1997, 44(10):963-970.

[13] Kocarev L. Chaos-based cryptography:A brief overview [J]. IEEE Circuits and Systems Magazine, 2002, 1(3):6-21.

[14] XIANG Tao, Wong Kwor-kwo, LIAO Xiao-feng. A novel symmetrical cryptosystem based on discretized two-dimensional chaotic map [J]. Physics Letters A, 2007, 364(3-4):252-258.

[15] Alvarez G, Montoya F, Romera M,etal. Cryptanalysis of a chaotic encryption system [J]. Physics Letters A, 2000, 276(1-4):191-196.

[16] 張 斌,金晨輝. 對迭代型混沌密碼的逆推壓縮攻擊[J]. 電子學(xué)報(bào), 2010, 38(1):129-134,140.

ZHANG Bin, JIN Chen-hui. Inversion and compression attacks to iterative chaotic ciphers [J]. Acta Electronica Sinica, 2010, 38(1):129-134,140. (in Chinese)

[17] 汪海明,李 明,金晨輝. 對XW混沌密碼算法的分割攻擊[J]. 計(jì)算機(jī)應(yīng)用研究, 2010, 27(7):2625-2628.

WANG Hai-ming, LI Ming, JIN Chen-hui. Divide-and-conquer attack on XW chaotic cipher [J]. Application Research of Computers, 2010, 27(7):2625-2628. (in Chinese)

[18] 尹汝明,袁 堅(jiān),山秀明,等. 混沌密碼系統(tǒng)弱密鑰隨機(jī)性分析[J]. 中國科學(xué):信息科學(xué), 2011, 41(7):777-788.

YIN Ru-ming, YUAN Jian, SHAN Xiu-ming,etal. Weak key analysis for chaotic cipher based on randomness properties [J]. Science in China:Information Sciences, 2011, 41(7):777-788. (in Chinese)

[19] 金晨輝. 一個(gè)基于混沌的分組密碼算法的分析[J]. 中國工程科學(xué), 2001, 3(6):75-80.

JIN Chen-hui. Analysis of a block cipher based on chaos [J]. Engineering Science, 2001, 3(6):75-80. (in Chinese)

[20] 金晨輝,高海英. 對兩個(gè)基于混沌的序列密碼算法的分析[J]. 電子學(xué)報(bào), 2004, 32(7):1066-1070.

JIN Chen-hui, GAO Hai-ying. Analysis of two stream ciphers based on chaos [J]. Acta Electronica Sinica, 2004, 32(7):1066-1070. (in Chinese)

[21] 金晨輝,楊 陽,祁傳達(dá). 對混沌序列密碼的相關(guān)密鑰攻擊[J]. 電子與信息學(xué)報(bào), 2006, 28(3):410-414.

JIN Chen-hui, YANG Yang, QI Chuan-da. A related-key attack on chaotic stream ciphers [J]. Journal of Electronics & Information Technology, 2006, 28(3):410-414. (in Chinese)

[22] 王麗燕,李永華,賈思齊,等. 一種基于復(fù)合混沌動(dòng)力系統(tǒng)的序列密碼算法[J]. 大連理工大學(xué)學(xué)報(bào), 2012, 52(5):730-735.

WANG Li-yan, LI Yong-hua, JIA Si-qi,etal. A stream cipher algorithm based on composite chaotic dynamical systems [J]. Journal of Dalian University of Technology, 2012, 52(5):730-735. (in Chinese)

[23] 孫小雁,張煥國,張茂勝,等. Logistic混沌擾動(dòng)三角形密碼體制[J]. 計(jì)算機(jī)應(yīng)用與軟件, 2014, 31(9):268-271.

SUN Xiao-yan, ZHANG Huan-guo, ZHANG Mao-sheng,etal. Triangular cryptosystem with Logistic chaos disturbance [J]. Computer Applications and Software, 2014, 31(9):268-271. (in Chinese)

[24] 王麗燕,許佳佳,李海燕. 基于兩個(gè)離散混沌動(dòng)力系統(tǒng)的序列密碼算法[J]. 大連理工大學(xué)學(xué)報(bào), 2014, 54(5):581-588.

WANG Li-yan, XU Jia-jia, LI Hai-yan. A stream cipher algorithm based on two discrete chaotic dynamical systems [J]. Journal of Dalian University of Technology, 2014, 54(5):581-588. (in Chinese)

[25] 張 順,高鐵杠. 基于類DNA編碼分組與替換的加密方案[J]. 電子與信息學(xué)報(bào), 2015, 37(1):150-157.

ZHANG Shun, GAO Tie-gang. Encryption based on DNA coding, codon grouping and substitution [J]. Journal of Electronics & Information Technology, 2015, 37(1):150-157. (in Chinese)

[26] Merah L, Ali-Pacha A, Naima H-S. Real-time cryptosystem based on synchronized chaotic system [J]. Nonlinear Dynamics, 2015, 82(1-2):877-890.

[27] 廖曉峰,肖 迪,陳 勇,等. 混沌密碼學(xué)原理及其應(yīng)用[M]. 北京:科學(xué)出版社, 2009.

LIAO Xiao-feng, XIAO Di, CHEN Yong,etal. Theory and Applications of Chaotic Cryptography [M]. Beijing:Science Press, 2009. (in Chinese)

Encryption algorithm based on contaminated two-dimensional chaotic dynamic system

WANG Li-yan*, LIU Yang

( College of Information Engineering, Dalian University, Dalian 116622, China )

An algorithm to construct a stream cipher system is presented by defining a contaminated dynamic system, in which the two-dimensional Henon dynamic system is contaminated by the two-dimensional Logistic dynamic system. By generating two columns of keys, this algorithm can effectively solve the problem that the output is not very sensitive to the change of the low bit of the input. The results of safety tests, such as computer simulation, runs test, correlation analysis, sensitivity analysis and balance test, etc. show the highly nonlinearity and sensitivity among ciphertext, plaintext and the key. Also the large key space of this algorithm can effectively prevent the statistical attacks, ciphertext only attacks and exhaustive attacks.

contaminated dynamic system; two-dimensional Henon dynamic system; two-dimensional Logistic dynamic system; stream cipher

2016-01-19；

2016-09-02．

國家自然科學(xué)基金資助項(xiàng)目(71072161)．

王麗燕*(1963-)，女，博士，教授，E-mail:wly1963@163.com；柳 揚(yáng)(1979-)，女，博士生，E-mail:lykx2001@163.com.

1000-8608(2016)06-0650-07

TN918

A

10.7511/dllgxb201606014