微分差分方程建模-320讲座2.pptVIP

基于不同目的的微分方程建模 北京理工大学 王宏洲 建模目的不同，处理方法不同 1、预测：已知当前情况，预测未来情况 x t 数值算法 一般的常微分方程(组)数值算法比较完备；复杂一些的实际问题也比较容易解决。 2、控制：已知当前情况，采取哪些措施才能实现设定的目标 如何调整 ，使得方程存在满足 的解，而且使得 最小或最大。 高速列车的节能控制：如何在确保准点达到的前提下，调整列车速度(牵引力)，使得能耗最低？ 3、稳定性态：波动稳定下来之后会如何？ 时间充分长会如何？事物最终的发展趋势。 比如，商品的价格与其价值的变化关系；食肉动物与草食性动物数量的变化规律；侵入人体的病菌与白血球的数量变化关系；投入一粒石子的池塘水面振幅变化规律。 xm t x 0 x0 xm/2 是已知的常数时，可以分析其几何性态；当 变化时，方程或方程组的几何性态会发生怎样的变化？ t x 事物发展的稳定与不稳定 时间 这些现象在现实中都有实用背景和研究价值 事物的某些特征 差分方程建模 北京理工大学 王宏洲 多数大型哺乳动物有繁殖周期 经济运行也具有一定的阶段性：农产品、股市 植物有固定的繁殖周期 每年多次繁殖的啮齿类动物，也有一定的繁殖间隔 差分方程模型的实际背景 实验室培养果蝇，不考虑自然灾害、疾病、资源突然短缺等因素的影响。另外假设： 第n个繁殖周期，果蝇的总量为x(n)； 每个繁殖周期，果蝇以一定的比例r产生下一代； 每个繁殖周期，果蝇有一定的比例d死亡； 不考虑果蝇的年龄结构、性别结构等。 相邻两个繁殖周期果蝇数量变化 x(n + 1) - x(n) = r x(n) – d x(n) x(n) = x(0)? (1 + r – d)n 1、指数增长模型 假设 果蝇的繁殖率r会随着总数x(n)增加而减少； 果蝇的死亡率d会随着总数x(n)增加而增加； 受生存资源限制，果蝇总数不会无限制增长，上限为N. 同样不考虑果蝇的年龄结构、性别结构等。 将r设为最简单的减函数r(x) = a – b x，将d设为最简单的增函数d(x) = p + qx，其中a、b、p、q均为非负常数。 x(n + 1) = [A – B x(n)] x(n) 2、阻滞增长模型 假设果蝇的生存环境始终保持稳定，而且 第n个繁殖周期，果蝇的总量为x(n)； 果蝇需要k个繁殖周期才能进入成年期； 每个繁殖周期，果蝇以一定的比例r产生下一代； 每个繁殖周期，果蝇有一定的比例d死亡； 不考虑果蝇的性别结构等。 从第n期到第n+1期新增个体数量，应由第n-k个周期时的种群总量决定，因为此后k个周期中新生的个体到第n个周期时均未成年。 x(n + 1) - x(n) = r (1 - d)k x(n - k) - d x(n). 3、考虑年龄因素 4、考虑突发因素 以田鼠为例，假设 第n年，鼠总量为x(n)； 每年鼠以一定的比例r产生下一代； 每年鼠有一定的比例d死亡； 城市每5年举行灭鼠运动杀灭90%，草场鼠数量超过2万只时人为灭杀90%； 不考虑果蝇的年龄结构、性别结构等。 x(n + 1) - x(n) = r x(n) – d x(n) –T(n)·0.9·x(n) 如果采取监控，发现数量达到一定程度就采取杀灭措施，可 以建立如下差分方程模型： x(n + 1) - x(n) = r x(n) – d x(n) –S(x(n))·0.9·x(n) 5、差分方程的稳定性 (1)离散化的阻滞增长模型 得到 下面准备讨论的是k??时， xk的发展趋势问题，即差分方程平衡点的稳定性。 1、求解平衡点 解方程 x = bx(1-x), 得到 x* = 1 – 1/b，x0 = 0 f `(x)=b (1-2x) f `(x* )=2 - b f `(x0 )=b 0 根据稳定性理论， |2 – b| 1，即1 b 3时x* 是稳定平衡点。 也就是说，0 r 2时，y* = N是稳定平衡点。 这与连续型模型中，r 取任何值N都是稳定不同 1 b 3时x* 是稳定的，即 xk? x*. 在 x-y 平面上分别做 y = bx(1-x)和 y = x 的图形 O 1 b 2 2 b 3 O 1/2 1 x y 1/2 图形分析表明： 当1 b 2时，x* 1/2 ，xk? x*的过程是单调的. 当2 b 3时，x* 1/2 ，xk? x*的过程呈现交替. 当b 3时，虽然 x*存在，但不稳定。 计算发现，当3 b 3.449时，虽然{xk}不收敛， {xk}有两个收敛子列： {x2k}和{x