M06稳定性模型导论.ppt

引理：设v(x)是连续的定正函数，则存在H0,使得对于每个Ck∈（０， H）, 曲面v(x)=Ck中有包围０点的封闭曲面Sk，且当　C1C2时，曲面S１被曲面S２包围 证明：设v(x)在球|x|≤a内定正的，令 H=min v(x), |x|=a, Ck∈（０， H）, 我们要证明在球|x|a内，曲面V(x)=Ck 中必存在包围原点的封闭曲面。如引理不真 ，则存在位于球|x|≤a内, 连接0点和球面|x|=a上的某点y的曲线x=x(s), 使得在整条曲线 上 v(x(s))-Ck≠0，但F(x(s))=v(x(s))-Ck是s的连续函数，且 v(0)-Ck0, v(y)-Ck0． 由连续函数性质，对于０＜Ｃk＜Ｈ,必定存在曲线上一点z , 使得v(z)-Ck=0. 矛盾。 6.5 食饵-捕食者模型 种群甲靠丰富的天然资源生存，种群乙(捕食者)靠捕食甲(食饵)为生，形成食饵-捕食者系统，如食用鱼和鲨鱼，美洲兔和山猫，落叶松和蚜虫。 模型的历史背景——一次世界大战期间地中海渔业的捕捞量下降(食用鱼和鲨鱼同时捕捞)，但是其中鲨鱼的比例却增加，为什么？ 食饵（甲）数量 x(t), 捕食者（乙）数量 y(t) 甲独立生存的增长率 r 乙使甲的增长率减小，减小量与 y成正比 乙独立生存的死亡率 d 甲使乙的死亡率减小，减小量与 x成正比 方程(1),(2) 无解析解 食饵-捕食者模型(Volterra) a ~捕食者掠取食饵能力 b ~食饵供养捕食者能力 Volterra模型的平衡点及其稳定性 平衡点 稳定性分析 P点稳定性不能用近似线性方程分析 特征值纯虚 P: 临界状态 鞍点 P′ 不稳定 t x(t) y(t) 0 20.0000 4.0000 0.1000 21.2406 3.9651 0.2000 22.5649 3.9405 0.3000 23.9763 3.9269 … … … 5.1000 9.6162 16.7235 5.2000 9.0173 16.2064 … … … 9.5000 18.4750 4.0447 9.6000 19.6136 3.9968 9.7000 20.8311 3.9587 用数学软件MATLAB求微分方程数值解 x~y 平面上的相轨线 计算结果（数值，图形） x(t), y(t)是周期函数，相图(x,y)是封闭曲线 观察，猜测 x(t), y(t)的周期约为9.6 xmax? 65.5, xmin ? 6, ymax ? 20.5, ymin ? 3.9 用数值积分可算出 x(t), y(t)一周期的平均值： x(t)的平均值约为25, y(t)的平均值约为10。 食饵-捕食者模型(Volterra) 消去dt 因 定正 时 对于每个C0, 相轨线是包围奇点P的封闭曲线, 即对应于周期解 用相轨线分析解的性质—方法一(较简单) 当C从0→∞时, 相轨线从中心奇点往外扩展, 充满第一象限 第六章 稳定性模型 6.1 捕鱼业的持续收获 6.2 军备竞赛 6.3 种群的相互竞争 6.4 种群的相互依存 6.5 食饵-捕食者模型 6.6 微分方程稳定性理论 稳定性模型 对象仍是动态过程，而建模目的是研究时间充分长以后过程的变化趋势 ——平衡状态是否稳定。 不求解微分方程，而是用微分方程稳定性理论研究平衡状态的稳定性。 6.1 捕鱼业的持续收获 再生资源（渔业、林业等）与非再生资源（矿业等） 再生资源应适度开发——在持续稳产前提下实现最大产量或最佳效益。 问题及 分析 在捕捞量稳定的条件下，如何控制捕捞使产量最大或效益最佳。 如果使捕捞量等于自然增长量，渔场鱼量将保持不变，则捕捞量稳定。 背景 产量模型 假设 无捕捞时鱼的自然增长服从 Logistic规律 单位时间捕捞量与渔场鱼量成正比 建模 捕捞情况下渔场鱼量满足 不需要求解x(t), 只需知道x(t)稳定的条件 r~固有增长率, N~最大鱼量 h(x)=Ex, E~捕捞强度 x(t) ~ 渔场鱼量 一阶微分方程的平衡点及其稳定性 一阶非线性（自治）方程 F(x)=0的根x0 ~微分方程的平衡点 设x(t)是方程的解，若从x0 某邻域的任一初值出发，都有 称x0是方程(1)的稳定平衡点 不求x(t), 判断x0稳定性的方法——直接法 (1)的近似线性方程 产量模型 平衡点 稳定性判断 x0 稳定, 可得到稳定产量 x1 稳定, 渔场干枯 E~捕捞强度 r~固有增长率 产量模型 在捕捞量稳定的条件下，控制捕捞强度使产量最大 图解法 P的横坐标 x0~平衡点 y=rx h P x0 y 0 y