微分方程建模学习.docVIP

微分方程建模 一般说来，微分方程建模的方法大概能够分为以下的几个步骤： 1.依据实质问题的要求确立要研究的量，包含自变量、未知函数、必需的参数等以及它 们各自的变化区间； 2.列方程。能够在合理假定的前提下，利用导数表示斜率、速度、变化率的实质意义， 依据一些基本定理（几何的、物理的、化学的或生物学的等等）或规律，找出未知函数的导 数（或微分）与有关各量之间的等量关系式，成立微分方程并确立定解条件（注：假如没有 现成的定理可供利用，也能够用微元剖析法与模拟近似法列出微分方程）； 3.解微分方程； 4.对模型的合用性作出评论，即用已知的数据查验微分方程的解能否与实质符合。若结 果与实质存在必定的差距，则还要对方程进行修正和调整，直到得出较满意的结果为止。 下边，我们就经过一些实例说明微分方程建模的详细步骤。 一.增添模型 在自然界和社会的经济活动中，很多量的变化都按照着一个基本的规律：任一单位时间 的增量都与该量自己当时的大小成正比。运用这一基本规律，就能够成立起各种各种的增添 模型。 1.马尔萨斯人口模型 严格地讲，议论人口问题所成立的模型应属于失散型模型。 但在人口基数很大的状况下， 忽然增添或减少的不过单一的个体或少量几个个体， 有关于全体数目而言，这类改变量是极 其细小的，所以，我们能够近似地假定人口随时间连续变化甚至是可微的。 这样，我们便可 以采纳微分方程的工具来研究这一问题。 最早研究人口问题的是英国的经济系家马尔萨斯 (Malthus)（1766—1834）。他依据百余 年的人口资料，经过专注研究，在 1798年发布的《人口论》中第一提出了人口增添模型。 他的基本假定是：任一单位时刻人口的增添量与当时的人口总数成正比， 且比率系数为常数。 于是，设t时刻的人口总数为y(t)，则单位时间人口的增添量即为 y(t t) y(t) t 依据基本假定，有 y(t t) y(t) y(t) （r为比率系数） r t 令t0，可得微分方程 dy r y （4.1） dt 这就是有名的马尔萨斯人口方程。 若假定t t0时的人口总数为 y0，则不难求得该方程的特 解为 y r(tt0) (4.2) y0e 即任一时刻的人口总数都按照指数规律向上增添。人们曾用这个公式对 1700—1961达二百 六十余年世界的人口资料进行了查验，发现计算结果与人口的实质状况居然是惊人的符合！ 放射性元素衰变模型 放射性元素的质量随时间的推移而渐渐减少（负增添），这类现象称为衰变。由物理学 定律知，放射性元素任一时刻的衰变速度与该时刻放射性元素的质量成正比。依据这一原理， 我们也能够经过微分方程研究放射性元素衰变的规律。 设放射性元素t时刻的质量m m(t)，则其衰变速度就是 dm，于是可得 dt dm (4.3) m dt 此中0是比率常数，可由该元素的半衰期（质量演变到一半所需要的时间）确立；前 置负号表示放射性元素的质量 m是随时刻t递减的。 假如在初始时刻（t 0）放射性元素的质量 m m0，则可求得该方程的特解为 m(t) m e t (4.4) 0 这说明放射性元素的质量也是随时刻按指数规律递减的。 为了能将求得的放射性元素衰变规律应用于实质， 还一定确立上式中的比率常数 。这 时，我们能够假定放射性元素的半衰期为 T，进而有 m0 m0 e T 2 解之，得 ln2 4.4）式又能够表示为 ，于是反应放射性元素衰变规律的（ T ln2 t m(t) m0 e T (4.5) 并由此可解得 T ln m0 (4.6) t m(t) ln2 它所反应的是放射性元素由初始时刻的质量 m0衰减到m(t)所需要的时间。 放射性元素的衰变规律常被考古、 地质方面的专家用于测定文物和地质的年月， 此中最 常用的是14C（碳-12的同位素）测定法。 这类方法的原理是：大气层在宇宙射线不断的轰 击下所产生的中子与氮气作用生成了拥有放射性的14C，并进一步氧化为二氧化碳，二氧化 碳被植物所汲取，而动物又以植物为食品，于是放射性碳就被带到了各种动植物的体。关于 拥有放射性的14C来说，无论是存在于空气中仍是生物体，它都在不断地演变。因为活着的 生物经过新代不断地摄入14C，因此使得生物体的14C与空气中的14C有相同的百分含量； 一旦生物死亡以后，跟着新代的停止，尸体的14C就会不断地演变而渐渐减少，所以依据14C 演变减少许的变化状况并利用（4.6）式，就能够判断生物死亡的时间。 下边，我们就来看一个运用14C测定法确立年月的详细实例： 1972年8月，出土了马王堆一号墓（注：出土时因墓中女尸历经千年而未腐以前惊动 世界）。经测定，出土的木炭标本中  14C的均匀原子演变速度为  29.78次/分，而新砍伐烧成 的木炭中  14C的均匀原子演变速