数学建模 微分方程模型.ppt

1、2 传染病模型,3 战争模型,4 最优捕鱼问题,1 微分方程模型,微 分 方 程 模 型,1 微分方程模型,一、微分方程模型的建模步骤在自然科学以及工程、经济、医学、体育、生物、 社会等学科中的许多系统，有时很难找到该系统有关 变量之间的直接关系函数表达式，但却容易找到 这些变量和它们的微小增量或变化率之间的关系式， 这时往往采用微分关系式来描述该系统即建立微 分方程模型 。我们以一个例子来说明建立微分方程模 型的基本步骤。,例1 某人的食量是10467（焦/天），其中5038 （焦/天）用于基本的新陈代谢（即自动消耗）。 在健身训练中，他所消耗的热量大约是69 （焦/公斤天）乘以他的体重（公

2、斤）。假设 以脂肪形式贮藏的热量100%地有效，而1公斤脂 肪含热量41868（焦）。 试研究此人的体重随时间变化的规律。,模型分析 在问题中并未出现“变化率”、“导数”这样的关键 词，但要寻找的是体重（记为W）关于时间t的 函数。如果我们把体重W看作是时间t的连续可 微函数，我们就能找到一个含有 的微分方程。,模型假设 1.以W(t)表示t时刻某人的体重，并设一天开始时 人的体重为W0。 2体重的变化是一个渐变的过程。因此可认为 W(t)是关于连续t而且充分光滑的。 3体重的变化等于输入与输出之差，其中输入 是指扣除了基本新陈代谢之后的净食量吸收； 输出就是进行健身训练时的消耗。,模型建立

3、问题中所涉及的时间仅仅是“每天”，由此， 对于“每天”体重的变化=输入-输出。由于考 虑的是体重随时间的变化情况，因此，可得 体重的变化/天=输入/天输出/天。代入具 体的数值，得输入/天 = 10467（焦/天）5038（焦/天） =5429（焦/天），输出/天 = 69（焦/公斤天）（公斤） = 69（焦/天）。,体重的变化/天=W/t（公斤/天）， 当t0时，它等于dW/dt。 考虑单位的匹配， 利用 “公斤/天=（焦/每天）/41868（焦/公斤）”, 可建立如下微分方程模型,模型求解 用变量分离法求解，模型方程等价于积分得,从而求得模型解就描述了此人的体重随时间变化的规律。,现在我们

4、再来考虑一下：此人的体重会达到平衡吗? 显然由W的表达式，当t时，体重有稳定值W 81 。 我们也可以直接由模型方程来回答这个问题。在平衡状态下， W是不发生变化的。所以 这就非常直接地给出了W平衡=81。所以，如果我们需要知道的仅仅是这个平 衡值，就不必去求解微分方程了！,至此，问题已基本上得以解决。 一般地，建立微分方程模型，其方法可归纳为： (1) 根据规律列方程。利用数学、力学、物理、 化学等学科中的定理或许多经过实践或实验检 验的规律和定律，如牛顿运动定律、物质放射 性的规律、曲线的切线性质等建立问题的微分 方程模型。,(3) 模拟近似法。在生物、经济等学科的实际问题中， 许多现象的

5、规律性不很清楚，即使有所了解也是极其 复杂的，常常用模拟近似的方法来建立微分方程模型、 建模时在不同的假设下去模拟实际的现象，这个过程 是近似的，用模拟近似法所建立的微分方程从数学上 去求解或分析解的性质，再去同实际情况对比，看这 个微分方程模型能否刻划、模拟、近似某些实际现象。 本章将结合例子讨论几个不同领域中微分方程模型的 建模方法。,2 传染病模型,问题,描述传染病的传播过程,分析受感染人数的变化规律,预报传染病高潮到来的时刻,预防传染病蔓延的手段,按照传播过程的一般规律，用机理分析方法建立模型,已感染人数 (病人) i(t),每个病人每天有效接触(足以使人致病)人数为,模型1,假设,若

6、有效接触的是病人，则不能使病人数增加,建模,?,模型2,区分已感染者(病人)和未感染者(健康人),假设,1）总人数N不变，病人和健康 人的 比例分别为,2）每个病人每天有效接触人数为, 且使接触的健康人致病,建模, 日 接触率,SI 模型,模型2,tm传染病高潮到来时刻, (日接触率) tm,病人可以治愈！,?,t=tm, di/dt 最大,模型3,传染病无免疫性病人治愈成为健康人，健康人可再次被感染,增加假设,SIS 模型,3）病人每天治愈的比例为, 日治愈率,建模, 日接触率,1/ 感染期, 一个感染期内每个病人的有效接触人数，称为接触数。,模型3,接触数 =1 阈值,感染期内有效接触感染

7、的健康者人数不超过病人数,模型2(SI模型)如何看作模型3(SIS模型)的特例,模型4,传染病有免疫性病人治愈后即移出感染系统，称移出者,SIR模型,假设,1）总人数N不变，病人、健康人和移出者的比例分别为,2）病人的日接触率 , 日治愈率, 接触数 = / ,建模,需建立 的两个方程,模型4,SIR模型,模型4,SIR模型,相轨线 的定义域,在D内作相轨线 的图形，进行分析,模型4,SIR模型,相轨线 及其分析,s(t)单调减相轨线的方向,P1: s01/ i(t)先升后降至0,P2: s01/ i(t)单调降至0,1/阈值,模型4,SIR模型,预防传染病蔓延的手段, (日接触率) 卫生水平

8、,(日治愈率) 医疗水平,传染病不蔓延的条件s01/, 的估计,降低 s0,提高 r0,提高阈值 1/,模型4,SIR模型,被传染人数的估计,记被传染人数比例, 小, s0 1,提高阈值1/降低被传染人数比例 x,s0 - 1/ = ,战争分类：正规战争，游击战争，混合战争,只考虑双方兵力多少和战斗力强弱,兵力因战斗及非战斗减员而减少，因增援而增加,战斗力与射击次数及命中率有关,建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例,第一次世界大战Lanchester提出预测战役结局的模型,3 战争模型,一般模型,每方战斗减员率取决于双方的兵力和战斗力,每方非战斗减员率与本方兵力成

9、正比,甲乙双方的增援率为u(t), v(t),f, g 取决于战争类型,x(t) 甲方兵力，y(t) 乙方兵力,模型假设,模型,正规战争模型,甲方战斗减员率只取决于乙方的兵力和战斗力,双方均以正规部队作战,忽略非战斗减员,假设没有增援,f(x, y)=ay, a 乙方每个士兵的杀伤率,a=ry py, ry 射击率， py 命中率,正规战争模型,为判断战争的结局，不求x(t), y(t)而在相平面上讨论 x 与 y 的关系,平方律 模型,游击战争模型,双方都用游击部队作战,甲方战斗减员率还随着甲方兵力的增加而增加,f(x, y)=cxy, c 乙方每个士兵的杀伤率,c = ry py ry射击

10、率 py 命中率,游击战争模型,线性律 模型,混合战争模型,甲方为游击部队，乙方为正规部队,乙方必须10倍于甲方的兵力,设 x0=100, rx/ry=1/2, px=0.1, sx=1(km2), sry=1(m2),再生资源（渔业、林业等）与非再生资源（矿业等）,再生资源应适度开发在持续稳产前提下实现最大产量或最佳效益。,问题及 分析,在捕捞量稳定的条件下，如何控制捕捞使产量最大或效益最佳。,如果使捕捞量等于自然增长量，渔场鱼量将保持不变，则捕捞量稳定。,背景,4 最优捕鱼问题,产量模型,假设,无捕捞时鱼的自然增长服从 Logistic规律,单位时间捕捞量与渔场鱼量成正比,建模,捕捞情况下

11、渔场鱼量满足,不需要求解x(t), 只需知道x(t)稳定的条件,r固有增长率, N最大鱼量,h(x)=Ex, E捕捞强度,x(t) 渔场鱼量,一阶微分方程的平衡点及其稳定性,一阶非线性（自治）方程,F(x)=0的根x0 微分方程的平衡点,不求x(t), 判断x0稳定性的方法直接法,(1)的近似线性方程,产量模型,稳定性判断,x0 稳定, 可得到稳定产量,x1 稳定, 渔场干枯,E捕捞强度,r固有增长率,产量模型,在捕捞量稳定的条件下，控制捕捞强度使产量最大,图解法,P的横坐标 x0平衡点,P的纵坐标 h产量,产量最大,控制渔场鱼量为最大鱼量的一半,效益模型,假设,鱼销售价格p,单位捕捞强度费用c,单位时间利润,在捕捞量稳定的条件下，控制捕捞强度使效益最大.,求E使R(E)最大,渔场鱼量,收入 T = ph(x) = pEx,支出 S = cE,捕捞过度,封闭式捕捞追求利润R(E)最大,开放式捕捞只求利润R(E) 0,R(E)=0时的捕捞强度(临界强度) Es=2ER,临界强度下的渔场鱼量,捕捞过度,令=0,