Microsoft Word - 薛?薇.doc

Size: px

Start display at page:

Download "Microsoft Word - 薛?薇.doc"

柿街屈
5 years ago
Views:

1 第 30 卷第 5 期 2015 年 10 月天津科技大学学报 Journal of Tianjin University of Science & Technology Vol. 30 No.5 Oct DOI: /j.issn 一个分数阶超混沌系统及其在图像加密中的应用薛薇, 肖慧, 徐进康, 贾红艳 ( 天津科技大学电子信息与自动化学院, 天津 ) 摘要 : 对一个同量分数阶超混沌 Qi 系统进行仿真分析, 绘制了该分数阶系统的相轨迹图 Poincaré 截面映射图 Lyapunov 指数谱图和分岔图, 仿真结果验证了该分数阶系统的混沌特性. 在此基础上, 将该分数阶超混沌 Qi 系统应用于数字图像加密, 并对密文的直方图相邻像素相关性密钥敏感性等进行分析. 结果表明 : 基于分数阶超混沌 Qi 系统的图像加密算法具有密钥空间大安全性高等特点. 关键词 : 分数阶 ; 超混沌 Qi 系统 ; 仿真分析 ; 图像加密中图分类号 :TN 文献标志码 :A 文章编号 : (2015) A Fractional-order Hyper-chaotic System and its Application in Image Encryption XUE Wei,XIAO Hui,XU Jinkang,JIA Hongyan (College of Electronic Information and Automation,TIanjin University of Science & Technology,Tianjin ,China) Abstract:A commensurate fractional-order hyper-chaotic Qi system was analyzed by means of numerical simulation including drawing its phase portraits,poincaré maps,lyapunov exponent spectra and bifurcation diagrams. Simulation results verified the chaotic characteristics of the fractional order system. After that the application of the fractional-order hyperchaotic Qi system to image encryption was investigated,and the histogram,adjacent pixels correlation,key sensitivity of the ciphertext were analyzed. Results show that the algorithm of image encryption based on fractional hyper-chaotic Qi system has a large key space and high security. Key words:fractional-order;hyper-chaotic Qi system;simulation analysis;image encryption 1963 年,Lorenz [1] 提出了著名的 Lorenz 混沌系统数学模型及混沌理论的基本特征. 之后, 许多学者相继提出一些新的混沌系统, 如 Chen 系统 [2] Lü 系统 [3] Liu 系统 [4] [5] Qi 系统等, 使得混沌动力学在非线性研究领域取得了重大的发展. 由于混沌系统具有内随机性不可预测性宽带特性遍历性以及对初始条件和系统参数的极端敏感性等, 因此, 混沌系统及其应用研究引起学者们极大的研究热情. 尤其是近年来, 随着分数阶微积分理论的发展和应用, 一些学者将分数阶微积分算子引入混沌系统中, 生成分数阶混沌系统并对其进行研究 [6 9]. 研究表明, 当混沌系统为分数阶时, 系统仍然能呈现出混沌状态, 且更能反映系统所呈现的实际物理现象 [10 11]. 因此, 研究分数阶混沌系统及其在各个领域的应用更具有普适性和实用价值. 混沌变换所具有的对参数和初值的敏感性等基本特性及其与密码学的天然关系早在 Shannon [12] 的经典文献中就已提到, 并据此提出了密码学中用于指导密码设计的 2 个基本原则 : 扩散和混乱. 现在混沌动力学正由数学和物理的基础理论研究逐步过渡到实际的工程应用领域, 并得到了很大发展, 例如混沌理论可用在保密通信图像加密等信息领域 [13 15]. 自 Matthews [16] 提出混沌加密思想以来, 有关混沌密码和混沌保密通信方案的研究逐步成为研究热收稿日期 : ; 修回日期 : 基金项目 : 国家自然科学基金青年基金资助项目 ( ) 作者简介 : 薛薇 (1963 ), 女, 河北河间人, 教授,xuewei@tust.edu.cn. 数字出版日期 : ; 数字出版网址 :

68 天津科技大学学报第 30 卷第 5 期点. 以往研究混沌加密技术大多数基于低维整数阶混沌系统, 虽然低维混沌系统由于形式简单而具有计算时间短的优点, 但其密钥空间小序列的复杂度不高, 导致密码系统的安全性不高. 而高维混沌系统尤其是超混沌系统, 由于具有 4 个及以上的状态变量, 因此密钥空间更大, 这些特点使得超混沌系统用于图像加密会提高系统的安全性.

2 68 天津科技大学学报第 30 卷第 5 期点. 以往研究混沌加密技术大多数基于低维整数阶混沌系统, 虽然低维混沌系统由于形式简单而具有计算时间短的优点, 但其密钥空间小序列的复杂度不高, 导致密码系统的安全性不高. 而高维混沌系统尤其是超混沌系统, 由于具有 4 个及以上的状态变量, 因此密钥空间更大, 这些特点使得超混沌系统用于图像加密会提高系统的安全性. 随着人们对分数阶混沌系统研究的深入, 发现分数阶混沌系统除具有整数阶混沌系统对初值敏感性伪随机性等优良特性外, 还具有分数阶混沌系统所特有的一些性质, 如能够反映系统的历史信息很强的历史记忆性等 [17]. 分数阶混沌系统增强了系统的非线性和复杂性, 可增大密码系统的密钥空间 [18]. 并且, 利用现有的整数阶混沌分析和预测方法并不能估计出分数阶混沌系统的分数阶次, 即现有的整数阶混沌系统分析和预测方法对分数阶混沌系统失效. 因此, 有必要探索研究基于分数阶混沌系统的数字图像加密方法. 考虑到 Qi 系统与 Chen 系统 Liu 系统等系统相比较, 数学模型在不增加复杂度的情况下, 既可以产生两翼的混沌吸引子, 又可以产生四翼的混沌吸引子, 其动力学特性更为复杂. 所以, 本文针对一个同量分数阶超混沌 Qi 系统及其在图像加密中的应用进行研究. 1 分数阶超混沌 Qi 系统的数值分析 Fig. 1 图 1 式 (1) 系统随变化的分岔图 Bifurcation diagram of system shown by euation (1) on x axis with parameter 1.1 相轨迹图当取参数 a = 15, b = 40, c = 5, d = 20, e = 5, f = 50, = 0.9 时, 分数阶超混沌 Qi 系统的相轨迹如图 2 所示. (a) y z 文献 [19] 中提出了超混沌 Qi 系统, 这个四维超混沌系统是通过对文献 [5] 中的混沌系统引入一个线性反馈控制器得到的. 本文在文献 [19] 的基础上, 分析同量分数阶超混沌 Qi 系统, 数学模型如下 : d x = a( y x) + eyz, dt d y = cx + dy xz w dt (1) d z = bz + xy dt d w = fy dt 其中 0<< 1. 当 a = 15, b = 40, c = 5, d = 20, e = 5, f = 50 时, 可绘制式 (1) 系统随阶次变化的分岔图如图 1 所示. 从图 1 中可以看出 [ 0.761,1] 时, 式 (1) 系统呈现混沌态, 这里选取 = 0.9 进行仿真分析和研究., 对 3.6 阶超混沌 Qi 系统 (b) x y z 图 2 式 (1) 系统的相轨迹图 Fig. 2 Phase portraits of the system shown by euation (1) 1.2 Lyapunov 指数及分岔图 Lyapunov 指数 (L) 可以表征系统运动特征, 它沿某一方向取值的正负和大小, 表示系统长时间在吸引子中相邻轨道沿该方向平均发散 (L>0) 或收敛 (L< 0) 的快慢程度. 对于超混沌系统必然存在 2 个或 2 个以上正的 Lyapunov 指数. 该分数阶超混沌系统随参数 a 变化的 Lyapunov 指数谱图和分岔图分别如图 3 和图 4 所示.

2015 年 10 月薛薇, 等 : 一个分数阶超混沌系统及其在图像加密中的应用 69 (a) a [0,60] 化时的分岔图与其对应的 Lyapunov 指数谱图一致. 此外, 系统的 Lyapunov 指数 L 4 始终小于 -10, 为了方便观察,y 轴截取为 [ 10,10], 如图 3(b) 所示. 当 a = 15 时, L 1 = 2.199, L 2 = 0.

3 2015 年 10 月薛薇, 等 : 一个分数阶超混沌系统及其在图像加密中的应用 69 (a) a [0,60] 化时的分岔图与其对应的 Lyapunov 指数谱图一致. 此外, 系统的 Lyapunov 指数 L 4 始终小于 -10, 为了方便观察,y 轴截取为 [ 10,10], 如图 3(b) 所示. 当 a = 15 时, L 1 = 2.199, L 2 = , L 3 = 0, L 4 = 44.64, 该系统是超混沌的. 随着参数 a 的改变, 式 (1) 系统平衡点的稳定性将会发生改变, 从而系统也将处于不同的状态.a 分别取时, 式 (1) 系统呈现混沌周期拟周期等不同的动力学特性, 相轨迹如图 5 所示. (b) a [0,60] ;y 轴截取 [ 10,10] 图 3 式 (1) 系统随参数 a 变化的 Lyapunov 指数谱图 Fig. 3 Lyapunov exponent diagrams of system shown by euation (1) with parameter a 由图 3 和图 4 可以看出, 式 (1) 系统随参数 a 变图 4 Fig. 4 式 (1) 系统随参数 a 变化在 x 方向上的分岔图 ( a [0,60] ) Bifurcation diagram of the system shown by euation (1) on x axis with parameter a( a [0,60] ) (a) a = 15 ( 混沌 ) (b) a = 16 ( 周期 ) (c) a = 35 ( 拟周期 ) 图 5 系统 (1) 在参数 a 取不同值时的相轨迹图 Fig. 5 Phase portraits of system (1) with different a 1.3 Poincaré 截面当 a = 15 时, 式 (1) 系统在 x = 25 和 y = 10 处两个截面上的 Poincaré 映射见图 6,Poincaré 截面是一些具有分形特征的密集点, 说明系统是混沌或超混沌的. (a) x=25 (b) y = 10 图 6 式 (1) 系统的 Poincaré 映射图 Fig. 6 Poincaré map of the system shown by euation (1)

70 天津科技大学学报第 30 卷第 5 期 2 分数阶超混沌 Qi 系统在图像加密中的应用 2.

因此, 采用一维时域乱序加密再三维空域加密的算法. 加密与解密流程见图 7. 初始值改为 x0 = y0 = z0 = w0 = 0.

图 9 错误密钥解密结果 Fig. 9 Decryption result of the wrong key 2.

可以看出 : 加密前各灰度级的像素分布波动较大, 波峰波谷间差值较大 ; 加密后各灰度级的像素分布比较均匀, 各灰度级的像素出现的频率基本相同,

7 Flow chart of encryption and decryption 最后得到彩色图像的加密结果如图 8(a) 所示, 解密结果如图 8(b) 所示.

4 70 天津科技大学学报第 30 卷第 5 期 2 分数阶超混沌 Qi 系统在图像加密中的应用 2.1 彩色图像的加密与解密彩色图像是由 RGB 三基色来描述的, 即图像在计算机中的存储是 m n 3 的矩阵形式, 而分数阶超混沌 Qi 系统是一个四维连续自治系统. 因此, 采用一维时域乱序加密再三维空域加密的算法. 加密与解密流程见图 7. 初始值改为 x0 = y0 = z0 = w0 = , 其他参数值不变时, 解密失败的图像. 可以看出 : 很微小的密钥变化也不能成功解密图像, 测试结果表明该加密系统对密钥具有较高的敏感性. 图 9 错误密钥解密结果 Fig. 9 Decryption result of the wrong key 直方图分析图 10 为加密前后 R 基色分量的灰度直方图, 描述图像中各灰度级的像素个数. 可以看出 : 加密前各灰度级的像素分布波动较大, 波峰波谷间差值较大 ; 加密后各灰度级的像素分布比较均匀, 各灰度级的像素出现的频率基本相同, 很好地掩盖了原图像的分布规律, 能够很好地抵抗统计分析的攻击. (a) 加密 (b) 解密图 7 加密与解密流程图 Fig. 7 Flow chart of encryption and decryption 最后得到彩色图像的加密结果如图 8(a) 所示, 解密结果如图 8(b) 所示. (a) 加密前 (a) 加密 (b) 解密图 8 原始图像的加密与解密结果 Fig. 8 Results of encryption and decryption 2.2 算法安全性分析密钥敏感性分析取系统参数 = 0.9, a = 15, b = 40, c = 5, d = 20, f = 50, 系统初始值 x0 = y0 = z0 = w0 = 0.1为密钥, 对图像进行加密. 当输入正确的密钥时, 可解密成功得到正确的解密图像如图 8(b) 所示, 图 9 为将 (b) 加密后图 10 图像加密前后 R 基色的直方图 Fig. 10 Histograms before and after image encryption 相邻像素的相关性分析分别从原始图像和加密后图像中任取一行一列像素点进行相邻像素相关系数计算, 求得图像在水平垂直相邻像素的相关系数见表 1. 由表 1 结果可知 : 此算法的加密信号具有较好的扩散和混淆. 原始

5 2015 年 10 月薛薇, 等 : 一个分数阶超混沌系统及其在图像加密中的应用 71 图像的相邻像素的相关系数接近 1, 具有高度相关性. 而加密后的相关性系数接近 0, 相邻像素已基本不相关. 因此, 该算法对统计分析的攻击具有更好的安全性. 表 1 相邻像素的相关系数 Tab. 1 The correlation of adjacent pixels 图像水平方向相关系数垂直方向原始图像 0.977, ,50 加密图像 0.004,485, ,917,6 3 结语本文分析了同量分数阶超混沌 Qi 系统的混沌特性, 计算并绘制了该分数阶超混沌系统的 Lyapunov 指数图分岔图 Poincaré 截面映射图等, 验证了该分数阶超混沌吸引子的存在性. 最后, 将该分数阶超混沌系统应用于彩色图像加密, 采用一维时域乱序加密三维空域再加密的算法, 利用所提出的分数阶超混沌 Qi 系统对数字图像进行加密, 并对该算法的安全性进行分析. 分析表明, 该算法具有密钥敏感性强密文分布均匀密文相邻像素的相关性低等特点, 说明了分数阶超混沌 Qi 系统在图像加密中的应用具有较好的前景. 参考文献 : [1] Lorenz E N. Deterministic non-periodic follow[j]. Journal of Atmospheric Sciences,1963,20: [2] Chen G R,Ueta T. Yet another chaotic attractor[j]. International Journal of Bifurcation and Chaos,1999, 9(7): [3] Lü J H,Chen G R. A new chaotic attractor coined[j]. International Journal of Bifurcation and Chaos,2002, 12(3): [4] Liu C X,Liu L,Liu T,et al. A new butterfly-shaped attractor of Lorenz-like system[j]. Chaos,Solitons and Fractals,2006,28(5): [5] Qi G Y,Chen G R,van Wyk M A,et al. A four-wing chaotic attractor generated from a new 3-D uadratic autonomous system[j]. Chaos,Solitons and Fractals, 2008,38(3): [6] Jia H Y,Chen Z Q,Xue W. Analysis and circuit implementation for the fractional-order Lorenz system[j]. Acta Physica Sinica,2013,62(14): [7] Huang L L,Xin F,Wang L Y. Circuit implementation and control of a new fractional-order hyperchaotic system[j]. Acta Physica Sinica,2011,60(1): [8] Grigorenko I,Grigorenko E. Chaotic dynamics of the fractional Lorenz system[j]. Physical Review Letters, 2003,91(3): [9] Li C P,Peng G J. Chaos in Chen s system with a fractional order[j]. Chaos,Solitons and Fractals,2004, 22(2): [10] Hifer R,Butzer P L,Westphal U,et al. Applications of Fractional Calculus in Physics[M]. New Jersey:World Scientific,2001. [11] Zhou P,Wei L J,Cheng X F. A novel fractional-order hyperchaotic system and its synchronization[j]. Chinese Physics B,2009,18(7): [12] Shannon C E. Communication theory of secrecy systems[j]. Bell System Technology Journal,1949,28(4): [13] Li W,Hao J H,Qi B. An encryption approach to chaos synchronization communications by using CPRNG [J]. Acta Physica Sinica,2008,57(3): [14] Peng Z P,Wang Ch H,Lin Y,et al. A novel fourdimensional multi-wing hyper-chaotic attractor and its application in image encryption[j]. Acta Physica Sinica, 2014,63(24): [15] Yu Ch Y,Zhu J L,Guo J Y. Four-dimensional hyperchaotic system and application research in signal encryption[j]. Electric Machines and Control,2012,16(3): [16] Matthews R. On the derivation of a chaotic encryption algorithm[j]. Cryptologia,1989,13(1): [17] Wang Q Y,Zhou S B. Image encryption algorithm based on fractional-order Chen chaotic system[j]. Journal of Computer Applications,2013,33(4): [18] Wu X,Wang H,Lu H. Modified generalized projective synchronization of a new fractional-order hyperchaotic system and its application in secure communication[j]. Nonlinear Analysis :Real World Applications,2012, 13(3): [19] Xue W,Qi G Y,Mu J J,et al. Hopf bifurcation analysis and circuit implementation for a novel four-wing hyperchaotic system[j]. Chinese Physics B,2013,22(8): 责任编辑 : 常涛

4 47 Lorenz [9-11]. M.. [16]... 1 [16] ẋ 1 = a( ) ẋ 2 = b + c (1) ẋ 3 = h + 1 a, b, c, h. a = 2, b = 14, c = 1.6, h = Lyapunov L 1 = , L

4 47 Lorenz [9-11]. M.. [16]... 1 [16] ẋ 1 = a( ) ẋ 2 = b + c (1) ẋ 3 = h + 1 a, b, c, h. a = 2, b = 14, c = 1.6, h = Lyapunov L 1 = , L 28 4 21 7 JOUNAL OF APPLIED SCIENCES Electronics and Information Engineering Vol. 28 No. 4 Jul. 21 DOI: 1.3969/j.issn.255-8297.21.4.13 1,2 1,3 1. 51641 2. 3439 3. 51641 Lyapunov. Lyapunov.... Lyapunov