CN111294198A - 一种基于混沌系统的自适应加密通讯方法 - Google Patents

一种基于混沌系统的自适应加密通讯方法 Download PDF

Info

Publication number
CN111294198A
CN111294198A CN202010251513.5A CN202010251513A CN111294198A CN 111294198 A CN111294198 A CN 111294198A CN 202010251513 A CN202010251513 A CN 202010251513A CN 111294198 A CN111294198 A CN 111294198A
Authority
CN
China
Prior art keywords
state
fractional order
chaotic system
receiving end
chaotic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010251513.5A
Other languages
English (en)
Other versions
CN111294198B (zh
Inventor
宋磊
刘玲亚
于文彬
王易因
雷军委
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202010251513.5A priority Critical patent/CN111294198B/zh
Publication of CN111294198A publication Critical patent/CN111294198A/zh
Application granted granted Critical
Publication of CN111294198B publication Critical patent/CN111294198B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L9/00Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
    • H04L9/001Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols using chaotic signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Security & Cryptography (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于混沌系统的自适应加密通讯方法,属于保密通讯领域。其通过在发射端构造分数阶微分型混沌系统,生成类随机的混沌状态序列对待加密信号进行调制并发送至接收端,然后在接收端根据模型参数构建结构相同的分数阶混沌系统,通过对比解调后的混沌状态得到状态误差量,再根据状态误差进行非线性积分构建滑模面信号,然后构建分数阶自适应估计器对系统不确定性进行估计。最后组成同步控制量实现接收端与发送端混沌系统的同步,从而实现最终加密信息的解密恢复。该方法利用了分数阶混沌系统对初始条件以及解算方法的极为敏感性,使得解密恢复十分困难,从而提高了加密通讯的安全性。

Description

一种基于混沌系统的自适应加密通讯方法
技术领域
本发明属于数据安全、加密解密与保密通讯领域,尤其是涉及基于分数阶自适应的混沌系统加密与解密通讯方法。
背景技术
保密通讯由于其在民用与军用领域都有着至关重要的应用,因而引起了世界范围内的广泛研究与关注。除了采用傅里叶变换等频域加密解密滤波恢复的方法外,采用混沌系统进行直接与间接进行加密掩盖的方法,在近年来也得到了深入的研究。其基本原理在于利用混沌系统对初值以及其它条件不同的敏感性以及混沌系统的类随机特性,使得有用信号完全掩盖在类随机的混沌信号中,难以被窃密方破解。随着分数阶微分的研究与计算机技术的发展,对分数阶的系统的解算方法越来越丰富与方便。而分数阶混沌系统相对于整数阶混沌系统来说,具有更加复杂的动态特性,而且其对解算方法要求非常严格,不同解算方法其计算量不同,而且解算结果也差别较大。因此分数阶混沌系统处理计算复杂的不足之处外,比普通混沌系统更适合应用于保密通讯中。基于以上背景原因,本发明提出了一种利用分数阶混沌与分数阶自适应同步来实现保密通讯的方法。
需要说明的是,在上述背景技术部分发明的信息仅用于加强对本发明的背景的理解,因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本发明的目的在于提供一种基于混沌系统的自适应加密通讯方法,进而至少在一定程度上克服整数阶混沌保密通讯解密解算过于简单与安全性不足的问题。
本发明提供了一种基于混沌系统的自适应加密通讯方法,包括以下步骤:
步骤S10:在发射端根据分数阶微分的计算机生成方法,设置发射端混沌系统的初始状态,构造分数阶混沌系统,得到发射端分数阶混沌的自由状态;
步骤S20:根据所述的发射端分数阶混沌系统状态,与需要进行加密传送的信息,进行信号混合调试与矩阵变换,然后将调试变换后的信号发送至接受端,同时将系统模型参数与调试变换参数作为密钥告知接收端混沌系统。
步骤S30:在远程接收端根据相应的系统模型参数,设置同步控制量的初始值,构造接收端的分数阶混沌系统,生成接受端分数阶混沌系统的状态。
步骤S40:在远程接收端接收发射端的发射状态数据,并根据相应的矩阵参数作为密钥进行解调,得到解调状态参数。并与接收端混沌系统状态进行对比,得到状态误差变量。
步骤S50:根据所述的误差状态,进行非线性积分,组成非线性滑模面信号;
步骤S60:根据所述的非线性滑模面信号,构造分数阶自适应估计器,得到系统不确定性的估计值;
步骤S70:根据所述的系统不确定性估计值与非线性滑模信号进行组合,形成最终的滑模自适应分数阶同步控制量,然后完成接收端混沌系统状态的解算,再根据调试参数对系统进行加密信息解密恢复。
在本发明的一种示例实施例中,根据分数阶微分的计算机生成方法,设置发射端混沌系统的初始状态,构造分数阶混沌系统,得到发射端分数阶混沌的自由状态包括:
zdα1(n)=c1(z2(n)-z1(n))+0.2z2(n)-0.05z3(n)+sin(z1(n)z3(n));
zdα2(n)=(c2-0.3)z1(n)-z2(n)-z1(n)z3(n)+cos(z2(n)+z3(n));
zdα3(n)=z1(n)z2(n)-(c3+0.5)z3(n)+z1(n)sin(z2(n));
其中z1、z2与z3为发射端三阶分数阶混沌系统的状态分别为,其初始状态为z1(1)、z2(1)、z3(1),而初始值设置为z1(1)=z10、z2(1)=z20、z3(1)=z30。z10、z20与z30为常值参数,c1、c2、c3为模型参数,其选取详见后文案例实施。z1(n)、z2(n)、z3(n)为发射端分数阶混沌系统的状态z1、z2与z3的第n个数据。zdα1、zdα2与zdα3分别为状态z1、z2与z3的α阶导数。α为分数阶的阶次,0<α<1,其详细设置见后文案例实施。
在本发明的一种示例实施例中,根据所述的发射端分数阶混沌系统状态,与需要进行加密传送的信息,进行信号混合调试与矩阵变换包括:
w1=z1+b1w;
w2=z2+b2w;
w3=z3+b3w;
Figure BDA0002435652090000041
其中w为待加密发送的信号,z1、z2、z3为发射端混沌系统的状态,w1、w2、w3为信号混合调试得到新的状态。b1、b2、b3为调试参数。
A为三阶可逆矩阵,
Figure BDA0002435652090000042
其逆矩阵记作A-1,满足A-1A=E,其中E为单位矩阵。m1、m2、m3为最终所求的发射状态。
在本发明的一种示例实施例中,根据相应的系统模型参数,设置同步控制量的初始值,构造接收端的分数阶混沌系统,生成接受端分数阶混沌系统的状态包括:
zdα4(n)=c1(z5(n)-z4(n))+0.2z5(n)-0.05z6(n)+u1(n);
zdα5(n)=(c2-0.3)z4(n)-z5(n)-z4(n)z6(n)+u2(n);
zdα6(n)=z4(n)z5(n)-(c3+0.5)z6(n)+u3(n);
其中u1、u2、u3为接收端分数阶混沌系统的同步控制量,u1、u2、u3的初始值为0,即u1(1)=0,u2(1)=0,u3(1)=0。z4、z5与z6为接收端三阶分数阶混沌系统的状态,其初始状态为z4(1)、z5(1)、z6(1),而初始值设置为z4(1)=z40、z5(1)=z50、z6(1)=z60,其中z40、z50与z60为常值参数,详细设计见后文案例实施。z4(n)、z5(n)、z6(n)为接收端分数阶混沌系统的状态z4、z5与z6的第n个数据。zdα4、zdα5与zdα6分别为状态z4、z5与z6的α阶导数。
在本发明的一种示例实施例中,根据相应的矩阵参数作为密钥进行解调,得到解调状态参数,并与接收端混沌系统状态进行对比,得到状态误差变量包括:
Figure BDA0002435652090000051
e1=z4-w1a
e2=z5-w2a
e3=z6-w3a
其中m1a、m2a、m3a为远程接收端接收发射状态数据,A-1为所述矩阵A的逆矩阵,w1a、w2a、w3a为解调状态参数,e1、e2、e3为相应的误差状态。
在本发明的一种示例实施例中,根据所述的误差状态,进行非线性积分,组成非线性滑模面信号包括:
Figure BDA0002435652090000052
Figure BDA0002435652090000053
Figure BDA0002435652090000054
f1s=∫f1dt;
f2s=∫f2dt;
f3s=∫f3dt;
Figure BDA0002435652090000061
Figure BDA0002435652090000062
Figure BDA0002435652090000063
其中e1、e2、e3为误差状态,f1、f2、f3为非线性误差量,f1s、f2s、f3s为非线性积分量,dt表示对时间信号进行积分,p1、p2、p3记为相应的非线性滑模面信号。其中d1、d2、d3、d4、d5、d6、d7、d8、d9、ε1、ε2、ε3以及j11、j12、j21、j22、j31、j32为常值参数,其详细设计见后文案例实施。
在本发明的一种示例实施例中,根据所述的非线性滑模面信号,构造分数阶自适应估计器,得到系统不确定性的估计值包括:
xdα1=-k11p1-k12u1-g1
g1=x1+p1
xdα2=-k21p2-k22u2-g2
g2=x2+p2
xdα3=-k31p3-k32u3-g3
g3=x3+p3
针对i=1,2,3,其中pi为非线性滑模面信号,xi为自适应估计器的分数阶状态,其初始值选为0。xdαi为xi的分数阶微分。ui为分数阶同步控制量,详细设计见下一步。gi即为所求的第i个子系统不确定性估计值。其中ki1与ki2为常值参数,其详细设计见后文案例实施。
在本发明的一种示例实施例中,根据所述的系统不确定性估计值与非线性滑模信号进行组合,计算最终的滑模自适应分数阶同步控制量,再根据调试参数对系统进行加密信息解密恢复包括:
u1=p1-kg1g1
u2=p2-kg2g2
u3=p3-kg3g3
w1b=w1a-z4
w2b=w2a-z5
w3b=w3a-z6
wa=wib/bi
其中g1、g2、g3为系统不确定性估计值,p1、p2、p3为非线性滑模信号,u1、u2、u3为最终的滑模自适应分数阶同步控制量,kg1、kg2、kg3为常值参数,其详细设计见后文案例实施。z4、z5与z6为接收端混沌系统状态,w1a、w2a与w3a为接收端解调状态,w1b、w2b、w2b为恢复状态,bi为解调参数且bi不全为0。
wa为最终的解密恢复的信号。从而通过上述方法,即可实现整个信息加密与解密恢复的全过程,实现数据的保密安全通讯。
有益效果
本发明方法提供的采用分数阶混沌系统进行加密与自适应解密的方法,其特点在于利用了分数阶混沌系统的复杂性,以及对解算方法的极为敏感性与复杂性,使得窃密方不仅在密钥、初始条件、以及解算方法上的任何差别,都将导致解密解算信号完全失真。同时在信号发送中,进行了调制与矩阵变换,不仅使得加密信号能够放缩掩盖在混沌信号中,而且使得同步容易实现,从而也增大了保密通讯的安全性。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种基于混沌系统的自适应加密通讯方法的设计实施流程图;
图2是本发明实施例所提供方法的混沌系统误差状态e1曲线(无单位);
图3是本发明实施例所提供方法的混沌系统误差状态e2曲线(无单位);
图4是本发明实施例所提供方法的混沌系统误差状态e3曲线(无单位);
图5是本发明实施例所提供方法的滑模自适应分数阶同步控制量u1曲线(无单位);
图6是本发明实施例所提供方法的滑模自适应分数阶同步控制量u2曲线(无单位);
图7是本发明实施例所提供方法的滑模自适应分数阶同步控制量u3曲线(无单位);
图8是本发明实施例所提供方法的解密恢复信息与原信息的对比图(无单位);
图9是本发明实施例所提供方法的接收端与发送端第一个混沌子系统的状态同步对比曲线(无单位);
图10是本发明实施例所提供方法的接收端与发送端第二个混沌子系统的状态同步对比曲线(无单位);
图11是本发明实施例所提供方法的接收端与发送端第三个混沌子系统的状态同步对比曲线(无单位)。
具体实施方式
现在将参考附图更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本发明将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中,提供许多具体细节从而给出对本发明的实施方式的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而省略所述特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知技术方案以避免喧宾夺主而使得本发明的各方面变得模糊。
本发明一种基于混沌系统的自适应加密通讯方法,其通过在发射端与接收端构建分数阶混沌系统,再对混沌系统与有用信号进行调制掩盖,再接收端利用矩阵逆变换进行初步解调,再利用误差反馈构造积分性非线性滑模面,以及利用滑模面信息构造分数阶不确定性自适应估计器,最终形成滑模自适应同步控制量,实现了发射端与接受端混沌信号的同步。该方法利用了分数阶混沌系统的解算复杂性以及对细微差别的敏感性,使得同步与解密十分困难,从而加大了保密通讯的安全性。
以下,将结合附图对本发明实例实施例中提及的一种基于混沌系统的自适应加密通讯方法进行解释以及说明。参考图1所示,一种基于混沌系统的自适应加密通讯方法可以包括以下步骤:
步骤S10:在发射端,根据分数阶微分的计算机生成方法,设置发射端混沌系统的初始状态,构造分数阶混沌系统,得到发射端分数阶混沌的自由状态;
具体的,首先设置分数阶的阶次α,0<α<1,其详细设置见后文案例实施。然后设置发射端三阶分数阶混沌系统的状态分别为z1、z2与z3,其初始状态为z1(1)、z2(1)、z3(1),而初始值设置为z1(1)=z10、z2(1)=z20、z3(1)=z30,其中z10、z20与z30为常值参数,详细设计见后文案例实施。
最后,按照下面的方程,构建发射端的三阶混沌系统。
zdα1(n)=c1(z2(n)-z1(n))+0.2z2(n)-0.05z3(n)+sin(z1(n)z3(n));
zdα2(n)=(c2-0.3)z1(n)-z2(n)-z1(n)z3(n)+cos(z2(n)+z3(n));
zdα3(n)=z1(n)z2(n)-(c3+0.5)z3(n)+z1(n)sin(z2(n));
其中c1、c2、c3为模型参数,其选取详见后文案例实施。z1(n)、z2(n)、z3(n)为发射端分数阶混沌系统的状态z1、z2与z3的第n个数据。zdα1、zdα2与zdα3分别为状态z1、z2与z3的α阶导数。其中发射端分数阶混沌系统的状态z1、z2与z3的更新计算方法按照一般预测递推的分数阶微分的求解方法即可,在此不在重复累述。
步骤S20:根据所述的发射端分数阶混沌系统状态,与需要进行加密传送的信息,进行信号混合调试与矩阵变换,然后将调试变换后的信号发送至接受端,同时将系统模型参数与调试变换参数作为密钥告知接收端混沌系统。
具体的,首先假设待加密发送的信号为w,选取发射端混沌系统的状态z1、z2、z3进行信号混合调试得到新的状态w1、w2、w3
其中
w1=z1+b1w;w2=z2+b2w;w3=z3+b3w;
其中b1、b2、b3为调试参数。
然后选取三阶可逆矩阵
Figure BDA0002435652090000111
其逆矩阵记作A-1,满足
A-1A=E,其中E为单位矩阵。然后进行如下的矩阵变换,得到发射状态,记作m1、m2、m3,其变换方式如下:
Figure BDA0002435652090000121
最后,将发射状态m1、m2、m3从发射端发送出去,有待接受端远程接收。同时将调试参数b1、b2、b3与矩阵A-1以及系统模型参数c1、c2、c3作为密钥告知接收方。
步骤S30:在远程接收端根据相应的系统模型参数,设置同步控制量的初始值,构造接收端的分数阶混沌系统,生成接受端分数阶混沌系统的状态。
具体的,设置接收端分数阶混沌系统的同步控制量u1、u2、u3的初始值为0,即u1(1)=0,u2(1)=0,u3(1)=0。
同时设置接收端三阶分数阶混沌系统的状态分别为z4、z5与z6,其初始状态为z4(1)、z5(1)、z6(1),而初始值设置为z4(1)=z40、z5(1)=z50、z6(1)=z60,其中z40、z50与z60为常值参数,详细设计见后文案例实施。
最后,按照下面的方程,构建接收端的三阶混沌系统。
zdα4(n)=c1(z5(n)-z4(n))+0.2z5(n)-0.05z6(n)+u1(n);
zdα5(n)=(c2-0.3)z4(n)-z5(n)-z4(n)z6(n)+u2(n);
zdα6(n)=z4(n)z5(n)-(c3+0.5)z6(n)+u3(n);
其中c1、c2、c3为模型参数,也是接收端获得的密钥。z4(n)、z5(n)、z6(n)为接收端分数阶混沌系统的状态z4、z5与z6的第n个数据。zdα4、zdα5与zdα6分别为状态z4、z5与z6的α阶导数。
步骤S40:在远程接收端接收发射端的发射状态数据,并根据相应的矩阵参数作为密钥进行解调,得到解调状态参数。并与接收端混沌系统状态进行对比,得到状态误差变量。
具体的,首先,在远程接收端接收发射状态数据,记作m1a、m2a、m3a
然后,根据矩阵A-1进行解调,得到解调状态,记作w1a、w2a、w3a,其解调过程如下:
Figure BDA0002435652090000131
最后,将解调状态与混沌系统状态进行比较,得到相应的误差状态,记作e1、e2、e3,其计算方式如下:
e1=z4-w1a;e2=z5-w2a;e3=z6-w3a
步骤S50:根据所述的误差状态,进行非线性积分,组成非线性滑模面信号;
具体的,针对所述的误差状态e1、e2、e3,进行非线性变换,得到如下,得到非线性误差量,分别记作f1、f2、f3,其计算方式如下:
Figure BDA0002435652090000132
Figure BDA0002435652090000133
Figure BDA0002435652090000134
其次,对所述的非线性误差量进行积分,得到非线性积分量,记作f1s、f2s、f3s,其计算方式如下:
f1s=∫f1dt;
f2s=∫f2dt;
f3s=∫f3dt;
其中dt表示对时间信号进行积分。
然后,对误差状态与非线性积分量进行组合,得到非线性滑模面信号,记作p1、p2、p3,其组合方式如下:
Figure BDA0002435652090000141
Figure BDA0002435652090000142
Figure BDA0002435652090000143
其中d1、d2、d3、d4、d5、d6、d7、d8、d9、ε1、ε2、ε3以及j11、j12、j21、j22、j31、j32为常值参数,其详细设计见后文案例实施。
步骤S60:根据所述的非线性滑模面信号,构造分数阶自适应估计器,得到不确定性的估计值;
具体的,首先根据所述的非线性滑模面信号p1,构造如下的分数阶自适应估计器,其中x1为自适应估计器的分数阶状态,其初始值选为0。xdα1为x1的分数阶微分。u1为分数阶同步控制量,详细设计见下一步。g1即为所求的第一个子系统不确定性估计值。其计算方式如下:
xdα1=-k11p1-k12u1-g1
g1=x1+p1
其中k11与k12为常值参数,其详细设计见后文案例实施。
其次,根据所述的非线性滑模面信号p2,构造如下的分数阶自适应估计器,其中x2为自适应估计器的分数阶状态,其初始值选为0。xdα2为x2的分数阶微分。u2为分数阶同步控制量,详细设计见下一步。g2即为所求的第二个子系统不确定性估计值。其计算方式如下:
xdα2=-k21p2-k22u2-g2
g2=x2+p2
其中k21与k22为常值参数,其详细设计见后文案例实施。
最后,根据所述的非线性滑模面信号p3,构造如下的分数阶自适应估计器,其中x3为自适应估计器的分数阶状态,其初始值选为0。xdα3为x3的分数阶微分。u3为分数阶同步控制量,详细设计见下一步。g3即为所求的第三个子系统不确定性估计值。其计算方式如下:
xdα3=-k31p3-k32u3-g3
g3=x3+p3
其中k31与k32为常值参数,其详细设计见后文案例实施。
步骤S70:根据所述的系统不确定性估计值与非线性滑模信号进行组合,形成最终的滑模自适应分数阶同步控制量,然后完成接收端混沌系统状态的解算,再根据调试参数对系统进行加密信息解密恢复。
具体的,首先根据所述的系统不确定性估计值g1、g2、g3与非线性滑模信号p1、p2、p3进行线性组合,得到最终的滑模自适应分数阶同步控制量u1、u2、u3如下:
u1=p1-kg1g1
u2=p2-kg2g2
u3=p3-kg3g3
其中kg1、kg2、kg3为常值参数,其详细设计见后文案例实施。
其次,根据上述同步控制量代入接收端混沌系统,解算得到混沌系统状态z4、z5与z6,以及接收端解调状态w1a、w2a与w3a,进行信号解密,得到恢复状态记作w1b、w2b、w2b,其解密方法如下:
w1b=w1a-z4;w2b=w2a-z5;w3b=w3a-z6
最后,选取bi满足,bi≠0,则wa=wib/bi
最终得到的wa记为解密恢复的信号。从而实现信号的混沌加密传输与解密恢复,最终案例仿真表明该方法具有很高的精度。
案例实施与计算机仿真模拟结果分析
为验证本发明例所提供方法的有效性,进行以下的案例仿真。
在步骤一中,设置z10=15、z20=4与z30=6为常值参数,c1=8、c2=150、c3=0.4,α=0.95。
在步骤二中,选取b1=0.1、b2=0、b3=0,设置待加密发送的信息w如下所示w=(0.4+3*rand)*sign(300t);选取
Figure BDA0002435652090000161
在步骤三中,设置u1(1)=0,u2(1)=0,u3(1)=0,z40=-23、z50=-24与z60=-26为常值参数。
在步骤四中,选取
Figure BDA0002435652090000171
得到相应的误差状态e1、e2与e3如图2、图3与图4所示。
在步骤五中,选取参数d1=-4、d2=-5、d3=-5、d4=-3.5、d5=-1、d6=-1、d7=-3、d8=-8、d9=-6、ε1=6、ε2=6、ε3=6以及j11=-40、j12=-20、j21=-35、j22=-25、j31=-45、j32=-15。在步骤六中,选取k11=1.2,k12=0.1,k21=0.8,k22=0.1,k31=1.4,k32=0.1。
在步骤七中,选取kg1=1、kg2=1、kg3=1,得到滑模自适应分数阶同步控制量u1、u2、u3如图5、图6与图7所示。最终得到的解密恢复信息与原信息的对比图如图8所示。而接收端与发送端混沌系统的状态同步对比曲线如图9、图10与图11所示。
有图9、10、11可以看出发送端混沌系统状态与接收端混沌系统状态可以快速实现同步,而图2、3、4的误差收敛状态也表明系统误差大约在0.1秒就能收敛到0附近,从而也能保证图8所示的解密信息能够恢复成原信息曲线,尽管放大后能够观察到两者有细微差别,但可以看出,0.1秒后,两者是几乎重合为一条曲线,因此该解密恢复的失真是非常小的,能满足工程应用的要求。同时上述案例也表明了分数阶微分的精密性,整个解算的精度非常高,仿真步长也比较小,否则迭代误差会导致系统同步失败而使得解密无法完成。从而也表明整个加密与解密方法安全性比较高,破解难度大,使得整个方法具有很高的实用价值。
本领域技术人员在考虑说明书及实践这里发明的发明后,将容易想到本发明的其他实施例。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由权利要求指出。

Claims (8)

1.一种基于混沌系统的自适应加密通讯方法,其特征在于,包括以下步骤:
步骤S10:在发射端根据分数阶微分的计算机生成方法,设置发射端混沌系统的初始状态,构造分数阶混沌系统,得到发射端分数阶混沌的自由状态;
步骤S20:根据所述的发射端分数阶混沌系统状态,与需要进行加密传送的信息,进行信号混合调试与矩阵变换,然后将调试变换后的信号发送至接受端,同时将系统模型参数与调试变换参数作为密钥告知接收端混沌系统。
步骤S30:在远程接收端根据相应的系统模型参数,设置同步控制量的初始值,构造接收端的分数阶混沌系统,生成接受端分数阶混沌系统的状态。
步骤S40:在远程接收端接收发射端的发射状态数据,并根据相应的矩阵参数作为密钥进行解调,得到解调状态参数。并与接收端混沌系统状态进行对比,得到状态误差变量。
步骤S50:根据所述的误差状态,进行非线性积分,组成非线性滑模面信号;
步骤S60:根据所述的非线性滑模面信号,构造分数阶自适应估计器,得到系统不确定性的估计值;
步骤S70:根据所述的系统不确定性估计值与非线性滑模信号进行组合,形成最终的滑模自适应分数阶同步控制量,然后完成接收端混沌系统状态的解算,再根据调试参数对系统进行加密信息解密恢复。
2.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据分数阶微分的计算机生成方法,设置发射端混沌系统的初始状态,构造分数阶混沌系统,得到发射端分数阶混沌的自由状态包括:
zdα1(n)=c1(z2(n)-z1(n))+0.2z2(n)-0.05z3(n)+sin(z1(n)z3(n));
zdα2(n)=(c2-0.3)z1(n)-z2(n)-z1(n)z3(n)+cos(z2(n)+z3(n));
zdα3(n)=z1(n)z2(n)-(c3+0.5)z3(n)+z1(n)sin(z2(n));
其中z1、z2与z3为发射端三阶分数阶混沌系统的状态分别为,其初始状态为z1(1)、z2(1)、z3(1),而初始值设置为z1(1)=z10、z2(1)=z20、z3(1)=z30。z10、z20与z30为常值参数,c1、c2、c3为模型参数。z1(n)、z2(n)、z3(n)为发射端分数阶混沌系统的状态z1、z2与z3的第n个数据。zdα1、zdα2与zdα3分别为状态z1、z2与z3的α阶导数。α为分数阶的阶次,0<α<1。
3.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据所述的发射端分数阶混沌系统状态,与需要进行加密传送的信息,进行信号混合调试与矩阵变换包括:
w1=z1+b1w;
w2=z2+b2w;
w3=z3+b3w;
Figure FDA0002435652080000031
其中w为待加密发送的信号,z1、z2、z3为发射端混沌系统的状态,w1、w2、w3为信号混合调试得到新的状态。b1、b2、b3为调试参数。
A为三阶可逆矩阵,其元素为常值,
Figure FDA0002435652080000032
其逆矩阵记作A-1,满足A-1A=E,其中E为单位矩阵。m1、m2、m3为最终所求的发射状态。
4.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据相应的系统模型参数,设置同步控制量的初始值,构造接收端的分数阶混沌系统,生成接受端分数阶混沌系统的状态包括:
zdα4(n)=c1(z5(n)-z4(n))+0.2z5(n)-0.05z6(n)+u1(n);
zdα5(n)=(c2-0.3)z4(n)-z5(n)-z4(n)z6(n)+u2(n);
zdα6(n)=z4(n)z5(n)-(c3+0.5)z6(n)+u3(n);
其中u1、u2、u3为接收端分数阶混沌系统的同步控制量,u1、u2、u3的初始值为0,即u1(1)=0,u2(1)=0,u3(1)=0。z4、z5与z6为接收端三阶分数阶混沌系统的状态,其初始状态为z4(1)、z5(1)、z6(1),而初始值设置为z4(1)=z40、z5(1)=z50、z6(1)=z60,其中z40、z50与z60为常值参数。z4(n)、z5(n)、z6(n)为接收端分数阶混沌系统的状态z4、z5与z6的第n个数据。zdα4、zdα5与zdα6分别为状态z4、z5与z6的α阶导数。
5.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据相应的矩阵参数作为密钥进行解调,得到解调状态参数,并与接收端混沌系统状态进行对比,得到状态误差变量包括:
Figure FDA0002435652080000041
e1=z4-w1a
e2=z5-w2a
e3=z6-w3a
其中m1a、m2a、m3a为远程接收端接收发射状态数据,A-1为所述矩阵A的逆矩阵,w1a、w2a、w3a为解调状态参数,e1、e2、e3为相应的误差状态。
6.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据所述的误差状态,进行非线性积分,组成非线性滑模面信号包括:
Figure FDA0002435652080000042
Figure FDA0002435652080000043
Figure FDA0002435652080000044
f1s=∫f1dt;
f2s=∫f2dt;
f3s=∫f3dt;
Figure FDA0002435652080000051
Figure FDA0002435652080000052
Figure FDA0002435652080000053
其中e1、e2、e3为误差状态,f1、f2、f3为非线性误差量,f1s、f2s、f3s为非线性积分量,dt表示对时间信号进行积分,p1、p2、p3记为相应的非线性滑模面信号。其中d1、d2、d3、d4、d5、d6、d7、d8、d9、ε1、ε2、ε3以及j11、j12、j21、j22、j31、j32为常值参数。
7.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据所述的非线性滑模面信号,构造分数阶自适应估计器,得到系统不确定性的估计值包括:
xdα1=-k11p1-k12u1-g1
g1=x1+p1
xdα2=-k21p2-k22u2-g2
g2=x2+p2
xdα3=-k31p3-k32u3-g3
g3=x3+p3
针对i=1,2,3,其中pi为非线性滑模面信号,xi为自适应估计器的分数阶状态,其初始值选为0。xdαi为xi的分数阶微分。ui为分数阶同步控制量,详细设计见下一步。gi即为所求的第i个子系统不确定性估计值。其中ki1与ki2为常值参数。
8.根据权利要求1所述的所述的一种基于混沌系统的自适应加密通讯方法,其特征在于,根据所述的系统不确定性估计值与非线性滑模信号进行组合,计算最终的滑模自适应分数阶同步控制量,再根据调试参数对系统进行加密信息解密恢复包括:
u1=p1-kg1g1
u2=p2-kg2g2
u3=p3-kg3g3
w1b=w1a-z4
w2b=w2a-z5
w3b=w3a-z6
wa=wib/bi
其中g1、g2、g3为系统不确定性估计值,p1、p2、p3为非线性滑模信号,u1、u2、u3为最终的滑模自适应分数阶同步控制量,kg1、kg2、kg3为常值参数,其详细设计见后文案例实施。z4、z5与z6为接收端混沌系统状态,w1a、w2a与w3a为接收端解调状态,w1b、w2b、w2b为恢复状态,bi为解调参数且bi不全为0。wa为最终的解密恢复的信号。
CN202010251513.5A 2020-04-01 2020-04-01 一种基于混沌系统的自适应加密通讯方法 Expired - Fee Related CN111294198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010251513.5A CN111294198B (zh) 2020-04-01 2020-04-01 一种基于混沌系统的自适应加密通讯方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010251513.5A CN111294198B (zh) 2020-04-01 2020-04-01 一种基于混沌系统的自适应加密通讯方法

Publications (2)

Publication Number Publication Date
CN111294198A true CN111294198A (zh) 2020-06-16
CN111294198B CN111294198B (zh) 2021-05-14

Family

ID=71030419

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010251513.5A Expired - Fee Related CN111294198B (zh) 2020-04-01 2020-04-01 一种基于混沌系统的自适应加密通讯方法

Country Status (1)

Country Link
CN (1) CN111294198B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112887094A (zh) * 2021-02-26 2021-06-01 湖南科技大学 一种基于单片机的分数阶混沌滑模同步保密通信方法
CN113114453A (zh) * 2021-04-16 2021-07-13 安徽大学 一种基于滑膜控制器的复杂网络保密通信方法
CN113541917A (zh) * 2021-03-16 2021-10-22 新疆大学 一种有限时间修正函数投影同步及参数调制的控制方法
CN114884648A (zh) * 2021-02-05 2022-08-09 中国人民解放军海军航空大学 一种输入饱和不确定分数阶混沌系统的保密通信方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN205453724U (zh) * 2016-03-01 2016-08-10 哈尔滨理工大学 一种多功能通用混沌电路系统
CN107086916A (zh) * 2017-05-23 2017-08-22 西安理工大学 一种基于分数阶自适应滑模控制的混沌系统同步方法
CN107359980A (zh) * 2017-07-31 2017-11-17 哈尔滨理工大学 一种六维分数阶超混沌系统及混沌信号发生器设计
CN108242995A (zh) * 2018-02-26 2018-07-03 长春工业大学 一种基于分段函数法多涡卷混沌吸引子的实现方法
CN108833075A (zh) * 2018-06-21 2018-11-16 东北大学 一种基于非奇异终端滑模控制器的二阶混沌投影同步方法
CN108964871A (zh) * 2018-07-06 2018-12-07 烟台大学 一种基于双陈氏混沌系统与终端滑模的三通道安全通讯方法
CN109104387A (zh) * 2018-08-23 2018-12-28 上海交通大学 一种水声OFDM通信中基于Dual-Zadoff-Chu序列的信号同步方法
US20190156066A1 (en) * 2016-04-07 2019-05-23 The Johns Hopkins University System and method for physical one-way function authentication via chaotic integrated photonic resonators
CN110568759A (zh) * 2019-09-26 2019-12-13 南京理工大学 分数阶混沌系统的鲁棒同步控制方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN205453724U (zh) * 2016-03-01 2016-08-10 哈尔滨理工大学 一种多功能通用混沌电路系统
US20190156066A1 (en) * 2016-04-07 2019-05-23 The Johns Hopkins University System and method for physical one-way function authentication via chaotic integrated photonic resonators
CN107086916A (zh) * 2017-05-23 2017-08-22 西安理工大学 一种基于分数阶自适应滑模控制的混沌系统同步方法
CN107359980A (zh) * 2017-07-31 2017-11-17 哈尔滨理工大学 一种六维分数阶超混沌系统及混沌信号发生器设计
CN108242995A (zh) * 2018-02-26 2018-07-03 长春工业大学 一种基于分段函数法多涡卷混沌吸引子的实现方法
CN108833075A (zh) * 2018-06-21 2018-11-16 东北大学 一种基于非奇异终端滑模控制器的二阶混沌投影同步方法
CN108964871A (zh) * 2018-07-06 2018-12-07 烟台大学 一种基于双陈氏混沌系统与终端滑模的三通道安全通讯方法
CN109104387A (zh) * 2018-08-23 2018-12-28 上海交通大学 一种水声OFDM通信中基于Dual-Zadoff-Chu序列的信号同步方法
CN110568759A (zh) * 2019-09-26 2019-12-13 南京理工大学 分数阶混沌系统的鲁棒同步控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
KARIMA RABAH ECT.: "Fractional Adaptive Sliding Mode Control Laws for Fractional Order Chaotic Systems Synchronization", 《17TH INTERNATIONAL CONFERENCE ON SCIENCES AND TECHNIQUES OF AUTOMATIC CONTROL & COMPUTER ENGINEERING - STA"2016》 *
LICAI LIU, CHUANHONG DU ECT.: "An Circuit Implementation for a Different Fractional-order Chaotic Switching An Circuit Implementation for a Different Fractional-order Chaotic Switching", 《IEEE》 *
SONG,JIANYING YANG,SHIYUN XU: "Chaos synchronization for a class of nonlinear oscillators with fractional order", 《NONLINEAR ANALYSIS》 *
仲启龙,邵永晖,郑永爱: "分数阶混沌系统的主动滑模同步", 《动力学与控制学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114884648A (zh) * 2021-02-05 2022-08-09 中国人民解放军海军航空大学 一种输入饱和不确定分数阶混沌系统的保密通信方法
CN114884648B (zh) * 2021-02-05 2024-01-26 中国人民解放军海军航空大学 一种输入饱和不确定分数阶混沌系统的保密通信方法
CN112887094A (zh) * 2021-02-26 2021-06-01 湖南科技大学 一种基于单片机的分数阶混沌滑模同步保密通信方法
CN112887094B (zh) * 2021-02-26 2022-06-21 湖南科技大学 一种基于单片机的分数阶混沌滑模同步保密通信方法
CN113541917A (zh) * 2021-03-16 2021-10-22 新疆大学 一种有限时间修正函数投影同步及参数调制的控制方法
CN113541917B (zh) * 2021-03-16 2023-05-05 新疆大学 一种有限时间修正函数投影同步及参数调制的控制方法
CN113114453A (zh) * 2021-04-16 2021-07-13 安徽大学 一种基于滑膜控制器的复杂网络保密通信方法

Also Published As

Publication number Publication date
CN111294198B (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
CN111294198B (zh) 一种基于混沌系统的自适应加密通讯方法
Portmann et al. Cryptographic security of quantum key distribution
CN107196763B (zh) Sm2算法协同签名及解密方法、装置与系统
Cherrier et al. Observers-based synchronization and input recovery for a class of nonlinear chaotic models
Tyagi et al. Universal hashing for information-theoretic security
CN111277417B (zh) 一种基于国网安全技术架构的电子签章实现方法
CN107359983B (zh) 基于无非线性项混沌信号掩盖的信号加密传输与恢复方法
Chang et al. Privacy-preserving federated learning via functional encryption, revisited
Gan et al. Design and ARM‐embedded implementation of a chaotic map‐based multicast scheme for multiuser speech wireless communication
Sheu et al. A two-channel secure communication using fractional chaotic systems
Wang et al. Secure compressive sensing of images based on combined chaotic DWT sparse basis and chaotic DCT measurement matrix
CN110012443A (zh) 一种全同态的数据加密聚合方法及其系统
Wu et al. Privacy-Preserving Serverless Federated Learning Scheme for Internet of Things
Shen et al. A novel digital image encryption algorithm based on orbit variation of phase diagram
CN111641500A (zh) 一种用于无人机无线视频传输安全的加解密方法
CN110149201B (zh) 基于误差掩盖与混沌同步的保密通信方法
CN108650074B (zh) 一种基于参数辨识的单通道混沌系统加密通讯方法
CN106656461A (zh) 一种信号量化情形下的混沌神经网络保密通信方法
CN108964871B (zh) 一种基于双陈氏混沌系统与终端滑模的三通道安全通讯方法
Alsaedi Colored image encryption and decryption using multi-chaos 2D quadratic strange attractors and matrix transformations
Han et al. A novel physical layer key generation method based on WGAN-GP adversarial autoencoder
Zhou et al. A survey of security aggregation
CN112751661B (zh) 一种基于同态加密的工业现场设备隐私数据保护方法
CN114337849A (zh) 基于互信息量估计神经网络的物理层保密方法及系统
WO2022185328A1 (en) System and method for identity-based key agreement for secure communication

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210514