非线性动力学定性理论方法

合集下载

定性与定量相结合的分析方法引文分析法

《信息素养与信息检索》精品资源共享课
4.2 定性分析法
4.2.3 归纳推理法
作用
提供假说。证明假说和理论。确定假说的支持度。理论择优。对事件未来情况进行预测。各种管理决策。
《信息素养与信息检索》精品资源共享课
4.3 定性与定量相结合的分析方法
4.3.1 层次分析法
含义
将一个复杂的多目标决策问题作为一个系统，将目标分解为多个目标或准则，进而分解为多指标（或准则、约束）的若干层次，通过定性指标模糊量化方法算出层次单排序（权数）和总排序，以作为目标（多指标）、多方案优化决策的系统方法。
《信息素养与信息检索》精品资源共享课
4.2 定性分析法
4.2.3 归纳推理法
局限性
它只能得出不充分可靠的结论。它未必把握住事物的本质。概括是静态的、过去的，难以概括它的发展和未来。
《信息素养与信息检索》精品资源共享课
4.2 定性分析法
4.2.3 归纳推理法
局限性
它只能得出不充分可靠的结论。它未必把握住事物的本质。概括是静态的、过去的，难以概括它的发展和未来。
《信息素养与信息检索》精品资源共享课
4.1 定量分析方法
4.1.1 时间序列法
含义
按照时间的顺序把随机事件变化发展的过程记录下来就构成了一个时间序列。对时间序列进行观察、研究，找寻它变化发展的规律，预测它将来的走势。挖掘事物随时间变化规律的方法
《信息素养与信息检索》精品资源共享课
4.1 定量分析方法
含义
通过研究系统内部诸多因素形成的各种反馈环，同时搜集与系统行为有关数据和信息，采用计算机仿真技术来对大系统、巨系统进行长期预测的方法。系统性、动态性、仿真性

非线性动力学入门-西安交通大学教师个人主页

. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . .
另一方面梁的轴向应变的表达式也会因变形大小的不同而采用不同的表达式比如小变形时应变而当考虑大变形时可能采用的应变表达式就是进而得到的梁的振动方程将会是一个含有高度非线性项的偏微分方程组
非线性动力学入门
张新华
西安交通大学工程力学系 2011 年 07 月
前言
─1687 年，牛顿(Isaac Newton, 1643 ～ 1727)发表了《自然哲学之数学原理》(Mathematical Principles of Natural Philosophy)，标志着经典力学(亦即牛顿力学)的正式诞生。牛顿力学主要研究自由质点系的宏观运动规律。 ─1788 年，拉格朗日(Joseph Louis Lagrange, 1736 ～ 1813)发表了分析力学教程(Analytical Mechanics)，标志着拉格朗日力学的诞生。Lagrange 力学属于分析力学的主要内容之一，在位形空间中研究带有约束的质点系动力学。 ─1833 年，哈密尔顿(William Rowen Hamilton, 1805 ～ 1865)对 Lagrange 力学进行了改造，引进了相空间(2n 维空间)，对系统内在的对称性(辛对称， Symplectic)进行了刻画。狭义上的哈密尔顿力学只适用于保守系统，而广义的哈密尔顿力学在适用于非保守系统。哈密尔顿力学也属于分析力学的主要组成部分。在此后发展起来的量子力学中 Hamilton 力学发挥着巨大的作用。目前在天体力学、计算 Hamilton 力学，量子力学，甚至弹性力学(即所谓的辛弹性力学)中哈密尔顿力学依然发挥着重要作用。 ─1927 年，Birkhoff(George David Birkhoff, 1844 ～ 1944)发表了“动力系统”(Dynamical Systems)，标志着 Birkhoff 动力学的正式问世。Birkhoff 动力学建立了研究非完整力学的框架。 ─1892 ～ 1899，彭加莱(Henri Poincaré, 1854 ～ 1912)发表了三卷本的“天体力学中的新方法”(New Methods of Celestial Mechanics)，系统性地提出了研究动力学系统的定性方法，即几何方法。经典力学的目标之一就是设法求得系统的解析解，而 Poincaré意识到对于大多数非线性系统而言，求其解析解是不可能的，而必须发展新的研究方法。他超越了他的时代，极富远见地预测到了非线性系统混沌现象(系统的解对初始条件具有极端敏感依赖性)的存在。更为重要的是，Poincaré开创了研究非线性动力系统的几何方法，当之无愧地被誉为非线性科学之父，其影响是划时代的。 ─1892 年，李亚普诺夫(Aleksandr Mikhailovich Lyapunov, 1857 ～ 1918)在他的博士论文“运动稳定性的一般问题”(General problem of the stability of motion )中，系统地探讨了非线性动力学系统的稳定性问题。他提出了两种研究稳定性的方法：李亚普诺夫第一方法(间接方法)和李亚普诺夫第二方法(直接方法)。他从代数角度出发，对动力学系统的研究开创了一个崭新的领域。彭加莱与李亚普诺夫，前者从几何角度，后者从代数角度，开拓了非线性科学的研究疆域和研究手段。 ─1963 年，Lorenz(Edward Norton Lorenz, 1917 ～ 2008)发表了“确定性非周期流”(Deterministic Nonperiodic Flow)的论文，认为大气系统的性态对初值极为敏感，从而导致准确的长期天气预报是不可能的。该文标志着人类首次借助于计算机发现了混沌(Chaos)现象的存在。 ─1757 年，欧拉(Leonhard Euler, 1707 ～ 1783)发表了压杆稳定性的论文，首次探讨了力学系统的分岔现象。作为分岔理论重要分支的突变理论(Catastrophe Theory)则主要由法国数学家托姆(René Thom, 1923 ～ 2002)于上个世纪 60 年代创立，由齐曼(Christopher Zeeman，1925 ～)在 70 年代大力推广普及。 ─1834 年，英国的罗素(John Scott Russell, 1808 ～ 1882)骑着马在 Union 运河上散步时，发现了现在称之为孤立波(又称作孤波，Solitary wave)的 i

非线性动力学(中科大课件)

习题1：求G的显示表达式
Байду номын сангаас
有些系统，其相空间不是线性空间，这类系统必定是非线性的。例如：刚体动力学方程，相空间为SO3 流形。称为运动(Kinetic)非线性。 F的非线性称为力学非线性。
14

2018/12/23
b) 什么是特征非线性现象？
例：
考虑映射
设
f(x)>x
n
则xn+1> xn
2018/12/23 5
19世纪 Hamilton 方程
引进广义动量
Hamilton
函数
2018/12/23 6
Lagrange,
Hamilton方程均与Newton方程等价力学的基本问题、基本方法未变结论：二百多年来牛顿力学无实质性进展
2018/12/23 7
19世纪末、20世纪初发生的对于牛顿力学的三大变革
涨落项
18
宏观描述

在长时间后，绝大多数集体运动模式
由于耗散而衰减掉，可以考虑剩余运动模式。
设长时间后只剩下 x1 ,…, xm 的运动，则在 t→∞ 时有

Haken称之为随动原理(slaving principle). 代入前一方程，消去 xj , j=m+1,…,n 得

xj , j=1,…,m 称作序参量（基本力学量）绝热消除法：中心流形→ 惯性流形
19
2018/12/23
FTP Complexity
IP address：

202.38.83.243 User name： guest Password： guest
20
2018/12/23

数学建模中对非线性动力系统模型的认识和体会

数学建模中对非线性动力系统模型的认识和体会真实动力系统几乎总是含有各种各样的非线性因素，诸如机械系统中的间隙、干摩擦，结构系统中的材料弹塑性和黏弹性、构件大变形，控制系统中的元器件饱和特性、控制策略非线性等等。

非线性动力学理论的研究和发展已经经历了一个多世纪，在新世纪之初，为了使非线性动力学理论得到更好的发展，非常有必要回顾一下非线性动力学研究和发展的历史。

非线性动力学理论的发展大致经历了三个阶段。

第一个阶段是从1881年到1920年前后，第二阶段从20世纪20年代到70年代，第三阶段从20世纪70年代至今。

人们对于非线性系统的动力学问题的研究可以追溯到1673年Huygens对单摆大幅摆动非等时性的观察。

第一阶段的主要进展是动力系统的定性理论，其标志性成果是法国科学家Poincare从1881年到1886年期间发表的系列论文“微分方程定义的积分曲线”，俄罗斯科学家Liapunov 从1882年到1892年期间完成的博士论文“运动稳定性通论”，以及美国科学家Birkhoff在1927年出版的著作“动力系统"。

第二阶段的主要进展是提出了一系列求解非线性振动问题的定量方法，代表人物有俄罗斯科学家Krylov、Bogliubov，乌克兰科学家Mitropolsky，美国科学家Nayfeh等等。

他们系统地发展了各种摄动方法和渐近方法，解决了力学和程科学中的许多问题。

在这个阶段中抽象提炼出了若干著名的数学模型，如Duffing方程、vander Pol方程、Mathieu方程等，至今仍被人们用以研究非线性系统动力学现象的本质特征。

从20世纪60~70年代开始，原来独立发展的分岔理论汇入非线性动力学研究的主流当中，混沌现象的发现更为非线性动力学的研究注入了活力，分岔、混沌的研究成为非线性动力学理论新的研究热点。

俄罗斯科学家Arnold和美国科学家Smale等数学家和力学家相继对非线性系统的分岔理论和混沌动力学进行了奠基性和深入的研究，Lorenz和Ueda等物理学家则在实验和数值模拟中获得了重要发现。

动态系统理论解读非线性动力学行为

动态系统理论解读非线性动力学行为动态系统理论是研究系统随时间变化的数学理论，在物理学、生物学、经济学等领域广泛应用。

非线性动力学是动态系统理论的一个重要分支，研究的是非线性系统的行为。

非线性动力学行为指的是系统中存在非线性因素导致的复杂行为，这些行为通常无法通过简单的线性理论来解释。

非线性动力学行为的研究领域包括混沌理论、奇异吸引子、分岔现象等。

混沌理论是非线性动力学行为的重要组成部分。

混沌现象指的是一个看似没有规律的、极为敏感的动态行为，它对初始条件极为敏感，微小的初始条件差异可能会导致系统最后的行为完全不同。

混沌现象的典型例子是著名的“蝴蝶效应”，即一个蝴蝶在巴西扇动翅膀可能最终引起美国得克萨斯州发生龙卷风的现象。

奇异吸引子是一种特殊的吸引子，它具有分岔结构。

吸引子是动态系统中一组确定的状态，而奇异吸引子则是一种分维度小于系统自身维度的吸引子。

奇异吸引子的特点是具有分形结构，即在不同尺度上具有相似的形状。

分岔现象是非线性动力学中的一个重要现象，它表示系统参数改变时出现的定性变化。

在分岔现象中，随着参数的改变，系统从一个稳定状态转变为多个稳定状态或不稳定状态。

这种转变可以是突然的、跳跃的或连续的，而且是可逆的。

非线性动力学行为的研究对于理解现实世界中复杂系统的行为模式具有重要意义。

在物理学中，非线性动力学行为可以帮助解释天体运动、流体力学等现象。

在生物学中，非线性动力学行为可以解释生物系统中的自组织、自适应等特性。

在经济学中，非线性动力学行为可以用来解释经济周期、市场波动等现象。

非线性动力学行为的研究方法包括数学建模、理论分析和计算机模拟等。

数学建模是非线性动力学研究的基础，可以将系统的动力学行为用方程或规则来描述。

理论分析通过数学方法对系统的动态行为进行解析，寻找系统的稳定状态和边界条件等。

计算机模拟则可以通过计算机程序对系统进行模拟，观察系统的行为变化。

然而，非线性动力学行为的研究也面临着一些挑战。

非线性振动

非线性振动的研究包括理论分析方法和数值分析方法。

其中理论分析方法有是沿着两个方向发展，第一是定性方法，第二是定量方法，也称为解析法。

定性方法是对方程解的存在性、唯一性、周期性和稳定性等的研究；定量方法是对方程解的具体表达形式、数量大小和解的数目等的研究。

数值方法目前已广泛用于计算非线性振动系统，是一种求解非线性方程的有效方法。

本文在查询相关文献的基础上，对非线性振动理论的分析方法最新研究成果做简要概括和分析比较。

1、平均法平均法是求解非线性振动最常见和最实用的近似方法之一。

其基本思想是设待解微分方程与派生方程具有相同形式的解，只是振幅和相位随时间缓慢变化。

将振幅和相位的导数用一个周期的平均值替代，得到平均化方程，求解平均化方程，得到振幅和相位的表达式，从而求解出原方程的近似解析解。

1.1利用平均法分析多自由度非线性振动平均法主要是用在单自由度非线性振动的分析中，是一种求近似解的方法，虽然精度较低，但可避免繁琐的中间运算，具有便于应用的突出优点。

将其推广的到多自由度系统，导出了平均化方程，由此能够得到多自由度非线性振动的幅频特性。

1.2用改进平均法求解自由衰减振动用平均法求解自由衰减振动方程时，无论是线性阻尼还是平方阻尼，在阻尼常量很小的情况下，平均法解均有较高的精度。

但随阻尼常量的增加，阻尼对振动周期的影响已不能忽略，此时平均法解的结果与实际振动情况有了明显的偏离，需要改进。

改进平均法是将待解微分方程的圆频率与派生方程圆频率的差异函数表示为阻尼系数的多项式。

2、FFT多谐波平衡法分析非线性系统非线性动力系统的响应可能含有几个主导频率，且有可能与激振频率不成倍数关系。

现有的单一谐波法和多谐波法仅限于系统响应主导频率为激振频率的非线性系统，因此在某些情况下使用单一谐波法或多谐波法研究非线性系统动力学特性是不可靠的，而基于快速傅立叶变换（FFT）和主导频率的 FFT 多谐波平衡法能够依据所有的主导频率构筑多谐波平衡方程，因此其解析解精确度高，并能广泛适用于单倍周期、多倍周期、与初始条件有关的多解性及拟周期响应等典型的非线性特征响应。

为什么要学些火灾动力学？

为什么要学些火灾动力学?前不久，消防圈有一种论调为：“指挥员源于战斗员，亘古不变的道理”，言下之意，只有干过,才有发言权。

这种思想似乎很普遍，听上去似乎很有道理。

本来嘛，救人靠技能，技能靠实战,有实战经验的似乎比没有实战经验的强.然而，2015年8月12日的天津港大爆炸表明,救援除了依赖经验，还需要有必要的知识和战略的眼光。

说起军事指挥家，近代中国无人超过“曾胡左李”。

曾国藩、胡林翼、左宗棠、李鸿章,这四人无一当过战斗员,但当时人都毫不怀疑他们的指挥能力.即使是解放战争期间，林彪和粟裕，都是黄埔军校出身，典型的学院派，也不是战斗员出身，丝毫不能妨害他们的指挥能力。

今天，高层次的指挥员大多是学院派，美国的消防队长也要求有本科BS学位，所以指挥这个岗位,恰恰不是战斗员所能够胜任的,因为他需要理科知识（Fire Science)和工科素养(Fire Protection Engineering），才能够认识火场的复杂现象，才能发出更好的指令,取得更有效的战术效果。

战斗员自救是自己的事，指挥员的命令可以不听（不能影响生命安全）,因此指挥员是否干过战斗员，并不是很重要的一件事。

我们听到有些声音，火灾动力学对灭火工作没什么用。

但其实这门知识博大精深，涉及很多方面，是非常有用的.Drysdale 教授阅读并引用了800篇文章，才写出了一本《火灾动力学》教材，无人不说有用。

（京东有售,见下图，点击’阅读原文’，查看火灾动力学试读链接。

）图：火灾动力学第三版你说无用，就是对消防没有入门的标志。

其实，在理论的知识和应用的知识之间，存在沟通问题。

不能很好地应用和传递知识,相当于没有知识。

美国消防队员需要学习火灾动力学知识，源于一项技术革命。

Nomex纤维，是一种特殊的阻燃材料,刚出现的时候，非常昂贵,只能用在高度风险的场合,如赛车手的特殊防护服。

随着成本的下降，现代消防队员也开始普及Nomex纤维制成的防护服，并携带自己的自助呼吸器（SCBA）。

系统工程导论第3讲

第三讲自然物理系统及其复杂性系统的层次来看，自然物理系统处于系统底层，可以看作是由非智能的机械物体构成的，满足特定的自然物理规律，其复杂性的根源在“非线性”。

在系统科学中，迄今真正成熟的主要是线性系统理论。

系统科学重点研究的是非线性系统，是处理非线性问题的一种方法论。

同时，相比于静态系统，动态系统是系统科学讨论的重点对象。

一、线性系统1 线性特性和线性系统一般地说，能够用线性数学模型描述的系统，称为线性系统。

线性系统的基本特性，即输出响应特性、状态响应特性、状态转移特性等，都满足叠加原理。

具体讲，令f代表某种数学操作，如关系、变换、运算、方程或其他，x为数学操作的对象，f(x)表示对x施行操作f的结果。

若f(x)满足以下两个条件(1) 加和性：f(x1+x2)＝f(x1)+f(x2)(2) 齐次性：f(kx)＝kf(x)即f(ax1+bx2z)＝af(x1)+bf(x2)，就称操作f为线性的，满足叠加原理。

2、线性系统的动态行为描述描述连续动态系统的数学模型是微分方程。

如果状态变量只是时间的函数，与空间分布无关，则称集中参数系统，用常微分方程描述；状态变量同时依赖于时间和空间分布的是分布参数系统，须用偏微分方程描述。

线性连续动态系统的数学模型为线性常微分方程，可以使用一元高阶方程，也可以使用多元一阶联立方程组，这两种形式是等价的。

一般形式如下：nnn n n nn x a x a x x a x a x ++=++=11'1111'1记系数矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=nn n n a a a a A 1111，状态向量()Tn x x x X ,,,21 =则上述动态方程可表示为：AX X ='A 不随时间变化→常系数方程→时不变线性系统A 随时间变化 →变系数方程→时变线性系统A 包含了线性系统一切行为特性的信息。

3、线性动态系统的求解分析线性系统的常用方法是求它的通解：∑=k tk V e C X k λ其中，λ为特征方程0=-I A λ的特征根；V k 是一组线性无关的特征向量；C k 是由初值决定的常系数。

非平衡非线性化学动力学侯中怀

第12章非平衡非线性化学动力学侯中怀 hzhlj@中国科学技术大学化学物理系合肥 230026非线性化学动力学的研究对象，是化学体系在远离平衡条件下，由体系中非线性过程的作用，自发形成的宏观尺度上的各种复杂的时空有序结构，包括多重定态，化学振荡，图灵斑图，化学波和化学混沌等[1-3]。

这些现象都是非平衡条件下大量分子的集体行为，因此非线性化学动力学的研究，属于物理化学和非平衡统计物理的交叉领域。

随着20世纪50年代BZ化学反应体系中各类非线性化学现象的实验发现，非线性化学动力学的研究便成为物理化学研究中的一个新的生长点。

20世纪70年代，以普里高津(Prigogine)为首的比利时布鲁塞尔学派提出了著名的“耗散结构”理论[4,5]，奠定了非线性化学现象的热力学基础。

过去20年，计算机技术和非线性科学的发展，使得人们能从理论上再现实验上观测到的各种非线性现象，以深入了解非线性化学现象的动力学机制，从而进一步推动非线性化学动力学在实际体系中的应用。

近年来，随着化学研究的对象向生命和纳米等复杂体系的深入，非平衡、非线性和复杂性之间的相互作用，目前是非线性化学动力学研究的一个主要发展方向。

在生命和表面催化等体系中，实验上已发现大量的非线性动力学行为，如细胞体系内的钙振荡及钙波[6]，生理时钟振荡[7]，单晶表面催化过程中的化学振荡、螺旋波、化学混沌等[8,9]。

研究表明，这些非线性化学动力学行为，对生命体系的功能和催化过程的活性与选择性等，起着非常重要的作用；要深入理解这些作用的机制，必须考虑到实际体系中的各种复杂性因素，包括噪声和无序等随机因素，环境和体系以及体系内部的复杂相互作用等。

本章中，我们将对非线性化学动力学的基本内容和研究进展作一简单概述。

为使内容具有相对完整性，第一节主要介绍非线性化学动力学的基本概念和研究方法。

在第二节和第三节，将重点介绍近年来复杂体系非线性化学动力学的一些研究结果，主要包括环境噪声、空间和拓扑无序、介观反应体系内涨落对非线性化学动力学的调控作用等。

数学中的动力系统理论及其应用研究

数学中的动力系统理论及其应用研究数学中的动力系统理论是一门非常重要的学科，它研究的是不断变化的系统及其变化规律。

这门学科有着广泛的应用，如生物学，物理学，经济学，力学等各领域都有其应用。

动力系统学说的发展历程长达两个世纪，最早的研究者是拉格朗日和欧拉，其后，哈密顿、普朗克、庞加莱、尼古拉斯·罗伯特∙比尔、基尔霍夫等数学家都取得了丰硕的成果。

一、动力系统的基本概念动力系统指的是一个具有演化规律的系统，通常用一个向量场表示，向量场用来表示物体的速度和方向，从而得到物体运动的轨道。

用数学方程式表示一个动力系统通常有两种方式，一种是演化方程式，即在以时间为自变量的空间中描述变化；另一种是离散映射，即在以整数为自变量的空间中描述变化。

我们称系统动力学模型研究的是与时间变化有关的系统的动力学过程。

二、动力系统的重要性动力学是研究非线性系统的一个重要领域，非线性系统是指一类物理或化学系统，其各种成分之间的关系为非线性的，所以同样的初始条件可能会导致不同的结果。

动力学的研究有助于人们认识复杂的现象，分析现象的发展趋势，预测未来的变化过程。

三、动力系统的应用1.生物学领域动力系统在生物学领域有着广泛的应用，例如生物钟模型的研究。

生物钟是指生物体内内在的自然节律，非常重要的生理功能，它的研究是生物学领域的热点之一。

利用动力学模型，可以对生物钟的内部规律进行研究和解析。

2.物理学领域在物理学领域，动力学模型可以用来描述各种系统的运动规律和相互作用。

例如，研究天体运动、宇宙演化和分子运动都需要使用动力学模型。

3.经济学领域经济学是有关经济及其相关现象的一门社会科学，其中的动力学模型可以用来分析经济体系中不断发生变化的现象，包括商业、金融等各个方面。

四、动力系统的研究方法1.定性方法定性方法是动力学研究中的基本方法，它着重于从数学上分析动力系统的本质，研究系统的稳定性和不稳定性等基本性质。

2.定量方法定量方法是通过使用精确的数学工具来研究动力系统，例如微积分、差分方程、拓扑学、几何等过程。

Brusselator模型的定性分析及数值模拟

第３卷第２期１
２００７年３月
江西师范大学学报（自然科学版）
ＪＵＲＡＬＯＦＪＧＲＭＡＬＵＮＩＲＩＹ（ＴＡＬＳＩＮＣ）ＯＮＩＡＮＸＩＮＯＶＥＳＴＮＡＵＲＣＥＥ
ＶＯ．１ＮＯ２１３．
Ｍａｒ．２７Ｏ０
这样我们就证明了是系统4的不变区域即表明对于满足条件5以及边界条件6或7的系统4只要初始条件uo髫及边界条件g戈或日戈位于内则4的解就永远保持在中即对于第一象限中的任何点p都可以取到足够大的使得p所以对于系统4只要其初值及边界位于第一象限则整体解在第一象限中存在
维普资讯
ｆｕｄ＝Ａ＋Ｍ一Ｂ＋），ｄ／ｔ（１Ｍ
【／ｔ：ＢＭ — ｌ．ｄｄＺ２
ｒ、，
一
从动力学角度我们来分析系统（）的解的性质．先，２只有唯一的平衡解（ＢＡ）这是因为方程２首（）Ａ，／，（）２的两式相加得到Ｍ：，Ｍ：代人第二式得到＝ＢＡ．将／其次，讨论这一平衡解的稳定性．为此考虑（）Ｊｃｂａ阵：２的ａｏｉｎ矩
间元素起自催化作用，且两种元素呈周期性变化，其结果不同于一般的化学反应最终达到一个平稳，而是两种元素长时间停留在一个时间上振动（ｓｉａｉｎ，间上又有扩散的过程＿．一反应称为Ｂｌｓｖｏｃｌｔ）空ｌｏ２这Ｊｅｏｏ— ｕＺａｏｉｋ反应，ｈｂｔｓｖｎ具体分下面４步进行．
２
关键词：反应扩散；极限环；显式差分中图分类号：２１０２．４文献标识码：Ａ

常微分方程-定性理论

混沌理论主要研究的是非线性常微分方程的复杂动态行为，包括对初值敏感依赖、长期行为的不可预测性和奇怪吸引子等特性。
分岔理论
当系统的参数发生变化时，常微分方程的解可能会发生突然变化，这种现象被称为分岔。分岔理论研究了分岔的产生条件和分岔的类型。
数值解法
对于无法解析求解的常微分方程，定性理论还研究了各种数值解法，如欧拉法、龙格-库塔法等，以近似求解其解。
稳定性是指系统在平衡点附近的动态行为，如果一个平衡点是稳定的，那么当系统受到小的扰动时，它会回到平衡点；反之，如果平衡点是不稳定的，那么系统会远离平衡点。
线性常微分方程
线性常微分方程是指形式为dy/dt = ay + b的常微分方程，其中a和b是常数，y是未知函数。
VS
线性常微分方程的解可以通过求解线性代数方程得到，其解的性质可以通过特征值和特征向量来描述。
定性理论通过分析微分方程的解轨线在相空间中的行为，来理解和预测系统的动态行为。它为解决实际问题提供了重要的数学工具和理论基础。
研究目的和意义
研究目的
常微分方程定性理论的研究目的是揭示微分方程解的内在性质和规律，理解解的动态行为，并应用于解决实际问题。
研究意义
定性理论在数学、物理学、工程学、经济学等领域有广泛的应用，是解决实际问题的重要工具。通过研究常微分方程定性理论，可以深入理解自然现象和社会现象的动态演化过程，为解决实际问题提供重要的数学方法和理论基础。
人口动态
常微分方程可以用来描述人口的变化规律，例如 Malthus模型。
行为科学
常微分方程可以用来描述人类行为的变化规律，例如心理动力学模型。
05 结论与展望
研究结论
稳定性理论

CCEBC_EEAC方法的定性分析_廖浩辉

818 中国科学E辑工程科学材料科学2004, 34(7): 818~831CCEBC/EEAC方法的定性分析*廖浩辉**唐云***(清华大学数学科学系, 北京 100084)摘要 CCEBC/EEAC 方法是对电力系统暂态稳定性作量化分析的一个很有效的方法. 对CCEBC/EEAC方法进行了定性分析, 说明从几何观点出发, CCEBC/EEAC方法中的CCCOI-RM可以看作是对系统模型的变量在权向量空间上的投影, 由此得到广义P-轨迹. 由于电力系统的暂态过程可以近似地看作分时间段简单Hamilton系统, 为了给出CCEBC/EEAC方法的定性分析, 比较了两机无穷大系统的势能与CCEBC/EEAC能量. 计算结果表明两者相似. 弄清这两者之间的定性关系, 对于从数学上论证CCEBC/EEAC方法能够用于判定系统的暂态稳定性将起到重要作用. 通过定性分析, 也能更好地了解 CCEBC/EEAC 方法基于量化分析所揭示的一些重要现象, 如多摆失稳和孤立稳定域.关键词CCEBC/EEAC 分时间段简单Hamilton系统暂态稳定性CCEBC/EEAC(complementary-cluster energy-barrier criterion/extended equal area criterion)方法是电力系统暂态稳定性分析的一个有效方法, 其商品化软件已得到重要应用[1]. 该方法是薛禹胜[2]给出的等面积法则(EAC)对多机系统的一种推广. CCEBC/EEAC方法的基本思想最初在文献[3]中以EEAC方法的形式提出, 它通过保稳线性变换把多机问题化为两机问题, 然后应用EAC方法计算临界清除时刻和稳定裕度. EEAC方法依赖于电力系统的经典模型[3], 为了把EEAC方法应用到电力系统精确模型的暂态稳定性分析, 薛禹胜在文献[2]中提出了CCEBC 方法, 并用于分析和推广EEAC方法, 得到了CCEBC/EEAC方法, 而原来的EEAC方法可看作CCEBC/EEAC方法的一个特例, 称之为SEEAC方法.CCEBC/EEAC 方法的推动力是计算机速度的不断提高. 在传统的时域分析中, 用来判定系统是否趋于稳定的数值积分所需要的时间是比较长的, 这是由于2003-05-22收稿, 2003-10-27收修改稿*国家重点基础研究专项经费(批准号: G1998020309)和国家自然科学基金(批准号: 10272059)资助项目**E-mail:liaohaohui00@***E-mail: ytang@第7期　廖浩辉等: CCEBC/EEAC方法的定性分析　819820　中国科学E辑工程科学材料科学第34卷　第7期　廖浩辉等: CCEBC/EEAC方法的定性分析　821822　中国科学E辑工程科学材料科学第34卷　第7期　廖浩辉等: CCEBC/EEAC方法的定性分析　823824　中国科学E辑工程科学材料科学第34卷　第7期　廖浩辉等: CCEBC/EEAC方法的定性分析　825图1 系统当τ =3.7时的实际势能与CCEBC/EEAC能量的比较1示实际势能, 2示CCEBC/EEAC能量图 2 广义P-δ轨迹(τ =3.7)量之间的联系, 是建立引言中第(Ⅱ)~(Ⅴ)的数学理论的基础. 其中(Ⅴ)是最重要的问题, 因为电力系统工程关心在特定时间内是否可以判定系统的暂态稳定性, 而不需要花费较长的积分时间.3CCEBC/EEAC揭示了电力系统中的重要现象以上的分析指出了对CCEBC/EEAC方法进行严格化的困难. 不管是薛禹胜的工程观点, 还是本文采用的几何观点, 都未对CCEBC/EEAC方法进行严格化. 但是, 采用几何观点可以对CCEBC/EEAC方法所揭示的多摆失稳和孤立稳定域826　中国科学E辑工程科学材料科学第34卷　图3 系统当τ = 3.9时的实际势能与CCEBC/EEAC能量的比较1示实际势能, 2示CCEBC/EEAC能量图4 广义P-δ轨迹(τ =3.9)(ISD, isolated stable domain)这些在电力系统可能发生的重要现象[2]给出定性解释, 为认识电力系统中暂态能量函数(TEF)方法与CCEBC/EEAC方法的联系提供了一个途径.系统(6)的一个摆次是指势能V(δ)极大的两个相邻时刻所决定的时间区间. 如果系统(6)在出现多个摆次后失稳, 则称系统发生多摆失稳. 多摆失稳可以看作是n>个自由度Hamilton系统中的一个自然的现象. 这是因为系统的轨迹保持(1)Hamilton量, 由势能与动量的相互转换关系, 系统会出现多摆现象. 因此, 从数学的角度来说, 多摆失稳并不是很特别的现象. 但是, 通过对CCEBC/EEAC与暂第7期　廖浩辉等: CCEBC/EEAC方法的定性分析　827图6 系统出现ISD和NARI时, 用二分法可能求到错误的cct1~4表示二分法的划分对于不同的清除时刻τ, 系统的暂态稳定性可从图7~9观察到. 从这些图容易明白ISD产生的原因在于失稳轨迹可以离开故障后系统的稳定域后在有限时刻内又重新进入.对于实际的电力系统而言, 由于故障轨迹可以离开后又进入故障后系统的稳定域内, 所以也和分时间段常微分方程一样, 会出现ISD. 在计算稳定域的临界边界时, 因为ISD的出现, 搜索方法可能会误求错误的临界点, 从而导致NARI[2]. NARI可以认为是稳定域与ISD在参数空间中的一个凸包. 这容易通过图6来描述.ISD的发现意味着在设计程序时必须注意搜寻算法的有效性. 时域方法通常通过二分法来搜寻T cct, 而利用二分法所求得的T cct不一定是系统真正的临界故障清除时刻(见图6), 因此必须注意搜寻步长. TEF方法通过故障轨迹第一次离开稳定域来判定T cct, 因此在理论上不会出现NARI现象. 而CCEBC/EEAC方法通过排除掉首摆和多摆失稳来确定出正确的T cct, 理论上避开了NARI. 在实际计算当中, 由于缺乏比较, 还难以下定论.图7 系统(12)当τ = 5.00, 5.52, 5.56时的轨迹图8 系统(12)当τ = 8.60, 8.73, 8.80时的轨迹图9 系统(12)当τ = 18.40, 18.48, 18.50时的轨迹致谢作者感谢薛禹胜教授的热情帮助以及审稿人的批评指正.参考文献1 薛禹胜. 滑步与发散, 运动系统与一般动力学系统. 电力系统自动化, 2003, 27(7): 17~212 薛禹胜.　运动稳定性量化理论:　非自治非线性多刚体系统的稳定性分析.　南京:　江苏科学技术出版　社,　19993 Xue Y, Van Cutsem T, Ribbens-Pavella M. A simple direct method for fast transient stability assessment oflarge power systems. IEEE Trans on Power Systems, 1988, 3(2): 400~4114 Zhang Y, Wehenkel L, Rousseaux P, et al. SIME: A hybrid approach to fast transient stability assessmentand contingency selection. Int J Electrical Power and Energy Systems, 1997, 19(3): 195~2085 Zhang Y, Wehenkel L, Pavella M. SIME: A comprehensive approach to fast transient stability assessment.Trans of IEE Japan, 1998, 118-B(2): 127~1326 Ernst D, Ruiz-Vega D, Pavella M, et al. A unified approach to transient stability contingency filtering,ranking and assessment. IEEE Trans on Power Systems, 2001, 16(3): 435~4437 邱志鹏,　邹云,　薛禹胜.　非线性系统同步稳定分析的互补簇簇际能量壁垒准则的数学描述.　非线性动力学学报, 2002, 9(3-4): 120~1318 邹云, 邱志鹏, 薛禹胜. 非自治动力学系统的同步稳定性. 应用数学学报, 2001, 24(1): 155~1579 Tang Y. Synchronous stability of multi-machine systems. Proceedings of the 3rd WSEAS Symposium onMath Methods and Comput Tech in Electrical Engineering, WSEAS Press, 2001. 5181~518410 薛禹胜. EEAC与直接法的机理: (一)受扰程度函数. 电力系统自动化, 2001, 25(11): 6~1111 薛禹胜. EEAC与直接法的机理: (二)壁垒点与观察点. 电力系统自动化, 2001, 25(12): 1~712 薛禹胜. EEAC与直接法的机理: (三)定性判稳与定量分析. 电力系统自动化, 2001, 25(13): 1~513 薛禹胜. EEAC与直接法的机理: (四)回顾与瞻望. 电力系统自动化, 2001, 25(14): 1~614 廖浩辉, 唐云. 电力系统中CCEBC/EEAC方法的数学描述. 电力系统自动化, 2003, 27(7): 12~1615 Univ California, Berkeley. Bifurcation and chaos in power systems: A survey, Aug 1992. EPRT TR-100834,Project 8010-1016 Chiang H D, Chu C C. Theoretical foundation of the BCU method for direct stability analysis ofnetwork-reduction power system models with small transfer conductances. IEEE Trans on Circuits and Systems—I: Fundamental Theory and Applications, 1995, 42(5): 252~26517 Chiang H D, Chu C C, Cauley G. Direct stability analysis of electric power systems using energy functions:Theory, applications, and perspective. Proceedings of the IEEE, 1995, 83(11): 1497~152918 Chiang H D. The BCU method for direct stability analysis of electric power systems: Theory andapplications. In: Chow J, Kokotovic P V, Thomas R J, eds. System Control Theory for Power Systems, 1995, 64: 39~9419 Chiang H D. Power system stability. In: John G W, ed. Wiley Enclyclopedia of Electrical and ElectronicEngineering. New York: John Wiley and Sons Inc, 1999。

分岔理论

⎧ x = 0和x = ± − μ ⎨ x=0 ⎩
chaos
chaos
μ <0 μ >0
(8-26) (8-27)
而对应本征值则为
λ = μ + 3x
2
如图8-14所示
chaos
chaos
图8-14 叉型分岔——亚临界情况
对于图8-12，令 μ < μ c 时，平衡态的一个分支( x = 0) 是稳定的；然而当 μ = μ c 时，这一支就变得不稳定了；一旦当 μ > μ c有新的平衡分支解 ( x = ± μ ) 又变成稳定的了，这种情况被称为超临界分岔。
chaos
chaos
Euler直杆弯曲满足下列非线性微分方程及边值 ⎧•• ⎪θ + μ sin θ = 0 ⎨• (8-1) • ⎪θ (0 ) = θ (1) = 0 ⎩ 方程(8-1)所表示的是一个本征值问题。当θ << 1 时，其对应的线性方程是
chaos
⎧•• ⎪θ + μθ = 0 ⎨• • ⎪θ (0 ) = θ (1) = 0 ⎩
∂f <0 ∂x
(8-13)
chaos
时平衡态稳定，而在
∂f >0 ∂x 时平衡态是不稳定的。
(8-14)
故对于
chaos
∂f =0 ∂x
(8-15)
正为由稳定变为不稳定的临界点。这个点我们即称之为实分岔点。再考虑到在分岔点处状态变量和参数的关系不唯一，进一步在实分岔点还有 ∂f (8-16) =0
chaos
chaos
是不稳定点
chaos
图8-11 x ~ α (t ) 注意到在图中只有稳定解才是按顺序相连的，如实线所示。而图上的虚线反映解在分支点的邻域，状态可以发生突变。

非线性动力学之一瞥——Lorenz系统

非线性动力学非线性系统之一瞥——Lorenz系统2013-01-300 前言非线性系统动力学线性系统是状态变量和输出变量对于所有可能的输入变量和初始状态都满足叠加原理的系统；非线性系统就是这些量不满足叠加原理的系统。

非线性系统在日常生活和自然界中不胜枚举，也远远多于线性系统。

非线性动力学是研究非线性系统的各种运动状态的定性和定量变化规律，尤其是系统的长时期行为。

研究的对象主要有分叉、混沌和孤立子等。

洛伦兹方程洛伦兹方程是美国气象学家洛伦兹在模拟天气这一非周期性现象时确定，这个方程的三个变量分别模拟温度、湿度和压力。

可以得出结论，初期微小的差别随着时间推移差别会越来越大，洛伦兹基于此提出长期的天气预报是不可能的。

这也被视为研究非线性混沌理论的开始，所以洛伦兹系统在研究非线性系统中具有举足轻重的地位。

本文借助洛伦兹系统对非线性进行简单的介绍。

洛伦兹方程如下。

方程中，、和都为实参数。

实参不同，系统的奇点及数目也是不同的。

1 奇点和稳定性奇点洛伦兹系统含有三个实参数，当参数变化，奇点的数目可能不同。

首先，一定是系统的奇点。

时，当时，系统仅有一个奇点；当时，系统还有另外两个奇点。

下面仅解时的两个非原点奇点。

令方程第一式得，第三式可得，将两式代入第二式得即，。

奇点稳定性判别下面根据Liapunov稳定性判别方法，找出系统在原点处大范围渐进稳定的条件，取Liapunov函数。

考虑，的情况。

则有将洛伦兹方程代入上式，可得变换为二次型，系数矩阵为已知，，则系数矩阵负定的条件是。

所以该系统是大范围渐进稳定的条件是，前提是，。

Liapunov函数V总是存在的，只要构造出合适的Liapunov函数，就可以通过Liapunov稳定性定理直接判断奇点的稳定性，而不需要求解非线性方程组。

有的Liapunov函数不易构造，则可以通过奇点处导算子的特征值来判断：若所有的特征值实部都小于0，则方程组在该奇点是局部渐进稳定的；若特征值实部至少有一个为正，该奇点是不稳定的。

第五章随机平均法非线性随机动力学教案 .docx

El。；］=『/卩(4，%，卩1，卩2)呦1勿2如眇2
□ 015
(a)
Scห้องสมุดไป่ตู้
(b)
图5. 2-1例4.2-1系统（a）平稳响应。S]=S?=二0.3 ,其它参数与图4.2-1中相同。（a）系统总能量的概率密度；（b）位移均方值。
—随机平均法结果；---第一种准则等效非线性系统方法结果；••数字模拟结果。
5.1
随机平均原理是随机平均法的严格数学基础。本书中所发展的拟HamUton随机平均法用到两种形式随机平均原理，一是被经常引用的Stratonovich-Khasminskii极限定理，该定理乃由Stratonovich⑹基于物理考虑提出，然后Khasminskii⑺为该定理提供了严格的数学提法与证明，Papanicolao与Kohle严贝II对该定理作了改进与引申。该定理的数学提法较为一般，此处仅限于本书用到的特殊形式。
立，此时(5.1-2)化为
% (兀)=”(X, 0 +Dki笃;"Sji(兀,0bij(x)=Qd禺kgt'gjid,叽
若进一步假定/.、緞不显含丫，则(5.1-7)化为
叫(兀)= /(*) +Dki气“gji (*) bij(x) = 2Dklgik(x)gjl(x)
第一式中的第二项即为Wong-Zakai修正项。
再设以H为Hamilton函数的Hamilton系统为不可积，即H是与（522）相应的HamUton系统的唯一独立首次积分。弓【入变
H=H(QQ(5.2-3)
应用It6微分公式(261),可由(5.2-2)导得HamUton过程H(f)所满足的It6随机微分方程
以(524)代替(5.2-2)中关于P}的方程，并在(522)的其余方程及(5.2-4)中，按(5.2-3)以H代替円。这组新方程形同(5.1-9),的)，AU)，…，几⑺为快变过程，而H(f)为慢变过程。根据Khasminskii定理在etO时，在訂量级时间区间上，刊了)弱收敛于一维扩散过程。仍以丹⑴表示这一极限扩散过程，则支配该过程的平均It6随机微分方程形为

双星形成的非线性动力学机制和定性分析理论

一
ＬＫｌｉｅｎ就提出双星因 “ 转裂变 ” 造成．是ｖ旋而但
Ｄｕｉｎ等Ｊｒｅｓ已经证明裂变在理论上行不通
目前，星形成的最成功的理论之一，由双是Ｂ等提出的碎裂理论．１８ｃ从９６年起，他们对此展开ｊ长期的一系列研究＿．假定双星在原始星，６其Ｊ体塌缩过程中．自身引力作用下形成计算机模在拟表明最初的球形星云，迅速旋转中首先塌缩为在
当Ｊａｓ等式＞（＇（成立时．ｅｎ不／）引力不稳定，且星云将塌缩我们模型的基础是恒并星起源于星云，主要由氢和氮等离子体组成的星而云服从非线性的磁流体动力学方程．它们的一般形
收璃日期２０ —０ —２０１４９基金项目：南省自然科学基盎资助项目（０１，）０１）云２（４０２ｔ０／％作者简彳：一方（９６）只，教授．要从事理论物理研究ｒ张１４一，副主
双星形成的非线性动力学机制和定－性分析理论 ’
张一方，李艳梅
（云南大学物理系，１云南昆明６０９；５０１２云南大学数学系，云南昆明６０９）５０１
摘要：把星云的复杂的磁流体动力学方程，似简化为只与速度相关的方程，近则得到一个可解的简单模型应用非线性方程的定性分析理论．获得双星演化的非线性动力学机制一定条件下，在二维空间中的某些奇点捕应于双星．而在别的条件下，奇点相垃于单星因此．各种星的形成是星云演化的结果美键词：双星；流体动力学；非线性方程；化；演星云