基于宽参数双涡卷保守混沌系统的数字信号发生器技术方案

本实用新型专利技术提供了一种基于宽参数双涡卷保守混沌系统的数字信号发生器，通过在FPGA开发板上设计的混沌微分方程计算电路模块可复现一种很大系统参数范围的双涡卷四维保守混沌系统，且该系统具有很大范围的宽参数特性和很强的空间遍历性。并且通过在FPGA数字电路结构的有效改进，极大的节省了资源并且提高了效率。对于求解连续混沌系统的离散化输出时，设计了一种新的基于龙格库塔法的组合逻辑电路和状态机时序，可以在五个时钟周期完成一次迭代，有效加快迭代速率。有效加快迭代速率。有效加快迭代速率。

【技术实现步骤摘要】
基于宽参数双涡卷保守混沌系统的数字信号发生器


[0001]本技术属于混沌系统领域，涉及基于宽参数双涡卷保守混沌系统的数字信号发生器。

技术介绍

[0002]探索混沌系统非线性原理，如何构造结构相对简单、动力学特性更复杂的混沌系统有助于非线性系统研究。目前主流的连续混沌系统有普通的三维连续混沌系统和四维以上的超混沌连续系统，其中对于保守连续混沌系统的研究和应用相对较少。但保守混沌系统往往具有更强的空间遍历性，且已经创立的连续混沌系统大多为单涡卷，而双涡卷或双翅膀混沌系统具有更为复杂的拓扑结构和动力学特性，因此也更适用于信息安全领域。尤其是可以将混沌信号发生器所产生的混沌数字信号，通过PCIE和光纤等快速通信协议输出到需要的应用装置，其应用包含但不限于信息安全、图像加密、机密运算和区块链等领域。
[0003]现在越来越多的研究者选择将混沌系统放在FPGA上去实现，因FPGA相对于普通的冯诺依曼结构的计算机，具有体积更小且没有串行运行速度慢的缺点，即可以做到并行运算、多通道同时输入和输出，极大加快了混沌信号产生速率。而如何设计好实现混沌系统的数字电路，在FPGA中进行系统布局，并因不同的电路结构对于输出的结果影响巨大，如何选择合适的电路结构来平衡硬件资源和输出速率等问题有待商榷。在数字电路的设计中，混沌系统微分方程的最大逻辑组合结构的计算长度，决定了FPGA板内电路模块的最大时钟频率，所以如何用最短的逻辑组合结构的计算长度争取最快的时钟频率而又不会丧失复现混沌系统的准确性，如何用FPGA更节省资源且更快速地复现保守混沌系统，对于混沌系统的应用也是一个可进行扩展性研究的问题。

技术实现思路

[0004]本技术的目的在于针对
技术介绍
中存在的问题，提供一种节约资源且输出速率高的宽参数双涡卷保守混沌系统的信号发生器。
[0005]本技术针对宽参数双涡卷保守混沌系统的信号发生器，设计的一种新的四维双涡卷保守混沌系统方程如下:
[0006][0007]式中，x、y、z、u表示系统的状态变量，a、b、c、d、e表示系统参数，具体地a=
‑
4，b=8，c=
‑
8，d=
‑
8，e=8。本技术使用Matlab2021b软件对系统进行动态特性的分析，并初步采用四阶龙格库塔法对本混沌系统进行离散化求解，系统时间步长h取值为1/64。
[0008]混沌系统动力学特性分析
[0009]李雅普诺夫指数
[0010]李雅普诺夫指数(Lyapunov Exponent，LE)可以分析混沌系统的混乱程度， LE的计算公式如下：
[0011][0012]系统初值分别取x0=1，y0=1，z0=1，u0=1，各系统变量分别为a=
‑
4，b=8，c=
‑
8，d=
‑
8，e=8，步长取1/64，时间长度t为1000，可求得LE1=0.26164433487438，LE2=0.007998255633539，LE3=
ꢀ‑
0.006780134019510，LE4=
‑
0.262862456917076，LE图谱如图4所示：
[0013]由图4可见，
∂
x(t)/
ꢀ∂
x(0)曲线表示x维的LE有一个明显大于0的状态，说明该混沌系统会进入无法预测的运动状态，因而可判定该系统属于混沌系统。
[0014]此外，表征了相空间系统体积的增长率或减小率，计算公式如下：
[0015][0016]通过公式（1）可以计算出公式（2）：
[0017][0018]根据上式可证明该系统的变化是有界的，倘若很明显的小于0时，说明系统空间体积在不断减小，即为耗散型。而此四维混沌系统值非常接近于0，则说明系统的相空间内体积增长和减小的速度在整体上可近似等价为0，因而该系统并非耗散型更不是增长型，而是保守型混沌系统。
[0019]稳定性分析
[0020]为了进一步验证该系统是否为保守混沌系统，求得系统的平衡状态方程为:
[0021][0022]当系统初值分别取x0=1，y0=1，z0=1，u0=1，各系统变量分别为a=
‑
4，b=8，c=
‑
8，d=
‑
8，e=8，通过解方程组（4）可得出三个实数平衡点，依次为R1=（0，0，1，0），R2=（0，0，0，0），R3=（0，0，
‑
1，0），该系统的雅可比矩阵为：
[0023][0024]将系统(1)在平衡点R1处线性化，所对应的Jacobi矩阵为：
[0025][0026]将系统(1)在平衡点R2处线性化，所对应的Jacobi矩阵为：
[0027][0028]将系统(1)在平衡点R3处线性化，所对应的Jacobi矩阵为：
[0029][0030]求特征根和特征向量的公式为：
[0031][0032]根据上式可求得 R1和R3所对应的四个特征值，分别为（8.8098i，
‑
8.8098i，5.1369i，
‑
5.1369i），R2所对应的四个特征值分别为（
‑
7.6693，5.9007i，
‑
5.9007i，7.6693），三个平衡点对应特征值在复平面的分布如图5所示。从图5可知，平衡点R2的特征根分别有两个纯虚数和两个纯实数，且为一正一负，即可判定该平衡点存在一个非稳定的特征根（正实数）及其对应的特征向量，也存在一个稳定的特征根（负实数）及其对应特征向量。因而在相平面上，解的运动轨线将会有一个既发散又存在收敛的趋势，此时可认为该平衡点是系统的一个鞍点。而平衡点R1和R3所对应的特征根在复平面上的分布与平衡点R2的稍有不同，由图5可见，这两个平衡点对应特征根的分布都在虚轴上，即实部为0，因而系统的解轨迹将会是一个以原点为圆心的圆，故可以判定这两个平衡点是系统的一组中心型平衡点，进一步证明了该系统的相空间运动既没有耗散也没有增大的趋势，且这种特点大多存在于保守系统中，因此可以判定该系统为保守混沌系统。
[0033]系统的相图分析
[0034]从相图上可以宏观观察系统特性，在绘制相图的实验中，状态变量初值分别取x0=1，y0=1，z0=1，u0=1，各系统变量分别为a=
‑
4，b=8，c=
‑
8，d=
‑
8，e=8，步长h取值为1/64，时间长度t为600，即每个相图是由38400个坐标点绘制而成，本系统经仿真得到的所有二维平面相图如图6和图7所示。
[0035]从图6和图7可以看出，本技术设计的保守混沌系统呈现出双涡卷混沌吸引子的现象，为了测试该系统的遍历性，本技术将经典的Lorenz普通耗散型系统的三维相图和该技术设计的保守混沌系统的三维相图进行对比如图8所示。从图8可以看出，保守混沌系统本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于宽参数双涡卷保守混沌系统的数字信号发生器，其特征在于，包括：输入模块，所述输入模块用于初始值的输入；S
n
寄存器单元，所述S
n
寄存器单元与输入模块连接，S
n
寄存器单元用于数据的存储；F（Sn）混沌微分方程计算电路模块，所述F（Sn）混沌微分方程计算电路模块与S
n
寄存器单元相连接，F（Sn）混沌微分方程计算电路模块用于对S
n
寄存器单元所传输的数据进行组合逻辑计算；四阶龙格库塔法的四维数值寄存器单元，所述四维数值寄存器单元与F（Sn）混沌微分方程计算电路模块相连接，四维数值寄存器单元存储组合用于存储龙格库塔法的逻辑组合计算电路所计算的数据，四维数值寄存器单元还可向F（Sn）混沌微分方程计算电路模块传送数据；S
n+1
寄存器单元，所述S n+1
寄存器单元分别与四维数值寄存器单元和S
n
寄存器单元相连接，S n+1
寄存器单元用于接收与四维数值寄存器单元所传输的数据，并向S
n
寄存器单元传送数据；状态机选择器，所述状态机选择器分别与S
n
寄存器单元、四维数值寄存器单元和S n+1
寄存器单元相连接，状态机选择器用于输出时序信号；输出模块，所述输出模块与S n+1
寄存器单元相连接。2.根据权利要求1所述的基于宽参数双涡卷保守混沌系统的数字信号发生器，其特征在于，所述F（Sn）混沌微分方程计算电路模块包括x输入寄存器、y输入寄存器、z输入寄存器和u输入寄存器，具体地：所述x输入寄存器连接有第一移位器和第一取反器，所述第一移位器连接有
∂
y输出寄存器，所述第一取反器连接有可执行加1操作的第一加法器，所述第一加法器连接有第二移位器，所述第二移位器连接有第二加法器，所述第二加法器连接有三选一选择器和
∂
u输出寄存器；所述z输入寄存器连接有第三移位器、第四移位器和第二取反器，所述第三移位器连接有可执行加
‑
8操作的第三加法器，所述第三加法器与三选一选择器相连接，所述第四移位器连接有可执行加8操作的第四...

【专利技术属性】
技术研发人员：杜娟，王晓彤，李守亮，李浩然，
申请(专利权)人：兰州大学，
类型：新型
国别省市：