《稳定性模型1.ppt》由会员分享,可在线阅读,更多相关《稳定性模型1.ppt(44页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、稳定性模型稳定性模型1 18.1 捕鱼业的持续收获捕鱼业的持续收获 再生资源(渔业、林业等)与再生资源(渔业、林业等)与非再生资源(矿业等)非再生资源(矿业等)再生资源应适度开发再生资源应适度开发在持续稳在持续稳产前提下实现最大产量或最佳效益。产前提下实现最大产量或最佳效益。问题问题及及 分析分析 在在捕捞量稳定捕捞量稳定的条件下,如何控的条件下,如何控制捕捞使产量最大或效益最佳。制捕捞使产量最大或效益最佳。如果使捕捞量等于自然增长量,如果使捕捞量等于自然增长量,渔渔场鱼量将保持不变场鱼量将保持不变,则捕捞量稳定。,则捕捞量稳定。背景背景产量模型产量模型假设假设 无捕捞时鱼的自然增长服从无捕捞
2、时鱼的自然增长服从 Logistic规律规律 单位时间捕捞量与渔场鱼量成正比单位时间捕捞量与渔场鱼量成正比建模建模 捕捞情况下捕捞情况下渔场鱼量满足渔场鱼量满足 不需要求解不需要求解x(t),只需知道只需知道x(t)稳定稳定的条件的条件r固有增长率固有增长率,N最大鱼量最大鱼量h(x)=Ex,E捕捞强度捕捞强度x(t)渔场鱼量渔场鱼量一阶微分方程的平衡点及其稳定性一阶微分方程的平衡点及其稳定性一阶非线性(自治)方程一阶非线性(自治)方程F(x)=0的根的根x0 微分方程的微分方程的平衡点平衡点设设x(t)是方程的解,若从是方程的解,若从x0 某邻域的任一初值出发,某邻域的任一初值出发,都有都有
3、称称x0是方程是方程(1)的的稳定平衡点稳定平衡点不求不求x(t),判断判断x0稳定性的方法稳定性的方法直接法直接法(1)的近似线性方程的近似线性方程产量模型产量模型平衡点平衡点稳定性判断稳定性判断x0 稳定稳定,可得到稳定产量可得到稳定产量x1 稳定稳定,渔场干枯渔场干枯E捕捞强度捕捞强度r固有增长率固有增长率产量模型产量模型在捕捞量稳定的条件下,在捕捞量稳定的条件下,控制捕捞强度使产量最大控制捕捞强度使产量最大图解法图解法P的横坐标的横坐标 x0平衡点平衡点y=rxhPx0y0y=h(x)=ExxNy=f(x)P的纵坐标的纵坐标 h产量产量产量最大产量最大f 与与h交点交点Phmx0*=N
4、/2P*y=E*x控制渔场鱼量为最大鱼量的一半控制渔场鱼量为最大鱼量的一半效益模型效益模型假设假设 鱼销售价格鱼销售价格p 单位捕捞强度费用单位捕捞强度费用c 单位时间利润单位时间利润在捕捞量稳定的条件下,控制捕捞在捕捞量稳定的条件下,控制捕捞强度使效益最大强度使效益最大.稳定平衡点稳定平衡点求求E使使R(E)最大最大渔场渔场鱼量鱼量收入收入 T=ph(x)=pEx支出支出 S=cEEsS(E)T(E)0rE捕捞捕捞过度过度 封闭式捕捞封闭式捕捞追求利润追求利润R(E)最大最大 开放式捕捞开放式捕捞只求利润只求利润R(E)0R(E)=0时的捕捞强度时的捕捞强度(临界强度临界强度)Es=2ER临
5、界强度下的渔场鱼量临界强度下的渔场鱼量捕捞过度捕捞过度ERE*令令=08.3 种群的相互竞争种群的相互竞争 一个自然环境中有两个种群生存,它们之间的一个自然环境中有两个种群生存,它们之间的关系:相互竞争;相互依存;弱肉强食。关系:相互竞争;相互依存;弱肉强食。当两个种群为争夺同一食物来源和生存空间相当两个种群为争夺同一食物来源和生存空间相互竞争时,常见的结局是,竞争力弱的灭绝,互竞争时,常见的结局是,竞争力弱的灭绝,竞争力强的达到环境容许的最大容量。竞争力强的达到环境容许的最大容量。建立数学模型描述两个种群相互竞争的过程,建立数学模型描述两个种群相互竞争的过程,分析产生这种结局的条件。分析产生
6、这种结局的条件。模型假设模型假设 有甲乙两个种群,它们独自生存有甲乙两个种群,它们独自生存时数量变化均服从时数量变化均服从Logistic规律规律;两种群在一起生存时,乙对甲增长的阻滞作两种群在一起生存时,乙对甲增长的阻滞作用与乙的数量成正比用与乙的数量成正比;甲对乙有同样的作用。甲对乙有同样的作用。对于消耗甲的资源而言,对于消耗甲的资源而言,乙乙(相对于相对于N2)是甲是甲(相相对于对于N1)的的 1 倍。倍。对甲增长的阻滞对甲增长的阻滞作用,乙大于甲作用,乙大于甲乙的竞争力强乙的竞争力强模型模型模型模型分析分析(平衡点及其稳定性平衡点及其稳定性)(二阶二阶)非线性非线性(自治自治)方程方程
7、的平衡点及其稳定性的平衡点及其稳定性平衡点平衡点P0(x10,x20)代数方程代数方程的根的根若从若从P0某邻域的任一初值出发,都有某邻域的任一初值出发,都有称称P0是微分方程的是微分方程的稳定平衡点稳定平衡点模型模型判断判断P0(x10,x20)稳定稳定性的方法性的方法直接法直接法(1)的近似线性方程的近似线性方程平衡点平衡点 P0稳定稳定(对对2,1)p 0 且且 q 0平衡点平衡点 P0不稳定不稳定(对对2,1)p 0 或或 q 0仅当仅当 1,2 1时,时,P3才有意义才有意义模型模型平衡点稳平衡点稳定性分析定性分析平衡点平衡点 Pi 稳定条件:稳定条件:p 0 且且 q 0种群竞争模
8、型的平衡点及稳定性种群竞争模型的平衡点及稳定性不稳定不稳定平平 衡点衡点 21,11,P1,P2 是一个种群存活而另一灭绝的平衡点是一个种群存活而另一灭绝的平衡点P3 是两种群共存的平衡点是两种群共存的平衡点 11,21P1稳定的条件稳定的条件 11?11 21,11P1,P2都不都不(局部局部)稳定稳定0(3)11,21,21,21加上与加上与(4)相区别的相区别的 11 P2 稳定稳定 P3 稳定稳定P1全局稳定全局稳定结果解释结果解释对于消耗甲的资源而言,对于消耗甲的资源而言,乙乙(相对于相对于N2)是甲是甲(相对相对于于N1)的的 1 倍。倍。对甲增长的阻滞对甲增长的阻滞作用,乙小于甲
9、作用,乙小于甲乙的竞争力弱乙的竞争力弱 P1稳定的条件:稳定的条件:11 21 甲的竞争力强甲的竞争力强甲达到最大容量,乙灭绝甲达到最大容量,乙灭绝 P2稳定的条件:稳定的条件:11,21 P3稳定的条件:稳定的条件:11,21通常通常 1 1/2,P3稳定条件不满足稳定条件不满足8.4 种群的相互依存种群的相互依存甲乙两甲乙两种群的相互依存有三种形式种群的相互依存有三种形式1)甲可以独自生存,乙不能独自生存;甲甲可以独自生存,乙不能独自生存;甲乙一起生存时相互提供食物、促进增长。乙一起生存时相互提供食物、促进增长。2)甲乙均可以独自生存;甲乙一起生存甲乙均可以独自生存;甲乙一起生存 时时相互
10、提供食物、促进增长。相互提供食物、促进增长。3)甲乙均不能独自生存;甲乙一起生存时甲乙均不能独自生存;甲乙一起生存时相互提供食物、促进增长。相互提供食物、促进增长。模型模型假设假设 甲可以独自生存,数量变化服从甲可以独自生存,数量变化服从Logistic规律规律;甲乙一起生存时乙为甲提供食物、促进增长。甲乙一起生存时乙为甲提供食物、促进增长。乙不能独自生存;甲乙一起生存时甲为乙提乙不能独自生存;甲乙一起生存时甲为乙提供食物、促进增长;乙的增长又受到本身的供食物、促进增长;乙的增长又受到本身的阻滞作用阻滞作用(服从服从Logistic规律规律)。模型模型乙为甲提供食物乙为甲提供食物是甲消耗的是甲
11、消耗的 1 倍倍甲为乙提供食物甲为乙提供食物是乙消耗的是乙消耗的 2 倍倍种群依存模型的平衡点及稳定性种群依存模型的平衡点及稳定性P2是甲乙相互依存而共生的平衡点是甲乙相互依存而共生的平衡点稳定条件稳定条件不稳定不稳定平衡点平衡点平衡点平衡点P2稳定稳定性的相轨线性的相轨线0 11,1 21 P2稳定稳定 1 21 前提下前提下P2存在的必要条件存在的必要条件结果结果解释解释 21 甲必须为乙提供足够的食物甲必须为乙提供足够的食物甲为乙提供的食物是乙消耗的甲为乙提供的食物是乙消耗的 2 倍倍 11,1 21条件下使条件下使 1 21 成立成立 P2稳定条件:稳定条件:11,1 2 0P:临界状
12、态临界状态 q 0P 不稳定不稳定 tx(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)是封闭曲是
13、封闭曲线线观察,猜测观察,猜测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)消去消去dt用相轨线分析用相轨线分析 点稳定性点稳定性c 由初始条件确定由初始条件确定取指数取指数x0fmf(x)x0g(y)gmy0y0在相平面上讨论相轨线的图形在相平面上讨论相轨线的图形用相轨线分析用相轨线分析 点稳定性点稳定性相轨线相轨线时
14、无相轨线时无相轨线以下设以下设y2y1xQ3Q4qy1y2x1x2pyy0 xx0P0 x1x2Q1Q2Q1(x1,y0),Q2(x2,y0)Q3(x,y1),Q4(x,y2)相轨线相轨线退化为退化为P点点 存在x1x0 x2,使f(x1)=f(x2)=p存在y1y0y2,使g(y1)=g(y2)=q相轨线是封闭曲线族相轨线是封闭曲线族xQ3Q4f(x)xx0fm0g(y)gmy0y0相轨线相轨线P中心中心相轨线相轨线是封闭曲线是封闭曲线x(t),y(t)是周期函数是周期函数(周期记周期记 T)求求x(t),y(t)在一周期的平均值在一周期的平均值轨线轨线中心中心用相轨线分析用相轨线分析 点稳
15、定性点稳定性T2T3T4T1PT1 T2 T3 T4x(t)的的“相位相位”领先领先 y(t)模型解释模型解释初值初值相轨线的方向相轨线的方向模型解释模型解释r 食饵增长率食饵增长率d 捕食者死亡率捕食者死亡率b 食饵供养捕食者能力食饵供养捕食者能力捕食者捕食者 数量数量食饵食饵数量数量Pr/ad/ba 捕食者掠取食饵能力捕食者掠取食饵能力捕食者数量与捕食者数量与r成正比成正比,与与a成反比成反比食饵食饵数量与数量与d成正比成正比,与与b成反比成反比模型模型解释解释一次大战期间地中海渔业的捕捞量下降,一次大战期间地中海渔业的捕捞量下降,但是其中但是其中鲨鱼的比例却在增加,为什么?鲨鱼的比例却在
16、增加,为什么?rr-1,dd+1捕捞捕捞战时战时捕捞捕捞rr-2,dd+2,2 1xy食饵食饵(鱼鱼)减少,减少,捕食者捕食者(鲨鱼鲨鱼)增加增加自然环境自然环境 还表明:对还表明:对害虫害虫(食饵食饵)益虫益虫(捕食者捕食者)系统,系统,使用灭两种使用灭两种虫的虫的杀虫剂杀虫剂,会使害虫增加,益虫减少。会使害虫增加,益虫减少。食饵食饵-捕食者模型捕食者模型(Volterra)的缺点与改进的缺点与改进Volterra模型模型改写改写多数多数食饵食饵捕食者系统观察不到周期震荡捕食者系统观察不到周期震荡,而是趋向某个平衡状态而是趋向某个平衡状态,即存在稳定平衡点即存在稳定平衡点加加Logistic
17、项项有有稳定平衡点稳定平衡点 相轨线是封闭曲线,结构不稳定相轨线是封闭曲线,结构不稳定一旦离开某一旦离开某一条闭轨线,就进入另一条闭轨线,不恢复原状。一条闭轨线,就进入另一条闭轨线,不恢复原状。自然界存在的周期性平衡生态系统是结构稳定的,自然界存在的周期性平衡生态系统是结构稳定的,即偏离周期轨道后,内部制约使系统恢复原状。即偏离周期轨道后,内部制约使系统恢复原状。食饵食饵-捕食者模型捕食者模型(Volterra)的缺点与改进的缺点与改进r1=1,N1=20,1=0.1,w=0.2,r2=0.5,2=0.18相轨线趋向极限环相轨线趋向极限环结构稳定结构稳定 两种群模型的几种形式两种群模型的几种形
18、式 相互竞争相互竞争相互依存相互依存弱肉强食弱肉强食8.6 8.6 微分方程稳定性理论简介微分方程稳定性理论简介 在研究实际问题时,我们常常不能直接得出变量之间的关系,但却能容易得出包含变量导数在内的关系式,这就是微分方程.在现实社会中,又有许多变量是离散变化的,如人口数、生产周期与商品价格等,而且离散的运算具有可操作性,差分正是联系连续与离散变量的一座桥梁.不管是微分方程还是差分方程模型,有时无法得到其解析解(必要时,可以利用计算机求其数值解),既使得到其解析解,尚有未知参数需要估计 如果则称平衡点x0是稳定的.否则为不稳定称代数方程 f(x)=0 的实根x=x0为方程(6-1)的平衡点(或
19、奇点).它也是方程(6-1)的解.设 而在实际问题中,讨论问题的解的变化趋势很重要,因此,以下只对其平衡点的稳定性加以讨论.稳定性判别方法稳定性判别方法由于在讨论方程(6-1)的来代替.稳定性时,可用 易知 x0也是方程(6-2)的平衡点.(6-2)的通解为关于x0是否稳定有以下结论:若则x0是稳定的;若则x0是不稳定的.这个结论对这个结论对于于(6-1)也是也是成立的成立的.关于常微分方程组的平衡点及其稳定性,设代数方程组的实根x=x0,y=y0称为方程(6-3)的平衡点,记作P0(x0,y0).它也是方程(6-3)的解.如果则称平衡点P0是稳定的,否则为不稳定 下面给出判别平衡点P0是否稳定的判别准则.设 则当p0且q0时,平衡点P0是稳定的;当p0或q0时,平衡点P0是不稳定的.BACK结束结束