2010数模讲义.ppt

差分方程 定义 数列A={a0，a1，a2，……}的一阶差分是 例1 储蓄存单 考虑一开始价值为1000元的储蓄存单在月利率为1%的条件下的累积价值。 下面的数列表示逐月价值： 其一阶差分为 例1 储蓄存单 例1 储蓄存单 上面的方程可表示无穷多个代数方程，称为动力系统。动力系统能够描述从一个周期到下一个周期的变化。知道了该序列中的某一项，就可以知道通过差分方程算出紧接着它的下一项，但是不能直接算出任意特定的项。 差分方程组——Trafalgar战斗 在1805年的特拉法尔加(Trafalgar)战斗中，由拿破仑指挥的法国、西班牙海军联军和由海军上将纳尔逊指挥的英国海军作战。一开始，法西联军有33艘战舰，而英军有27艘战舰。在一次遭遇战中己方损失都是对方战舰的10%。分数值是有意义的，表示有一艘或多艘战舰不能全力以赴地参加战斗。 Trafalgar战斗 令n表示战斗过程中遭遇战的阶段并定义 Bn=第n阶段英军的战舰数 Fn=第n阶段法军的战舰数 于是在第n阶段的遭遇战后，各方剩余的战舰数为 Bn+1=Bn-0.1Fn, Fn+1=Fn-0.1Bn Trafalgar战斗 Trafalgar战斗 对于全部军力介入的情形，英军是全面失败的。 纳尔逊的分割并各个击破战略：拿破仑军队的33艘战舰基本上是分三个战斗编组（A=3, B=17, C=13）沿一条直线一字排开的。纳尔逊的战略是用13艘英军战舰去迎战战斗编组A（另外14艘备用）。战斗后留存下来的战舰加上备用的14艘战舰去迎战战斗编组B，最后所有剩下的战舰去迎战战斗编组C。 Trafalgar战斗 假设在三次战斗中，每一次战斗中每一方战舰损失都是对方战舰的5%。下面几幅图展示了每次战斗的数值解。 Trafalgar战斗A Trafalgar战斗B Trafalgar战斗B Trafalgar战斗C Trafalgar战斗 从图表可以看出，纳尔逊赢了。 利用分割并各个击败战略的模型的预测结果与历史上真正发生的战斗结果类似。只是法西联军没有参加第三次战斗，而是把约13艘战舰撤回到了法国。 纳尔逊的战略是光辉的，如果没有这样的战略，英军可能会损失掉他们的舰队。 微分方程（组）模型 一、人口模型 二、竞争捕猎模型 三、捕食者食饵模型 微分方模型简介 在研究实际问题时，常常会联系到某些变量的变化率或导数， 这样所得到变量之间的关系式就是微分方模型。微分方程模型反映的是变量之间的间接关系，因此，要得到直接关系，就得求微分方程。 求解微分方程有三种方法： 1 求精确解； 2 求数值解（近似解）； 3 定性理论方法。 人口模型 18世纪末马尔萨斯（Tomas Malthus, 1766-1834）发表了著作《人口原理》，从此激发起人们研究人口增长趋势的兴趣。马尔萨斯在他的这本书里提出了人口按指数增长的模型，并断言人口数量最终将超出食物增长所能提供的容纳能力。虽然马尔萨斯模型的假设忽略了人口增长中的一些重要因素，但是把这个模型用于以后的改进模型的基础是很有价值的。 马尔萨斯的基本观点 在不加抑制的情况下，人口以几何级数增长，而生活资料以算术级数增长。人口增长总是快于生活资料的增长。 当人口扩张到生活资料的极限时，就会出现饥饿、战争、疾病。 贫困是不可避免的，因此社会福利终归无用。 模型问题 假设我们知道在某个给定时刻的人口数量, 比如在 t=t0时刻为P0 ,我们感兴趣的是,预测在未来某个时刻的人口数量P 模型假设 考虑一些有关人口增长的因素. 两个最重要的因素是出生率和死亡率，另外人口增长还受到人口的迁移、生存空间的限制、可利用的资源等因素的影响。在我们的模型中忽略了这些其他因素，只考虑出生率和死亡率。 模型假设 假设出生率和死亡率均为常数，它们的差就是人口的增长率，设为k，这样我们得到 用瞬时变化率来逼近平均变化率，就得到下面的微分方程模型 模型的解 该方程是变量可分离的方程，解之得 对模型的检验 因为 ln(P/P0)=k(t-t0) ，所以如果作出ln(P/P0) 关于图 t-t0图形，就会得到一条经过原点斜率为 k 的直线。得出k 后就可以结合模型的解来预测人口数量。? 注意到，? 我们不得不认为这个模型从长期来讲是不合理的。 指数增长模型的应用及局限性 与19世纪以前欧洲一些地区人口统计数据吻合 适用于19世纪后迁往加拿大的欧洲移民后代 可用于短期人口增长预测 不符合19世纪后多数地区人口增长规律 不能预测较长期的人口增长过程 人口增长率 k 不是常数(逐渐下降) 反映有限增长的改进模型 人口增长到一定数量后，增长率下降的