1、第四节 微分方程稳定性理论简介这里简单介绍下面将要用到的有关内容：一、 一阶方程的平衡点及稳定性设有微分方程（1）()dxft右端不显含自变量 t，代数方程（2）()0f的实根 称为方程（1）的平衡点（或奇点） ，它也是方程（1）的解（奇解）0x如果从所有可能的初始条件出发，方程（1）的解 都满足()xt（3）0lim()tx则称平衡点 是稳定的（稳定性理论中称渐近稳定） ；否则，称 是不稳0 0x定的（不渐近稳定） 。判断平衡点 是否稳定通常有两种方法，利用定义即（3）式称间接法，不0x求方程（1）的解 ，因而不利用（3）式的方法称直接法，下面介绍直接法。()t将 在 做泰勒展开，只取一次项
2、，则方程（1）近似为：()fx0（4）0()dfxt（4）称为（1）的近似线性方程。 也是（4）的平衡点。关于平衡点0x的稳定性有如下的结论：0x若 ，则 是方程（1） 、 （4）的稳定的平衡点。0()f0x若 ，则 不是方程（1） 、 （4）的稳定的平衡点x对于方程（4）的稳定性很容易由定义（3）证明，因为（4）的一般解是0（5）0()fxttce其中 C 是由初始条件决定的常数。二、 二阶（平面）方程的平衡点和稳定性方程的一般形式可用两个一阶方程表示为（6）1122(),),dxtfxgt右端不显含 t，代数方程组（7）12(,)0fxg的实根 称为方程（6）的平衡点。记为012(,)x
3、012(,)Px如果从所有可能的初始条件出发，方程（6）的解 都满足2(,)t（8）10lim()tx20lim()tx则称平衡点 是稳定的（渐近稳定） ；否则，称 P0是不稳定的（不02,P渐近稳定） 。为了用直接法讨论方法方程（6）的平衡点的稳定性，先看线性常系数方程（9）1122()dxtabxt系数矩阵记作 12abA并假定 A 的行列式 det0于是原点 是方程（9）的唯一平衡点，它的稳定性由的特征方程0(,)PtI的根 （特征根）决定，上方程可以写成更加明确的形式：（10）2120()detpqabA将特征根记作 ，则12,（11）2(4)pq方程（9）的解一般有形式 （ ）或 （
4、12ttce1212()tce）12为任意实数。由定义（8） ，当 全为负数或有负的实部时 是稳,c 12, 0(,)P定的平衡点，反之，当 有一个为正数或有正的实部时 是不稳定的平12,衡点微分方程稳定性理论将平衡点分为结点、焦点、鞍点、中心等类型，完全由特征根 或相应的 取值决定，下表简明地给出了这些结果，表中最后12,pq一列指按照定义（8）式得下马看花关于稳定性的结论。表 1 由特征方程决定的平衡点的类型和稳定性12,pq平衡点类型 稳定性20,稳定结点 稳定12pq不稳定结点 不稳定0鞍点 不稳定1220,稳定退化结点 稳定pq不稳定退化结点 不稳定12, 2,稳定焦点 稳定0不稳定
5、焦点 不稳定12,pq中心 不稳定由上表可以看出，根据特征方程的系数 的正负很容易判断平衡点的稳定,pq性，准则如下：若（12）0,pq则平衡点稳定，若（13）或则平衡点不稳定以上是对线性方程（9）的平衡点 稳定性的结论，对于一般的非线0(,)P性方程（6） ，可以用近似线性方法判断其平衡点 的稳定性，在012(,)x点将 和 作泰勒展开，只取一次项，得（6）的近似线012(,)Px12(,)fx12(,)gx性方程(14)1 21 20000121122(),)(,), ,xxdtffxggt系数矩阵记作 12012(,)xPxfAg特征方程系数为,01212()(xPxpfdetqA显然，
6、 点对于方程（14）的稳定性由表 1 或准则（12） 、 （13）决012,)P定，而且已经证明了如下结论:若方程（14）的特征根不为零或实部不为零，则 点对于方程（6）012(,)Px的稳定性与对于近似方程（14）的稳定性相同。这样， 点对于方程（6）的稳定性也由准则（12） 、 （13）决定。012(,)Px第五节 种群的相互竞争与相互依存当某个自然环境中只有一种生物的群体（生态学上称为种群）生存时，人们常用 Logistic 模型来描述这个群数量的演变过程，即(1)(1)dxrtNx（t）是种群在时刻 t 的数量， 是固有增长率，N 是环境资源容许的种群r最大数量，在前面我们曾应用过这种
7、模型，由方程（1）可以直接得到， =N0x是稳定平衡点，即 t时 x(t)N，从模型本身的意义看这是明显的结果。如果一个自然环境中有两个或两个以上种群生存，那么它们之间就要存在着或是相互竞争，或是相互依存，或是弱肉强食（食饵与捕食者）的关系。这里将从稳定状态的角度分别讨论这些关系。一、种群的相互竞争当两个种群为了争夺有限的食物来源和生活空间而进行生存竞争时，最常见的结局是竞争力较弱的种群灭绝，竞争力较强的种群达到环境容许的最大数量。人们今天可以看到自然界长期演变成的这样的结局，例如一个小岛上虽然有四种燕子栖息，但是它们的食物来源各不相同，一种只在陆地上觅食，另两种分别在浅水的海滩上和离岸稍远的
8、海中捕鱼，第四种则飞越宽阔的海面到远方攫取海味，每一种燕子在它各自生存环境中的竞争力明显地强于其它几种，这里我们建立一个模型解释类似的现象，并分析产生这种结局的条件。模型建立 有甲乙两个种群，当它们独自在一个自然环境中生存时，数量的演变均遵从 Logistic 规律，记 是两个种群的数量， 是它们的12(),xt 12,r固有增长率，N 1、N 2是它们的最大容量， 于是对于种群甲有 1()dxrt其中因子 反映由于甲方有限资源的消耗导致的对它本身增长的阻滞1()作用， 可解释为相对于 N1而言单位数量的甲消耗的供养甲的食物量（设食1xN物总量为 1） 。当两个种群在同一自然环境中生存时，考察
9、由于乙消耗同一种有限资源对甲的增长产生的影响，可以合理地在因子 中再减去一项，该项与种群乙1()xN的数量 （相对于 N2而言）成正比，得到种群甲方增长的方程2x(2)112()dxrt这里 的意义是，单位数量乙（相对 N2而言）消耗的供养甲的食物量为单1位数量甲（相对 N1）消耗的供养甲的食物量的 倍。1类似地，甲的存在也影响了乙的增长，种群乙的方程应该是(3)2122()dxxrtN对 可作相应的解释。2在两种群的相互竞争中 、 是两个关键指标，从上面对它们的解释可知，12 1 表示在消耗供养甲的资源中，乙的消耗多于甲，因而对甲增长的阻滞作用乙大于甲，即乙的竞争力强于甲，对 1 可作相应的
10、理解。2一般地说， 与 之间没有确定的关系，但是可以把下面这种特殊情况作12为较常见的一类实际情况的典型代表，即两个种群在消耗资源中对甲增长的阻作用对乙增长的阻滞作用相同，具体地说就是，因为单位数量的甲和乙消耗的供养甲方食物量之比是 1： ，消耗的供养甲方食物量之比是 ：1，所谓阻2滞作用相同即 1： = ：1，所以这种特殊情形可以定量地表示为2=1 （4）1即 、 互为倒数，可以简单地理解为，如果一个乙消耗的食物是一个甲12的 = 倍，则一个甲消耗的食物是一个乙的 =1/ 。k 2k下面我们仍然讨论 、 相互独立的一般情况，而将条件（4）下对问题12的分析留给大家讨论。稳定性分析 为了研究两
11、个种群相互竞争的结局，即 t时 的12(),xt趋向，不必要解方程（2） 、 （3）,只需对它的平衡点进行稳定性分析。首先根据微分方程（2） 、 （3）解代数方程组（5）12121212(,)()0,xfxrNg得到 4 个平衡点：12123 41()()(,0)(,),(0)NPNPP因为仅当平衡点们于平面坐标系的第一象限时（ ）才有实际意义，12,0x所以对 而言要求 、 同时小于 1，或同时大于 1。3P12按照判断平衡点性的方法（见前面）计算 1212122211()()xxrxrfNNAg 12(),34ixPpfdetiqA将 4 个平衡点 p、q 的结果及稳定条件列入下表 *）表
12、 1 种群竞争模型的平衡点及稳定性平衡点 p q 稳定条件1(,0)PN12()r12()r12,2 2r11231()(),112()()r1221()r12,4(0,)P()2不稳定注：表中最后一列“稳定条件”除了要求 p0,q0 以外，还有其他原因，见下面的具体分析。为了便于对平衡点 P1、P 2、P 3的稳定条件进行分析，在相平面上讨论它们。在代数方程组（5）中记1212(,)0xxN12122(,)对于 、 的不同取值范围，直线 =0 和 =0 在相平面上的相对位置不12 同，下面给出它们的 4 种情况;并对这 4 种情况进行分析1、 。由表 1 知对于 有 0， 0， 稳定； 的稳
13、2,1(,)PNpq1P1定性还可以从 t时相轨线的趋向来分析，图 1 中 =0 和 =0 两条直线将相平面（ ）划分为 3 个区域：120,x12/N1N1xPO2N21/21S23S00（1） 2,稳定1P图 1 稳定12,1P(6)12:/0,/Sdxtdt(7)2(8)312:/,/xtt可以证明，不论轨线从哪个区域出发，t时都将趋向 P1（N 1，0） 。若轨线从 S1出发，由（6）可知随着 t 的增加轨线向右上方运动，必然进入 S2；若轨线从 S2出发，由（7）可知轨线向右下方运动，那么它或者趋向 点，1或者进入 S3，但是进入 S3是不可能的，因为，如果设轨线在某时刻 t1经直线
14、=0 进入 S3，则 (t1)/ =0,由方程（2）不难算出dxt212rdxtNt由（7） 、 （8）知 0, 故 ，表明 (t)在 t1达到极小值，2/xt21/0xt1x而这是不可能的，因为在 S2中 0,即 (t)一直是增加的；d若轨线从 S3出发，由（8）可知轨线向左下方运动，那么它或者趋向 点，1P或者进入 S2，而进入 S2后，根据上面的分析最终也将趋向 。1P综上分析可以画出轨线示意图（图 1） ，因为直线 =0 上 d =0，所以在x=0 上轨线方向垂直于 轴；在 =0 上 d =0,轨线方向平行于 轴。1x2x12、 ，类似的分析可知 稳定。12,(0,)PN1N12/N1
15、x2PO21/N2x00（2） 1,稳定2P图 2 稳定12,2P3、 ，由表 1 知对于 点 0， 0，故 稳定，对轨线趋12,3Ppq3势的分析见图 3。1N12/N1x3PO21/N2x00（3） 12,稳定3P图 3 稳定12,3P4、 ，由表 1 知对于 点 0，故 不稳定（鞍点） ，轨线或12,3Pq者趋向 ，或者趋向 ，由轨线的初始位置决定，示意图见图 4，在这种情况P下 和 都不能说是稳定的，正因为这样，所以 稳定（与初始条件无关）的12 1条件需要加上 ， 稳定的条件加上 。12212/N1x3PO21/N2x00（4） 12,不稳定3P1N图 4 不稳定12,3P结果解释
16、根据建模过程中 的含义，说明 、 、 点稳定在生态12,1P23上的意义。1、 ， 意味着在对供养甲的资源的竞争中乙弱于甲，2,1意味着在对供养乙的资源的竞争中甲强于乙，于是种群乙终灭绝，种群2甲趋向最大容量，即 趋向平衡点12(),xt1(,0)PN2、 ，情况与 1 正好的相反。12,3、 ，因为在竞争甲的资源中乙较弱，而在竞争乙的资源中甲较弱，于是可以达到一个双方共存的稳定的平衡状态 ，这是种群竞争中很少3P出现的情况。4、 ，请大家作出解释。12,生态学中有一个竞争排斥原理；若两个种群的单个成员消耗的资源差不多相同，而环境能承受的种群甲的最大容量比种群乙大，那么种群乙终将灭亡，用本节的
17、模型很容解释这个原理。将方程（2） 、 （3）改写为1121()Nxdxrt212()xxrdtN原理的两个条件相当于 12122,从这 3 个式子显然可得 ，这正是 稳定，即种群乙灭绝的条件。12,1P二、种群的相互依存自然界中处于同一环境下两个种群相互依存而共生的现象是很普遍的，植物可以独立生存。昆虫的的授粉作用又可以提高植物的增长率，而以花粉为食物的昆虫却不能离开植物单独存活，人类与人工饲养的牲畜之间也有类似的关系，这种共生现象可以描述如下。设种群甲可以独立存在，按 Logistic 规律增长，种群乙为甲提供食物，有助于甲的增长，类似于前面的方程（2） ，种群甲的数量演变规律可以写作（r
18、 1、N 1、N 2的意义同前）(9)112()dxxrtN前面的 号这里变成+号，表示乙不是消耗甲的资源而是为甲提供食物，1的含义是: 单位数量乙（相对于 N2）提供的供养甲的食物量为单位数量甲（相对于 N1）消耗的供养甲食物量的 倍。1种群乙没有甲的存在会灭亡，设其死亡率为 r2，则乙单独存在时有(10)22/dxtr甲为乙提供食物，于是（2）式右端应加上甲对乙增长的促进作用，有(11)1222/()xtrN显然仅当 时种群乙的数量才会增长，与此相同乙的增长又会受到12x自身的阻滞作用，所以 93）式右端还要添加 Logistic 项，方程变为(12)12222/()xdxtrN方程（9）
19、 、 （12）构成相互依存现象的数学模型，下面利用平衡点的稳定性分析，讨论时间足够长以后两个种群的变化趋向。类似于前面的作法将方程（9） 、 （12）的平衡点及其稳定性分析的结果列入表 2表 2 种群依存模型的平衡点及稳定性平衡点 pq稳定条件1(,0)PN12()r12()r212,1221()(),121()()r1212()12,3(0,)P12r12r不稳定显然， 点稳定才表明两个种群在同一环境里相互依存而共生，我们着重2P分析 稳定的条件。 。由 的表达式容易看出，要使平衡点 有实际意义，即位于相平面第一象2 2P限 ( )，必须满足下面两个条件中的一个：10,x12122:,A而由
20、表 2 中 点的 、 可知，仅在条件 下 才是稳定的（而在 下Ppq1A2P2A是鞍点，不稳定） ，图 5 画出了条件 下相轨线的示意图，其中2P1,12xN。直线 和 将相平面（ ）划分为 4 个1220120,x区域： ； ；12:/,/Sdxtdt212:/,/Sxdtt； 。从这 4 个区域中30,40,的正负不难看出其相轨线的趋向如图 5 所示。12/,/xtt1N12/1x2PO2x001S3S24S图 5 在条件 A1下 稳定的相轨线2P分析条件 的实际意义，其关键部分是 1，考虑到 的含义，这表1A2示种群甲要为乙提供足够的食物维持其生长，而 1 则是在 1 条件下为2位于相平面第一象限所必需的，当然这要求 很小（ 是必要条件） ，注2P 11意到 的含义，这实际上是对乙向甲提供食物加以限制，以防止甲的过份增长。1在种群依存模型（9） 、 （12）中如果平衡点 稳定，那么种群乙灭绝，1(,0)PN没有种群的共存，请大家分析导致 稳定的条件及在生态学上的意义。1(,0)评注 模型（9） 、 （12）是种群相互依存的一种类型，即种群甲可独立生存，而种群乙不能，依存模型还有其它类型，如两种群均能独立生存，及均不能独立生存的情况，这些情况的稳态结果如何，大家可以类似讨论。
