数学建模-微分方程-讲座.ppt

常微分方程理论在数学建模中的应用 北京理工大学数学系 王宏洲 1、微分方程的主要适用范围 2、微分方程模型的分析方法 4、建立微分方程模型的依据 5、案例－物体的运动、振动、受力形变 6、案例－生物的数量变化或密度变化 数学模型 ? 微分方程稳定性方法建模 北京理工大学数学系 王宏洲 事物发展的稳定与不稳定 微分、差分方程稳定性理论简介 常微分方程稳定性理论 差分方程稳定性理论 平衡点稳定的几何特征 一阶微分方程 2、二阶微分方程 二阶微分方程 二阶微分方程 二阶微分方程 生态模型－两种群竞争 考虑两个生物种群竞争同一种有限资源的问题。在自然条件下，适应环境能力弱的种群将趋于灭亡，适应能力强的种群将增长到环境允许的最大数量。种群竞争模型现在已经被广泛地应用到描述企业、国家等社会实体之间的竞争研究中。 下面通过建立模型来解释这种现象，并分析出现各种结局的条件。 2. 稳定性分析（竞争的结局） 2.1 求平衡点 令f(x1,x2)=g(x1,x2)=0，得到四个平衡点： P1(N1,0)，P2(0,N2)， P3(0,0)， 生态模型－食饵捕食者模型 平衡点稳定性分析 稳定性分析 模型结果分析 t 这些现象在现实中都有实用背景和研究价值 参考书：《常微分方程几何理论与分支问题》， 张锦炎，北京大学出版社 涉及“常微分方程定性理论”的其他书籍。 一、常微分方程稳定性理论 1、一阶微分方程 方程右端不显含t t x x0 稳定 不稳定 通常判断平衡点稳定性有两种方法，直接求解法和定性分析法。 定性分析法 1、若方程为线性，即 f(x) = ax + b，则 a 0稳定， a 0不稳定； 2、若方程为非线性，即 x `( t ) = f(x) ，考虑 f ` (x0) 。 f ` (x0) 0稳定， f ` (x0) 0不稳定。 则 a 0稳定，a 0不稳定. 由ax+b=0得到平衡点：x = -b/a 显然，如果 a 0，则方程的所有解在t趋于无穷时都会趋于平衡点，即平衡点是稳定的. 所以讨论二阶微分方程的稳定性往往就归结为对二维一阶方程组的讨论 求方程组的平衡点，即求解 首先将方程组线性化： 其系数矩阵为： 二阶微分方程的稳定性由 p 和 q 的正负决定。 p 0 且 q 0 时平衡点 P0 稳定； p 0 或 q 0 时平衡点 P0 不稳定. 生存空间论 战略空间论 能源竞争 企业间的市场竞争 1. 模型的建立 设同一环境中有甲、乙两个种群，x1(t)、x2(t)分别记t时刻甲、乙种群的数量；r1、r2为各自固有的增长率，N1、N2为各自环境最大容量。据此建立下面的模型： 其中?1，?2 是非常关键的指标，反映一个种群对另 一种群的竞争能力。 r1r2 -(r1+r2) P3 r1r2(1-?1)(1-?2)(1-?1?2)-1 [r1(1-?1)+ r2(1-?2)](1-?1?2)-1 P4 -r1r2(1-?1) -r1(1-?1)+ r2 P2 -r1r2(1-?2) r1-r2(1- ?2 ) P1 稳定条件 q p ?21 (?11) ?11 (?21) 不稳定 ?11 ?21 p0而且q0 2.2 平衡点的稳定性 根据前面的方法不能给出各个平衡点全部的稳定性条件。 下面对?1 和 ?2 分情况讨论平衡点的稳定性条件。 考虑转到相平面上，即在 x1 - x2 平面上研究 方程解沿着 t 增加所表现出的趋势。 x1’(t) = r1x1(1-x1/N1-?1x2/N2) x2’(t) = r2x2(1-?2x1/N1-x2/N2) 可知，在任意时刻， x1(t)和x2(t)是增是减由 ?=1-x1/N1-?1x2/N2 和 ?=1-?2x1/N1-x2/N2 决定。 1、 ?11，?21 S1 S2 S3 O N1/ ?2 N1 x1 x2 N2/ ?1 N2 ?=0 ?=0 这时 ?=0 和 ?=0 将相平 面分为三个区域： S1：x’10, x’20; S2：x’10, x’20; S3：x’10, x’20. t 增加时，所有解都将 趋于P1，所以 P1 是 稳定的。 O N1 x1 x2 N2 2、?11，?21，P2稳定 3、 ?11，?21，P3稳定 O N1 N1/ ?2 x1