微分方程模型 人口的预测.ppt

第5章 微分方程模型 §5.6 人口的预测 （Malthus模型与Logistic模型） § 5.6 人口的预测 § 5.6 人口的预测 新产品的推广 一、可分离变量微分方程 二、一阶线性微分方程 二、伯努利（Bernoulli）方程 五、常系数线性微分方程组解法步骤: 例1. 常数变易法 作变换 求解一阶非齐次线性微分方程 则 即 两边积分，得 故非齐次线性微分方方程的通解为 对应齐次方程通解 非齐次方程特解 把齐次方程通解中的常数变易为待定函数的方法. 解： 例1 伯努利方程的标准形式: 令 求出此方程通解后,换回原变量即得伯努利方程的通解. 除方程两边 , 得 解法: (线性方程) 的通解. 解: 令 则方程变形为 其通解为 将 代入, 得原方程通解: 例1 求方程 四、二阶线性微分方程 时, 称为非齐次方程; 时, 称为齐次方程. 齐次方程通解 非齐次方程特解 解的结构 线性无关） 特征根: (1) 当 时, 通解为 (2) 当 时, 通解为 (3) 当 时, 通解为 二阶常系数齐次线性微分方程 特征方程: ? 为特征方程的 k (＝0, 1, 2) 重根, 则设特解为 为特征方程的 k (＝0, 1 )重根, 则设特解为 二阶常系数非齐次线性微分方程 的通解. 解: 本题 特征方程为 其根为 对应齐次方程的通解为 设非齐次方程特解为 比较系数, 得 因此特解为 代入方程得 所求通解为 例1 的通解. 解: 特征方程为 其根为 对应齐次方程的通解为 比较系数, 得 因此特解为 代入方程: 所求通解为 为特征方程的单根 , 因此设非齐次方程特解为 例2 * * * 微分方程模型 在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题. 求出方程的解 ——求出未知函数的解析表达式 ——利用各种数值解法、数值软件（如Matlab）求近似解 不必求出方程的解 ——根据微分方程的理论研究某些性质，或它的变化趋势 为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。 本节将建立几个简单的单种群增长模型，以简略分析一下这方面的问题。一般生态系统的分析可以通过一些简单模型的复合来研究，大家若有兴趣可以根据生态系统的特征自行建立相应的模型。 美丽的大自然 种群的数量本应取离散值，但由于种群数量一般较大，为建立微分方程模型，可将种群数量看作连续变量，甚至允许它为可微变量，由此引起的误差将是十分微小的。 离散化为连续，方便研究 世界人口 年 1625 1830 1930 1960 1974 1987 1999 人口（亿） 5 10 20 30 40 50 60 中国人口 年 1908 1933 1953 1964 1982 1990 2000 人口（亿） 3 4.7 6 7.2 10.3 11.3 12.95 研究人口变化规律 控制人口过快增长 做出较准确的预报 建立人口数学模型 最简单的人口增长模型—常用公式 记今年人口为 x0, k 年后人口为 xk,，年增长率为 r 则 模型1 马尔萨斯（Malthus）模型——1798年提出 假设：单位时间内人口的增长量与当时的人口成正比。 （2） （1）的解为： 符号： 则 （1） 于是x（t）满足如下微分方程： r--人口增长率（常数） （可分离变量微分方程） 单位时间内人口的增长量 马尔萨斯模型的一个显著特点：种群数量翻一番所需的时间是固定的。 令种群数量翻一番所需的时间为T，则有： 故 （2） 当r>0时，表明人口将按指数规律无限增长，因此又称为人口指数模型。 与常用公式的一致 模型检验 用P164给出的近两个世纪的美国人口统计数据（以百万作单位），对模型作检验。 参数估计： r，x0可用已知数据利用线性最小二乘法进行估计 （2） （2）式两边取对数，得： （3） 以1790－1900年的数据拟合（3）式，用Matlab软件计算得：r＝0.2743/10年， 所有散点到曲线的距离平方和最小 Matlab计算示范 以1790-1900年共计12个数据为例进行拟合： t=[0：11]；