1。介绍
智能电网的定义是提供传统的电力网络和现代信息技术。对于智能电网,最基本的目标可以被解释为经济和清洁能源的成就一代在全球范围内(
1,
2]。传统电网目前正在配备智能设备和系统,一般来说,收集设备,分配控制器单元、数据管理、电力市场监控系统,远程终端单元,断路器,人机接口设备、网络管理工具和路由器,网络集中器、相量测量单元,智能电表,和保护继电器
3]。这些设备功能风险决策操作,比如能源截止,调试,和甩负荷中扮演至关重要的角色完成所需的操作以及网格中的监视和测量活动。这类设备的安装使电力系统容易受到外部网络威胁。如果双向通信阻塞由于网络攻击,可能会有不可预见的损害在网格中。因此,开发安全工具对网络威胁是必不可少的。在网格对网络威胁是一个战略预防措施来防止不可预见的损害时,双向沟通受阻。
智能电网的安全概念应该被认为在一个广泛的网格组成的各种类型的通信系统,如GSM网络,光纤线,WiMAX, rs - 232和rs - 485,电缆线路,射频线,电力系统通信线路(
4]。戒备森严的硬件和计算基础设施的发展,以确保所有这些硬件和系统工作可靠的智能电网内分别操作。
通过使用各种和这些复杂的基础设施,synchrophasor技术用于监视和控制电力系统实时状态通过synchrophasor跨大地理区域实时测量非常低的延迟。这些测量是由IEEE C37.118和IEC 61850 synchrophasor通信框架的安全是至关重要的任何错误的信息可能导致严重损害到物理设备。然而,几乎所有这些沟通框架没有内置的安全机制,并且限制通信只有本地网络。在这些通信框架,我们认为IEC 61850是一个完整的通信系统,解决电力系统组件的建模,抽象的服务,和通信协议
5),但只有IEC 61850-90-5通信框架使用内置的加密算法和密钥的定期刷新,这可以很容易地检测未经授权的修改(
6]。相反,IEC 61850-90-5通信框架有三次大的数据包大小由于元数据和在每个数据包携带完整的解码信息相比,IEEE C37.118.1和IEEE C37.118.2框架。
今天,各种先进的方法开始被用于提高电网的安全级别。网格安全的一个重要应用是连接每个网格组件在一个软件定义网络(SDN),其配置指定根据电力系统组件的漏洞来攻击(
7]。对称和非对称加密的隐私和身份验证方法也是目前使用在电网自动化系统(PLC)和RTU单位
8]。另一个网格中的安全应用程序是基于网格的控制组件,特别是控制pmu的数据包,对破坏数据由外部或内部攻击或故障原因用卡尔曼滤波方法(
9]。抵御网络攻击也通过使用安全VPN技术(
10]。所有这些方法实现非常复杂和紧密联系的规则来满足指定的安全目标。做所有这些都需要高计算成本,高内存和功耗。这可以防止这些安全系统和使他们的传播仍在有限数量的设备和区域。在这种背景下,混沌流密码可以被视为一种块密码关于低复杂度和低资源消耗(
11]。
混沌和混沌系统的问题是最复杂的动态行为的非线性系统。他们的科学领域,有助于解释非线性系统。Chaos-based工程应用,导致出现显著改善控制(
12,
13)、沟通(
14)、人工智能(
15,
16)和遗传算法(
17]近年来区域。他们也被用来作为随机数生成器(
18- - - - - -
20.)在加密等网络安全应用程序(
21- - - - - -
23和数据隐藏
24]。
在这项工作中,小说提出了微机与超混沌系统加密方法允许PMU设备相互通信在智能电网安全的方式。最好根据我们所知,这是第一次使用微机加密在智能电网环境下。提出的安全方法是基于加密引擎确保一个安全的通信层PMU和相量数据集中器(PDC)组件。为此,应用程序服务器软件,这与IEEE C37.118.1兼容PMU装置通过Wireshark实时平台,也是发达。评估结果数据交换的完整性,密码学的加密-解密的持续时间和强度。结果表明,快速和复杂的加密功能提出了微机加密方法可以很容易地用于确保安全通信在智能电网组件非常狭窄的pmu等时间间隔。
本文的组织结构如下:在部分
2介绍了现代电力系统分析方法。节
3,智能电网的安全方法。节
4非线性系统中使用该方法及其动力学分析。部分
5致力于随机数发生器的设计,提出了统计测试的部分
6。PMU数据加密的应用程序的细节和它的安全分析也提供了在这一节中。最后,本文的结论部分
7。
年代ec><年代ec id="sec2">
2。现代电力系统分析的方法
电力系统的协调和可靠性取决于建立一个和谐的力量和信息基础设施创建一个实时双向通信网络一代,传输和分配电能。双向通信技术,控制系统,智能电网和计算机处理是关键。先进的传感器、仪表、可编程继电器和自动化馈线开关的现代智能电网中使用的设备。电力系统健康运行取决于变电站操作数据的连续测量。在这种背景下,可以说,pmu被广泛认为是最有前途的一个领域的发展电力系统的实时监控。pmu的能力来计算同步相量电压和电流瞬间和更准确地鼓励他们一致的扩散在电力系统网络遍布世界各地。
<年代ecid="sec2.1">
2.1。相量计算和相量测量单元的原则
2.1.1。相量计算
相量计算电源总线在纯正弦信号给出
(1)米米l:mtext>
x米米l:mi>
t米米l:mi>
=米米l:mo>
X米米l:mi>
米米米l:mi>
因为米米l:mi>
ω米米l:mi>
t米米l:mi>
+米米l:mo>
φ米米l:mi>
,米米l:mo>
在哪里
ω米米l:mi>
和
φ米米l:mi>
的像频率信号弧度/秒和相位角弧度,分别
X米米l:mi>
米米米l:mi>
是信号的峰值振幅。图
1代表所有的部分纯正弦信号。
一个正弦信号(A)和其代表作为相量(b)。矢量的相位角是任意的,因为它取决于选择的轴
t米米l:mi>
=米米l:mo>
0米米l:mn>
。
相量表示只是纯正弦信号。在实践中,波形通常腐败与其他不同频率的信号。提取一个频率分量的信号,使用傅里叶变换。采样数据的提取、离散(DFT)或快速变换(FFT) (
25]。一个
x米米l:mi>
t米米l:mi>
正弦信号频率
k米米l:mi>
f米米l:mi>
0米米l:mn>
表示为傅里叶级数显示如下:
(2)米米l:mtext>
x米米l:mi>
t米米l:mi>
=米米l:mo>
一个米米l:mi>
k米米l:mi>
因为米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
f米米l:mi>
0米米l:mn>
t米米l:mi>
+米米l:mo>
b米米l:mi>
k米米l:mi>
罪米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
f米米l:mi>
0米米l:mn>
t米米l:mi>
=米米l:mo>
一个米米l:mi>
k米米l:mi>
2米米l:mn>
+米米l:mo>
b米米l:mi>
k米米l:mi>
2米米l:mn>
因为米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
f米米l:mi>
0米米l:mn>
t米米l:mi>
+米米l:mo>
φ米米l:mi>
在哪里米米l:mtext>
φ米米l:mi>
=米米l:mo>
反正切米米l:mi>
−米米l:mo>
b米米l:mi>
k米米l:mi>
一个米米l:mi>
k米米l:mi>
。米米l:mo>
在这里,傅里叶级数系数
一个米米l:mi>
k米米l:mi>
和
b米米l:mi>
k米米l:mi>
测量
因为米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
f米米l:mi>
0米米l:mn>
t米米l:mi>
和
罪米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
f米米l:mi>
0米米l:mn>
t米米l:mi>
出现在函数
x米米l:mi>
t米米l:mi>
。相量表示信号变成的
(3)米米l:mtext>
X米米l:mi>
k米米l:mi>
=米米l:mo>
1米米l:mn>
2米米l:mn>
一个米米l:mi>
k米米l:mi>
2米米l:mn>
+米米l:mo>
b米米l:mi>
k米米l:mi>
2米米l:mn>
e米米l:mi>
j米米l:mi>
φ米米l:mi>
。米米l:mo>
的复杂形式的方程(
3)是构成如下:
(4)米米l:mtext>
X米米l:mi>
k米米l:mi>
=米米l:mo>
1米米l:mn>
2米米l:mn>
一个米米l:mi>
k米米l:mi>
−米米l:mo>
j米米l:mi>
b米米l:mi>
k米米l:mi>
。米米l:mo>
使用傅里叶级数系数的关系和DFT,的相量表示
k米米l:mi>
th米米l:mtext>
谐波分量是由
(5)米米l:mtext>
X米米l:mi>
k米米l:mi>
=米米l:mo>
1米米l:mn>
2米米l:mn>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
n米米l:mi>
Δ米米l:mi>
T米米l:mi>
e米米l:mi>
−米米l:mo>
j米米l:mi>
2米米l:mn>
k米米l:mi>
n米米l:mi>
/米米l:mo>
N米米l:mi>
,米米l:mo>
(6)米米l:mtext>
X米米l:mi>
k米米l:mi>
=米米l:mo>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
n米米l:mi>
Δ米米l:mi>
T米米l:mi>
因为米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
n米米l:mi>
N米米l:mi>
−米米l:mo>
j米米l:mi>
罪米米l:mi>
2米米l:mn>
π米米l:mi>
k米米l:mi>
n米米l:mi>
N米米l:mi>
,米米l:mo>
x米米l:mi>
n米米l:mi>
Δ米米l:mi>
T米米l:mi>
=米米l:mo>
x米米l:mi>
n米米l:mi>
和米米l:mtext>
θ米米l:mi>
=米米l:mo>
2米米l:mn>
π米米l:mi>
N米米l:mi>
,米米l:mo>
在哪里
N是样品的数量在数据窗口中,
n样品数量,
x米米l:mi>
n米米l:mi>
是输入样本,
θ米米l:mi>
是指样本角度。然后,方程(
6)成为如下:
(7)米米l:mtext>
X米米l:mi>
k米米l:mi>
=米米l:mo>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
n米米l:mi>
因为米米l:mi>
k米米l:mi>
n米米l:mi>
θ米米l:mi>
−米米l:mo>
j米米l:mi>
罪米米l:mi>
k米米l:mi>
n米米l:mi>
θ米米l:mi>
。米米l:mo>
罪恶和余弦的部分
x米米l:mi>
n米米l:mi>
分离,分别吗
X米米l:mi>
k米米l:mi>
c米米l:mi>
和
X米米l:mi>
k米米l:mi>
年代米米l:mi>
的形式
(8)米米l:mtext>
X米米l:mi>
k米米l:mi>
=米米l:mo>
X米米l:mi>
k米米l:mi>
c米米l:mi>
−米米l:mo>
j米米l:mi>
X米米l:mi>
k米米l:mi>
年代米米l:mi>
,米米l:mo>
X米米l:mi>
k米米l:mi>
c米米l:mi>
=米米l:mo>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
n米米l:mi>
因为米米l:mi>
k米米l:mi>
n米米l:mi>
θ米米l:mi>
,米米l:mo>
X米米l:mi>
k米米l:mi>
年代米米l:mi>
=米米l:mo>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
n米米l:mi>
罪米米l:mi>
k米米l:mi>
n米米l:mi>
θ米米l:mi>
。米米l:mo>
最后,
X米米l:mi>
k米米l:mi>
用于表示大部分的相量的相量计算。在相量计算过程中,需要用新数据更新矢量估计。这是最简单的方式连续相量计算和称为nonrecursive相量更新方法。
(9)米米l:mtext>
X米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
=米米l:mo>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
n米米l:mi>
+米米l:mo>
1米米l:mn>
+米米l:mo>
r米米l:mi>
因为米米l:mi>
n米米l:mi>
θ米米l:mi>
−米米l:mo>
j米米l:mi>
罪米米l:mi>
n米米l:mi>
θ米米l:mi>
。米米l:mo>
在方程(
9),
r=−1,1,2,3,…;当
r=−1,
x米米l:mi>
0米米l:mn>
样品现在在右边,但当
r= 0,没有
x米米l:mi>
0米米l:mn>
样即使样本总数,也就是说,
N,是相同的。
见图
2相量计算1与样品
n= 0…
N2−1,相量计算样本
n= 1,2,…
N。
θ米米l:mi>
连续的样本之间的角度是基于周期的基本频率。在这里,
n是样本数量,
N在数据窗口中样本的数量。每个窗口的相量计算执行新的估计早不使用任何数据;这个算法是最简单的方法,称为“nonrecursive算法”。
(10)米米l:mtext>
X米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
=米米l:mo>
X米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
−米米l:mo>
1米米l:mn>
+米米l:mo>
2米米l:mn>
N米米l:mi>
∑米米l:mo>
n米米l:mi>
=米米l:mo>
0米米l:mn>
N米米l:mi>
−米米l:mo>
1米米l:mn>
x米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
−米米l:mo>
x米米l:mi>
r米米l:mi>
因为米米l:mi>
n米米l:mi>
θ米米l:mi>
−米米l:mo>
j米米l:mi>
罪米米l:mi>
n米米l:mi>
θ米米l:mi>
。米米l:mo>
Nonrecursive更新矢量的估计
N样本窗口。
在方程(
10),
r= 0、1、2、3;当
r表示目前的状态,(
r−1)表示过去的状态。如方程(
10),递归估计当前的输出,
X米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
,取决于前面的输出,
X米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
−米米l:mo>
1米米l:mn>
,当前的输入,
x米米l:mi>
N米米l:mi>
+米米l:mo>
r米米l:mi>
。递归相量估计是快而nonrecursive相量估计相量计算不是表现在每一个步骤。如果正弦波不是连续的,有一个小错误矢量更新。然后,在递归矢量估计,这个小错误积累导致更大的错误。
年代ec><年代ec id="sec2.1.2">
2.1.2。相量测量单元的原则
pmu安装了全球定位系统(GPS)接收器需要非常精确的测量电压和电流在定义的时间戳。这里,GPS接收器协调同步的电流和电压测量(
26]。pmu使用递归算法计算对称组件的电压和电流,如对称分量距离继电器(SCDR) [
27]。
全球定位系统(
28)是一个36卫星系统生产时间信号在地球表面。每次24卫星是用于生成时间信号。GPS接收机测量距离用无线电信号的传播时间。一微秒级精度是一个合适的范围测量工频电压和电流(
27]。
由于pmu的引入电力系统,电力系统设计师一直在寻找利用pmu的能力为了观察系统在一个更好的方法。今天,收集和合并信息的pmu是最好的输电线路故障检测方法(
28]。pmu的处理数据是用来改善稳定协调(
29日),电力系统状态估计(
30.],补救行动计划(
31日],和扰动监测[
32]。相反,线性化主要发生在状态向量,以及PMU测量,在直角坐标系中表示。
年代ec><年代ec id="sec2.1.3">
2.1.3。IEEE-14总线模型和用于PMU的放置算法
它可以表示,最重要的条件,提高电力系统的可靠性和稳定性取决于建立一个连续和精确的测量系统。基尔霍夫电流和电压法的基本原理是基于电气相量量。同样,可以假定电信号以光速移动,系统的电压和电流瞬时值有可能通过结合所有的电压和电流相量测量用一个精确的计时器。PMU装置,这是定制的相量测量装置在电力系统中,可以检测出电压和电流相量信息网格。PMU装置的最佳位置应考虑由于需要一个伟大的经济成本的安装PMU在每一个电力系统母线(
33,
34]。而不是每个母线安装pmu的电力系统,高动力和强烈的可变性的汇流条件选择安装。在这种背景下,首先,PMU装置的放置在IEEE-14公交系统已经建立的深度优先搜索方法(
35- - - - - -
38]。在这项研究中,一个IEEE-14母线系统,最优PMU位置安装在母线系统进行了研究。然后,提出了一种微机与超混沌系统加密方法作为一种新的安全方法PMU装置之一。
开源PSAT功率流分析工具用于安全通信的位置pmu的最合适的数量IEEE-14母线系统。根据深度优先搜索方法,PMU装置放置的数量是6,和最合适的汇流PMU放置汇流
1、4、6、8、10米米l:mn>
,米米l:mo>
和米米l:mtext>
14米米l:mn>
,分别。图
3显示了pmu的位置IEEE-14母线供电系统。
IEEE-14母线PMU内置模型与深度优先搜索算法。
由于这些设备的使用,它是确定16可以测量电流相量研究电力系统由14汇流和20行。列表可测相量的PMU-installed电力系统表
1,分别。
可测相量的PMU-installed IEEE-14母线系统。
| 1 - 2线 | 4 - 5行 | 6尺11寸线 | 9线 |
| 1 - 5行 | 4 - 7行 | 6 - 12行 | 9-14线 |
| 2 - 4线 | 4 - 9日线 | 6-13线 | - 11行 |
| 3 - 4线 | 5 - 6行 | 7 - 8线 | 13 - 14日线 |
同样,另一种IEEE测试系统可用于案例研究。在这样一个测试场景中,最重要的区别就是确定PMU的数量允许完整的监测系统,定义适当的PMU和PDC数据交换过程中的IEEE C37.118.1协议(
39]。例如,在IEEE-39, ieee - 118和ieee - 300是用作测试系统,只有将创建一个不同的数据交换方案放置14,39岁,和73年PMU装置,分别,这将提供完整的监测研究测试系统使用一些优化布置算法[
40]。图
4显示了数据包的结构符合IEEE C37.118.1 PMU装置。数据包结构包括设备ID、GPS时间戳,第二部分,用户定义的标记命令发送数据包内的通信协议标准,相量数据,数据,和一个16位CRC数据完整性控制器。相量数据由频率(
f)、电压大小(
米米米l:mi>
_米米l:mo>
V米米l:mi>
)、电压相角(
δ米米l:mi>
_米米l:mo>
V米米l:mi>
)、电流强度(
米米米l:mi>
_米米l:mo>
我米米l:mi>
)和电流相角(
δ米米l:mi>
_米米l:mo>
我米米l:mi>
母线的)。
IEEE C37.118.1兼容PMU包定义。
3所示。一个新颖的智能电网的安全方法
得到一个正确的相量计算研究IEEE-14母线系统取决于适当放置的PMU装置与足够的数量。一个正确的相量计算还取决于PMU的有效和适当的数据通信设备。在这项研究中,而不是获得时间和包延迟完成一次完整的相量计算,我们有,尤其集中在基于流的通信用于揭示加密一个PMU PDC设备之间数据通信的成功。这里,内置的PMU装置在IEEE-14母线系统通信利用IEEE C37.118.1同步相量与TCP / IP通信协议的支持。PDC单位负责收集所有pmu的数据网格使用IEEE C37.118.1协议。在这项研究中,Wireshark平台只负责接口数据存储库的PMU装置通过一个预先确定的IP和套接字地址。基于客户机的应用程序在PMU两侧和基于服务器的应用程序软件在PDC端一起处理了PMU的基于超混沌加密过程数据。
见图
5,Rasp3单位经营的PMU装置进行加密数字矢量分组数据在和谐与IEEE C37.118.1协议使用混沌加密引擎。而连接还是PMU IEEE-14总线系统上的单位除了Rasp3, PDC服务器执行微机加密使用Wireshark通过收集数据平台。加密和解密过程是由客户端和服务器使用提出了微机加密引擎。
整体PMU-based混沌加密方案。
见图
6,研究了实时数据采集PMU装置是通过IEEE C37.118.1兼容Wireshark的平台。应用服务器管理PMU和PDC双方与混沌加密引擎。
研究IEEE C37.118.1帧被Wireshark网络协议分析仪的软件。
PDC服务器负责确定合适的初始条件安全通信的混沌系统和其他参数。的沟通是进行会话1 Mbit随机数生成的混沌流密码系统使用。促进安全通信,PDC,首先生成随机数使用5 d超混沌系统和测试他们对nist - 800 - 22为随机性统计测试。然后,它将这些参数发送到PMU作为加密密钥的一部分。PMU使用参数来生成相同的序列使用相同的混沌系统。然后,PMU数据加密与密文流通过随机数进行测试。
年代ec><年代ec id="sec4">
4所示。利用非线性系统及其动态分析
在这项研究中,使用5 d超混沌洛伦兹系统(
41- - - - - -
43]。5 d超混沌系统及其动力学分析给出了在这一节中。平衡点、李雅普诺夫和分岔分析检查进行动力分析。同时,相图所示部分
5。5 d超混沌洛伦兹系统设置如下:
(11)米米l:mtext>
x米米l:mi>
˙米米l:mo>
=米米l:mo>
−米米l:mo>
σ米米l:mi>
x米米l:mi>
−米米l:mo>
y米米l:mi>
+米米l:mo>
w米米l:mi>
,米米l:mo>
y米米l:mi>
˙米米l:mo>
=米米l:mo>
r米米l:mi>
x米米l:mi>
−米米l:mo>
y米米l:mi>
−米米l:mo>
x米米l:mi>
z米米l:mi>
−米米l:mo>
w米米l:mi>
,米米l:mo>
z米米l:mi>
˙米米l:mo>
=米米l:mo>
−米米l:mo>
β米米l:mi>
z米米l:mi>
+米米l:mo>
x米米l:mi>
y米米l:mi>
,米米l:mo>
u米米l:mi>
˙米米l:mo>
=米米l:mo>
−米米l:mo>
x米米l:mi>
z米米l:mi>
+米米l:mo>
p米米l:mi>
w米米l:mi>
,米米l:mo>
w米米l:mi>
˙米米l:mo>
=米米l:mo>
问米米l:mi>
y米米l:mi>
,米米l:mo>
在哪里
w米米l:mi>
是第五个状态变量和
问是一个积极的实际参数。5 d系统(
11)有三个积极的,消极的,和一个零李雅普诺夫指数。
平衡的五维超混沌洛伦兹系统(
11)可以通过假设找到
x米米l:mi>
˙米米l:mo>
=米米l:mo>
0米米l:mn>
,
y米米l:mi>
˙米米l:mo>
=米米l:mo>
0米米l:mn>
,
z米米l:mi>
˙米米l:mo>
=米米l:mo>
0米米l:mn>
,
u米米l:mi>
˙米米l:mo>
=米米l:mo>
0米米l:mn>
,
w米米l:mi>
˙米米l:mo>
=米米l:mo>
0米米l:mn>
和解决以下方程:
(12)米米l:mtext>
−米米l:mo>
σ米米l:mi>
x米米l:mi>
−米米l:mo>
y米米l:mi>
+米米l:mo>
w米米l:mi>
=米米l:mo>
0米米l:mn>
r米米l:mi>
x米米l:mi>
−米米l:mo>
y米米l:mi>
−米米l:mo>
x米米l:mi>
z米米l:mi>
−米米l:mo>
w米米l:mi>
=米米l:mo>
0米米l:mn>
−米米l:mo>
β米米l:mi>
z米米l:mi>
+米米l:mo>
x米米l:mi>
y米米l:mi>
=米米l:mo>
0米米l:mn>
−米米l:mo>
x米米l:mi>
z米米l:mi>
+米米l:mo>
p米米l:mi>
w米米l:mi>
=米米l:mo>
0米米l:mn>
问米米l:mi>
y米米l:mi>
=米米l:mo>
0。米米l:mn>
。米米l:mo>
因此,这个系统只有零平衡点:
E米米l:mi>
0米米l:mn>
0,0,0,0,0米米l:mn>
。雅可比矩阵的线性化系统的超混沌系统在唯一的平衡
E米米l:mi>
0米米l:mn>
是
(13)米米l:mtext>
J米米l:mi>
=米米l:mo>
−米米l:mo>
σ米米l:mi>
σ米米l:mi>
0米米l:mn>
1米米l:mn>
0米米l:mn>
r米米l:mi>
−米米l:mo>
1米米l:mn>
0米米l:mn>
0米米l:mn>
−米米l:mo>
1米米l:mn>
0米米l:mn>
0米米l:mn>
−米米l:mo>
β米米l:mi>
0米米l:mn>
0米米l:mn>
0米米l:mn>
0米米l:mn>
0米米l:mn>
p米米l:mi>
0米米l:mn>
0米米l:mn>
问米米l:mi>
0米米l:mn>
0米米l:mn>
0米米l:mn>
。米米l:mo>
我们周围的动力学研究原来的平衡点
E0年代ub>当临界分岔参数的值作为
p= 2,
问= 8。因此,系统的特征方程变成如下:
(14)米米l:mtext>
λ米米l:mi>
−米米l:mo>
p米米l:mi>
λ米米l:mi>
+米米l:mo>
8米米l:mn>
3米米l:mn>
λ米米l:mi>
3米米l:mn>
+米米l:mo>
11米米l:mn>
λ米米l:mi>
2米米l:mn>
+米米l:mo>
问米米l:mi>
−米米l:mo>
270年米米l:mn>
λ米米l:mi>
+米米l:mo>
10米米l:mn>
问米米l:mi>
=米米l:mo>
0米米l:mn>
。米米l:mo>
通过使用狼算法(
44),5d超混沌系统的李雅普诺夫指数计算
0.4580米米l:mn>
,米米l:mo>
0.3371米米l:mn>
,米米l:mo>
0.0415米米l:mn>
,米米l:mo>
0.0000米米l:mn>
,米米l:mo>
和米米l:mtext>
12.5046米米l:mn>
。
的值
p
>米米l:mo>
0和
问
<米米l:mo>
2970年,唯一的平衡
E米米l:mi>
0米米l:mn>
由于的存在是一个不稳定的界定点特征值与积极的真实部分。
分岔图也检查来确定模拟的混乱。可视化的动态行为和收集5 d超混沌系统的解决方案通过使用MATLAB数值解算器和阶段的超混沌系统的肖像描绘在图
7显示了分岔图的
p−7.5和2.5之间的参数变化值通过保持其他参数固定的。5维超混沌洛伦兹系统的初始状态决定
x米米l:mi>
0米米l:mn>
=米米l:mo>
0米米l:mn>
,米米l:mo>
y米米l:mi>
0米米l:mn>
=米米l:mo>
−米米l:mo>
0.01米米l:mn>
,米米l:mo>
z米米l:mi>
0米米l:mn>
=米米l:mo>
9米米l:mn>
,米米l:mo>
u米米l:mi>
0米米l:mn>
=米米l:mo>
1米米l:mn>
,
w米米l:mi>
0米米l:mn>
=米米l:mo>
0米米l:mn>
。步长为0.005,混乱出现在两个地区参数
p米米l:mi>
小于
−米米l:mo>
1。5米米l:mn>
大于1。
分岔图
p米米l:mi>
参数之间的不同
−米米l:mo>
8米米l:mn>
和3。
在数据
8(一个)和
8 (b),还透露,系统分为超混沌的情况可以看出,它有着独特的特性在其李雅普诺夫指数。在这种背景下,分岔和李雅普诺夫指数图确认彼此。
5 d超混沌系统的李雅普诺夫指数
p参数之间的不同(一个)
−米米l:mo>
8米米l:mn>
3、(b)
−米米l:mo>
1。5米米l:mn>
和2.5。