[发明专利]一种基于数据时间序列的混沌检测方法在审
申请号: | 201810022002.9 | 申请日: | 2018-01-10 |
公开(公告)号: | CN108197712A | 公开(公告)日: | 2018-06-22 |
发明(设计)人: | 徐建新;朱秀乐;肖清泰;高勤 | 申请(专利权)人: | 昆明理工大学 |
主分类号: | G06N7/08 | 分类号: | G06N7/08 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 650093 云*** | 国省代码: | 云南;53 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及一种基于数据时间序列的混沌检测方法,属于基于数据时间序列检测的技术领域。本发明根据混沌时间序列的特点,将时间序列进行分组,分成Q(1≤Q≤N)个子集,其中N为时间序列的长度,对每个子集进行升序排序,并保留其相对应的索引序列,并求解每个子集的最大斜率,形成一个最大斜率序列。针对于斜率序列,将最大斜率序列在进行分组,划分为P(1≤P≤Q),通过统计学方法,对时间序列的最大斜率的行为进行描述,提供了一种新型的混沌检测方法,该方法首先根据渐进生长率K是否趋向于0或1从而检测出规则动力学和非规则动力学,再根据周期性指数λ与0的关系来检测出出周期和准周期序列。 | ||
搜索关键词: | 时间序列 最大斜率 检测 混沌 动力学 子集 混沌时间序列 索引序列 分组 统计学 非规则 准周期 渐进 求解 升序 排序 保留 | ||
【主权项】:
1.一种基于数据时间序列的混沌检测方法,其特征在于:包括如下步骤:STEP1.设接收到长度为N的离散时间序列x;STEP2.将STEP1中得到时间序列x进行分组,分成Q个子集xi,其每个子集的长为ni,使得每个子集包含于下一个子集,即
其中i=1,2,...Q,n1=n,n的数值为最大斜率的最小观察持续时间,且n<N,因此可以求得每个子集的长度为ni=(i‑1)p0+n,其中p0为积分步长,且p0<N‑n;STEP3.将时间序列x的每个子集序列xi进行升序排序,导出其索引序列Ii,j,该索引为排序后的子集序列对应于原子集序列中元素的位置索引,并将索引中元素相邻之间的最大差异为时间序列的最大斜率Si,即Si=max{Ii,j‑Ii,j+1},其中i=1,2,...,Q;j=1,2,...,nj,S是由Si组成的最大斜率序列;STEP4.由STEP3中计算出的最大斜率序列S,通过统计学方面标准差的知识对时间序列的特征进行研究,从而计算其最大斜率序列S的标准差,即
根据标准差从而推导出生长速率,即
最后由生长速率推导得到渐进生长率K(n),即
STEP5.由上述STEP4得出的渐进生长速率可以判断出K(n)是一个非负的指标,由此可以区分出来规则动力学和非规则动力学:K(n)→1时是非规则动力学;K(n)趋向于0或的等于0时是规则动力学;STEP6.在上述STEP5区分出的规则动力学中,区分周期和准周期,根据STEP3步骤将S划分为P个子集Vk,其子集长度为mk,使得
其中k=1,2,...P,且m1>2;STEP7.用指数拟合来定义
其中μ0、K0为常数,且μ0=‑K0,uk=mk/mp是归一化积分时间;STEP8.用STEP7中等式推导得出判断指标λ(n):
其中Γ(nk,n)=‑log(1‑μ(nk,n)),
分别为序列Γ,u的平均数;STEP9.根据λ(n)的值,得到最终的结果:当λ(n)=0,时间序列为周期;当λ(n)>0,时间序列为混沌的;当λ(n)<0,时间序列为准周期;在时间序列为周期和准周期的情况下,周期为T=SQ,SQ表示最大斜率序列S的第Q个数据。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于昆明理工大学,未经昆明理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201810022002.9/,转载请声明来源钻瓜专利网。