第五章稳定性模型

《第五章稳定性模型》由会员分享，可在线阅读，更多相关《第五章稳定性模型（39页珍藏版）》请在装配图网上搜索。

1、第六章第六章 稳定性模型稳定性模型6.1 微分方程稳定性理论简介微分方程稳定性理论简介6.2 军备竞赛军备竞赛6.3 种群的相互竞争种群的相互竞争6.4 种群的相互依存种群的相互依存6.5 种群的弱肉强食种群的弱肉强食稳定性模型稳定性模型 对象仍是动态过程，而建模目的是研究时对象仍是动态过程，而建模目的是研究时间充分长以后过程的变化趋势间充分长以后过程的变化趋势 平衡状平衡状态是否稳定。态是否稳定。不求解微分方程，而是用微分方程稳定性不求解微分方程，而是用微分方程稳定性理论研究平衡状态的稳定性。理论研究平衡状态的稳定性。6.2 军备竞赛军备竞赛 描述双方描述双方(国家或国家集团国家或国家集团)

2、军备竞赛过程军备竞赛过程 解释解释(预测预测)双方军备竞赛的结局双方军备竞赛的结局假设假设 1）由于相互不信任，一方军备越大，另一）由于相互不信任，一方军备越大，另一方军备增加越快；方军备增加越快；2）由于经济实力限制，一方军备越大，对）由于经济实力限制，一方军备越大，对自己军备增长的制约越大；自己军备增长的制约越大；3）由于相互敌视或领土争端，每一方都存）由于相互敌视或领土争端，每一方都存在增加军备的潜力。在增加军备的潜力。进一步进一步假设假设 1）2）的作用为线性；）的作用为线性；3）的作用为常数）的作用为常数目的目的gkyxtx)(建模建模军备竞赛的结局军备竞赛的结局微分方程的平衡点及其

3、稳定性微分方程的平衡点及其稳定性x(t)甲方军备数量，甲方军备数量，y(t)乙方军备数量乙方军备数量hylxty)(,本方经济实力的制约；本方经济实力的制约；k,l 对方对方军备数量的刺激；军备数量的刺激；g,h 本方本方军备竞赛的潜力。军备竞赛的潜力。t 时的时的x(t)，y(t)线性常系数线性常系数微分方程组微分方程组dycxtybyaxtx)()(的平衡点及其稳定性的平衡点及其稳定性平衡点平衡点P0(x0,y0)=(0,0)代数方程代数方程00dycxbyax的根的根若从若从P0某邻域的任一初值出发，都有某邻域的任一初值出发，都有,)(lim0 xtxt称称P0是微分方程的是微分方程的稳

4、定平衡点稳定平衡点,)(lim0ytyt记系数矩阵记系数矩阵dcbaA特征方程特征方程0)det(IAAqdapqpdet)(02特征根特征根2/)4(22,1qpp线性常系数线性常系数微分方程组微分方程组dycxtybyaxtx)()(的平衡点及其稳定性的平衡点及其稳定性特征根特征根2/)4(22,1qpp平衡点平衡点 P0(0,0)微分方程一般解形式微分方程一般解形式ttecec2121平衡点平衡点 P0(0,0)稳定稳定平衡点平衡点 P0(0,0)不稳定不稳定 1,2为负数或有负实部为负数或有负实部p 0 且且 q 0p 0 或或 q kl 下下 x(t),y(t)0,即友好邻国通过裁军

5、可达到永久和平。即友好邻国通过裁军可达到永久和平。hylxtygkyxtx)()(模型模型,本方经济实力的制约；本方经济实力的制约；k,l 对方对方军备数量的刺激；军备数量的刺激；g,h 本方本方军备竞赛的潜力。军备竞赛的潜力。3）若）若 g,h 不为零，即便双方一时和解，使某时不为零，即便双方一时和解，使某时x(t),y(t)很小，但因很小，但因 ，也会重整军备。，也会重整军备。0,0yx4）即使某时一方）即使某时一方(由于战败或协议由于战败或协议)军备大减军备大减,如如 x(t)=0，也会因也会因 使该方重整军备，使该方重整军备，gkyx 即存在互不信任即存在互不信任()或固有争端或固有争

6、端()的单方面的单方面裁军不会持久。裁军不会持久。0k0g模型的定性解释模型的定性解释,本方经济实力的制约；本方经济实力的制约；k,l 对方对方军备数量的刺激；军备数量的刺激；g,h 本方本方军备竞赛的潜力。军备竞赛的潜力。hylxtygkyxtx)()(模型模型6.3 种群的相互竞争种群的相互竞争 一个自然环境中有两个种群生存，它们之间一个自然环境中有两个种群生存，它们之间的关系：相互竞争；相互依存；弱肉强食。的关系：相互竞争；相互依存；弱肉强食。当两个种群为争夺同一食物来源和生存空间当两个种群为争夺同一食物来源和生存空间相互竞争时，常见的结局是，竞争力弱的灭绝，相互竞争时，常见的结局是，竞

7、争力弱的灭绝，竞争力强的达到环境容许的最大容量。竞争力强的达到环境容许的最大容量。建立数学模型描述两个种群相互竞争的过程，建立数学模型描述两个种群相互竞争的过程，分析产生这种结局的条件。分析产生这种结局的条件。221122221)(NxNxxrtx)1()(11111Nxxrtx111111)(Nxxrtx 模型假设模型假设 有甲乙两个种群，它们独自生存有甲乙两个种群，它们独自生存时数量变化均服从时数量变化均服从Logistic规律规律;)1()(22222Nxxrtx 两种群在一起生存时，乙对甲增长的阻滞作两种群在一起生存时，乙对甲增长的阻滞作用与乙的数量成正比用与乙的数量成正比;甲对乙有同

8、样的作用。甲对乙有同样的作用。对于消耗甲的资源而对于消耗甲的资源而言，乙言，乙(相对于相对于N2)是甲是甲(相对于相对于N1)的的 1 倍。倍。11对甲增长的阻滞对甲增长的阻滞作用，乙大于甲作用，乙大于甲乙的竞争力强乙的竞争力强模型模型221Nx模型模型分析分析221111111)(NxNxxrtx221122221)(NxNxxrtx的趋向时)(),(21txtxt(平衡点及其稳定性平衡点及其稳定性)(二阶二阶)非线性非线性(自治自治)方程方程),()(),()(212211xxgtxxxftx的平衡点及其稳定性的平衡点及其稳定性平衡点平衡点P0(x10,x20)代数方程代数方程0),(0)

9、,(2121xxgxxf的根的根若从若从P0某邻域的任一初值出发，都有某邻域的任一初值出发，都有,)(lim011xtxt称称P0是微分方程的是微分方程的稳定平衡点稳定平衡点,)(lim022xtxt模型模型判断判断P0(x10,x20)稳定稳定性的方法性的方法直接法直接法(1)的近似线性方程的近似线性方程)1(),()(),()(212211xxgtxxxftx)2()(,()(,()()(,()(,()(0220201011020120220201011020112121xxxxgxxxxgtxxxxxfxxxxftxxxxx02121PxxxxggffAAqgfpqpPxxdet)(00

10、212平衡点平衡点 P0稳定稳定(对对2,1)p 0 且且 q 0平衡点平衡点 P0不稳定不稳定(对对2,1)p 0 或或 q 0),0(),0,(2211NPNP平衡点：01),(01),(221122221221111121NxNxxrxxgNxNxxrxxf221111111)(NxNxxrtx221122221)(NxNxxrtx仅当仅当 1,2 1时，时，P3才有意义才有意义模型模型)0,0(,1)1(,1)1(4212221113PNNP2211221222211122111121212121NxNxrNxrNxrNxNxrggffAxxxx平衡点稳平衡点稳定性分析定性分析4,3,

11、2,1,det,)(21iAqgfpipipxx2211222212211111211),(1),(NxNxxrxxgNxNxxrxxf平衡点平衡点 Pi 稳定条件：稳定条件：p 0 且且 q 0种群竞争模型的平衡点及稳定性种群竞争模型的平衡点及稳定性不稳定不稳定平平 衡点衡点)0,(11Np)1(221rrpq)1(221 rr),0(22Np211)1(rr)1(121 rr2122211131)1(,1)1(NNp2121211)1)(1(rr)0,0(4p)(21rr 21rr2122111)1()1(rr 21,11,P1,P2 是一个种群存活而另一灭绝的平衡点是一个种群存活而另一灭

12、绝的平衡点P3 是两种群共存的平衡点是两种群共存的平衡点 11,21P1稳定的条件稳定的条件 11?11 21,11P1,P2都不都不(局部局部)稳定稳定1x2x12/N21/N1N2N03P00(3)11,21,21,21加上与加上与(4)相区别的相区别的 11 P2 稳定稳定 P3 稳定稳定P1全局稳定全局稳定结果解释结果解释对于消耗甲的资源而言，对于消耗甲的资源而言，乙乙(相对于相对于N2)是甲是甲(相对相对于于N1)的的 1 倍。倍。11对甲增长的阻滞对甲增长的阻滞作用，乙小于甲作用，乙小于甲乙的竞争力弱乙的竞争力弱 P1稳定的条件：稳定的条件：11 21 甲的竞争力强甲的竞争力强甲达

13、到最大容量，乙灭绝甲达到最大容量，乙灭绝 P2稳定的条件：稳定的条件：11,21 P3稳定的条件：稳定的条件：11,21通常通常 1 1/2，P3稳定条件不满足稳定条件不满足6.4 种群的相互依存种群的相互依存甲乙两甲乙两种群的相互依存有三种形式种群的相互依存有三种形式1)甲可以独自生存，乙不能独自生存；甲甲可以独自生存，乙不能独自生存；甲乙一起生存时相互提供食物、促进增长。乙一起生存时相互提供食物、促进增长。2)甲乙均可以独自生存；甲乙一起生存甲乙均可以独自生存；甲乙一起生存 时相互提供食物、促进增长。时相互提供食物、促进增长。3)甲乙均不能独自生存；甲乙一起生存甲乙均不能独自生存；甲乙一起

14、生存时相互提供食物、促进增长。时相互提供食物、促进增长。1111111)(Nxxrtx 模型模型假设假设 甲可以独自生存，数量变化服从甲可以独自生存，数量变化服从Logistic规律规律;甲乙一起生存时乙为甲提供食物、促进增长。甲乙一起生存时乙为甲提供食物、促进增长。乙不能独自生存；甲乙一起生存时甲为乙乙不能独自生存；甲乙一起生存时甲为乙提供食物、促进增长；乙的增长又受到本身提供食物、促进增长；乙的增长又受到本身的阻滞作用的阻滞作用(服从服从Logistic规律规律)。模型模型乙为甲提供食物乙为甲提供食物是甲消耗的是甲消耗的 1 倍倍221Nx甲为乙提供食物甲为乙提供食物是乙消耗的是乙消耗的

15、2 倍倍1)(222xrtx 1122221)(Nxxrtx221122221)(NxNxxrtx种群依存模型的平衡点及稳定性种群依存模型的平衡点及稳定性P2是甲乙相互依存而共生的平衡点是甲乙相互依存而共生的平衡点稳定条件稳定条件不稳定不稳定1,12121,1,12121平衡点平衡点pq)0,(11NP)1(221rr)1(221rr)0,0(3P21rr 21rr2122211121)1(,1)1(NNP2121211)1)(1(rr2122111)1()1(rr平衡点平衡点P2稳定稳定性的相轨线性的相轨线2122211121)1(,1)1(NNP),(1)(2111221111111xxx

16、rNxNxxrtx),(1)(212222112222xxxrNxNxxrtx.0,0:;0,0:;0,0:;0,0:214213212211xxSxxSxxSxxS1x2x021/N1N1S2S3S2P004S 11,1 21 P2稳定稳定 1 21 前提下前提下P2存在的必要条件存在的必要条件结果结果解释解释2122211121)1(,1)1(NNP 21 甲必须为乙提供足够的食物甲必须为乙提供足够的食物甲为乙提供的食物是乙消耗的甲为乙提供的食物是乙消耗的 2 倍倍 11,1 21条件下使条件下使 1 21 成立成立 P2稳定条件：稳定条件：11,1 2 0P:临界状态临界状态 q 0P

17、不稳定不稳定 0/0abrbadAP)0,0(Ptx(t)y(t)020.00004.00000.100021.24063.96510.200022.56493.94050.300023.97633.92695.10009.616216.72355.2000 9.017316.20649.500018.47504.04479.600019.61363.99689.700020.83113.9587用数学软件用数学软件MATLAB求求微分方程数值解微分方程数值解xy 平面上的相轨线平面上的相轨线计算结果（数值，图形）计算结果（数值，图形）x(t),y(t)是周期函数，相图是周期函数，相图(x,y

18、)是封闭曲线是封闭曲线xayrtx)()(ybxdty)()(观察，猜测观察，猜测x(t),y(t)的周期约为的周期约为9.6xmax 65.5,xmin 6,ymax 20.5,ymin 3.9用数值积分可算出用数值积分可算出 x(t),y(t)一周期的平均值：一周期的平均值：x(t)的平均值约为的平均值约为25,y(t)的平均值约为的平均值约为10。食饵食饵-捕食者模型捕食者模型(Volterra)()(bxdyayrxdydxdyyayrdxxbxd 消去消去dt1lnlncayyrbxxdybxdtyxayrtx)()()()(ceyexayrbxd)(用相轨线分析用相轨线分析 点稳定

19、性点稳定性)/,/(arbdPc 由初始条件确定由初始条件确定取指数取指数x0fmf(x)x0g(y)gmy0y0arygygggm/,)(,0)()0(00,0)()0(ffcygxf)()(ceyexayrbxd)(在相平面上讨论相轨线的图形在相平面上讨论相轨线的图形用相轨线分析用相轨线分析 点稳定性点稳定性)/,/(arbdP相轨线相轨线)(xf)(ygbdxfxfm/,)(00mmgfc 时无相轨线时无相轨线以下设以下设mmgfc y2y1xQ3Q4qy1y2x1x2pyy0 xx0P0 x1x2Q1Q2Q1(x1,y0),Q2(x2,y0)Q3(x,y1),Q4(x,y2)mmgfc

20、 00,yyxx相轨线相轨线退化为退化为P点点mmgfc 0yy令mfpxf)(mpgc 设 存在x1x0 x2,使f(x1)=f(x2)=pmgyg)(,21xxx考察pxf)(存在y1y0y2,使g(y1)=g(y2)=qmpgygxf)()(mgqyg)(内任意点是,21xxx相轨线是封闭曲线族相轨线是封闭曲线族xQ3Q4f(x)xx0fm0g(y)gmy0y0cygxf)()(相轨线相轨线P中心中心相轨线相轨线是封闭曲线是封闭曲线x(t),y(t)是周期函数是周期函数(周期记周期记 T)求求x(t),y(t)在一周期的平均值在一周期的平均值yx,ybxdty)()()(1)(dyybt

21、xTdttxTx0)(1xayrtx)()(arybdxyxP/,/:),(0000轨线轨线中心中心00,yyxxbdx/ary/用相轨线分析用相轨线分析 点稳定性点稳定性)/,/(arbdPTdtdyybT0)(11)0(ln)(ln(1bdTbyTyT02040608010012005101520253000yx，00yx，00yx00yx0PT2T3T4T1P024681012020406080100120 x(t)y(t)T1 T2 T3 T4x(t)的的“相位相位”领先领先 y(t)xayrtx)()(ybxdty)()()()(:1tytxT)()(:2tytxT)()(:3tyt

22、xT)()(:4tytxT模型解释模型解释)/,/(arbdP),(000yxP初值初值相轨线的方向相轨线的方向模型解释模型解释r 食饵增长率食饵增长率d 捕食者死亡率捕食者死亡率b 食饵供养捕食者能力食饵供养捕食者能力ary 捕食者捕食者 数量数量bdx 食饵食饵数量数量020406080100120051015202530P)/,/(arbdPr/ad/ba 捕食者掠取食饵能力捕食者掠取食饵能力捕食者数量与捕食者数量与r成正比成正比,与与a成反比成反比食饵食饵数量与数量与d成正比成正比,与与b成反比成反比模型模型解释解释一次大战期间地中海渔业的捕捞量下降，一次大战期间地中海渔业的捕捞量下降

23、，但是其中但是其中鲨鱼的比例却在增加，为什么？鲨鱼的比例却在增加，为什么？rr-1,dd+1捕捞捕捞战时战时捕捞捕捞rr-2,dd+2,2 1yyxx11,),(111yxP),(222yxP),(yxPxy食饵食饵(鱼鱼)减少，减少，捕食者捕食者(鲨鱼鲨鱼)增加增加自然环境自然环境arybdx/,/),(yxP1212,yyxx1PP 21PP 还表明：对还表明：对害虫害虫(食饵食饵)益虫益虫(捕食者捕食者)系统，系统，使用灭两种使用灭两种虫的虫的杀虫剂杀虫剂,会使害虫增加，益虫减少。会使害虫增加，益虫减少。1PP 食饵食饵-捕食者模型捕食者模型(Volterra)的缺点与改进的缺点与改进2

24、211111)(Nxxrtx1122221)(Nxxrtx221111111)(NxNxxrtx221122221)(NxNxxrtxxayrtx)()(ybxdty)()(Volterra模型模型改写改写多数多数食饵食饵捕食者系统观察不到周期震荡捕食者系统观察不到周期震荡,而是趋向某个平衡状态而是趋向某个平衡状态,即存在稳定平衡点即存在稳定平衡点加加Logistic项项有有稳定平衡点稳定平衡点 相轨线是封闭曲线，结构不稳定相轨线是封闭曲线，结构不稳定一旦离开某一旦离开某一条闭轨线，就进入另一条闭轨线，不恢复原状。一条闭轨线，就进入另一条闭轨线，不恢复原状。自然界存在的周期性平衡生态系统是结构

25、稳定的，自然界存在的周期性平衡生态系统是结构稳定的，即偏离周期轨道后，内部制约使系统恢复原状。即偏离周期轨道后，内部制约使系统恢复原状。食饵食饵-捕食者模型捕食者模型(Volterra)的缺点与改进的缺点与改进1211111111)(wxxNxxrtx11222211)(wxxxrtxr1=1,N1=20,1=0.1,w=0.2,r2=0.5,2=0.18相轨线趋向极限环相轨线趋向极限环051015200102030结构稳定结构稳定 两种群模型的几种形式两种群模型的几种形式 相互竞争相互竞争221111111)(NxNxxrtx221122221)(NxNxxrtx相互依存相互依存2211111111)(NxNxxrtx221122221)(NxNxxrtx弱肉强食弱肉强食221111111)(NxNxxrtx221122221)(NxNxxrtx