分叉计算的若干问题.pptVIP

分岔计算的若干问题 武际可 (WU Jike) 北京大学力学与工程科学系 Department of Mechanics and Engineering Sciences Peking University §1 前言§2 弧长法与分岔点的判断§3 动力系统分岔的定义§4 非自治动力系统的等价分类§5 同宿(异宿)轨道的寻求§6 求解动力系统的伪弧长法  §1 前言 著名理论物理学家、诺贝尔奖金获奖者费音曼（Richard Phillips Feynmann,1918－1988）在他的《物理学讲义》（Lectures on Physics）中的一段话：“我们已经写下水流的方程，用实验的方法，我们发现了一系列用于理解的近似概念――涡街、湍流尾迹、边界层。如果我们有稍许不同形式的类似方程，而且我们没有法子对这方程作实验，我想以哪怕是原始的、可怀疑的、含混的方式去求解这一方程，以求确定数量性态或数量结果。例如，我们的方程是对于太阳，把它看为无太阳黑子、表面无硬块的结构、没有凸凹不平的氢气球。尽管作了这样的简化，我们仍然还未找到求解的方法。…。人类智慧觉醒的下一世纪，也许可能产生了解方程定量含义的新方法。” 空穴的萌生 假设 ，即弹性球的内半径趋于零。并且考虑，即在某个p作用下空穴半径从零突然增加到非零。在这样的条件下解有关方程，就会得到空穴萌生的条件。我们对这个问题进行了分析求解，并且对相应的分析表达式进行了数值计算，对于不同的值和不同的v值，得到值的结果表示在图上。从图上可以看出，当 时，趋于一个确定的值。这种情况正好是空穴从半径为零突然产生的临界情形。 §2 弧长法与分岔点的判断 用伪弧长法追踪非线性方程组的解曲线：对于有限维静力状态方程，若记 ，忽略状态变量和参变量的差别，则原来方程改变为： 对(1)微分得 这里 因而 若干算例利用弧长法可以没有困难地计算材料软化、塑性流动曲屈后行为等一系列以前认为困难的问题。 §3 动力系统分岔的定义 §5 同宿(异宿)轨道的寻求 §6 求解动力系统的伪弧长法 对于非自治系统 我们最近证明了如下的结果：在非自治系统（10）中，若矩阵A(t) 和它的逆矩阵都连续且有界，它的特征值的实部有 个恒大于一 个正实数，有 个恒小于一个负实数 .而且 它的全部特征值满足如下条件：如果存在正常数 ，满足： 使对于充分大 的 或者 则（10）拓扑等价于(9)。 的一个解是 ，对这个解做小扰动 代入原方程得 把这两个方程相减，并且对右端的 线性化，就会得到一个关 于 的线性非自治系统: 在有了关于线性非自治系统等价类的结果以后，我们还可以讨论 动力系统的一个具体的解是否是分岔解。设动力系统 的分岔定义，如果解 有分岔点，我们就可以说 是一个分岔解。 仿照定义4可以给出关于 的一个以T为周期的闭轨 如果我们能够判断在相应的闭轨(12)上有一个（或多个）奇点 即 则这个闭轨就是同宿（或异宿）轨。 一般认为，同宿和异宿轨道的出现是动力系统从规则运动转化为混沌的通道。所以确定动力系统在什么条件下出现同宿和异宿轨道，对于讨论系统向混沌运动的转化具有十分重要的意义。 设我们已经追踪到动力系统 令 为闭轨上的点， 为奇点，把(13)在 展开，就得到 对上式补充一个方程 在上面 个方程中， 是非零向量，所以系数行列 式必然满足： 再代入(14)就可以得到 于是就可以判断 是否在闭轨上，从而确定已知的闭轨是否是同 宿(异宿)轨。 上面的求解过程是在参数 给定的条件下来求解的。所得到 的奇点,不一定在闭轨上，也不一定唯一。但是我们可以大致了解 奇点与闭轨的位置随 变化情况从而确定闭轨上奇点对应的 下面接着讨论在 变化的条件下确定同宿异宿轨的方法，令 为闭轨上的点， 对应于同宿(异宿)轨上的奇点,把 (13)在点 展开，就得到 对上式求解就可以得到闭轨上的 对上式补充两个方程 在上面 个方程中， 是非零向量，所以 系数行列式必然满足： 闭轨上的 再代入（12）就可以得到 于是就可以判断 是否在闭轨上，从而确定是否达到了同宿 (异宿)轨。 把上式与追踪闭轨的弧长方程联立求解就可以得到对应于 这里 文献[1]已经详细讨论了对于求解静力问题的弧长法。至于对动力问题，即带时间的微分方程的求解，有时引进弧长可以克服许多困难，使本来难于求解的问题变得容易。现在就来介绍我们在[9]中发展的弧长法。 考虑一个常微分方程组的初值问题