一种具有恒李亚普诺夫指数谱的混沌电路的制作方法

本发明涉及一种混沌电路，特别涉及一种具有恒李亚普诺夫指数谱的混沌电路。

背景技术：

1963年，混沌之父e.n.lorenz发现了第一个混沌系统，它是大气对流问题的简化模型，并成为了后人研究混沌理论的出发点和基石。之后，s.smale提出了21世纪的18个著名的数学问题，其中第14个问题就是关于lorenz系统的研究，可见lorenz系统的科学意义和研究价值是非常重要的。随着混沌理论研究的发展，混沌的应用也在各领域迅速扩大，目前，混沌科学的应用研究已经发展到如何有效地利用混沌，动力系统的混沌性在通信技术和信号处理方面的应用研究更是如火如荼，其中混沌保密通信和混沌电路研究就是这样一个方兴未艾的前沿领域。显而易见，混沌应用离不开混沌系统的设计，而对于混沌保密通信而言，混沌电路的设计是混沌应用的先决条件。

电路系统中对于混沌现象的研究一直独树一帜，非线性动态电路是动力学系统的一个分支，利用电路来产生、处理复杂的混沌信号就成了非线性科学研究中的一个热点。由于混沌电路与其对应的数学模型具有很好的吻合性，使混沌电路能够方便地模拟各种非线性混沌系统，并能重现各种复杂的非线性现象，因此在混沌的理论探索和应用研究中非线性电路充当着一个非常重要的角色。非线性电路理论为非线性元器件的应用提供了理论基础，应用非线性元器件能构造出产生人们所需要的混沌信号的电路。从电路的设计这个角度来考虑，应该在混沌理论分析的基础上，将状态方程中的各种理论参数数值通过某种对应关系的变换，如：比例变换、微分-积分转换，加法减法变换等，最后将理论参数值转换成实际所需要的电路参数值。这种用理论指导电路设计的方法，是证实用电路实验产生混沌吸引子具有可行性的关键技术所在。利用这种方法获得的电路参数具有较高的准确性，可进一步用于指导硬件电路的设计与实验。而现在的混沌信号产生电路均不能实现将状态方程中的各种理论参数数值通过变换转换成实际所需要的电路参数值。

技术实现要素：

为了解决上述技术问题，本发明提供一种结构简单、成本低的具有恒李亚普诺夫指数谱的混沌电路。

本发明解决上述问题的技术方案是：一种具有恒李亚普诺夫指数谱的混沌电路，包括信号源模块、线性模块、分段线性函数模块和负电阻模块，信号源模块的输出端与线性模块、负电阻模块相连，分段线性函数模块的输出端与线性模块相连，负电阻模块与线性模块相连。

上述具有恒李亚普诺夫指数谱的混沌电路，所述分段线性函数模块包括第一至第五运算放大器、第一至第十三电阻，所述第一运算放大器的反相输入端经第一电阻后接地，第二电阻跨接于第一运算放大器的反相输入端和输出端之间，第一运算放大器的同相输入端与第四运算放大器的输出端相连，第四运算放大器的反相输入端与第四运算放大器的输出端相连，第四运算放大器的同相输入端经第五电阻后与第三运算放大器的输出端相连，第三运算放大器的同相输入端接地，第三运算放大器的反相输入端经第三电阻后与第一运算放大器的输出端相连，所述第四电阻跨接于第三运算放大器的反相输入端和输出端之间，所述第二运算放大器的同相输入端接地，第二运算放大器的反相输入端经第七电阻后与第一运算放大器的同相输入端相连，第二运算放大器的输出端经第九电阻、第十电阻后与第三运算放大器的反相输入端相连，所述第八电阻跨接于第二运算放大器的反相输入端与输出端之间，所述第十一电阻的一端接地，另一端连接在第九电阻与第十电阻之间，所述第五运算放大器的同相输入端与第四运算放大器的同相输入端相连，第五运算放大器的反相输入端经第六电阻后接地，所述第十二电阻跨接于第五运算放大器的反相输入端与输出端之间，所述第十三电阻跨接于第五运算放大器的同相输入端与输出端之间。

上述具有恒李亚普诺夫指数谱的混沌电路，所述线性模块包括一个电感和一个电容，信号源模块的输出端与电容的一端相连，电容的另一端分别与第五运算放大器的同相输入端、电感的一端相连，电感的另一端与负电阻模块相连。

上述具有恒李亚普诺夫指数谱的混沌电路，所述负电阻模块包括第六运算放大器、第十四电阻、第十五电阻和第十六电阻，所述第六运算放大器的反相输入端与电感的另一端相连，第六运算放大器的同相输入端经第十六电阻后与信号源模块的输出端相连，所述第十四电阻跨接于第六运算放大器的反相输入端与输出端之间，所述第十五电阻跨接于第六运算放大器的同相输入端与输出端之间。

上述具有恒李亚普诺夫指数谱的混沌电路，所述第一电阻、第七电阻的阻值为10kω，第二电阻的阻值为58.2kω，第三电阻、第四电阻、第十电阻的阻值为100kω，第五电阻、第六电阻、第十二电阻、第十三电阻的阻值为0.4kω，第八电阻的阻值为130kω，第九电阻、第十一电阻、第十五电阻、第十六电阻的阻值为1kω，第十四电阻的阻值为33ω。

所述混沌电路方程如下：

α＝0.9,β＝0.11；

g(x-sint)＝m1(x-sint)+0.5(m0-m1)[|x-sint+x1|-|x-sint-x1|]；

m0＝-0.33,m1＝4.9,x1＝1.1。

本发明的有益效果在于：本发明包括线性模块、分段线性函数模块和负电阻模块，通过在线性模块中加入非线性元器件，将状态方程中的各种理论参数数值通过变换转换成实际所需要的电路参数值，获得的电路参数的准确性高，证实了用电路实验产生混沌吸引子具有可行性，整个电路具有结构简单、成本低的优点；并且本发明对混沌电路的状态方程进行了各种理论研究，得到了李亚普诺夫指数、分岔图、混沌相图等，证实了能产生一种新的具有恒李亚普诺夫指数谱的混沌电路。

附图说明

图1为本发明的结构框图。

图2为本发明的电路图。

图3为本发明的simulink系统仿真示意图。

图4为本发明的simulink系统仿真相图。

图5为本发明的x的时间序列波形图。

图6为本发明的李亚普诺夫指数谱图。

图7为本发明的分岔图。

图8为本发明的采用pspice16.3仿真的结果图。

图9为本发明的实物x的时间序列波形图。

图10为本发明的实物混沌吸引子相图。

图11为本发明的实物混沌信号频谱图。

具体实施方式

下面结合附图和实施例对本发明作进一步的说明。

如图1所示，一种具有恒李亚普诺夫指数谱的混沌电路，包括信号源模块1、线性模块、分段线性函数模块2和负电阻模块3，信号源模块1的输出端与线性模块、负电阻模块3相连，分段线性函数模块2的输出端与线性模块相连，负电阻模块3与线性模块相连。

如图2所示，所述分段线性函数模块2包括第一至第五运算放大器、第一至第十三电阻，第一运算放大器u1的反相输入端经第一电阻r1后接地，第二电阻r2跨接于第一运算放大器u1的反相输入端和输出端之间，第一运算放大器u1的同相输入端与第四运算放大器u4的输出端相连，第四运算放大器u4的反相输入端与第四运算放大器u4的输出端相连，第四运算放大器u4的同相输入端经第五电阻r5后与第三运算放大器u3的输出端相连，第三运算放大器u3的同相输入端接地，第三运算放大器u3的反相输入端经第三电阻r3后与第一运算放大器u1的输出端相连，所述第四电阻r4跨接于第三运算放大器u3的反相输入端和输出端之间，所述第二运算放大器u2的同相输入端接地，第二运算放大器u2的反相输入端经第七电阻r7后与第一运算放大器u1的同相输入端相连，第二运算放大器u2的输出端经第九电阻r9、第十电阻r10后与第三运算放大器u3的反相输入端相连，所述第八电阻r8跨接于第二运算放大器u2的反相输入端与输出端之间，所述第十一电阻r11的一端接地，另一端连接在第九电阻r9与第十电阻r10之间，所述第五运算放大器u5的同相输入端与第四运算放大器u4的同相输入端相连，第五运算放大器u5的反相输入端经第六电阻r6后接地，所述第十二电阻r12跨接于第五运算放大器u5的反相输入端与输出端之间，所述第十三电阻r13跨接于第五运算放大器u5的同相输入端与输出端之间。

取n＝2，分段线性函数表达式如下：

g(x-sint)＝m1(x-sint)+0.5(m0-m1)[|x-sint+x1|-|x-sint-x1|]；其中参数为：m0＝-0.33,m1＝4.9,x1＝1.1，x-sint为x与正弦信号之差。

为便于电路实现，分段线性函数表达式为：

rf(vc2-vc1)＝g1(vc2-vc1)+0.5(g0-g1)(|vc2-vc1+e1|-|vc2-vc1-e1|)；其中g0＝-0.33ms,g1＝4.9ms,e1＝1.1,e1为转折点电压值。

分段线性函数模块2的各电阻参数的理论计算值为：

r为负电阻的阻值28ω，r1、r2、r7、r8、r9、r11分别表示第一、第二、第七、第八、第九、第十一电阻r11的阻值。

如图2所示，所述线性模块包括一个电感l和一个电容c，信号源模块1的输出端与电容c的一端相连，电容c的另一端分别与第五运算放大器u5的同相输入端、电感l的一端相连，电感l的另一端与负电阻模块3相连。

如图2所示，所述负电阻模块3包括第六运算放大器u6、第十四电阻r14、第十五电阻r15和第十六电阻r16，所述第六运算放大器u6的反相输入端与电感l的另一端相连，第六运算放大器u6的同相输入端经第十六电阻r16后与信号源模块1的输出端相连，所述第十四电阻r14跨接于第六运算放大器u6的反相输入端与输出端之间，所述第十五电阻r15跨接于第六运算放大器u6的同相输入端与输出端之间。

所述第一电阻r1、第七电阻r7的阻值为10kω，第二电阻r2的阻值为58.2kω，第三电阻r3、第四电阻r4、第十电阻r10的阻值为100kω，第五电阻r5、第六电阻r6、第十二电阻r12、第十三电阻r13的阻值为0.4kω，第八电阻r8的阻值为130kω，第九电阻r9、第十一电阻r11、第十五电阻r15、第十六电阻r16的阻值为1kω，第十四电阻r14的阻值为33ω。

混沌电路方程如下：

α＝0.9,β＝0.11；

g(x-sint)＝m1(x-sint)+0.5(m0-m1)[|x-sint+x1|-|x-sint-x1|]；

m0＝-0.33,m1＝4.9,x1＝1.1。

如图3-图5所示，此时，可设李亚普诺夫指数工具箱内参数分别为：

finaltime:1000,step:0.01,updatalyapunov:10,initialconditions:000,no.oflinearizedodes:9得电路的李亚普诺夫指数为:0.084267,0，-4.3843(2.0192)，由此可知该电路是混沌状态的。如图6所示，从李亚普诺夫指数谱来看，该系统具有恒李亚普诺夫指数谱的特性。因为该电路的时间t与系统没有关系，这是由于在计算时把时间t的微分等于1。所以随着时间的变化李亚普诺夫指数谱是恒定的。

参数变化的分岔分析：

如图7所示，当固定参数β＝0.11m0＝-0.33,m1＝4.9,x1＝1.1时，α∈(0,1)之间时系统的分岔图如图7中(a)所示；当固定参数α＝0.9m0＝-0.33,m1＝4.9,x1＝1.1时，β∈(0,1)之间时系统的分岔图如图7中(b)所示；当固定参数α＝0.9,β＝0.11m0＝-0.33,m1＝4.9,x1＝1.1时，正弦函数的频率f∈(0,1)之间时系统的分岔图如图7中(c)所示。

如图8所示，pspice16.3仿真和matlab的仿真结果一致，图8中(a)为pspice仿真时x的时间序列图，(b)为pspice仿真时的相图。

如图9、图10、图11所示，根据原理图进行实验，选取图中所有运放为tl082cp，采用双电源±15v，函数信号发生器为tfg6060，输入正弦波幅值为3v，频率为3500hz，示波器采用泰克的tds1002b。得到结论：实物实验结果与理论仿真结果一致，说明了该混沌信号的可实现性和存在性。