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

PAGE/NUMPAGES第四节微分方程稳定性理论简介这里简单介绍下面将要用到的有关内容：一阶方程的平衡点及稳定性设有微分方程SKIPIF1<0（1）右端不显含自变量t，代数方程SKIPIF1<0（2）的实根SKIPIF1<0称为方程（1）的平衡点（或奇点），它也是方程（1）的解（奇解）如果从所有可能的初始条件出发，方程（1）的解SKIPIF1<0都满足SKIPIF1<0（3）则称平衡点SKIPIF1<0是稳定的（稳定性理论中称渐近稳定）；否则，称SKIPIF1<0是不稳定的（不渐近稳定）。判断平衡点SKIPIF1<0是否稳定通常有两种方法，利用定义即（3）式称间接法，不求方程（1）的解SKIPIF1<0，因而不利用（3）式的方法称直接法，下面介绍直接法。将SKIPIF1<0在SKIPIF1<0做泰勒展开，只取一次项，则方程（1）近似为：SKIPIF1<0（4）（4）称为（1）的近似线性方程。SKIPIF1<0也是（4）的平衡点。关于平衡点SKIPIF1<0的稳定性有如下的结论：若SKIPIF1<0，则SKIPIF1<0是方程（1）、（4）的稳定的平衡点。若SKIPIF1<0，则SKIPIF1<0不是方程（1）、（4）的稳定的平衡点SKIPIF1<0对于方程（4）的稳定性很容易由定义（3）证明，因为（4）的一般解是SKIPIF1<0（5）其中C是由初始条件决定的常数。二阶（平面）方程的平衡点和稳定性方程的一般形式可用两个一阶方程表示为SKIPIF1<0（6）右端不显含t，代数方程组SKIPIF1<0（7）的实根SKIPIF1<0称为方程（6）的平衡点。记为SKIPIF1<0如果从所有可能的初始条件出发，方程（6）的解SKIPIF1<0都满足SKIPIF1<0SKIPIF1<0（8）则称平衡点SKIPIF1<0是稳定的（渐近稳定）；否则，称P0是不稳定的（不渐近稳定）。为了用直接法讨论方法方程（6）的平衡点的稳定性，先看线性常系数方程SKIPIF1<0（9）系数矩阵记作SKIPIF1<0并假定A的行列式SKIPIF1<0于是原点SKIPIF1<0是方程（9）的唯一平衡点，它的稳定性由的特征方程SKIPIF1<0的根SKIPIF1<0（特征根）决定，上方程可以写成更加明确的形式：SKIPIF1<0（10）将特征根记作SKIPIF1<0，则SKIPIF1<0（11）方程（9）的解一般有形式SKIPIF1<0（SKIPIF1<0）或SKIPIF1<0（SKIPIF1<0）SKIPIF1<0为任意实数。由定义（8），当SKIPIF1<0全为负数或有负的实部时SKIPIF1<0是稳定的平衡点，反之，当SKIPIF1<0有一个为正数或有正的实部时SKIPIF1<0是不稳定的平衡点微分方程稳定性理论将平衡点分为结点、焦点、鞍点、中心等类型，完全由特征根SKIPIF1<0或相应的SKIPIF1<0取值决定，下表简明地给出了这些结果，表中最后一列指按照定义（8）式得下马看花关于稳定性的结论。表1由特征方程决定的平衡点的类型和稳定性SKIPIF1<0SKIPIF1<0SKIPIF1<0平衡点类型稳定性SKIPIF1<0SKIPIF1<0稳定结点稳定SKIPIF1<0SKIPIF1<0不稳定结点不稳定SKIPIF1<0SKIPIF1<0鞍点不稳定SKIPIF1<0SKIPIF1<0稳定退化结点稳定SKIPIF1<0SKIPIF1<0不稳定退化结点不稳定SKIPIF1<0SKIPIF1<0稳定焦点稳定SKIPIF1<0SKIPIF1<0不稳定焦点不稳定SKIPIF1<0SKIPIF1<0中心不稳定由上表可以看出，根据特征方程的系数SKIPIF1<0的正负很容易判断平衡点的稳定性，准则如下：若SKIPIF1<0（12）则平衡点稳定，若SKIPIF1<0SKIPIF1<0（13）则平衡点不稳定以上是对线性方程（9）的平衡点SKIPIF1<0稳定性的结论，对于一般的非线性方程（6），可以用近似线性方法判断其平衡点SKIPIF1<0的稳定性，在SKIPIF1<0点将SKIPIF1<0和SKIPIF1<0作泰勒展开，只取一次项，得（6）的近似线性方程SKIPIF1<0(14)系数矩阵记作SKIPIF1<0特征方程系数为SKIPIF1<0,SKIPIF1<0显然，SKIPIF1<0点对于方程（14）的稳定性由表1或准则（12）、（13）决定，而且已经证明了如下结论:若方程（14）的特征根不为零或实部不为零，则SKIPIF1<0点对于方程（6）的稳定性与对于近似方程（14）的稳定性相同。这样，SKIPIF1<0点对于方程（6）的稳定性也由准则（12）、（13）决定。第五节种群的相互竞争与相互依存当某个自然环境中只有一种生物的群体（生态学上称为种群）生存时，人们常用Logistic模型来描述这个群数量的演变过程，即SKIPIF1<0(1)x（t）是种群在时刻t的数量，SKIPIF1<0是固有增长率，N是环境资源容许的种群最大数量，在前面我们曾应用过这种模型，由方程（1）可以直接得到，SKIPIF1<0=N是稳定平衡点，即t→∞时x(t)→N，从模型本身的意义看这是明显的结果。如果一个自然环境中有两个或两个以上种群生存，那么它们之间就要存在着或是相互竞争，或是相互依存，或是弱肉强食（食饵与捕食者）的关系。这里将从稳定状态的角度分别讨论这些关系。一、种群的相互竞争当两个种群为了争夺有限的食物来源和生活空间而进行生存竞争时，最常见的结局是竞争力较弱的种群灭绝，竞争力较强的种群达到环境容许的最大数量。人们今天可以看到自然界长期演变成的这样的结局，例如一个小岛上虽然有四种燕子栖息，但是它们的食物来源各不相同，一种只在陆地上觅食，另两种分别在浅水的海滩上和离岸稍远的海中捕鱼，第四种则飞越宽阔的海面到远方攫取海味，每一种燕子在它各自生存环境中的竞争力明显地强于其它几种，这里我们建立一个模型解释类似的现象，并分析产生这种结局的条件。模型建立有甲乙两个种群，当它们独自在一个自然环境中生存时，数量的演变均遵从Logistic规律，记SKIPIF1<0是两个种群的数量，SKIPIF1<0是它们的固有增长率，N1、N2是它们的最大容量，于是对于种群甲有SKIPIF1<0其中因子SKIPIF1<0反映由于甲方有限资源的消耗导致的对它本身增长的阻滞作用，SKIPIF1<0可解释为相对于N1而言单位数量的甲消耗的供养甲的食物量（设食物总量为1）。当两个种群在同一自然环境中生存时，考察由于乙消耗同一种有限资源对甲的增长产生的影响，可以合理地在因子SKIPIF1<0中再减去一项，该项与种群乙的数量SKIPIF1<0（相对于N2而言）成正比，得到种群甲方增长的方程SKIPIF1<0(2)这里SKIPIF1<0的意义是，单位数量乙（相对N2而言）消耗的供养甲的食物量为单位数量甲（相对N1）消耗的供养甲的食物量的SKIPIF1<0倍。类似地，甲的存在也影响了乙的增长，种群乙的方程应该是SKIPIF1<0(3)对SKIPIF1<0可作相应的解释。在两种群的相互竞争中SKIPIF1<0、SKIPIF1<0是两个关键指标，从上面对它们的解释可知，SKIPIF1<0＞1表示在消耗供养甲的资源中，乙的消耗多于甲，因而对甲增长的阻滞作用乙大于甲，即乙的竞争力强于甲，对SKIPIF1<0＞1可作相应的理解。一般地说，SKIPIF1<0与SKIPIF1<0之间没有确定的关系，但是可以把下面这种特殊情况作为较常见的一类实际情况的典型代表，即两个种群在消耗资源中对甲增长的阻作用对乙增长的阻滞作用相同，具体地说就是，因为单位数量的甲和乙消耗的供养甲方食物量之比是1：SKIPIF1<0，消耗的供养甲方食物量之比是SKIPIF1<0：1，所谓阻滞作用相同即1：SKIPIF1<0=SKIPIF1<0：1，所以这种特殊情形可以定量地表示为SKIPIF1<0SKIPIF1<0=1（4）即SKIPIF1<0、SKIPIF1<0互为倒数，可以简单地理解为，如果一个乙消耗的食物是一个甲的SKIPIF1<0=SKIPIF1<0倍，则一个甲消耗的食物是一个乙的SKIPIF1<0=1/SKIPIF1<0。下面我们仍然讨论SKIPIF1<0、SKIPIF1<0相互独立的一般情况，而将条件（4）下对问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 的分析留给大家讨论。稳定性分析为了研究两个种群相互竞争的结局，即t→∞时SKIPIF1<0的趋向，不必要解方程（2）、（3）,只需对它的平衡点进行稳定性分析。首先根据微分方程（2）、（3）解代数方程组SKIPIF1<0（5）得到4个平衡点：SKIPIF1<0因为仅当平衡点们于平面坐标系的第一象限时（SKIPIF1<0）才有实际意义，所以对SKIPIF1<0而言要求SKIPIF1<0、SKIPIF1<0同时小于1，或同时大于1。按照判断平衡点性的方法（见前面）计算SKIPIF1<0SKIPIF1<0SKIPIF1<0将4个平衡点p、q的结果及稳定条件列入下表*）表1种群竞争模型的平衡点及稳定性平衡点pq稳定条件SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0不稳定注：表中最后一列“稳定条件”除了要求p>0,q>0以外，还有其他原因，见下面的具体分析。为了便于对平衡点P1、P2、P3的稳定条件进行分析，在相平面上讨论它们。在代数方程组（5）中记SKIPIF1<0SKIPIF1<0对于SKIPIF1<0、SKIPIF1<0的不同取值范围，直线SKIPIF1<0=0和SKIPIF1<0=0在相平面上的相对位置不同，下面给出它们的4种情况;并对这4种情况进行分析1、SKIPIF1<0。由表1知对于SKIPIF1<0有SKIPIF1<0＞0，SKIPIF1<0＜0，SKIPIF1<0稳定；SKIPIF1<0的稳定性还可以从t→∞时相轨线的趋向来分析，图1中SKIPIF1<0=0和SKIPIF1<0=0两条直线将相平面（SKIPIF1<0）划分为3个区域：SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0OSKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0（1）SKIPIF1<0SKIPIF1<0稳定图1SKIPIF1<0SKIPIF1<0稳定SKIPIF1<0(6)SKIPIF1<0(7)SKIPIF1<0(8)可以证明，不论轨线从哪个区域出发，t→∞时都将趋向P1（N1，0）。若轨线从S1出发，由（6）可知随着t的增加轨线向右上方运动，必然进入S2；若轨线从S2出发，由（7）可知轨线向右下方运动，那么它或者趋向SKIPIF1<0点，或者进入S3，但是进入S3是不可能的，因为，如果设轨线在某时刻t1经直线SKIPIF1<0=0进入S3，则SKIPIF1<0SKIPIF1<0(t1)/SKIPIF1<0=0,由方程（2）不难算出SKIPIF1<0由（7）、（8）知SKIPIF1<0＜0,故SKIPIF1<0，表明SKIPIF1<0(t)在t1达到极小值，而这是不可能的，因为在S2中SKIPIF1<0＞0,即SKIPIF1<0(t)一直是增加的；若轨线从S3出发，由（8）可知轨线向左下方运动，那么它或者趋向SKIPIF1<0点，或者进入S2，而进入S2后，根据上面的分析最终也将趋向SKIPIF1<0。综上分析可以画出轨线示意图（图1），因为直线SKIPIF1<0=0上dSKIPIF1<0=0，所以在SKIPIF1<0=0上轨线方向垂直于SKIPIF1<0轴；在SKIPIF1<0=0上dSKIPIF1<0=0,轨线方向平行于SKIPIF1<0轴。2、SKIPIF1<0，类似的分析可知SKIPIF1<0稳定。SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0OSKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0（2）SKIPIF1<0SKIPIF1<0稳定图2SKIPIF1<0SKIPIF1<0稳定3、SKIPIF1<0，由表1知对于SKIPIF1<0点SKIPIF1<0＞0，SKIPIF1<0＞0，故SKIPIF1<0稳定，对轨线趋势的分析见图3。SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0OSKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0（3）SKIPIF1<0SKIPIF1<0稳定图3SKIPIF1<0SKIPIF1<0稳定4、SKIPIF1<0，由表1知对于SKIPIF1<0点SKIPIF1<0＜0，故SKIPIF1<0不稳定（鞍点），轨线或者趋向SKIPIF1<0，或者趋向SKIPIF1<0，由轨线的初始位置决定，示意图见图4，在这种情况下SKIPIF1<0和SKIPIF1<0都不能说是稳定的，正因为这样，所以SKIPIF1<0稳定（与初始条件无关）的条件需要加上SKIPIF1<0，SKIPIF1<0稳定的条件加上SKIPIF1<0。SKIPIF1<0SKIPIF1<0SKIPIF1<0OSKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0（4）SKIPIF1<0SKIPIF1<0不稳定SKIPIF1<0图4SKIPIF1<0SKIPIF1<0不稳定结果解释根据建模过程中SKIPIF1<0的含义，说明SKIPIF1<0、SKIPIF1<0、SKIPIF1<0点稳定在生态上的意义。1、SKIPIF1<0，SKIPIF1<0意味着在对供养甲的资源的竞争中乙弱于甲，SKIPIF1<0意味着在对供养乙的资源的竞争中甲强于乙，于是种群乙终灭绝，种群甲趋向最大容量，即SKIPIF1<0趋向平衡点SKIPIF1<02、SKIPIF1<0，情况与1正好的相反。3、SKIPIF1<0，因为在竞争甲的资源中乙较弱，而在竞争乙的资源中甲较弱，于是可以达到一个双方共存的稳定的平衡状态SKIPIF1<0，这是种群竞争中很少出现的情况。4、SKIPIF1<0，请大家作出解释。生态学中有一个竞争排斥原理；若两个种群的单个成员消耗的资源差不多相同，而环境能承受的种群甲的最大容量比种群乙大，那么种群乙终将灭亡，用本节的模型很容解释这个原理。将方程（2）、（3）改写为SKIPIF1<0SKIPIF1<0原理的两个条件相当于SKIPIF1<0从这3个式子显然可得SKIPIF1<0，这正是SKIPIF1<0稳定，即种群乙灭绝的条件。二、种群的相互依存自然界中处于同一环境下两个种群相互依存而共生的现象是很普遍的，植物可以独立生存。昆虫的的授粉作用又可以提高植物的增长率，而以花粉为食物的昆虫却不能离开植物单独存活，人类与人工饲养的牲畜之间也有类似的关系，这种共生现象可以描述如下。设种群甲可以独立存在，按Logistic规律增长，种群乙为甲提供食物，有助于甲的增长，类似于前面的方程（2），种群甲的数量演变规律可以写作（r1、N1、N2的意义同前）SKIPIF1<0(9)SKIPIF1<0前面的SKIPIF1<0号这里变成+号，表示乙不是消耗甲的资源而是为甲提供食物，SKIPIF1<0的含义是:单位数量乙（相对于N2）提供的供养甲的食物量为单位数量甲（相对于N1）消耗的供养甲食物量的SKIPIF1<0倍。种群乙没有甲的存在会灭亡，设其死亡率为r2，则乙单独存在时有SKIPIF1<0(10)甲为乙提供食物，于是（2）式右端应加上甲对乙增长的促进作用，有SKIPIF1<0(11)显然仅当SKIPIF1<0时种群乙的数量才会增长，与此相同乙的增长又会受到自身的阻滞作用，所以93）式右端还要添加Logistic项，方程变为SKIPIF1<0(12)方程（9）、（12）构成相互依存现象的数学模型，下面利用平衡点的稳定性分析，讨论时间足够长以后两个种群的变化趋向。类似于前面的作法将方程（9）、（12）的平衡点及其稳定性分析的结果列入表2表2种群依存模型的平衡点及稳定性平衡点SKIPIF1<0SKIPIF1<0稳定条件SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0不稳定显然，SKIPIF1<0点稳定才表明两个种群在同一环境里相互依存而共生，我们着重分析SKIPIF1<0稳定的条件。。由SKIPIF1<0的表达式容易看出，要使平衡点SKIPIF1<0有实际意义，即位于相平面第一象限(SKIPIF1<0)，必须满足下面两个条件中的一个：SKIPIF1<0而由表2中SKIPIF1<0点的SKIPIF1<0、SKIPIF1<0可知，仅在条件SKIPIF1<0下SKIPIF1<0才是稳定的（而在SKIPIF1<0下SKIPIF1<0是鞍点，不稳定），图5画出了条件SKIPIF1<0下相轨线的示意图，其中SKIPIF1<0,SKIPIF1<0。直线SKIPIF1<0和SKIPIF1<0将相平面（SKIPIF1<0）划分为4个区域：SKIPIF1<0；SKIPIF1<0；SKIPIF1<0；SKIPIF1<0。从这4个区域中SKIPIF1<0SKIPIF1<0的正负不难看出其相轨线的趋向如图5所示。SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0OSKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0图5在条件A1下SKIPIF1<0稳定的相轨线分析条件SKIPIF1<0的实际意义，其关键部分是SKIPIF1<0＞1，考虑到SKIPIF1<0的含义，这表示种群甲要为乙提供足够的食物维持其生长，而SKIPIF1<0<1则是在SKIPIF1<0＞1条件下为SKIPIF1<0位于相平面第一象限所必需的，当然这要求SKIPIF1<0很小（SKIPIF1<0是必要条件），注意到SKIPIF1<0的含义，这实际上是对乙向甲提供食物加以限制，以防止甲的过份增长。在种群依存模型（9）、（12）中如果平衡点SKIPIF1<0稳定，那么种群乙灭绝，没有种群的共存，请大家分析导致SKIPIF1<0稳定的条件及在生态学上的意义。评注模型（9）、（12）是种群相互依存的一种类型，即种群甲可独立生存，而种群乙不能，依存模型还有其它类型，如两种群均能独立生存，及均不能独立生存的情况，这些情况的稳态结果如何，大家可以类似讨论。友情提示： 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 范本是经验性极强的领域，本范文无法思考和涵盖全面，供参考!最好找专业人士起草或审核后使用。