第1部分微分方程简单模型_精品文档PPT格式课件下载.ppt
《第1部分微分方程简单模型_精品文档PPT格式课件下载.ppt》由会员分享,可在线阅读,更多相关《第1部分微分方程简单模型_精品文档PPT格式课件下载.ppt(22页珍藏版)》请在冰豆网上搜索。
方程,并得出理想单摆运动的周期公式。
从图从图3-1中不难看出,小球所受的合力为中不难看出,小球所受的合力为mgsin,根据根据牛顿第二定律牛顿第二定律可得:
可得:
从而得出两阶微分方程:
(3.1)这是理想单摆应这是理想单摆应满足的运动方程满足的运动方程(3.13.1)是一个两阶非线性方程,不是一个两阶非线性方程,不易求解。
当易求解。
当很小时,很小时,sin,此时,此时,可考察(可考察(3.13.1)的近似线性方程:
)的近似线性方程:
(3.2)由此即可得出由此即可得出(3.23.2)的解为)的解为:
(t)=0cost其中其中当当时时,(t)=0故有故有MQPmg图图3-1(3.13.1)的的近似方程近似方程例例求平面上求平面上过点点(1,3)且每点切且每点切线斜率斜率为横坐横坐标2倍的曲倍的曲线方程方程.解解:
设所求的曲所求的曲线方程方程为由由导数的几何意数的几何意义,应有有即即又由条件又由条件:
曲曲线过(1,3),即即于是得于是得故所求的曲故所求的曲线方程方程为:
导弹追踪问题导弹追踪问题设位于坐标原点的甲舰向位于设位于坐标原点的甲舰向位于x轴上点轴上点A(1,0)处的乙舰处的乙舰发射导弹,导弹头始终对准乙舰如果乙舰以最大的速度发射导弹,导弹头始终对准乙舰如果乙舰以最大的速度v0(常数常数)沿平行于沿平行于y轴的直线行驶,导弹的速度是轴的直线行驶,导弹的速度是5v0,求导求导弹运行的曲线方程乙舰行驶多远时,导弹将它击中?
弹运行的曲线方程乙舰行驶多远时,导弹将它击中?
(解析法)(解析法)由
(1),
(2)消去t,整理得模型:
马尔萨斯(马尔萨斯(MalthusMalthus)模型模型马尔萨斯在分析人口出生与死亡情况的资料后马尔萨斯在分析人口出生与死亡情况的资料后发现,人口净增长率发现,人口净增长率rr基本上是一常数,(基本上是一常数,(rr=bb-dd,bb为出生率,为出生率,dd为死亡率),因而提出了著名的为死亡率),因而提出了著名的人口指数增长模型人口指数增长模型。
分析与建模:
人口的净增长率是一个常数,也就是单位时间内人口增长量与当时人口数成正比。
设t时刻人口数为N(t),t=t0时,N(t0)=N0,则这个方程的解为:
马尔萨斯模型的一个显著特点马尔萨斯模型的一个显著特点:
种群数量翻种群数量翻一番所需的时间是固定的一番所需的时间是固定的。
令种群数量翻一番所需的时间为T,则有:
故故即Malthus模型模型检验模型检验比较历年的人口统计资料,可发现人口增长的实际情况比较历年的人口统计资料,可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符,例如,与马尔萨斯模型的预报结果基本相符,例如,1961年世界人年世界人口数为口数为30.6(即(即3.06109),人口增长率约为),人口增长率约为2%,人口数大,人口数大约每约每35年增加一倍。
检查年增加一倍。
检查1700年至年至1961的的260年人口实际数年人口实际数量,发现两者几乎完全一致,且按马氏模型计算,人口数量量,发现两者几乎完全一致,且按马氏模型计算,人口数量每每34.6年增加一倍,两者也几乎相同。
年增加一倍,两者也几乎相同。
模型预测模型预测假如人口数真能保持每假如人口数真能保持每34.6年增加一倍,那么人口数将年增加一倍,那么人口数将以几何级数的方式增长。
例如,到以几何级数的方式增长。
例如,到2510年,人口达年,人口达21014个,个,即使海洋全部变成陆地,每人也只有即使海洋全部变成陆地,每人也只有9.3平方英尺的活动范围,平方英尺的活动范围,而到而到2670年,人口达年,人口达361015个,只好一个人站在另一人的肩个,只好一个人站在另一人的肩上排成二层了。
上排成二层了。
故故马尔萨斯模型是不完善的。
马尔萨斯模型是不完善的。
几何级数的增长MalthusMalthus模型模型实际上只有在群体总数实际上只有在群体总数不太大时才合理,到总数增大时,不太大时才合理,到总数增大时,生物群体的各成员之间由于有限的生物群体的各成员之间由于有限的生存空间,有限的自然资源及食物生存空间,有限的自然资源及食物等原因,就可能发生生存竞争等现等原因,就可能发生生存竞争等现象。
象。
所以所以MalthusMalthus模型假设的人口模型假设的人口净净增长率不可能始终保持常数,增长率不可能始终保持常数,它应当与人口数量有关。
它应当与人口数量有关。
LogisticLogistic模型模型人口净增长率应当与人口数量有关,即:
人口净增长率应当与人口数量有关,即:
r=r(N)从而有从而有:
(1)r(NN)是未知函数,但根是未知函数,但根据实际背景,它无法用据实际背景,它无法用拟合方法来求拟合方法来求。
为了得出一个有实际意义的为了得出一个有实际意义的模型,我们不妨采用一下工模型,我们不妨采用一下工程师原则。
工程师们在建立程师原则。
工程师们在建立实际问题的数学模型时,总实际问题的数学模型时,总是采用尽可能简单的方法。
是采用尽可能简单的方法。
r(N)最简单的形式是常数,此最简单的形式是常数,此时得到的就是马尔萨斯模型。
时得到的就是马尔萨斯模型。
对马尔萨斯模型的最简单的改对马尔萨斯模型的最简单的改进就是引进一次项(竞争项)进就是引进一次项(竞争项)此时得到微分方程:
此时得到微分方程:
或或
(2)
(2)被称为被称为LogisticLogistic模型或生物总数增长的统计筹算律,是由荷兰数学生模型或生物总数增长的统计筹算律,是由荷兰数学生物学家弗赫斯特(物学家弗赫斯特(VerhulstVerhulst)首先提出的。
一次项系数是负的,因为当种群数首先提出的。
一次项系数是负的,因为当种群数量很大时,会对自身增大产生抑制性,故一次项又被称为竞争项。
量很大时,会对自身增大产生抑制性,故一次项又被称为竞争项。
(22)可改写成:
可改写成:
(3)(3)式还有另一解释,由于空间和资源都是有限的,不可能供养无限增式还有另一解释,由于空间和资源都是有限的,不可能供养无限增长的种群个体,当种群数量过多时,由于人均资源占有率的下降及环境恶长的种群个体,当种群数量过多时,由于人均资源占有率的下降及环境恶化、疾病增多等原因,出生率将降低而死亡率却会提高。
设环境能供养的化、疾病增多等原因,出生率将降低而死亡率却会提高。
设环境能供养的种群数量的上界为种群数量的上界为K(近似地将近似地将K看成常数),看成常数),N表示当前的种群数量,表示当前的种群数量,K-N恰为环境还能供养的种群数量,(恰为环境还能供养的种群数量,(3)指出,种群增长率与两者的乘积成)指出,种群增长率与两者的乘积成正比,正好符合统计规律,得到了实验结果的支持,这就是(正比,正好符合统计规律,得到了实验结果的支持,这就是(3)也被称)也被称为统计筹算律的原因。
为统计筹算律的原因。
对对(33)分离变量:
分离变量:
两边积分并整理得:
令令N(0)=N0,求得:
求得:
故故(33)的满足初始条件的满足初始条件N(0)=N0的解为:
的解为:
(4)易见:
易见:
N(0)=N0,N(t)的图形请看右图的图形请看右图模型检验模型检验用用LogisticLogistic模型来描述种群增长的规律效果如何呢?
模型来描述种群增长的规律效果如何呢?
19451945年克朗比克(年克朗比克(CrombicCrombic)做了一个人工饲养小谷虫的实验,数做了一个人工饲养小谷虫的实验,数学生物学家高斯(学生物学家高斯(EEFFGaussGauss)也做了一个原生物草履虫实也做了一个原生物草履虫实验,实验结果都和验,实验结果都和LogisticLogistic曲线十分吻合。
曲线十分吻合。
大量实验资料表明用大量实验资料表明用LogisticLogistic模型来描述种群的增长,效模型来描述种群的增长,效果还是相当不错的。
例如,高斯果还是相当不错的。
例如,高斯把把5只草履虫放进一个盛有只草履虫放进一个盛有0.5cm3营养液的小试管,他发现,开始时草履虫以每天营养液的小试管,他发现,开始时草履虫以每天230.9%的速率增长,此后增长速度不断减慢,到第五天达到最大量的速率增长,此后增长速度不断减慢,到第五天达到最大量375个,实验数据与个,实验数据与r=2.309,a=0.006157,N(0)=5的的LogisticLogistic曲线:
曲线:
几乎完全吻合,见右图几乎完全吻合,见右图MalthusMalthus模型和模型和LogisticLogistic模型的总结模型的总结MalthusMalthus模型和模型和LogisticLogistic模型模型均为对微分方程(均为对微分方程
(1)所)所作的模拟近似方程。
前一模型假设了种群增长率作的模拟近似方程。
前一模型假设了种群增长率r为一常数,为一常数,(r被称为该种群的内禀增长率)。
后一模型则假设环境只被称为该种群的内禀增长率)。
后一模型则假设环境只能供养一定数量的种群,从而引入了一个竞争项。
能供养一定数量的种群,从而引入了一个竞争项。
用模拟近似法建立微分方程来研究实际问题时必须对用模拟近似法建立微分方程来研究实际问题时必须对求得的解进行检验,看其是否与实际情况相符或基本相符。
求得的解进行检验,看其是否与实际情况相符或基本相符。
相符性越好则模拟得越好,否则就得找出不相符的主要原相符性越好则模拟得越好,否则就得找出不相符的主要原因,对模型进行修改。
因,对模型进行修改。
MalthusMalthus模型与模型与LogisticLogistic模型虽然都是为了研究种群数模型虽然都是为了研究种群数量的增长情况而建立的,但它们也可用来研究其他实际问题,量的增长情况而建立的,但它们也可用来研究其他实际问题,只要这些实际问题的数学模型有相同的微分方程即可。
只要这些实际问题的数学模型有相同的微分方程即可。
以前,美国原子能委员会把浓缩的放射性废料装入密封的圆桶里,然后仍到水深为300英尺的海里。
1问题(这是一场笔墨官司)问题(这是一场笔墨官司):
生态学家和科学家提出生态学家和科学家提出:
圆桶是否会在运输过圆桶是否会在运输过程中破裂而造成放射性污染?
程中破裂而造成放射性污染?
美国原子能委员会:
不会破裂(用实验证明)。
又有几位工程师提出:
圆桶扔到海洋中时是否圆桶扔到海洋中时是否会因与海底碰撞而破裂?
会因与海底碰撞而破裂?
决不会。
放射性核废料处理问题放射性核废料处理问题圆桶与海底的碰撞时的速度会不会超过圆桶与海底的碰撞时的速度会不会超过4040英英尺尺/秒?
秒?
若圆桶与海底碰撞时的速度超过若圆桶与海底碰撞时的速度超过4040英尺英尺/秒时,秒时,就会因碰撞而破裂。
就会因碰撞而破裂。
这几位工程师通过大量的实验证明:
通过建立数学模型来解决这一问题。
一些参数及假设:
假设圆筒下沉时,所受海水的阻力与其速度成正比,即受力分析:
xyGfo2建模与求解建模与求解根据牛顿第二定理可解得