Duffing的报告

1、关于 Duffing 方程的报告1 Duffing 方程的基本概念混沌系统对微弱信号具有极强的敏感性同时对噪声具有极大的抑制能力，它的这种 性质证明了混沌系统具有可应用于小信号检测的潜力，从检测过程中分析混沌运动发生 的间歇性。Duffing 方程是一个在混沌系统小信号检测中被广泛使用的一个典型的非线性 方程，即存在于噪声中的信号可以被 Duffing 振子通过从混沌运动状态到周期振荡状态 的改变测试出来。Duffing 方程是描述共振现象、调和振动、次调和振动、拟周期振动、概周期振动、 奇异吸引子和混沌现象 (或随机过程 )的简单数学模型。因此，在非线性振动理论中研究 Duffing 方程具

2、有重要的意义。它的标准形式为 : 其中， > 0为阻尼系数。 g(x) 是含有三次方项的非线性函数， f(x,t) 为一周期函数。 Duffing 方程通常作如下分类：1) 假设 g(x) 满足超线性条件 :则称 Duffing 方程是超线性的 ;2) 假设 g(x) 满足次线性条件 :则称 Duffing 方程是次线性的 ;3) 假设 g(x) 满足半线性条件 :则称 Duffing 方程式半线性的。若将 Duffing 方程规范化，有以下四种基本类型：其中，类型 1 为硬特性 Duffing 方程；类型 2 和 4成为软特性 Duffing 方程；类型 3 称 为日本型，日本学者上田

3、研究较多，并发现了日本吸引子，也称为Ueda吸引子；美国科学家对类型 4的Duffing 方程进行了深入的研究，因此类型 4也称为型 Duffing 方程。Duffing 方程系统是一个典型的非线性振动系统， 尽管是从简单物理模型中得出来的 非线性振动模型，但是其模型具有代表性。工程实际中的许多非线性振动问题的数学模 型都可以转化为该方程来研究，如船的横摇运动、结构振动、化学键的破坏等，横向波 动方程的轴向张力扰动模型，转子轴承的动力学方程也与 Duffing 系统基本相似，另外 Duffing 系统也非常广泛地被应用到实际工程中， 例如尖锐碰摩转子的故障检测、 微弱周 期信号检测、电力系统周

4、期振荡分析、周期电路系统的模拟与控制等。关于 Duffing 系 统还有许多问题尚未彻底研究清楚，如 Duffing 方程的分数谐波振动、超谐波振动、组 合振动等等，而且研究结果中规律性的成果可以推广到其他类似系统。因此从某种角度 来说，对非线性 Duffing 系统的研究是研究许多复杂动力学系统的基础。2 典型的 Duffing 方程的具体形式典型的 Duffing 方程的具体形式为：其中 x(t) 和 x(t) 为状态变量， rcos(t) 为周期驱动力， M为阻尼比， -x(t)+x 3 (t) 为 非线性恢复力。一般为了便于对系统的分析，同时达到降阶的目的，将系统描述如下：在时域中，

5、用 x(t)-t 或 x(t)-t 关系曲线描述状态变量 x(t) 或 x(t) 随时间 t 的变化 规律称为时序图。以 x(t) 为横轴和 x(t) 为纵轴所构成的平面称为相平面。在相平面上做 x(t)-x(t) 关系曲线，表征系统状态的相点 (x(t) ，x(t) )随时间变化的轨迹称为相迹图。3 混沌理论的特征和论点混沌理论的特征 目前人们把混沌看成是一种无周期的有序。它包括如下特征： 混沌的内在确定性：它表面上看起来跟噪声很相似，但是又完全不同于噪声，在系 统方程完全确定的情况下，不用添加随机参数，系统却会表现出随机的、任意的，无周 期的、不可预测又永不重复的行为。混沌的分形性质：比如

6、 Lorenz 吸引子， Henon 吸引子都是混沌吸引子，它们都具 有分形的结构；混沌的标度不变性：混沌现象是没有周期的，它是一种无序中的有序。再由经分岔 而走向混沌的通道中，还遵从 Feigenbaum 常数系。混沌系统对初始值极其敏感：就像“蝴蝶效应”一样，只要系统的初始值有微小的 变化，就会使系统的结果出现巨大的改变。由此可知对系统长期行为的预测是不能进行 的，是不准确的。混沌理论的主要论点耗散结构：混沌系统被人们认为是具有耗散结构的，其中存在不同的子系统，它们 都是不稳定的，互相不成比例也不相称，有时候达到某种状态时会出现混沌，此状态成 为临界点或临界状态。有时候也会趋向另一个层次的

7、耗散结构。蝴蝶效应：这其实指的就是系统的初值敏感性。比如，巴西的一只蝴蝶拍打翅膀， 在一段时间以后，就会造成美国某个州的飓风。这就是着名的蝴蝶效应，大气就是一个 非线性的系统，有时候会出现混沌。奇异的引子：通俗的讲，吸引子只存在于耗散系统当中，由于有能量的消耗，随着 时间的推移，相体积就会越来越小，维数也会越来越小。最终达到一个极限状态，这就 是吸引子，而不会回到原状态了。而混沌吸引子不仅对初始状态高度敏感，而且还具有 分数维结构。它在整体上是稳定的，而在局部上有时是不稳定的，无法预测它的未来状 态，所以法国物理学家儒勒和泰肯首次把它叫做奇异吸引子。回馈机能： 1984 年 Prigogine

8、 提出，系统的过去决定着它未来的走向。在系统的 运行过程中系统的结构可由非线性方程表示出由吸引子产生出结果，并且回馈 至系统成为新的输入。 系统就这样反复的进行并产生新的结构， 所以此系统是不可逆的。非线性：就像前面提到的一样，理想的线性系统是不存在的，都是经过非线性系统 的线性化近似得来的，世界上的系统都是非线性系统，世界万物之间都存在着复杂的联 系，是不可分割的。混沌系统也是非线性的。混沌的检测方法检测混沌现象的方法有很多种，比如自关联函数、 KS 熵、宽频带功率谱、波形图与 相轨图、庞加莱映射图等，其中最主要的就是以下几种：功率谱：系统运动的功率谱是数字信号处理的主要内容之一，可以通过信

9、号处理系 统得到，也可以通过 A/D 转换，把模拟信号转换成数字信号以后， 通过计算机运算得出。 在功率谱图中周期运动通常表现为尖锋，当系统出现混沌时，功率谱图会中出现" 噪声背景"和宽锋。这也是研究系统从分岔走向混沌的重要方法。相空间重构：相空间重构技术是检测混沌的简单而有效的办法。它的主要原理就是 区别噪声和混沌相轨线的区别。混沌是一个成体稳定和局部不稳定的混合体，它对初值 相当的敏感。而且混沌运动在相轨线中表现出不规则的螺旋型曲线，这是它和噪声的重 大区别。这种方法也存在很大的局限性，当系统超过 2 维的时候，你将很难把它同噪声 区别开来。李雅普诺夫指数和维数：在非线

10、性动力学的长期研究中，普遍使用李雅普诺夫指数 和维数来判断系统的状态。 因为，混沌吸引子是整体上稳定的， 而在局部又是不稳定的， 是伸长与折叠的结果。所以，根据李雅普诺夫指数的定义，如果混沌吸引子形成局部不 稳定的身长，则至少有一个李雅普诺夫指数为正，只有混沌吸引子才有这样的特征，因 此也成为定量判断混沌吸引子的依据。与指数类似，根据李雅普诺夫维数的计算公式也 可以定量的判断混沌是否存在。众所周知，混沌吸引子是具有分数维结构的，它的分数 维维数反映了它与其他吸引子相区别的特征。4 典型混沌系统和混沌同步的简介典型混沌系统的介绍 混沌从表述形式上大体包括两大类：以微分方程表述的时间连续函数和以状

11、态方程 表述的时间离散函数。时间离散系统多用于扩频通信，而时间连续函数多见于保密通信 之中。介于本文主要考虑连续系统在保密通信之中的应用，这里就重点介绍连续时间混 沌系统中的典型模型： Lorenz 系统、蔡氏电路、统一混沌系统。Lorenz 系统混沌的最早实例是由美国麻省理工学院的气象学家洛伦兹在1963年研究大气运动时描述的。他提出了着名的 Lorenz 方程组：。x a y x。y bx xz y。 xy cz(4-1) 这是一个三阶常微分方程组。它以无限平板间流体热对流运动的简化模型为基础，由于 它的变量不显含时间 t ，一般称作自治方程。式中 x表示对流强度， y表示向上流和向下 流

12、在单位元之间的温度差， z 表示垂直方向温度分布的非线性强度， -xz 和 xy 为非线性 项， b 是瑞利数，它表示引起对流和湍流的驱动因素 ( 如贝纳对流上下板的温度差 T) v a 和抑制对流因素 ( 如(Prandtl) 数粘性)之比，是系统 (4-1) 的主要控制参数。 k 是普 朗特数(v 和 k分别为分子粘性系数和热传导系数 )，c 代表与对流纵横比有关的外形比， 且 a和 c为无量纲常数。在参数范围为 b a (a c 3)/(a c 1)时，Lorenz 系统均处于 混沌态。在混沌区域内选择系统参数 a=10, b=28 ，c=8/3 ，取系统的初始状态为 x(0), y(0

13、), z(0)=10, 10, 10 ，此时，系统为一混沌系统，系统的三维吸引子如图所示，二维吸 引子如图所示，图所示分别为分量 x、y 随时间 t 的变化情况。图 Lorenz 系统的吸引子 图 分量 x 随时间 t 的变化情况 图 Lorenz 系统的 x-y 相图 总体上， Lorenz 吸引子由左右两个环套而成，每个环绕着一个不动点，它实际上是 一条双螺旋的曲线，就像以十分灵巧的方式交织起来的一对蝴蝶的翅膀。这个吸引子中 的环和螺线有无穷的深度，它们之间可以无限靠近，但永远不会相交，仅占据有限的空 间，具有无穷嵌套的复杂结构。例如，随着时间的演化，每一个环都靠得很近的无穷多 层，每层上

14、都密密麻麻的排列着无穷多个螺线，它代表系统的相点在右侧转几圈后又跳 到左侧转几圈，运动轨道无法预测什么时候从这一侧过渡到另一侧，并且它所绕各自中 心的方式和圈数也是个明显的随机数。这就是混沌状态。在诸多用于混沌研究的非线性电路系统中，有一类是分段线性的非线性系统，它们 共同的特点是易于进行数学分析和物理实现。其中蔡氏电路最为着名。在许多文献中都 以蔡氏电路为基础研究了混沌现象。它的优点在于极为简单的系统就产生了极为复杂的 动力学行为，人们从不同角度研究蔡氏电路以及它的变形电路，发现了几乎所有目前已 知的动力学现象，如倍周期分叉、周期运动、阵发性混沌、奇异非混沌、混沌等。蔡氏 电路及其变形电路己

15、经成为研究动力学行为的一个重要模型，而且蔡氏电路的应用研究 已经遍及诸多领域，如保密通信、手写特征识别、音乐的产生，以及利用蔡氏电路组成 的时空系统进行轨迹识别等。蔡氏电路是一个三维自治振荡系统，如图 (a) 所示，由四个线性元件电感 L,电阻 R，电容 C1，C2和非线性电阻 N 组成。蔡氏电路可由以下状态方程描述：dvc11G vC2vC1g(vc1)dtC1dvC21G vC1vC2iLdtC2diL1vC2dtL(4-2)其中 G=1/R， g(vC1)是蔡氏二极管的分段线性 v i 特性，如图 (b) ，表示g(vC1) m0vC112(m1m0)( vC1BPvC1BP )(4-3

16、)现简化电路，令t G ,x vC1 ,yC2BPvC2,ziLBPBPGC2C1C2RL2,aR m1,bR m0，我们可得到该电路的无量纲方程：y x f(x)yz(4-4)其中f (x) bx12(ab) x 1(4-5)图蔡氏电路及其非线性电阻特性a)蔡氏电路(b) NR 伏安特性分析该混沌电路模型有如下特点：(1) 方程(4-4) 关于状态空间原点是对称的，即式 (4-4) 中的(x, y, z)用(-x,-y,-z) 代替时，方程保持不变。且在状态空间的三个子空间中各有一个平衡点，分别记为： p k,0, k ,p 0,0,0 ,pk,0,k ，其中ba b 1 。A a,b,c(

17、2) 在上述的子空间中，方程均属线性，取：ac10令xx,y,z ,K k,0,k , 则系统的状态方程可写成：A,a x KxDxA, ,b xx D0A, ,a x KxD(4-6) 蔡氏二极管的伏安特性与参数 m0,m1,BP 的选取有关，其中 m0 , m1的大小直接关系到蔡氏 电路的解行为和吸引子的结构。 为简便起见，将蔡氏电路中的吸引子分成两类： 当 m0 m1 称为第一种蔡氏吸引子；当 m0 m1 称为第二种蔡氏吸引子。这里只仿真第一类混沌吸引 子。取 10, 14.87,a 1.27,b 0.68时，初始值为 x0 0.513, y0 0.067,z0 0.385 其吸引子如图

18、，所示。图蔡氏电路第一类混沌吸引子图蔡氏电路的 (x,y) 相图1963 年，美国气象学家 . Lorenz 在对流实验中发现了第一个混沌吸引子， Lorenz 吸引子为混沌研究提供了一个重要模型。七十年代以来掀起了一股揭示混沌现象，研究 混沌理论的热潮。 1999 年，陈关荣教授在研究混沌反控制过程中发现了一个与 Lorenz 类似的混沌吸引子 -Chen 吸引子 24 。2002年，吕金虎、陈关荣等又提出了一个新的混 沌系统 - 统一混沌系统 25-27 ，这一系统连接了 Lorenz 吸引子和 Chen吸引子。 统一混沌系统的数学模型为：x2510x2x1。x22835x129 1 x2 x1x3。x3x1x28x3 /3(4-7)式(4-7) 中 为系统参数， 当0,1时系统均为混沌态。 当 0,0.8) 时，统一系统属于广义 Lorenz 系统；当 (0.8,1时，统一系统属于广义 Chen系统； 0.8属于广义 Lu 系 统。现以 0.46 为例，给出统一混沌系统吸引子的数值仿真结果，如图，所示。图