基于勒贝格测度和wigner变换混沌判定方法.docx

基于勒贝格测度和wigner变换混沌判定方 【摘要】本文提出利用Wigner-Ville分布的良好 频聚性，将混沌吸引子特性充分的展现出来，并结合勒贝格 测度，对混沌状进行识别。避免了以往的混沌识别方法大多 数依赖相空间重构、需要大量的实时数据等缺点。最后通过 duffing系统验证了这种方法的有效性 【关键词】混沌判定;Wigner-Ville分布；勒贝格测度 0引言 目前混沌判定的方法大致分为定性和定量两类，定性 有：庞加莱截面法、频闪法、功率谱法、主分量分析、小波 变换法、自相关函数法等等；定量的分析的方法有李亚普诺 夫指数法、关联维数、K炳、替代数据法、参差诊断、Hust 指数、“搅拌诊断”等。其中交通流混沌现象的判断主要集 中在采用替代数据法、G-P算法及其改进的Lyapunov指数法、 Poincare截面法、功率谱法等方法。 这些方法存在很多问题，比如需要选取具有量大、时间 长的要求，不能满足交通系统短时性的要求，从而造像了交 通信号控制的滞后性。 Wigner-V订le分布具有良好的频聚性，在其频谱图上能 够充分的表现出信号的频率和能量的分布。因此在信号时频 分析领域得到了广泛的应用。混沌时间序列是由一个低维的 具有非线性和确定性的动态系统产生的外表像随机信号但 并非是随机信号的的时间序列。而Wigner-Ville良好的频 聚性能够将混沌时间序列的内部层次很好的展现出来，通过 勒贝格测度对时间序列的混沌性进行定量的分析。本文将通 过Duffing混沌动力学系统为例子，具体的阐明这种混沌判 定方法。 1 Duffing系统的动力学特性 Duffing振子在非线性动力学系统的研究中占有重要的 地位。Duffing振子方程形式简单，却能呈现出丰富的非线 性动态特性。这是由于Duffing方程等号右边加上外加强迫 项后，系统的本征频率与外加周期强迫项的频率相互作用的 结果。 本文选取Holmes型Duffing振子作为混沌识别的系统 模型，其一般形式为： +k■-x3+x5二 Y cos ( o t) (1) 设fd为系统的临界值(阀值)，fd/k的解析值为一个常 数，实验证明了阻尼比k的取值范围为0. 2-0.5,在这里选 取k=0.5,上式的状态方程形式如下： = 3 y ■= co [-ky+x?+x?+rcos(3 t) ] (2) 式中，Y为周期策动力的幅值，k为阻尼比，-xB+xH 为非线性恢复力。Yeos (t)为周期策动力，系统频率为 ?=lrad/so取阻尼比20.5。在k固定的情况下，系统状态 随Y的变化出现有规律变化:历经同宿轨道?分叉?混沌轨 道?临界周期轨道?大尺度周期轨道状态。其各个状态的时 域波形及相平面轨迹。分析析系统时域波形及相平面轨迹变 化可知： 当Y=0时，系统相平面鞍点为(0, 0),焦点为 (±1, 0)。点(X, ■)将最终停在两焦点之一。 当Y0时，系统表现出复杂的动力学形态，当 V较小时，相轨迹表现为Poincare映射意义下的吸引子， 相点围绕一个焦点或另一个焦点做周期震荡。当Y超过一 定阀值Y c时(Y c的大小可由Melnikov方法求出)。随着 Y的增大，系统经历同宿轨道?周期分叉直至混沌状态。这 一过程随着V的变化非常迅速。V在很长时间内，系统都 处于混沌状态。直到大于另一阀值Yd时，系统进入大尺度 周期状态。此时相轨即将焦点?鞍点团团围住，其对应的 Poincare映射亦为不动点。 2基于Wigner变换和勒贝格测度的混沌时间序列的识 别方法 本文选取-x?+x■为恢复力项的Duffing振子方程 ■ +kH-x3+x5= Y cos (st) (3) 通过调整策动力的幅值Y使得系统分别处于混沌状态 与大周期态采用四阶Runge-Kutta算法求出方程的解，对求 解出来的x (t)做Wigner变换记为W (t, w )o其幅频如 图lo (I)混沌态下的幅频图(II)周期态的幅频图 图1 Duffing系统混沌态和周期态的Wigner分布 图中横坐标为频率，纵坐标为幅值。从图2.1可以很明 显的看出混沌态与大周期态之间有很大的不同。混沌状态 下，其wigner分布整体区域“平整”，但它具有一定的规 律性，正如前文混沌性质所说的有一定的内在规律。而在大 周期状态下，其wigner分布图显示出一个明显的波峰，而 其余部分相对平整。这充分说明了此时duffing系统进入到 了一个统一的周期状态中。即大周期状态。但是这些频谱图 无法更确切的描述与判定这两种状态。这就需要利用勒贝格 测度对对其进行更详细的研究。 在这里我们取： L (kO) ={? 的勒贝格测度 I I |W(3, t) | |2k0} (4) 图2勒贝格测度 根据勒贝格测度互不相交的可列个区间长