理学微分方程建模.pptx

《理学微分方程建模.pptx》由会员分享，可在线阅读，更多相关《理学微分方程建模.pptx（93页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、例1 （理想单摆运动）建立理想单摆运动满足的微分方程，并得出理想单摆运动的周期公式。从图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页/共93页例2 我方巡逻艇发现敌方潜水艇。与此同时敌方潜水艇也发现了我方巡逻艇，并迅速下潜逃逸。

2、设两艇间距离为6060哩，潜水艇最大航速为3030节而巡逻艇最大航速为6060节，问巡逻艇应如何追赶潜水艇。这一问题属于对策问题，较为复杂。讨论以下简单情形：敌潜艇发现自己目标已暴露后，立即下潜，并沿着直 线方向全速逃逸，逃逸方向我方不知。设巡逻艇在A处发现位于B处的潜水艇，取极坐标，以B为极点，BA为极轴，设巡逻艇追赶路径在此极坐标下的方程为r=r()，见图3-2。BAA1drdsd图3-2由题意，故ds=2dr图3-2可看出，第2页/共93页故有:即:（3.3）解为：（3.4）先使自己到极点的距离等于潜艇到极点的距离,然后按（3.4）对数螺线航行，即可追上潜艇。追赶方法如下：第3页/共93

3、页例3 一个半径为Rcm的半球形容器内开始时盛满了水，但由于其底部一个面积为Scm2的小孔在t=0时刻被打开，水被不断放出。问：容器中的水被放完总共需要多少时间？解:以容器的底部O点为 原点，取坐标系如图3.3所示。令h(t)为t时刻容器中水的高度，现建立h(t)满足的微分方程。设水从小孔流出的速度为v(t)，由力学定律，在不计水的内部磨擦力和表面张力的假定下，有：因体积守衡，又可得：易见：故有：即：这是可分离变量的一阶微分方程，得 RxySO图图3-3hr第4页/共93页例4 一根长度为l的金属杆被水平地夹在两端垂直的支架上，一端的温度恒为T1，另一端温度恒为T2，（T1、T2为常数，T1

4、T2）。金属杆横截面积为A，截面的边界长度为B，它完全暴露在空气中，空气温度为T3，（T3钍234-24天-钋234-6/5分-铀234-257亿年-钍230-8万年-镭226-1600年-氡222-19/5天-钋218-3分-铅214-27分-钋214-铅210-20年-铋210-5天-钋210-138天-铅206（一种非放射性物质）注：时间均为半衰期 (2)地壳里几乎所有的岩石中均含有微量的铀。一方面，铀系中的各种放射性物质均在不断衰减，而另一方面，铀又不断地衰减，补充着其后继元素。各种放射性物质（除铀以外）在岩石中处于放射性平衡中。根据世界各地抽样测量的资料，地壳中的铀在铀系中所占平均重

5、量比约为百万分之2.7（一般含量极微）。各地采集的岩石中铀的含量差异很大，但从未发现含量高于23%的。第17页/共93页简化假定：本问题建模是为了鉴定几幅不超过300年的古画，为了使模型尽可能简单，可作如下假设：(1)由于镭的半衰期为1600年，经过300年左右，应用微分方程方法不难计算出白铅中的镭至少还有原量的90%，故可以假定，每克白铅中的镭在每分钟里的分解数是一个常数。(2)铅210210的衰变为：铅铅210T=22年年钋钋210铅铅206T=138天天若画为真品，颜料应有300年左右或300年以上的历史，容易证明：每克白铅中钋210的分解数等于铅210的分解数（相差极微，已无法区别）。

6、可用前者代替后者，因钋的半衰期较短，易于测量。第18页/共93页建模：(1)记提炼白铅的时刻为t=0，当时每克白铅中铅210的分子数为y0，由于提炼前岩石中的铀系是处于放射性平衡的，故铀与铅的单位时间分解数相同。可以推算出当时每克白铅中铅210每分钟分解数不能大于30000个。若则（个）这些铀约重 （克）即每克白铅约含0.040.04克铀，含量为4%4%以上确定了每克白铅中铅分解数的上界，若画上的铅分解数大于该值，说明画是赝品；但若是小于不能断定画一定是真品。第19页/共93页 (2)设t时刻1克白铅中铅210含量为y(t)，而镭的单位时间分解数为r（常数），则y(t)满足微分方程：由此解得：

7、故：画中每克白铅所含铅210目前的分解数y(t)及目前镭的分解数r均可用仪器测出，从而可求出y0的近似值，并利用（1）判断这样的分解数是否合理。第20页/共93页Carnegie-MellonCarnegie-Mellon大学的科学家们利用上述模型对部分有疑问的油画作了鉴定，测得数据如下（见表3-13-1）。油画名称210210分解数（个/分）镭226226分解数（个/分）1 1、在埃牟斯的门徒 8.5 8.50.80.82 2、濯足12.612.60.260.263 3、看乐谱的女人10.310.30.30.34 4、演奏曼陀琳的女人8.28.20.170.175 5、花边织工1.51.51

8、.41.46 6、笑女5.25.26.06.0计算y0（个/分）98050980501571301571301273401273401022501022501274.81274.8-10181-10181表3-1 3-1 对“在埃牟斯的门徒”，y y0 09805098050（个/每克每分钟），它必定是一幅伪造品。类似可以判定（2 2），（3 3），（4 4）也是赝品。而（5 5）和（6 6）都不会是几十年内伪制品，因为放射性物质已处于接近平衡的状态，这样的平衡不可能发生在十九世纪和二十世纪的任何作品中。判定结果：第21页/共93页 利用放射原理，还可以对其他文物的年代进行测定。利用放射原理，

9、还可以对其他文物的年代进行测定。例如对有机物（动、植物）遗体，考古学上目前流行的测例如对有机物（动、植物）遗体，考古学上目前流行的测定方法是放射性碳定方法是放射性碳1414测定法，这种方法具有较高的精确度，测定法，这种方法具有较高的精确度，其基本原理是：由于大气层受到宇宙线的连续照射，空气其基本原理是：由于大气层受到宇宙线的连续照射，空气中含有微量的中微子，它们和空气中的氮结合，形成放射中含有微量的中微子，它们和空气中的氮结合，形成放射性碳性碳1414（C C1414）。有机物存活时，它们通过新陈代谢与外界）。有机物存活时，它们通过新陈代谢与外界进行物质交换，使体内的进行物质交换，使体内的C

10、C1414处于放射性平衡中。一旦有机处于放射性平衡中。一旦有机物死亡，新陈代谢终止，放射性平衡即被破坏。因而，通物死亡，新陈代谢终止，放射性平衡即被破坏。因而，通过对比测定，可以估计出它们生存的年代。例如，过对比测定，可以估计出它们生存的年代。例如，19501950年年在巴比伦发现一根刻有在巴比伦发现一根刻有HammurabiHammurabi王朝字样的木炭，经测定，王朝字样的木炭，经测定，其其C C1414衰减数为衰减数为4.094.09个个/每克每分钟，而新砍伐烧成的木炭中每克每分钟，而新砍伐烧成的木炭中C C1414衰减数为衰减数为6.686.68个个/每克每分钟，每克每分钟，C C14

11、14的半衰期为的半衰期为55685568年，年，由此可以推算出该王朝约存在于由此可以推算出该王朝约存在于3900-40003900-4000年前。年前。第22页/共93页例例6 6 新产品的推广新产品的推广 经济学家和社会学家一直很关心新产品的推销速度问题。怎样建立一个数学模型来描述它，并由此析出一些有用的结果以指导生产呢？以下是第二次世界大战后日本家电业界建立的电饭包销售模型。设需求量有一个上界，并记此上界为K，记t时刻已销售出的电饭包数量为x(t)，则尚未使用的人数大致为Kx(t)，于是由统计筹算律：记比例系数为k k，则x(t)满足：此方程即LogisticLogistic模型，解为：还

12、有两个奇解:x=0和和x=K 对x(t)求一阶、两阶导数：第23页/共93页x(t)0，即x(t)单调增加。令x(t0)=0，有当tt0时，x(t)单调减小。在销出量小于最大需求量的一半时，销售速度是不断增大的，销出量达到最大需求量的一半时，该产品最为畅销，接着销售速度将开始下降。所以初期应采取小批量生产并加以广告宣传；从有20%20%用户到有80%80%用户这段时期，应该大批量生产；后期则应适时转产，这样做可以取得较高的经济效果。第24页/共93页3.33.3 为什么要用三级火箭来发射人造卫星为什么要用三级火箭来发射人造卫星构造数学模型，以说明为什么不能用一级火箭而必须用多级火箭来发射人造卫

13、星？为什么一般都采用三级火箭系统？1 1、为什么不能用一级火箭发射人造卫星?（1 1）卫星能在轨道上运动的最低速度 假设：（i i）卫星轨道为过地球中心的某一平面上的圆，卫星 在此轨道上作匀速圆周运动。（iiii）地球是固定于空间中的均匀球体，其它星球对卫 星的引力忽略不计。分析：根据牛顿第三定律，地球对卫星的引力为:在地面有:得:k=gR2 R R为地球半径，约为64006400公里 故引力:假设(ii)第25页/共93页dmm-dmvu-v假设(i)卫星所受到的引力也就是它作匀速圆周运动的向心力故又有:从而:设g=9.81g=9.81米/秒2 2，得:卫星离地面高度 (公里)卫星速度 (公

14、里/秒)100100200200400400600600800800100010007.807.807.697.697.587.587.477.477.377.377.867.86（2 2）火箭推进力及速度的分析 假设：火箭重力及空气阻力均不计 分析：记火箭在时刻t的质量和速度分别为m(t)和(t)有：记火箭喷出的气体相对于火箭的速度为u（常数），由动量守恒定理：0 0和m m0 0一定的情况下，火箭速度(t)(t)由喷发速度u u及质量比决定。故：由此解得：(3.11)第26页/共93页（2 2）火箭推进力及速度的分析 现将火箭卫星系统的质量分成三部分：（i）mP（有效负载，如卫星）（ii）

15、mF（燃料质量）（iii）mS（结构质量如外壳、燃料容器及推进器）。最终质量为mP+mS，初始速度为0，所以末速度：根据目前的技术条件和燃料性能，u只能达到3公里/秒，即使发射空壳火箭，其末速度也不超过6.6公里/秒。目前根本不可能用一级火箭发射人造卫星火箭推进力在加速整个火箭时，其实际效益越来越低。如果将结构质量在燃料燃烧过程中不断减少，那么末速度能达到要求吗？第27页/共93页2 2、理想火箭模型 假设：记结构质量mS在mS+mF中占的比例为，假设火箭能随时抛弃无用的结构，结构质量与燃料质量以与（1-）的比例同时减少。建模:由 得到：解得：理想火箭与一级火箭最大的区别在于，当火箭燃料耗尽时

16、，结构质量也逐渐抛尽，它的最终质量为mP，所以最终速度为：只要m0足够大，我们可以使卫星达到我们希望它具有的任意速度。考虑到空气阻力和重力等因素，估计（按比例的粗略估计）发射卫星要使=10.5公里/秒才行，则可推算出m0/mp约为51,即发射一吨重的卫星大约需要50吨重的理想火箭 第28页/共93页3 3、理想过程的实际逼近多级火箭卫星系统 记火箭级数为n，当第i级火箭的燃料烧尽时，第i+1级火箭立即自动点火，并抛弃已经无用的第i级火箭。用mi表示第i级火箭的质量，mP表示有效负载。先作如下假设：（i）设各级火箭具有相同的,即i级火箭中mi为结构质量，（1-）mi为燃料质量。（ii）设燃烧级初

17、始质量与其负载质量之比保持不变，并记比值为k k。考虑二级火箭：由3.11式，当第一级火箭燃烧完时，其末速度为：当第二级火箭燃尽时，末速度为：该假设有点强加的味道，先权作讨论的方便吧第29页/共93页又由假设（ii），m2=kmP，m1=k(m2+mP)，代入上式，仍设u=3公里/秒，且为了计算方便，近似取=0.1，则可得：要使2=10.5公里/秒，则应使:即k11.2，而:类似地，可以推算出三级火箭：在同样假设下:要使3=10.5公里/秒，则(k+1)/(0.1k+1)3.21，k3.25，而（m1+m2+m3+mP）/mP77。三级火箭比二级火箭几乎节省了一半 是否三级火箭就是最省呢？最简

18、单的方法就是对四级、五级等火箭进行讨论。第30页/共93页考虑N N级火箭：记n级火箭的总质量（包含有效负载mP）为m0，在相同的假设下可以计算出相应的m0/mP的值，见表3-2n（级数）（级数）1 2 3 4 5 （理想）（理想）火箭质量（吨）火箭质量（吨）/149 77 65 60 50表3-2由于工艺的复杂性及每节火箭都需配备一个推进器，所以使用四级或四级以上火箭是不合算的，三级火箭提供了一个最好的方案。当然若燃料的价钱很便宜而推进器的价钱很贵切且制作工艺非常复杂的话，也可选择二级火箭。第31页/共93页4 4、火箭结构的优化设计 3 3中已经能说过假设(ii)(ii)有点强加的味道；现

19、去掉该假设，在各级火箭具有相同的粗糙假设下，来讨论火箭结构的最优设计。W1=m1+mn+mP W2=m2+mn+mPWn=mn+mPWn+1=mP记应用（3.113.11）可求得末速度：记则又问题化为，在n一定的条件下，求使k1 k2kn最小 解条件极值问题：或等价地求解无约束极值问题：可以解出最优结构设计应满足：火箭结构优化设计讨论中我们得到与假设（ii）相符的结果，这说明前面的讨论都是有效的！第32页/共93页3.43.4 药物在体内的分布药物在体内的分布 何为房室系统？在用微分方程研究实际问题时，人们常常采用一种叫“房室系统”的观点来考察问题。根据研究对象的特征或研究的不同精度要求，我们

20、把研究对象看成一个整体（单房室系统）或将其剖分成若干个相互存在着某种联系的部分（多房室系统）。房室具有以下特征：它由考察对象均匀分布而成，房室中考察对象的数量或浓度（密度）的变化率与外部环境有关，这种关系被称为“交换”且交换满足着总量守衡。在本节中，我们将用房室系统的方法来研究药物在体内的分布。在下一节中，我们将用多房室系统的方法来研究另一问题。交换环境内部单房室系统均匀分布第33页/共93页 药物的分解与排泄（输出）速率通常被认为是与药物当前的浓度成正比的，即：药物分布的单房室模型 单房室模型是最简单的模型，它假设：体内药物在任一时刻都是均匀分布的，设t时刻体内药物的总量为x(t)；系统处于

21、一种动态平衡中，即成立着关系式：药物的输入规律与给药的方式有关。下面，我们来研究一下在几种常见的给药方式下体内药体的变化规律。机体环境药物总量图3-8 假设药物均匀分布第34页/共93页情况1 1 快速静脉注射机体环境只输出不输入房室其解为：药物的浓度：与放射性物质类似，医学上将血浆药物浓度衰减一半所需的时间称为药物的血浆半衰期：负增长率的Malthus模型 在快速静脉注射时，总量为D的药物在瞬间被注入体内。设机体的体积为V，则我们可以近似地将系统看成初始总量为D，浓度为D/V，只输出不输入的房室，即系统可看成近似地满足微分方程：（3.12）第35页/共93页情况2 2 恒速静脉点滴 机体环境

22、恒定速率输入房室药物似恒速点滴方式进入体内，即:则体内药物总量满足：（x(0)=0）（3.13）这是一个一阶常系数线性方程，其解为：或易见：称为稳态血药浓度 对于多次点滴，设点滴时间为T1，两次点滴之间的间隔时间设为T2，则在第一次点滴结束时病人体内的药物浓度可由上式得出。其后T2时间内为情况1 1。故：（第一次）0tT1 T1tT1+T2 类似可讨论以后各次点滴时的情况，区别只在初值上的不同。第二次点滴起，患者 体内的初始药物浓度不为零。第36页/共93页情况3 3 口服药或肌注 y(t)x(t)K1yK1x环境机体外部药物 口服药或肌肉注射时，药物的吸收方式与点滴时不同，药物虽然瞬间进入了

23、体内，但它一般都集中与身体的某一部位，靠其表面与肌体接触而逐步被吸收。设药物被吸收的速率与存量药物的数量成正比，记比例系数为K1，即若记t时刻残留药物量为y(t)，则y满足：D为口服或肌注药物总量 因而：所以：解得：从而药物浓度：第37页/共93页 图3-93-9给出了上述三种情况下体内血药浓度的变化曲线。容易看出，快速静脉注射能使血药浓度立即达到峰值，常用于急救等紧急情况；口服、肌注与点滴也有一定的差异，主要表现在血药浓度的峰值出现在不同的时刻，血药的有效浓度保持时间也不尽相同。图3-9 我们已求得三种常见给药方式下的血药浓度C C(t t)，当然也容易求得血药浓度的峰值及出现峰值的时间，因

24、而，也不难根据不同疾病的治疗要求找出最佳治疗方案。第38页/共93页 新药品、新疫苗在临床应用前必须经过较长时间的基础研究、小量试制、新药品、新疫苗在临床应用前必须经过较长时间的基础研究、小量试制、中间试验、专业机构评审及临床研究。当一种新药品、新疫苗研制出来后，中间试验、专业机构评审及临床研究。当一种新药品、新疫苗研制出来后，研究人员必须用大量实验搞清它是否真的有用，如何使用才能发挥最大效用，研究人员必须用大量实验搞清它是否真的有用，如何使用才能发挥最大效用，提供给医生治病时参考。在实验中研究人员要测定模型中的各种参数，搞清提供给医生治病时参考。在实验中研究人员要测定模型中的各种参数，搞清血

25、药浓度的变化规律，根据疾病的特点找出最佳治疗方案（包括给药方式、血药浓度的变化规律，根据疾病的特点找出最佳治疗方案（包括给药方式、最佳剂量、给药间隔时间及给药次数等），这些研究与试验据估计最少也需最佳剂量、给药间隔时间及给药次数等），这些研究与试验据估计最少也需要数年时间。在要数年时间。在20032003年春夏之交的年春夏之交的SARSSARS（非典）流行期内，有些人希望医药（非典）流行期内，有些人希望医药部门能赶快拿出一种能治疗部门能赶快拿出一种能治疗SARSSARS的良药或预防的良药或预防SARSSARS的有效疫苗来，但这只能的有效疫苗来，但这只能是一种空想。是一种空想。SARSSARS的

26、突如其来，形成了的突如其来，形成了“外行不懂、内行陌生外行不懂、内行陌生”的情况。国的情况。国内权威机构一度曾认为这是内权威机构一度曾认为这是“衣原体衣原体”引起的肺炎，可以用抗生素控制和治引起的肺炎，可以用抗生素控制和治疗。但事实上，抗生素类药物对疗。但事实上，抗生素类药物对SARSSARS的控制与治疗丝毫不起作用。以钟南山的控制与治疗丝毫不起作用。以钟南山院士为首的广东省专家并不迷信权威，坚持认为院士为首的广东省专家并不迷信权威，坚持认为SARSSARS是病毒感染引起的肺炎，是病毒感染引起的肺炎，两个月后（两个月后（4 4月月1616日），世界卫生组织正式确认日），世界卫生组织正式确认SA

27、RSSARS是冠状病毒的一个变种引是冠状病毒的一个变种引起的非典型性肺炎（注：这种确认并非是由权威机构定义的，而是经对猩猩起的非典型性肺炎（注：这种确认并非是由权威机构定义的，而是经对猩猩的多次实验证实的）。发现病原体尚且如此不易，要攻克难关，找到治疗、的多次实验证实的）。发现病原体尚且如此不易，要攻克难关，找到治疗、预防的办法当然就更困难了，企图几个月解决问题注定只能是一种不切实际预防的办法当然就更困难了，企图几个月解决问题注定只能是一种不切实际的幻想。的幻想。第39页/共93页 上述研究是将机体看成一个均匀分布的同质单元，故被称单房室模型，但机体事实上并不是这样。药物进入血液，通过血液循环

28、药物被带到身体的各个部位，又通过交换进入各个器官。因此，要建立更接近实际情况的数学模型就必须正视机体部位之间的差异及相互之间的关联关系，这就需要多房室系统模型。IIIk12k21两房室系统图3-10 图3-103-10表示的是一种常见的两房室模型，其间的k k1212表示由室I I渗透到室IIII的变化率前的系数，而k k2121则表示由室IIII返回室I I的变化率前的系数，它们刻划了两室间的内在联系，其值应当用实验测定，使之尽可能地接近实际情况。当差异较大的部分较多时，可以类似建立多房室系统，即N N房室系统第40页/共93页hy3.53.5 传染病模型传染病模型 传染病是人类的大敌，通过

29、疾病传播过程中若干重要因素之间的联系建立微分方程加以讨论，研究传染病流行的规律并找出控制疾病流行的方法显然是一件十分有意义的工作。在本节中，我们将主要用多房室系统的观点来看待传染病的流行，并建立起相应的多房室模型。医生们发现，在一个民族或地区，当某种传染病流传时，波及到的总人数大体上保持为一个常数。即既非所有人都会得病也非毫无规律，两次流行（同种疾病）的波及人数不会相差太大。如何解释这一现象呢？试用建模方法来加以证明。问题的提出：第41页/共93页 设某地区共有n+1人，最初时刻共有i人得病，t时刻已感染（infective）的病人数为i(t)，假定每一已感染者在单位时间内将疾病传播给k个人（

30、k称为该疾病的传染强度），且设此疾病既不导致死亡也不会康复模型1 此模型即MalthusMalthus模型，它大体上反映了传染病流行初期的病人增长情况，在医学上有一定的参考价值，但随着时间的推移，将越来越偏离实际情况。已感染者与尚未感染者之间存在着明显的区别，有必要将人群划分成已感染者与尚未感染的易感染，对每一类中的个体则不加任何区分，来建立两房室系统。则可导出：故可得：（3.15）第42页/共93页模型2 记t时刻的病人数与易感染人数（susceptible）分别为i(t)与s(t)，初始时刻的病人数为 i。根据病人不死也不会康复的假设及（竞争项）统计筹算律，其中：解得：（3.17）可得：（

31、3.16）统计结果显示，(3.17)(3.17)预报结果比(3.15)(3.15)更接近实际情况。医学上称曲线 为传染病曲线，并称 最大值时刻t1为此传染病的流行高峰。令：得：此值与传染病的实际高峰期非常接近，可用作医学上的预报公式。模型2 2仍有不足之处，它无法解释医生们发现的现象，且当时间趋与无穷时，模型预测最终所有人都得病，与实际情况不符。为了使模型更精确，有必要再将人群细分，建立多房室系统第43页/共93页infectiverecoveredsusceptiblekl（3.18）l 称为传染病恢复系数 求解过程如下：对（3）式求导，由（1）、（2）得：解得：记：则：将人群划分为三类（见

32、右图）：易感染者、已感染者和已恢复者（recovered）。分别记t时刻的三类人数为s(t)、i(t)和r(t)，则可建立下面的三房室模型：模型3第44页/共93页infectiverecoveredsusceptiblekl 由(1)(1)式可得：从而解得：积分得：（3.19）不难验证，当t+时，r(t)趋向于一个常数，从而可以解释医生们发现的现象。为揭示产生上述现象的原因（3.183.18）中的第（1 1）式改写成：其中 通常是一个与疾病种类有关的较大的常数。下面对 进行讨论，请参见右图如果 ，则有 ，此疾病在该地区根本流行不起来。如果 ，则开始时 ，i(t)单增。但在i(t)增加的同时，

33、伴随地有s(t)单减。当s(t)减少到小于等于 时，i(t)开始减小，直至此疾病在该地区消失。鉴于在本模型中的作用，被医生们称为此疾病在该地区的阀值。的引入解释了为什么此疾病没有波及到该地区的所有人。图3-14 第45页/共93页综上所述，模型3 3指出了传染病的以下特征：（1 1）当人群中有人得了某种传染病时，此疾病并不一定流传，仅当易受感染的人数与超过阀值时，疾病才会流传起来。（2 2）疾病并非因缺少易感染者而停止传播，相反，是因为缺少传播者才停止传播的，否则将导致所有人得病。（3 3）种群不可能因为某种传染病而绝灭。模型检验：医疗机构一般依据r(t)来统计疾病的波及人数，从广义上理解，r

34、(t)为t时刻已就医而被隔离的人数，是康复还是死亡对模型并无影响。及：注意到：可得：（3.20）第46页/共93页 通常情况下，传染病波及的人数占总人数的百分比不会太大，故 一般是小量。利用泰勒公式展开取前三项，有：代入（3.203.20）得近似方程：积分得：其中：这里双曲正切函数 ：而：对r(t)求导：（3.21）第47页/共93页曲线 在医学上被称为疾病传染曲线。图3-14给出了（3.21）式曲线的图形，可用医疗单位每天实际登录数进行比较拟合得最优曲线。图3-14（a）第48页/共93页3.63.6 糖尿病的诊断糖尿病的诊断 糖尿病是一种新陈代谢疾病，它是由胰岛素缺乏引起的新陈代谢紊乱造成

35、的。糖尿病的诊断是通过葡萄糖容量测试（GTTGTT）来检查的，较严重的糖尿病医生不难发现，较为困难的是轻微糖尿病的诊断。轻微糖尿病诊断时的主要困难在于医生们对葡萄糖容许剂量的标准看法不一。例如，美国罗得岛的一位内科医生看了一份GTTGTT测试的报告后认为病人患有糖尿病，而另一位医生则认为此人测试结果应属正常。为进一步诊断，这份检测报告被送到波士顿，当地专家看了报告后则认为此人患有垂体肿瘤。二十世纪6060年代中期，北爱尔兰马由医院的医生RosevearRosevear和MolnarMolnar以及美国明尼苏达大学的AckemanAckeman和GatewoodGatewood博士研究了血糖循环

36、系统，建立了一个简单的数学模型，为轻微糖尿病的诊断提供了较为可靠的依据。模型假设模型假设根据生物、医学等原理，作如下假设：根据生物、医学等原理，作如下假设：(1)(1)葡萄糖是所有细胞和组织的能量来源，在新陈代谢葡萄糖是所有细胞和组织的能量来源，在新陈代谢 中起着十分重要的作用。每个人都有自己最适当的中起着十分重要的作用。每个人都有自己最适当的 血糖浓度，当体内的血糖浓度过渡偏离这一浓度血糖浓度，当体内的血糖浓度过渡偏离这一浓度 时，将导致疾病甚至死亡。时，将导致疾病甚至死亡。(2)(2)血糖浓度是处于一个自我调节系统之中的，它受到血糖浓度是处于一个自我调节系统之中的，它受到 生理激素和其他代

37、谢物的影响和控制，这些代谢物生理激素和其他代谢物的影响和控制，这些代谢物 包括胰岛素、高血糖素、肾上腺素、糖皮质激素、包括胰岛素、高血糖素、肾上腺素、糖皮质激素、生长激素、甲状腺素等，统称为内分泌激素。生长激素、甲状腺素等，统称为内分泌激素。(3)(3)内分泌激素中对血糖起主要影响的是胰岛素，葡萄内分泌激素中对血糖起主要影响的是胰岛素，葡萄 糖只有在胰岛素的作用下才能在细胞内进行大量的糖只有在胰岛素的作用下才能在细胞内进行大量的 生化反应，降低血糖浓度。此外，高血糖素能将体生化反应，降低血糖浓度。此外，高血糖素能将体 内过量的糖转化为糖元储存于肝脏中，从而降低血内过量的糖转化为糖元储存于肝脏中

38、，从而降低血 糖的浓度。糖的浓度。第49页/共93页模型用一、两个参数来区分正常人与轻微病人（测量若干次），根据上述假设，建模时将研究对象集中于两个浓度：葡萄糖浓度和激素浓度。以G表示血糖浓度，以 H表示内分泌激素的浓度。根据上述假设血糖浓度的变化规律依赖于体内现有的血糖浓度及内分泌激素的浓度，记这一依赖关系为函 数F(G,H)。而内分泌激素浓度的变化规律同样依赖于体内现有的血糖 浓度以及内分泌激素的浓度，记其依赖关系为函 数F(G,H)，故有：=(G,H)+J(t)=(G,H)（3.19）其中J(t)为被检测者在开始检测后服下的一定数量的葡萄糖。病人在检测前必须禁食，故可设检测前病人血糖浓度

39、及内分泌激素的浓度均已处于平衡状态 第50页/共93页即可令 t=0时 G=G0,H=H0且 F1(G0,H0 )=0F2(G0,H0 )=0从而有 在测试过程中 G,H 均为变量，而我们关心的却只是它们的改变量，故令g=G G0,h=H H0,在（3.19）中将 展开，得到其中 、是g 和h 的高阶无穷小量。第51页/共93页 很小时（即检测者至多为轻微病人时），为求解方 便，我们考察不包含它们的近似方程组 方程组（3.20 3.20）是一个非线性方程组，较难求解。当 、首先，我们来确定右端各项的符号。从图 中可看出，当J(t)=0 J(t)=0 时,若g g 0 0 且 h=0h=0，则此

40、人血糖浓度高于正常值，内分泌激素将促使组织吸收葡萄糖，并将其存储进肝脏，此时有 0，从而应有：0 第52页/共93页其激素浓度将增加以抑制血糖浓度的增高，因而又有：其激素浓度将增加以抑制血糖浓度的增高，因而又有：:0 0反之，当J(t)=0而g=0且h0 时，此人激素浓度高于正常值，血糖浓度及激素浓度均将减少，从而必有将方程组（3.20）改写成其中 均为正常数。第53页/共93页（3.21）是关于 g、h的一阶常系数微分方程组，因激素浓度不易测得，对前式再次求导化为：由于由于故故 或或 (3.22),令则则(3.22)可简写成可简写成 (3.23)其中其中,第54页/共93页设在设在t=0 时

41、患者开始被测试，他需在很短时间内喝下一定数量时患者开始被测试，他需在很短时间内喝下一定数量 的外加葡萄糖水，如忽略这一小段时间，此后方程可写成的外加葡萄糖水，如忽略这一小段时间，此后方程可写成(3.24)（注：要考虑这一小段时间的影响可利 用Dirac的函数）(3.24)式具有正系数，且式具有正系数，且 当当t趋于无穷时趋于无穷时g趋于趋于0，（体内的葡萄，（体内的葡萄 糖浓度将逐渐趋于平衡值），不难证糖浓度将逐渐趋于平衡值），不难证 明明G将趋于将趋于g（t）的解有三种形式，取决于的解有三种形式，取决于 的符号。的符号。0 0时可得时可得 (1)当其中其中 ，所以，所以 (3.25)(3.2

42、5)式中含有5个参数，即 、A、和，用下述方法可以确定它们的值。在外加葡萄糖水喝入前患者血糖浓度应为 （检查前患者是禁食的）,可先作一次测试将其测得。第55页/共93页进而，取 t=(i=1、2、3、4)各测一次，将测得的值代入（3.25），得到一个方程组，由此可解得相应的参数值。一般，为了使测得的结果更准确，可略多测几次，如 测5-6次，再根据最小平方误差来求参数，即求解 min min 解出所需的参数当 0时可类似加以讨论。实际计算时不难发现，G的微小误差会引 起的很大偏差，故任一包含的诊断标准都将是不可靠的。同时也可发 现G对 并不十分敏感（计算结果与实际值相差较小），故可用 的测试结果

43、作为GTT检测值来判断此人是否真的患有轻微的糖尿病。为了判断上的方便，一般利用所谓自然周 期T作为判别标准根据人们的生活习惯，两餐之间的间隔时间大体 为4小时。临床应用显示，在T 4（小时）时一般表示为正常情况，当T 明显大于4小时时一般表示此人的确患有轻微的糖尿病。由于内分泌激素浓度不易测量，在上面的建模过程中对各种不同激素未加以一一区别，即对其采用了集中参数法。这样做虽大大简化了模型，但也在一定程度上影响了模型的应用效果。临床应用时发现，在患者饮下葡萄糖水大 约3-5小时后，测得的数据有一定的偏差，其原因可能是内分泌激素的作用造成的，因而，要得到更精确的结果，当然要考虑到内分泌激素浓度的变

44、化，建立更精确的模型。罗德岛医院已找到一种测量内分泌浓度的方法，相信在此基础上一定可以设计出诊断轻微糖尿病的更好方法。第56页/共93页3.73.7 稳定性问题稳定性问题 在研究许多实际问题时，人们最为关心的也许并非系统与时间有关的变化状态，而是系统最终的发展趋势。例如，在研究某频危种群时，虽然我们也想了解它当前或今后的数量，但我们更为关心的却是它最终是否会绝灭，用什么办法可以拯救这一种群，使之免于绝种等等问题。要解决这类问题，需要用到微分方程或微分方程组的稳定性理论。在下两节，我们将研究几个与稳定性有关的问题。第57页/共93页一般的微分方程或微分方程组可以写成：定义 称微分方程或微分方程组

45、 为自治系统或动力系统。（3.28）若方程或方程组f f(x x)=0)=0有解X Xo o，X=XX=Xo o显然满足（显然满足（3.28）。）。称点X Xo o为微分方程或微分方程组（3.28)3.28)的平衡点或奇点。第58页/共93页例7 7 本章第2 2节中的LogisticLogistic模型 共有两个平衡点：N=0和N=K，分别对应微分方程的两两个特殊解。前者为No=0时的解而后者为No=K时的解。当NoK时，则位于N=K的上方。从图3-17中不难看出，若No0，积分曲线在N轴上的投影曲线（称为轨线）将趋于K。这说明，平衡点N=0和N=K有着极大的区别。图3-17 定义1 1 自

46、治系统 的相空间是指以（x1,xn）为坐标 的空间Rn。特别，当n=2时，称相空间为相平面。空间Rn的点集(x1,xn)|xi=xi(t)满足(3.28)，i=1,n称为系统的轨线，所有轨线在相空间的分布图称为相图。第59页/共93页定义2 2 设x x0 0是（3.283.28）的平衡点，称：（1 1）x x0 0是稳定的，如果对于任意的00，存在一个00，只要|x x(0)-(0)-x x0 0|，就有|x x(t t)-)-x x0 0|对所有的t t都成立。（2 2）x x0 0是渐近稳定的，如果它是稳定的且 。微分方程平衡点的稳定性除了几何方法，还可以通过解析方法来讨论，所用工具为以

47、下一些定理。（3 3）x x0 0是不稳定的，如果（1 1）不成立。根据这一定义，LogisticLogistic方程的平衡点N=KN=K是稳定的且为渐近稳定的，而平衡点N=0N=0则是不稳定的。第60页/共93页解析方法定理1 1 设x xo o是微分方程 的平衡点：若 ,则x xo o是渐近稳定的若 ,则x xo o是渐近不稳定的证 由泰勒公式，当x x与x xo o充分接近时，有：由于x xo o是平衡点，故f f(x xo o)=0)=0。若 ，则当x x 0)0，从而x x单增；当x x x xo o时，又有f f(x x)0)00，可能出现以下情形：若q0，120。当p0时，零点不

48、稳定；当p p00时，零点稳定 若q q00，1 12 200时，零点不 稳定 当p0时，零点稳定（2）0，零点稳定若a=0，有零点为中心的周期解 综上所述：仅当p0时，（3.30）零点才是渐近稳定的；当p=0且q0时（3.30）有周期解，零点是稳定的中心（非渐近稳定）；在其他情况下，零点均为不稳定的。非线性方程组（3.29）平衡点稳定性讨论可以证明有下面定理成立:定理2 若（3.30）的零点是渐近稳定的，则（3.29）的平衡点 也是渐近稳定的；若（3.30）的零点是不稳定的，则（3.29）的平衡点也是不稳定的。第64页/共93页3.83.8 捕食系统的捕食系统的VolterraVolterr

49、a方程方程 问题背景：意大利生物学家DAncona曾致力于鱼类种群相互制约关系的研究，在研究过程中他无意中发现了一些第一次世界大战期间地中海沿岸港口捕获的几种鱼类占捕获总量百分比的资料，从这些资料中他发现各种软骨掠肉鱼，如鲨鱼、鳐鱼等我们称之为捕食者（或食肉鱼）的一些不是很理想的鱼类占总渔获量的百分比。在 19141923年期间，意大利阜姆港收购的鱼中食肉鱼所占的比例有明显的增加：年代1914191419151915191619161917191719181918百分比11.911.921.421.422.122.121.221.236.436.4年代191919191920192019211

50、9211922192219231923百分比27.327.316.016.015.915.914.814.810.710.7 他知道，捕获的各种鱼的比例近似地反映了地中海里各种鱼类的比例。战争期间捕鱼量大幅下降，但捕获量的下降为什么会导致鲨鱼、鳐鱼等食肉鱼比例的上升，即对捕食者有利而不是对食饵有利呢？他百思不得其解，无法解释这一现象，就去求教当时著名的意大利数学家V.VolterraV.Volterra，希望他能建立一个数学模型研究这一问题。第65页/共93页 Volterra将鱼划分为两类。一类为食用鱼（食饵），数量记为x1(t)，另一类为食肉鱼（捕食者），数量记为x2(t)，并建立双房室系