数学建模 6.6 微分方程稳定性理论简介.DOCVIP

6 . 6 微分方程稳定性理论简介一阶方程的平衡点及稳定性设有微分方程(1) 右端不含字变量t ，称为自治方程. 代数方程f(x) = 0 (2) 的实根x = x0 称为方程(1) 的平衡点( 或奇点). 它也是(1) 的解( 奇解). 如果存在某个邻域，使方程(1) 的解x(t) 从这个邻域内的某个x(0) 出发，满足(3) 则称平衡点x0 是稳定的( 稳定性理论中称渐进稳定); 否则，称x0 是不稳定的( 不渐进稳定). 判断平衡点x0 是否稳定通常有两种方法. 利用定义即(3) 式称间接法. 不求方程(1) 的解x(t) ，因而不利用(3) 式的方法称直接法. 下面介绍直接法. 将f(x) 在x0 点作Taylor 展开，只取一次项，方程(1) 近似为, (4) (4) 称为(1) 的近似线性方程，x0 也是方程(4) 的平衡点. 关于x0 点稳定性有如下结论：若f (x0) 0 ，则x0 对于方程(4) 和(1) 都是稳定的; 若f (x0) 0 ，则x0 对于方程(4) 和(1) 都是不稳定的. x0 点对方程(4) 稳定性很容易由定义(3) 证明：记f (x0) = a ，则(4) 的一般解为x(t) = ceat + x0 (5) 其中常数c 由初始条件确定，显然，a 0 时(3) 式成立. 二阶方程的平衡点和稳定性二阶方程可用两个一阶方程表为(6) 右端不显含t ，是自治方程. 代数方程组(7) 的实根x1 = x10, x2 = x20 称为方程(6) 的平衡点，记作P0(x10, x20). 如果存在某个邻域，使方程(6) 的解x 1 (t), x2(t) 从这个邻域内的某个(x1(0), x2(0)) 出发，满足(8) 则称平衡点P0 是稳定的( 渐进稳定); 否则，称P0 是不稳定的( 不渐进稳定). 先看线性常系数方程(9) ( 非齐次方程组，可用平移的方法(x1 = u1 + c1, x2 = u2 + c2) 化为齐次方程组) 系数矩阵记作(10) 为研究方程(9) 的唯一平衡点P0(0, 0) 的稳定性，假定A的行列式detA ( 0 (11) P0(0, 0) 的稳定性由(9) 的特征方程det(A ( (I) = 0 (12) 的根(( 特征根) 决定. 方程(12) 可以写成更加明晰的形式(13) 将特征根记作(1, (2 ，则(14) 方程(9) 的一般解具有形式或，c1, c2 为任意常数. ( 注：课本p 199误为) 按照稳定性的定义(8) 式可知，当(1, (2 均为负数或均有负实部时P0(0, 0) 是稳定平衡点; 而当(1, (2 有一个为正数或有正实部时P0(0, 0) 是不稳定平衡点. 在条件(11) 下(1, (2 均不为零. 按上述理论可得根据特征方程的系数p, q 的正负来判断平衡点稳定性的准则：若p 0, q 0 ，则平衡点稳定; 若p 0, q 0 ，则平衡点不稳定. 对一般的非线性方程(6) ，仍可在平衡点作一次Taylor 展开，得常系数的近似线性方程来讨论. 6 . 2 军备竞赛两个国家或国家集团之间由于相互不信任和各种矛盾的存在、发展而不断增加自己的军事力量，防御对方可能发动的战争. 本节介绍L. F. Richardson 1939年提出的一个模型. 模型的假设和构成为了方便起见，用军备表示军事力量的总和，如兵力、装备、军事预算等. 甲乙双方在时刻t 的军备分别记作x(t) 和y(t) ，假设它们的变化只取决于下面3个因素：1 . 由于相互不信任及矛盾的发展，一方军备越大，另一方军备增加得越快; 2. 由于各方本身经济实力的限制，任一方军备越大，对军备增长的制约作用越大; 3 . 由于相互敌视或领土争端，每一方都存在着增加军备的固有潜力. 进一步假定前两个因素的影响是线性的，第3个因素的影响是常数，那么x(t) 和y(t)的变化过程可用微分方程组(1) 表示，其中的系数均大于或等于零. k, l 是对方军备刺激程度的度量; ( , (是己方经济实力制约程度的度量; g, h 是己方军备竞赛的固有潜力. 如果我们感兴趣的是军备竞赛的结局由什么因素决定，而不关心竞赛的过程，那么只需用微分方程稳定性理论讨论时间充分长以后x ( t ) , y ( t )的变化趋势，即方程(1)的平衡点的稳定情况. 令(1)式右端等于零，容易算出平衡点( x0, y0 )为(2) 方程(1)的系数矩阵为于是按照判断平衡点稳定性的方法计算(见6 . 6节(9)～(13)式)