一类微分方程模稳定性的数值模拟 毕业论文.docVIP

一类微分方程模型稳定性的数值模拟 摘要：数值模拟也叫计算机模拟。数值模拟20世纪以来，随着大量的边缘科学诸如电磁流体力学、化学流体力学、动力气象学、海洋动力学、地下水动力学等等的产生和发展，以来，国际数值模拟研究取得显著发展，它已不仅是一种现代化的实验手段，而且已发展为具有独立特征的学科分支。在以计算机为基础的技术带动下将的变化要考虑 （1） 右端不显含t，代数方程组 （2） 的实根称为方程（1）的平衡点，记为。如果从所有可能的初始条件出发，方程（1）的解都满足 （3） 则称平衡点是稳定的（渐近稳定）；否则，称P0是不稳定的（不渐近稳定）。为了用直接法讨论方法方程（1）的平衡点的稳定性，先看线性常系数方程： （4） 系数矩阵记作： 并假定A的行列式。于是原点是方程（4）的唯一平衡点，它的稳定性由的特征方程：的根（特征根）决定，上方程可以写成更加明确的形式： （5） 将特征根记作，则： （6） 方程（4）的解一般有形式（）或（），为任意实数。由定义（3），当全为负数或有负的实部时是稳定的平衡点，反之，当有一个为正数或有正的实部时是不稳定的平衡点。微分方程稳定性理论将平衡点分为结点、焦点、鞍点、中心等类型，完全由特征根或相应的取值决定，下表简明地给出了这些结果，表中最后一列指按照定义（3）式得下列关于稳定性的结论。 表1 ：由特征方程决定的平衡点的类型和稳定性 平衡点类型 稳定性 稳定结点 稳定 不稳定结点 不稳定 鞍点 不稳定 稳定退化结点 稳定 不稳定退化结点 不稳定 稳定焦点 稳定 不稳定焦点 不稳定 中心 不稳定 由上表可以看出，根据特征方程的系数的正负很容易判断平衡点的稳定性，准则如下：若A1：则平衡点稳定。若A2： 则平衡点不稳定。以上是对线性方程（4）的平衡点稳定性的结论，对于一般的非线性方程（1），可以用近似线性方法判断其平衡点的稳定性，在点将和作泰勒展开，只取一次项，得（1）的近似线性方程 (7) 系数矩阵记作： 特征方程系数：, 。显然，点对于方程（7）的稳定性由表1或准则A1、A2决定，而且已经证明了如下结论: 若方程（7）的特征根不为零或实部不为零，则点对于方程（1）的稳定性与对于近似方程（7）的稳定性相同。这样，点对于方程（1）的稳定性也由准则A1、A2决定。 2.2 微分方程与极限环及其稳定性 设系统 （8） 具有闭轨线C。假如在C充分小邻域中，除C之外，轨线全不是闭轨线，且这些非闭轨线当t→＋∞或t→－∞时趋近于闭轨线C，则说闭轨线C是孤立的，并称之为(15)的一个极限环.极限环C将相平面分成两个区域：内域和外域。 如果极限环C的内域的靠近C的轨线当t→+∞(-∞)时盘旋地趋近于C，则称C是内稳定(内不稳定的)；如果在极限环C的外域的靠近C的轨线当t→+∞(－∞)时盘旋地趋近于C，侧称C是外稳定的(外不稳定的)；如果当t→＋∞(－∞)时，C的内部及外部靠近C的轨线都盘旋地趋近于C，则称C是稳定的(不稳定的),如果当t→＋∞(－∞)时，C的内外部的稳定性相反，则称C为半稳定的。 3 一类微分方程模型稳定性的数值模拟 3.1 模型: 相互竞争的种群模型 1、模型假设及说明: （1）、假设有甲乙两个种群，当它们独自在一个自然环境中生存时，数量的演变均遵从Logistic规律，记,是两个种群的数量，，是它们的固有增长率，,是它们的最大容量。 （2）、当甲、乙两个种群在同一自然环境中生存时，乙对甲增长的阻滞作用与乙的数量成正比；甲对乙有同样的作用。 （3）、：单位数量乙（相对于而言）消耗的供养甲的食物量为单位数量甲（相对于而言）消耗的供养甲的食物量的倍。表示的意义与之相反。 （4）、可解释为相对于而言单位数量的甲消耗的供养甲的食物量（设食物总量为1）。可做类似解释。 2、模型的建立: 3、稳定性分析：利用平衡点的稳定性分析，