分叉计算的若干问题.pptxVIP

§1 前言§2 弧长法与分岔点的判断§3 动力系统分岔的定义§4 非自治动力系统的等价分类§5 同宿(异宿)轨道的寻求§6 求解动力系统的伪弧长法 §1 前言著名理论物理学家、诺贝尔奖金获奖者费音曼（Richard Phillips Feynmann,1918－1988）在他的《物理学讲义》（Lectures on Physics）中的一段话：“我们已经写下水流的方程，用实验的方法，我们发现了一系列用于理解的近似概念――涡街、湍流尾迹、边界层。如果我们有稍许不同形式的类似方程，而且我们没有法子对这方程作实验，我想以哪怕是原始的、可怀疑的、含混的方式去求解这一方程，以求确定数量性态或数量结果。例如，我们的方程是对于太阳，把它看为无太阳黑子、表面无硬块的结构、没有凸凹不平的氢气球。尽管作了这样的简化，我们仍然还未找到求解的方法。…。人类智慧觉醒的下一世纪，也许可能产生了解方程定量含义的新方法。”Richard Phillips Feynmann1918－1988分岔是非常普遍的现象。分岔是系统两种性质上不同的状态转变点。掌握系统的分岔点对把握系统的性质和行为非常重要。分岔问题是实质非线性问题。分岔的研究目前引起所有的可以精确化的科学的共同兴趣。 中青海省桥头发电厂六期扩建工程施工现场1998年9月19日上午发生一起重大伤亡事故，一座正在修建中的冷却塔在修至20米高处时，用于浇筑混凝土的钢模板突然全部坍塌脱落，致使正在施工的数十名工人全部从空中跌落下来，造成4人死亡，52人受伤。目前，事故原因正在调查中。图为事故现场。 中新社深圳2000年十一月二十八日电：正在施工中的深圳盐坝高速公路起点高架引桥昨晚九时左右突然坍塌，在桥面作业的六十九名工人随桥面滚落坠下，截至今天凌晨一时三十分，所有人员均被救出，三十三人被送医院急救，其中十人重伤，二十三人轻伤。图为桥梁倒塌现场(陈文 摄) 。初步认定事故是因技术问题造成建桥支架失稳引发的。空穴的萌生 假设 ，即弹性球的内半径趋于零。并且考虑，即在某个p作用下空穴半径从零突然增加到非零。在这样的条件下解有关方程，就会得到空穴萌生的条件。我们对这个问题进行了分析求解，并且对相应的分析表达式进行了数值计算，对于不同的值和不同的v值，得到值的结果表示在图上。从图上可以看出，当 时，趋于一个确定的值。这种情况正好是空穴从半径为零突然产生的临界情形。ba§2 弧长法与分岔点的判断1. 武际可，苏先樾. 弹性系统的稳定性, 科学出版社,1994.2. 周 鲲. 动力系统Hopf分岔的数值计算及应用,硕士论文.3. 武际可, 黄克服. 分岔的数值方法, 待出版. 用伪弧长法追踪非线性方程组的解曲线：对于有限维静力状态方程，若记 ，忽略状态变量和参变量的差别，则原来方程改变为： 而且已知有一点 是满足（1）的。对(1)微分得若记若令这里表示划去该列. 显然, 这个辅助参数,并定义 引进弧长 空间中一维曲线的弧长由 来确定。因而 是(1)式解曲线 的切向量, 是其单位切向量. 于是(1)式的求解可以化归为求解类似地, 可将上式转化为常微分方程的Cauchy问题： (2)数值求解(2)式的最简单方式是Euler预报法，在求解过程中每前进一步或若干步后，可以做Newton型迭代修正，联合使用这两种方法，就形成了追踪解流形的一种有效算法，即预报—修正法过程的延续算法。考虑动力系统 是它的平衡解即 对给定的动力系统作微小扰动 u 有把上面的两个式子相减就可以得到如果u有非零解，就说 是动力系统的分岔点矩阵A,B的特征值具有如下关系： 它将A的特征值所在复平面上的虚轴映射成了B的特征值所在复平面上的特征圆. 判断和确定解曲线上的分岔点:为了寻求一种适于高维问题的算法来判断和确定分岔点，我们可以引进一个变换：因此, 在初始状态(稳定状态), 的特征值在复平面 的左半平面上, 而对应的 的特征值均在复平面的单位圆内. 这样,通过变换将寻求系统的分岔点的问题转化为用乘幂法（文献[1]）求矩阵 的具最大模的共轭复特征值的问题. 如果我们搜索出具有最大的点并能进一步判断 即对应 ，则无疑 模 对应的点 就是静分岔点；否之则为Hopf分岔点。由于在实际计算中很难恰有 因此我们需要借助于变号法来搜索解曲线上的分岔点。1. 武际可，苏先樾. 弹性系统的稳定性, 科学出版社,1994.(3).分岔方向的寻求：对于静分岔方向的寻求，这里仅给出简单分叉，即当时的数值分析方法。根据定理5.5.1(文献[1]),可以把求分岔方向化为线性特征值问题其中，是的单位正交基底。的两个单位向量 ，为了