《微分方程建模案例2》由会员分享，可在线阅读，更多相关《微分方程建模案例2（20页珍藏版）》请在金锄头文库上搜索。
1、为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。本节将建立几个简单的单种群增长模型，以简略分析一下这方面的问题。,种群的数量本应取离散值，但由于种群数量一般较大，为建立微分方程模型，可将种群数量看作连续变量,由此引起的误差将是十分微小的。,Malthus模型与Logistic模型,模型1 马尔萨斯（Malthus）模型,马尔萨斯在分析人口出生与死亡情况的资料后发现，人口净增长率r基本上是一常数，（r=b-d,b为出生率，d为死亡率）， 既：,马尔萨斯模型的一个显著特点：种群数量翻一番所需的时间是固定的。,Malthus模型实际上只有在群体总数不太大时才
2、合理，到总数增大时，生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。,所以Malthus模型假设的人口净增长率不可能始终保持常数，它应当与人口数量有关。,模型2 Logistic模型,人口净增长率应当与人口数量有关，即： r=r(N),r(N)是未知函数，但根据实际背景，它无法用拟合方法来求 。,为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。,r(N)最简单的形式是常数，此时得到的就是马尔萨斯模型。对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）,(3.9)式还有另一解释
3、，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为K（近似地将K看成常数），N表示当前的种群数量，K-N恰为环境还能供养的种群数量，（3.9）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（3.9）也被称为统计筹算律的原因。,图3-5,对（3.9）分离变量：,两边积分并整理得：,令N(0)=N0，求得：,N(t)的图形请看图3.5,大量实验资料表明用Logistic模型来描述种群的增长，效果还是相当不错的。例如，高斯把5只
4、草履虫放进一个盛有0.5cm3营养液的小试管，他发现，开始时草履虫以每天230.9%的速率增长，此后增长速度不断减慢，到第五天达到最大量375个，实验数据与r=2.309，a=0.006157，N(0)=5的Logistic曲线： 几乎完全吻合，见图3.6。,图3-6,Malthus模型和Logistic模型的总结,Malthus模型和Logistic模型均为对微分方程（3.7）所作的模拟近似方程。前一模型假设了种群增长率r为一常数，（r被称为该种群的内禀增长率）。后一模型则假设环境只能供养一定数量的种群，从而引入了一个竞争项。,用模拟近似法建立微分方程来研究实际问题时必须对求得的解进行检验，
5、看其是否与实际情况相符或基本相符。相符性越好则模拟得越好，否则就得找出不相符的主要原因，对模型进行修改。,Malthus模型与Logistic模型虽然都是为了研究种群数量的增长情况而建立的，但它们也可用来研究其他实际问题，只要这些实际问题的数学模型有相同的微分方程即可。,例5 赝品的鉴定,在第二次世界大战比利时解放以后，荷兰野战军保安机关开始搜捕纳粹同谋犯。他们从一家曾向纳粹德国出卖过艺术品的公司中发现线索，于1945年5月29日以通敌罪逮捕了三流画家范梅格伦（HAVanmeegren），此人曾将17世纪荷兰名画家扬弗米尔（Jan Veermeer）的油画“捉奸”等卖给纳粹德国戈林的中间人。可
6、是，范梅格伦在同年7月12日在牢里宣称：他从未把“捉奸”卖给戈林，而且他还说，这一幅画和众所周知的油画“在埃牟斯的门徒”以及其他四幅冒充弗米尔的油画和两幅德胡斯（17世纪荷兰画家）的油画，都是他自己的作品，这件事在当时震惊了全世界，为了证明自己是一个伪造者，他在监狱里开始伪造弗米尔的油画“耶稣在门徒们中间”，当这项工作接近完成时，范梅格伦获悉自己的通敌罪已被改为伪造罪，因此他拒绝将这幅画变陈，以免留下罪证。,为了审理这一案件，法庭组织了一个由著名化学家、物理学家和艺术史学家组成的国际专门小组查究这一事件。他们用X射线检验画布上是否曾经有过别的画。此外，他们分析了油彩中的拌料（色粉），检验油画中
7、有没有历经岁月的迹象。科学家们终于在其中的几幅画中发现了现代颜料钴兰的痕迹，还在几幅画中检验出了20世纪初才发明的酚醛类人工树脂。根据这些证据，范梅格伦于1947年10月12日被宣告犯有伪造罪，被判刑一年。可是他在监狱中只待了两个多月就因心脏病发作，于1947年12月30日死去。,然而，事情到此并未结束，许多人还是不肯相信著名的“在埃牟斯的门徒”是范梅格伦伪造的。事实上，在此之前这幅画已经被文物鉴定家认定为真迹，并以17万美元的高价被伦布兰特学会买下。专家小组对于怀疑者的回答是：由于范梅格伦曾因他在艺术界中没有地位而十分懊恼，他下决心绘制“在埃牟斯的门徒”，来证明他高于三流画家。当创造出这样的
8、杰作后，他的志气消退了。而且，当他看到这幅“在埃牟斯的门徒”多么容易卖掉以后，他在炮制后来的伪制品时就不太用心了 。这种解释不能使怀疑者感到满意，他们要求完全科学地、确定地证明“在埃牟斯的门徒”的确是一个伪造品。这一问题一直拖了20年，直到1967年，才被卡内基梅伦（Carnegie-Mellon）大学的科学家们 基本上解决。,原理与模型,测定油画和其他岩石类材料的年龄的关键是本世纪初发现的放射性现象。,放射性现象：著名物理学家卢瑟夫在本世纪初发现，某些“放射性”元素的原子是不稳定的，并且在已知的一段时间内，有一定比例的原子自然蜕变而形成新元素的原子，且物质的放射性与所存在的物质的原子数成正比
9、。,用N(t)表示时间t时存在的原子数,则：,用来计算半衰期T:,其解为:,与本问题相关的其他知识:,(1)艺术家们应用白铅作为颜料之一，已达两千年以上。白铅中含有微量的放射铅210，白铅是从铅矿中提炼出来的，而铅又属于铀系，其演变简图如下（删去了许多中间环节）,(3)从铅矿中提炼铅时，铅210与铅206一起被作为铅留下，而其余物质则有9095%被留在矿渣里，因而打破了原有的放射性平衡。,铀238-45亿年-钍234-24天-钋234-6/5分-铀234-257亿年-钍230-8万年-镭226-1600年-氡222-19/5天-钋218-3分-铅214-27分-钋214-铅210-20年-铋2
10、10-5天-钋210-138天-铅206（一种非放射性物质） 注：时间均为半衰期,(2)地壳里几乎所有的岩石中均含有微量的铀。一方面，铀系中的各种放射性物质均在不断衰减，而另一方面，铀又不断地衰减，补充着其后继元素。各种放射性物质（除铀以外）在岩石中处于放射性平衡中。根据世界各地抽样测量的资料，地壳中的铀在铀系中所占平均重量比约为百万分之2.7（一般含量极微）。各地采集的岩石中铀的含量差异很大，但从未发现含量高于23%的。,简化假定：,本问题建模是为了鉴定几幅不超过300年的古画，为了使模型尽可能简单，可作如下假设：,(1)由于镭的半衰期为1600年，经过300年左右，应用微分方程方法不难计算
11、出白铅中的镭至少还有原量的90%，故可以假定，每克白铅中的镭在每分钟里的分解数是一个常数。,建模：,(1)记提炼白铅的时刻为t=0，当时每克白铅中铅210的分子数为y0，由于提炼前岩石中的铀系是处于放射性平衡的，故铀与铅的单位时间分解数相同。可以推算出当时每克白铅中铅210每分钟分解数不能大于30000个。,以上确定了每克白铅中铅分解数的上界，若画上的铅分解数大于该值，说明画是赝品；但若是小于不能断定画一定是真品。,(2)设t时刻1克白铅中铅210含量为y(t)，而镭的单位时间分解数为r（常数），则y(t)满足微分方程：,由此解得：,故：,画中每克白铅所含铅210目前的分解数y(t)及目前镭的
12、分解数r均可用仪器测出，从而可求出y0的近似值，并利用（1）判断这样的分解数是否合理。,Carnegie-Mellon大学的科学家们利用上述模型对部分有疑问的油画作了鉴定，测得数据如下（见表3-1）。,例6 新产品的推广,经济学家和社会学家一直很关心新产品的推销速度问题。怎样建立一个数学模型来描述它，并由此析出一些有用的结果以指导生产呢？以下是第二次世界大战后日本家电业界建立的电饭包销售模型。,设需求量有一个上界，并记此上界为K，记t时刻已销售出的电饭包数量为x(t)，则尚未使用的人数大致为Kx(t)，于是由统计筹算律：,记比例系数为k，则x(t)满足：,此方程即Logistic模型，解为：,对x(t)求一阶、两阶导数：,x(t)0，即x(t)单调增加。,令x(t0)=0，有,当tt0时，x(t)单调减小。,在销出量小于最大需求量的一半时，销售速度是不断增大的，销出量达到最大需求量的一半时，该产品最为畅销，接着销售速度将开始下降。,所以初期应采取小批量生产并加以广告宣传；从有20%用户到有80%用户这段时期，应该大批量生产；后期则应适时转产，这样做可以取得较高的经济效果。,知识回顾Knowledge Review,
