最新微分方程1PPT课件

《最新微分方程1PPT课件》由会员分享，可在线阅读，更多相关《最新微分方程1PPT课件（197页珍藏版）》请在金锄头文库上搜索。

1、微分方程1建议: 模型要详略得当 在用微分方程解决实际问题的过程中一定要意识到实际问题是十分复杂的，微分方程只能是在一定程度上对问题的一种近似描述，只要结果的误差在一定范围内即可.任何模型都不可能把影响问题的所有因素都反映在微分方程中，或者要求所得结果十分精确.一个好的微分方程模型是在实际问题的精确性和数学处理的可能性之间的一个平衡.初值问题初值问题 的解为的解为 初值问题初值问题 的解为的解为 雪球融化问题雪球融化问题设雪球在融化时体积的变化率与表面积成比例，且融化过程中它始终为球体，该雪球在开始时的半径为6cm ，经过2小时后，其半径缩小为3cm。求雪球的体积随时间变化的关系。解:设t时刻

2、雪球的体积为 ，表面积为 ，球体与表面积的关系为 由题得引入新常数 再利用题中的条件得分离变量积分得方程得通解为再利用条件 确定出常数C和r代入关系式得 t的取值在 之间。 有一段时间,美国原子能委员会(现为核管理委员会)是这样处理浓缩放射性废物的,他们把这些废物装入密封性能很好的圆桶中,然后扔到水深300英尺的海里。 这种做法是否会造成放射性污染,很自然地引起了生态学家及社会各界的关注。原子能委员会一再保证,圆桶非常坚固,决不会破漏,这种做法是绝对安全的。然而一些工程师们却对此表示怀疑,他们认为圆桶在和海底相撞时有可能发生破裂。而原子能委员会有专家们则仍然坚持自己的看法。于是,双方展开了一场

3、笔墨官司。 究竟谁的意见正确呢?看来只能让事实说话了。问题的关键在于圆桶到底能承受多大速度的碰撞,圆桶和海底碰撞时的速度有多大?放射性废物的处理问题放射性废物的处理问题 大量破坏性实验,发现圆桶在40英尺秒的冲撞下会发生破裂,剩下的问题就是计算圆桶沉入300英尺深的海底时,其末速度究竟有多大了。 美国原子能委员会使用的是55加仑的圆桶,装满放射性废物时的圆桶重量为W527.436磅,而在海水中受到的浮力B470.327磅。此外,下沉时圆桶还要受到海水的阻力,阻力Dv,其中C为常数。工程师们做了大量实验,测得C0.08。现在,取一个垂直向下的坐标,并以海平面为坐标原点(0)。于是,根据牛顿第二定

4、律建立圆桶下沉时应满足方程 质量质量加速度加速度= =重力重力- -浮力浮力- -摩擦阻力摩擦阻力 模型及其解oymgBD困难：无法知道下沉到海底的时间积分和代入初始条件得：最后再用数值计算可以得到水深300时的速度大小。 借助数值方法求出v(300)的近似值。计算结果表明, v(300)45.1英尺秒40英尺秒。 工程师们的猜测是正确的,他们打赢了这场官司。现在,美国原子能委员会已改变了他们处理放射性废物的方法,并明确规定禁止将放射性废物抛入海中。 ： 有高为1米的半球形容器, 水从它的底部小孔流出, 小孔横截面积为1平方厘米(如图). 开始时容器内盛满了水, 求水从小孔流出过程中容器里水面

5、的高度h(水面与孔口中心间的距离)随时间t的变化规律.例例 解：解： 由力学知识得,水从孔口流出的流量为流量系数孔口截面面积重力加速度小孔流水问题设在微小的时间间隔水面的高度由h 降至 ,比较(1)和(2)得:即为未知函数的微分方程.可分离变量所求规律为值得进一步探讨的问题: 漏斗型的容器 由于水的张力的原因，每次水都无法全部留尽，总会剩一小部分在容器中。如何才能让水尽可能少的留在容器中？我们知道，水与容器接触的面积越大，留在容器中的水就越多先讨论一下漏斗的模型。 y 容器的位置容器的位置 可否将容器倾斜，使上部的面积大于下部的面积，使水流的速度更快？倾斜角度？容器的运动状态容器的运动状态 容

6、器的运动状态对流水的速度是肯定会造成影响的，考虑极限的状态，如果容器以大于等于当地重力加速度的加速度竖直向下运动，那么，容器里的水就不会流出。容器以不同的方式运动时对水的流出时间有多少影响？有没有一种运动状态能加快水流的速度呢？ 涡流的影响涡流的影响 涡流对水流的速度是有一定影响的。拿一个水桶反复做这样的试验：首先将桶装满水，记录水面的高度，然后拔出塞住孔口的塞子，让水自然从桶破了的孔中流出，测量流出的时间，然后反复从同一高度作相同的试验，最后求出水自然流尽所需时间的平均值；然后从同一高度作相同的试验，不同的是用一根棍子绕同一方向在水中搅动，使其产生涡流，然后重复上面的步骤。最后发现通过两种方

7、法测得的水流尽所需时间的平均值有较大的差距，于是猜想有无涡流或许对水流的速度也是有一定影响的。 疾病的传播疾病的传播 SARS的预测 SARS是21世纪第一个在世界范围内传播的传染病。SARS从2002年11月份开始在我国和世界范围内流行，到2003年6月23日为止，世界卫生组织(WHO)报道的SARS患者已经达到了8459人，其中802人死亡。 中国是SARS流行的重灾区，到2003年6月23日为止的SARS患者为5326人，其中347人死亡。给人民生活和国民经济的发展带来了巨大的影响。 SARS是由一种冠状病毒引起的传染性很强的呼吸道传染病，它主要通过近距离空气飞沫以及接触病人呼吸道分泌物

8、和密切接触进行传播，也可能通过病人飞沫污染物、如通过手、衣物、食物、水或环境等途径传播。 SARS潜伏期一般为2-11天，在潜伏期无感染 。SARS患者的主要症状有：发热（体温38以上）为首发症状，多为高热，并可持续1-2周以上，可伴有寒战或其他症状，包括头痛、全身酸痛和不适、乏力，部分病人在早期也会有轻度的呼吸道症状(如咳嗽、咽痛等）。 SARS患者治愈后不会再被感染。假设：假设： 1）单位时间感染的人数与现有的感染者成比例；单位时间感染的人数与现有的感染者成比例； 2）单位时间内治愈的人数与现有的感染者成比单位时间内治愈的人数与现有的感染者成比 例；例； 3）单位时间内死亡的感染者人数与现

9、有的感染者 成比例； 4）SARS患者治愈恢复后不再被感染； 5）各类人口的自然死亡可以忽略； 6）忽略迁移的影响。 令I(t)是第t天时SARS感染者的数量， b(t)为感染率，d(t)为死亡率，c(t)为治愈率模型根据当时SARS数据是按天公布这个特点,离散化得 I(t+1)=I(t)+ r(t)I(t) 只要知道开始时SARS的感染人数和函数r(t)，就可以利用该模型进行预测。r(t)的估计 r(t)=I(t+1)-I(t)/I(t) 利用实际数据计算，再进行曲线拟合服药问题服药问题 医生给病人开处方时必须注明两点:服药的剂量和服药的时间间隔.超剂量的药品会对身体产生严重不良后果,甚至死

10、亡,而剂量不足,则不能达到治病的目的.已知患者服药后,随时间推移,药品在体内逐渐被吸收,发生生化反应,也就是体内药品的浓度逐渐降低.药品浓度降低的速率与体内当时药品的浓度成正比.当服药量为A、服药间隔为T,试分析体内药的浓度随时间的变化规律.体内药的浓度随时间的变化规律井深的估算井深的估算假如你站在井边且身上带着一只具有秒表功假如你站在井边且身上带着一只具有秒表功 能的计算器，你也许会出于好奇心想用扔下能的计算器，你也许会出于好奇心想用扔下 一块石头听回声的方法来估计井的深度，一块石头听回声的方法来估计井的深度， 假定你能准确地测定时间，你又怎样来推算假定你能准确地测定时间，你又怎样来推算 井

11、的深度呢，请你分析一下这一问题。井的深度呢，请你分析一下这一问题。我有一只具有秒我有一只具有秒 表功能的计算器。表功能的计算器。方法一方法一假定空气阻力不计，可以直接利用自由落体运动的公式假定空气阻力不计，可以直接利用自由落体运动的公式来计算。例如，来计算。例如， 设设t=4=4秒，秒，g=9.81=9.81米米/ /秒秒2 2，则可求得，则可求得h78.578.5米。米。我学过微积分，我可以做我学过微积分，我可以做 得更好，呵呵。得更好，呵呵。 除去地球吸引力外，对石块下落影响最大的当属空气阻力。根据流体力学知识，此时可设空气阻力正比于石块下落的速度，阻力系数K为常数，因而，由牛顿第二定律可

12、得： 令k=K/m，解得 代入初始条件 v(0)=0，得c=g/k，故有 再积分一次，得： 若设k=0.05并仍设 t=4秒，则可求 得h73.6米。 听到回声再按跑表，计算得到的时间中包含了 反应时间 进一步深入考虑进一步深入考虑不妨设平均反应时间 为0.1秒 ，假如仍 设t=4秒，扣除反应时间后应 为3.9秒，代入 式，求得h69.9米。 多测几次，取平均值再一步深入考虑再一步深入考虑代入初始条 件h(0)=0，得到计算水井高度的公式： 将e-kt用泰勒公式展开并 令k 0+ ，即可得出前面不考虑空气阻力时的结果。还应考虑回声传回来所需要的时间。为此，令石块下落 的真正时间 为t1，声音传

13、回来的时间记 为t2，还得解一个方程组： 这一方程组是非线性的，求解不太容易，为了估算井深竟要去解一个非线性主程组似乎不合情理 相对于石块速度，声音速度要快得多，我们可 用方法二先求一次 h，令t2=h/340，校正t，求石块下落时间 t1t-t2将t1代入式再算一次，得出井深的近似值。例如， 若h=69.9米，则 t20.21秒，故 t13.69秒，求得 h62.3米。 速度V运动，方向永远指向P点,求M点的运动在轴上有一点P以常速度a沿着轴追线问题追线问题平面上另有一点M,它以常正向移动;在轨迹.解: :首先我们建立点首先我们建立点M M运动时所满足的微分方程模型运动时所满足的微分方程模型

14、标,根据条件有 图3.1标,根据条件有 （3.1.7）以记点M在时刻t的坐标，以X记图3.1点P在时刻t的横坐标， 表示P点在t=0的横坐（3.1.6）把（3.1.6）代入（3.1.8），并记（3.1.8）得:上式两边关于作为自变量，把求导得（3.1.9）即由（3.1.9）和（3.1.10）得到M的追线方程 于是（3.1.11）变为 又由得:（3.1.10）（3.1.11）令则讨论由图3.1知,在点M未追上点P之前,点P的横坐标总大于点M 的横坐标,即当时所以积分上式得:即由此得 （3.1.12）从 得解,即点M沿Ox轴移动.假设开始追逐是,点 P 和 M 同在一条平行于轴的直线上,并记它们的

15、位置为及显然有,由此得,从而得:由上式得:（3.1.13）（3.1.14）(3.1.13)减去(3.1.14)得:即为了使点M有可能追上点P,我们假设（3.1.15）此时,由(3.1.15)得到追线方程为:由于初始点在追线上,即当时,因此得:从而得追线方程:当时,就得到相遇点的坐标是追上所需的时间是悬链线问题悬链线问题有一绳索悬挂在A和B两点(不一定是在同一水平线),如图3.2所示.设绳索是均匀的,柔软的,仅受绳本身的重量作用,它弯曲如图中的形状,试确定该绳索在平衡状态时的形状.解解: 设C是其最低点,选取坐标系如图中所示,且轴通过C点.A AB BC CO O图图3.23.2考虑绳索在最低点

16、C与点之间的一段,这一段在下面三个力的作用下平衡:(1)在点P的张力T,方向沿着P点的切线方向;(2)在点C的水平张力H;(3)CP段的垂直的重量,记为,设它作用在某一点Q处,不一定是CP的中心,见图3.3由于平衡关系,这些力在轴(水平)方向的代数和为0,在轴(垂直)方向的代数和也必须为0.T TQ QC CH H图图3.33.3现将张力分解为两个分力：水平方向分，垂直方向分力为此时，在力为轴方向向左而向右；在轴方向，向下而向上，按平衡关系有：两式相除，并利用关系式（为点的切线斜率）得：是在最低点处的张力，是常数，但依赖于，将上式两边对微分得其中表示在水平方向上， 每增加单位距离时，段弧所增加

17、的重量设绳索的密度为，则有其中表示从点算起的弧长，我们需要求出，因为（3.1.1）或又由于故从而方程（3.1.16）化为：（3.1.1）记为绳索最低点到坐标原点的距离，则有：（3.1.1）是一个不显含自变量的方程，令则方程（3.1.1）化为分离变量，积分得：即式中把初始条件代入（3.1.1）上式得：（3.1.1）上式两端同时乘以，故（3.1.1）变为得：两式相加得：即积分上式，得：把初始条件代入上式得为简单起见，假设此时，从而得绳索的方程：上式表示的曲线叫做悬链线此时，绳索在最低点与点之间的一段弧长是：Volterra 捕食捕食-被捕食模型被捕食模型 设有捕食种群和食饵种群生活在同一小环境中,

18、建立微分方程组来研究两种群个体数量随时间的变化趋势.设 t 时刻食饵和捕食者的数量或密度分别为假设个体不区分大小, 而且没有个体向环境输入或从环境输出, 当环境中不存在捕食者时, 食饵种群的增长规律用下述Logistic方程来描述上式左端表示被捕食者的相对增长率; 右端的常数称为内禀增长率,为环境的容纳量, 由 (4.1.6) 可以看出,(4.1.6)因此当时,种群规模增长, 时, 种群规模减小.反映了环境能保证食饵个体数量变化时最合适的容量, 把(4.1.6) 改写形式是其出生率 减去死亡率(4.1.7)其中项反映了以下事实: 即在容纳量一定的条件下,的增大, 将使每一个体平均的生活条件降低

19、, 从而影响种群的相对增长率, 因此 或称为密度制约项.由于捕食者的存在, 将使食饵的增长率减少, 设单位总量成正比, 注意到 t 时刻有y(t) 个捕食者, 它们在时间内每个捕食者吃掉的食饵数量与该时刻食饵的单位时间内吃掉食饵的总数量应为为常数,对于捕食种群, 当不存在食饵种群时, 仍用Logistic于是(4.1.7) 变为方程来描述增长规律, 即当存在食饵种群时, 被捕食者吃掉的食饵将转化为能量去生育后代, 设转化系数为则捕食种群的增长规律为其中式中项反映了捕食者仅以食饵 为生.这样我们得到一个Volterra 捕食-食饵系统(4.1.10)池水含盐量问题的推广 现在有3个水池，盛有不含

20、盐的水各 1000 ，从t=0时刻起以每分钟4 /分的速率向第一个水池内注入盐水，同时又以每分钟3 /分的速率从该水池中流出搅拌均匀的盐水。第一个水池中流出的盐水进入第2个水池，同时再以2(1-cost) /分的速率向第2个水池中加入盐水，第2个水池中搅拌均匀的盐水又以2 /分的速率流入第3个水池，同时再以2(1-sint) /分的速率向第3个水池中加入盐水，第3个水池中搅拌均匀的盐水又以1 /分的速率出。已知输入盐水的含盐率为每立方米1千克，用计算机仿真这3个水池里盐水的变化过程，计算各水池中水的体积、含盐量、含盐率的变化。池水含盐量推广问题的示意图池水含盐量推广问题的示意图池水含盐量推广问

21、题的微分方程模型池水含盐量推广问题的微分方程模型 令 分别是3个池子中t时刻的含盐量, 分别是3个池子中t时刻水的体积,则restart;sys1:= diff(w1(t),t)=1, diff(w2(t),t)=1+2*(1-cos(t*Pi/180), diff(w3(t),t)=1+2*(1-sin(t*Pi/180);dsolve(sys1,w1(0)=1000,w2(0)=1000,w3(0)=1000, w1(t),w2(t),w3(t);assign(%);w1(t);plot(w1(t),w2(t),w3(t),t=0.200);dsolve(diff(x1(t),t)=16-

22、x1(t)/w1(t), x1(0)=0,x1(t);assign(%), x1(t);dsolve( diff(x2(t),t)=3*x1(t)/w1(t)+8*(1-cos(t*Pi/180)-2*x2(t)/w2(t),x2(0)=0,x2(t);assign(%),x2(t);dsolve(diff(x3(t),t)=2*x2(t)/w2(t)+8*(1-sin(t*Pi/180)-x3(t)/w3(t),x3(0)=0,x3(t);assign(%);x3(t);plot(x1(t),t=0.200);水的变化水的变化池1-红，池2-黄，池3-绿池水含盐量推广问题的计算流程池水含盐量

23、推广问题的计算流程 推广问题比较复杂，理论分析很困难，数值计算时难度相当，只是计算量有所增加输入初始值，取定步长在每一个循环内 计算各池中水的体积 含盐量和含盐率输出需要的结果数值计算程序数值计算程序:Euler :Euler 折线法折线法(Maple)(Maple)restart: h:=1.; t1:=h:s11:=0.: w11:=1000.: r11:=s11/w11: s21:=0.: w21:=1000.: r21:=s21/w21: s31:=0.: w31:=1000.: r31:=s31/w31: for i from 2 to 200 do ti:=i*h; s1i:=s1

24、i-1+4*h-3*h*r1i-1; w1i:=w1i-1+h; r1i:=s1i/w1i; s2i:=s2i-1+r1i-1*3*h+2*(1-cos(ti *3.14/180)*h-2*h*r2i-1; w2i:=w2i-1+h+2*(1-cos(ti *3.14/180)*h; r2i:=s2i/w2i; s3i:=s3i-1+r2i-1*2*h+2*(1-sin(ti *3.14/180)*h-h*r3i-1; w3i:=w3i-1+h+2*(1-sin(ti*3.14/180)*h; r3i:=s3i/w3i;end do:p1:=seq(tn,s1n,n=1.200): p2:=s

25、eq(tn,s2n,n=1.200):p3:=seq(tn,s3n,n=1.200): q1:=seq(tn,r1n,n=1.200): q2:=seq(tn,r2n,n=1.200): q3:=seq(tn,r3n,n=1.200):plot(p1,p2,p3);plot(q1,q2,q3);含盐率的变化含盐率的变化池1-红，池2-绿，池3-黄微分方程组应用举例微分方程组应用举例 微分方程组在工程技术中的应用是非常广泛的，涉及的领域有机械、电工技术、通讯、医学等许多领域。下面给出几个例子，用以阐述关于微分方程组的实际应用及其建模思想。1 1、 两自由度的振动问题两自由度的振动问题 常系数线性

26、方程组在工程技术与科学研究中有很多应用，不少问题都归结为它的求解问题。例如两个自由度的振动问题。为了消除不需例如两个自由度的振动问题。为了消除不需要的振动，常常在振动系统中设置减振器。要的振动，常常在振动系统中设置减振器。其中的受力情况。是原机械部件的质量；是减振器的质量；和是两个弹簧，它们的弹性系数（或称为刚度）也分别用和是减速器表示；（假定阻力与速度成正比）的阻尼系数；是强迫力；和分别表示和距它们的平衡位置的位移。体和下面来建立这个系统的运动方程。先分别考虑物律，有方程(1) 物体的受力情况假定弹簧和都满Hook定律。当物体位移时，物体同时位移同时位移这时，弹簧变形（拉长或压缩）的长度为。

27、因此，这时弹簧的弹性力是（力的方向与位移方向相反），由牛顿第二定(2) 物体 的受力情况1) 沿位移方向的外力：2) 阻尼力： （方向与速度方向相反）；3) 这时物体 受到两个弹簧的作用：弹簧 的弹性力 ；弹簧 的弹性力 。由牛顿第二定律，有因此，上述运动系统满足微分方程组上述方程组在变换 及 之下，就变成了一个常系数线性非齐次方程组 (4.6.1)为了简单，只考虑无阻尼自由振动的情形，即 。于是，(4.6.1) 变成 (4.6.2) 它的特征方程为 设 ，则 (4.6.3) 可写成解之，得容易看出， ，因此，可令 ，这时因此，方程 (4.6.3) 最后可写成故 (4.6.3) 的特征根全是纯

28、虚根。(4.6.1) 解的第一个函数可写成形如部件的运动频率。这个事实可以用来防止机器或这是两个简谐运动的叠合。每一个简谐运动的角频（即与）均与减振器的参数与有关。因此，调整与可以改变原设备与外力发生共振现象，以及减轻外力的干扰等。 在研究传染病的传播、种群生态、环境的污染、药物在人体内的分布等问题中，经常把所研究的事物看成由有限个部分组成的系统，而每个部分称为一个仓室仓室(compartment)。它具有以下特点：(1)每个仓室有固定的容量，内含每个时刻都均匀分布着的物质（或能量）。(2)各个仓室间以及仓室与外部环境间均可进行物质（或能量）交换，并服从物质（或能量）守恒定律。这样的系统称为仓

29、室系统仓室系统。 下面根据人体内胆固醇的吸收、合成、排泄等机理来建立胆固醇流动仓室模型。2 2 、胆固醇流动的仓室模型、胆固醇流动的仓室模型 人体内血浆中的胆固醇既可从食物中吸取，也可以从体内通过肝脏等器官合成，且血浆中过量的胆固醇通过肝脏、肠等排出体外。把具有吸收和排泄胆固醇功能的器官组成的系统称仓室 ，如血浆、肝脏、肠道和皮肤等都属于仓室 。把与胆固醇有关的其他器官构成的系统称为仓室 ，如动脉壁等属于仓室 ，即在仓室 中，胆固醇既不能从外部吸收，也不能向外部排放，两个仓室都可以合成胆固醇，且两个仓室的胆固醇可以相互交换如图4.4所示图图4.44.4图 4.4 给出了两个仓室中胆固醇的转移率

30、及 ，其中 代表单位时间每克胆固醇从仓室 流动到仓室 的量， 代表单位时间每克胆固醇从仓室 流动到仓室 的量。其他参数意义为： 代表单位时间内每克胆固醇从仓室 排除的量。 代表单位时间通过从食物中吸收而输进仓室 的胆固醇的量。 是代表单位时间通过人体合成而流入仓室 的胆固的量。 (4.6.4)设和分别代表在时刻，仓室和中胆固醇的总量。通过上述分析和假设，可建立下面关于胆固醇的仓室模型为求解方程组 (4.6.4)，将 (4.6.4) 写成向量形式 ，即这里 。方程组 (4.6.4) 对应的齐次方程组的特征方程为其判别式为因此 ，故系数矩阵 有两个负实部特征根的通解为其相应的特征向量为因此，是方程

31、的基解矩阵，又因为是方程的解。故方程直接计算 得这样可得方程组 (4.6.4) 的通解为3 3、人造卫星的轨迹方程、人造卫星的轨迹方程我们知道，人造卫星在最后一段时间运载火箭熄灭之后，即进入它的轨道，轨道的形状因发射角度和发射速度的不同，而分别出现椭圆、抛物线或双曲线。下面来讨论这些问题。 地球与人造卫星是相互吸引的，但因二者的质量相差很大，因此，可假设地球是不动的，又因人造卫星的体积与地球相比是很小的，故可把它看作质点。为简单起见，我们不考虑太阳、月亮和其它星球的作用，并略去空气阻力。在上面的假设下，可把问题归结为： 从地球表面上一点 ，以倾角 ，初速度 射出一质量为 的物体，如图4.5所示

32、，求此物体运动的轨道方程。 图4.5 根据万有引力定律，地球对卫星的引力大小为过发射点和地心的直线作轴，轴与发射方向所成的平面为平面，平面通过地心，取垂直于轴且过地心的直线为轴，取开始发射时间为，经过时间后，卫星位于点，下面建立和所满足的方程。其方向指向地心，其中 是引力系数， ， 是地球质量， ， 是地球与卫星间的距离。如图4.6所示，这个引力在 轴方向上的分力分别为卫星在 轴上所获得的分加速度分别为 和 。由牛顿第二定律，得到卫星的运动方程为： 轴上的分量分别为当时，卫星在地表面以倾角，初速度射出，所以，在时，。卫星的初速度在因此，初始条件为 （4.6.6）下面利用首次积分法来求方程组 (

33、4.6.5) 满足初始条件(4.6.6)的解。将方程(4.6.5)两端消去 后，以 乘第一个方程，以 乘第二个方程。然后相减，得因为故有对上式两边积分，得首次积分为再以 乘方程组 (4.6.5) 中第一个方程，以 乘第二个方程，然后两式相加，得由于及从而得积分得上式另一首次积分于是，原方程组 (4.6.5) 降为较低阶的方程组 将它们代入 (4.6.7) 得作极坐标变换，并求导 两式联立消去 ，得 这就是卫星运动轨道的极坐标参数方程。若将参这里得到一个仅含有一个未知数得一阶微分方程，若由此解出，代入(4.6.9)中第一式，便可确定，由此得到数消去，便得出卫星运动轨道的极坐标方程。利用分离变量求

34、该方程的解得为此，由 (4.6.9) 的第一式求得, 并代入 (4.6.10) 得，则上式化为整理得令知或(4.6.11)这就是所求的卫星运动轨道的极坐标方程，其中有三个任意常数 (或)，它们可由初始条件 (4.6.6) 确定。注意到当时，因此，。且由(4.6.6)及(4.6.8) 把它们代入 (4.6.9) 及 (4.6.11) 得 (4.6.12)我们知道，(4.6.11)式是圆锥曲线的极坐标方程当 时，轨道是圆；当 时，轨道是椭圆；当 时，轨道是抛物线；当 时，轨道是双曲线 下面进一步讨论卫星发射的初速度与卫星轨道形状的关系因为故上故将(4.6.12)式中的代入上式，并整理得(4.6.1

35、3)注意到 (4.6.12) 式及式可化为因此，当 时，得由此得卫星的轨道是一个圆.即把地球半径，质量及引力常数的具体数值代入上式，并计算得即称为第一宇宙速度第一宇宙速度，此时当时，由(4.6.13)得即所求速度是第一宇宙速度的倍，即，称为第二宇宙速度第二宇宙速度，它的轨道是抛物线。是双曲线（一支）。当时，因，由 (4.6.13) 式可知这表明，当初速度小于第二宇宙速度时，卫星轨道是一个椭圆。当时，由(4.6.13)式可得因此，当初速度大于第二宇宙速度时，它的轨道由电路和机械装置组装成的永磁体扩音器模型如图4.7所示。4 4、扩音器振动模型、扩音器振动模型图4.7 一个时变电源电压 驱动着一个

36、音圈能换器，能换器的转化系数为 ，通过能换器使扬声器的振动膜发生振动。由音圈组成的能换器，本质上是一个在永磁场内的自由运动。当变化的电流通过音圈时，音圈在永磁体的磁力和电流产生的磁力相互作用下进行运动。用 代表能换器与扬声器相互作用力， 是能换器电阻， 是能换器的感应系数.压定律，可得的微分方程组质量为的扬声器的振动是一个具有阻尼的弹簧系统振动。阻尼系数为，弹簧的弹性系数为，能换器转化系数与有关变量间的相互关系为这里是音圈两端的电压降；是音圈位移； 代表音圈的速度，由牛顿第二定律及回路电(4.6.14)令，把 (4.6.14) 化为一阶微分方程组(4.6.15)方程组(4.6.15)的向量形式

37、为，这里设 则矩阵 的特征方程为特征根为 ，其相应特征向量分别为因此，可求得齐次线形方程组的基解矩阵为为求(4.6.15)的通解，需求出(4.6.15)的一个特写成解，我们用待定系数法来求此特解。把因此，方程组 (4.6.15) 有特解将 代入方程 得比较上述方程两边和的系数得把 (4.6.16) 的消去得由此可求得和分别为的通解为方程这里 是任意常向量。进一步，可求得原方程中音圈的位移 和驱动能换器的电流 分别为5、人体含铅量问题、人体含铅量问题铅是如何进入体内的铅是如何进入体内的 ? ?消化道消化道 通常是铅吸收的主要途径，儿童由这一途径吸收的铅所占比例约为85。铅在肠道内通过主动转运和被

38、动扩散两种方式由小肠吸收人血液。主动转运占吸收总量的80以上。主动转运依赖铅与肠粘膜上一种转运蛋白质结合，由转运蛋白作载体将铅转运人血液。被动扩散则是铅由肠腔向血液自然扩散。肠腔中铅浓度越高、血液中铅浓度越低，被动扩散的量就越大。呼吸道呼吸道 空气中含铅的灰尘经鼻孔吸人呼吸道后，一部分被鼻毛、气管纤毛和支气管纤毛、细支气管纤毛阻挡，最后以痰液的形式返回口腔，儿童将其咽人消化道，再由消化道吸收入血。另一部分特别微小的铅尘到达肺泡，沉积在肺泡里，然后由吞噬细胞等吸收，进入血液。皮肤皮肤 当我们接触有机铅的时候，铅会被皮肤直接吸收。 铅中毒会引起神经、消化、循环系统紊乱，表现为贫血、腹痛、高血压等。

39、血铅高组的儿童的总智商、操作智商、语言智商分别比低血铅组落后14分、14分和13分，而每升血液中的铅浓度上升100微克，儿童的身高将降低13厘米。 低剂量的铅接触可以对人体的红细胞、肾脏、免疫系统、骨髓和中枢神经的功能产生不良影响，而所有这些影响发生前都可能没有明显的临床症状。铅中毒会影响婴幼儿最初站立、行走和说话的年龄，也可能引起孩子注意力涣散、记忆力减退、理解力降低及学习困难等。 儿童体内的铅水平可分从血、骨、齿、尿、发检测到，儿童血铅水平分为5级。 级：血铅值低于10微克分升，身体处于相对安全状态； 级：血铅值为10微克19微克分升，属于轻度铅中毒。影响造血、神经传导和认知能力，儿童易出

40、现头昏、烦躁、注意力涣散、多动、厌食、腹胀、轻度贫血； 级：血铅值达到20微克44微克分升，为中度铅中毒。引起缺钙、缺锌、缺铁、免疫力低下，运动不协调，视力和听力受损，学习困难、智力下降，生长发育迟缓，贫血、腹绞痛、食癖、反应迟钝等；级：血铅值达到45微克69微克分升，就是重度铅中毒。可出现性格改变、易激怒、攻击性行为、运动失调、贫血、腹绞痛、高血压、心律失常和痴呆等； 级：当血铅值大于70微克分升时，为极重度铅中毒，可导致脏器损害、铅性脑病、瘫痪、昏迷等。 问题和模型问题和模型 铅是一种人体所必须的微量元素，但体内铅含量过多时就会引起铅中毒.铅是一种重金属元素，通过食物、饮料、空气进入人体，

41、经过呼吸和消化系统后进入血液，再经过血液循环慢慢进入人体和骨头中.铅可以经过人体的排泄系统、通出出汗、剪头发、剪指甲排出体外. 我们根据铅在人体内的变化情况将人体分为血液、组织、骨头3个仓室，铅在这三个仓室中的转化关系如图所示.x1(t)表示t时刻血液中的含铅量，x2(t)a02*x3(t).表示t时刻组织中的含铅量，x3(t)表示t时刻骨头中的含铅量.假设在单位时间内从环境经过消化、吸收系统进入血液的铅为L，从血液进入组织、骨头的铅分别为a31*x1(t)和a21*x1(t)，从组织、骨头再进入血液的铅分别为a13*x2(t)、a12*x2(t)，血液和骨头向外界排出的铅分别为a01*x1(

42、t)和示意图示意图 (1)根据前面的假设建立人体血液、组织和骨头中含铅量所满足的微分方程组,并求其通解. (2)取时间单位为天，对一个志愿者在一种环境中的生活情况进行测定得a21=0.011，a12=0.012, a31=0.0039，a13=0.000035,a01=0.021, a02=0.016,L=49.3. 假设该志愿者开始时体内的含铅量为0，求他体内血液、组织和骨头中含铅量随时间变化的关系，画出这些解曲线的图，并求当x+时这些函数的极限. (4)假设该志愿者在此环境中生活了365天后搬到了一个无铅的环境中去(不再有外界的铅进入体内，即L=0)，再讨论他体内血液、组织和骨头中含铅量随

43、时间变化的关系，并画出0t1460时这些含铅量曲线的图形.简化图简化图模型模型模型 x = Ax + b, 数据数据Lead Transfer Coefficients Lead Transfer Coefficients (Rabinowitz, et al.)Units: days-1a a2121 = 0.011 = 0.011a a1212 = 0.012 = 0.012from blood to tissue and backa a3131 = 0.0039 = 0.0039a a1313 = 0.000035 = 0.000035from blood to bone and bac

44、ka a0101 = 0.021 = 0.021a a0202 = 0.016 = 0.016excretion from blood and tissue时间分为两段: 0,365和365,1460#定义常数定义常数a01:= 0.021; a02:=0.016; a12:=0.012; a13:=0.000035; a21:=0.011; a31:=0.0039; L:=49.3;#定义矩阵和向量定义矩阵和向量A := matrix(3,3,-(a01+a21+a31), a12,a13,a21, -(a02+a12), 0, a31, 0, -a13);b :=matrix(3,1,L,

45、 0, 0);#定义方程定义方程equn1:=diff(x1(t),t)=-(a01+a21+a31)*x1(t)+a12*x2(t)+a13*x3(t)+L;equn2:=diff(x2(t),t)=a21*x1(t)-(a02+a12)*x2(t);equn3:=diff(x3(t),t)=a31*x1(t)-a13*x3(t);# #解初始值问题解初始值问题dsolve(equn1,equn2,equn3,x1(0)=0,x2(0)=0,x3(0)=0,x1(t),x2(t),x3(t);结果太复杂结果太复杂, ,利用迭加原理、特征值和特利用迭加原理、特征值和特征向量法征向量法# #非齐

46、次方程的特解非齐次方程的特解with(linalg): with(plots):xe := evalm(-(inverse(A)&*b);#特征值和特征向量特征值和特征向量eigenvals(A);eigenvects(A);lambda:=-.3061796847e-4, -.1980315266e-1,-.4410122938e-1;v1:=.1124436946e-1, .442226656e-2, 10.00746814;v2:=-.5975248926, -.8018660787, .1178839056;v3:=-.8256380196, .5640574390, .7307156

47、328e-1;#特征向量组成矩阵特征向量组成矩阵 P:=augment(v1,v2,v3);# #得到解的表达式得到解的表达式x1:=c1*P1,1*exp(lambda1*t)+c2*P1,2*exp(lambda2*t)+c3*P1,3*exp(lambda3*t)+xe1,1;x2 :=c1*P2,1*exp(lambda1*t)+c2*P2,2*exp(lambda2*t)+c3*P2,3*exp(lambda3*t)+xe2,1;x3 :=c1*P3,1*exp(lambda1*t)+c2*P3,2*exp(lambda2*t)+c3*P3,3*exp(lambda3*t)+xe3,

48、1;#在解的表达式用在解的表达式用t=0t=0代入代入x10:=simplify(subs(t=0,x1);x20:=simplify(subs(t=0,x2);x30:=simplify(subs(t=0,x3);#求解常数求解常数solve(x10=0,x20=0,x30=0,c1,c2,c3);assign(%);#作图作图plot(x1, x2, x3, t=0.365, color=red,blue,green, thickness=2); plot1:=%:#重新显示解重新显示解x1;x2;x3;#求极限求极限limit(x1,t=infinity);limit(x2,t=infi

49、nity);limit(x3,t=infinity);#得到得到365365天后解的表达式天后解的表达式xx1:=cc1*P1,1*exp(lambda1*t)+cc2*P1,2*exp(lambda2*t)+cc3*P1,3*exp(lambda3*t);xx2 :=cc1*P2,1*exp(lambda1*t)+cc2*P2,2*exp(lambda2*t)+cc3*P2,3*exp(lambda3*t);xx3 :=cc1*P3,1*exp(lambda1*t)+cc2*P3,2*exp(lambda2*t)+cc3*P3,3*exp(lambda3*t);# 计算解在计算解在t=365

50、t=365时的值时的值 x1365:=simplify(subs(t=365,x1);x2365:=simplify(subs(t=365,x2);x3365:=simplify(subs(t=365,x3);xx1365:=simplify(subs(t=365,xx1);xx2365:=simplify(subs(t=365,xx2);xx3365:=simplify(subs(t=365,xx3);#求解常数求解常数solve(xx1365=x1365,xx2365=x2365,xx3365=x3365, cc1,cc2,cc3); assign(%);#重新显示解重新显示解xx1;xx

51、2;xx3;#求极限求极限limit(xx1,t=infinity);limit(xx2,t=infinity);limit(xx3,t=infinity);#作图作图plot(xx1, xx2, xx3, t=365.1460, color=red,blue,green, thickness=2); plot2:=%:#将两个图画在一个图上将两个图画在一个图上display(plot1,plot2);6、血液中的酒精浓度的变化简单情形简单情形 一种啤酒的酒精度是4%，一个人单位时间内排除的酒精与当时体内酒精含量成正比，他在晚会上瞬间喝下2瓶啤酒，问他身体中的酒精含量怎样变化？若啤酒是持续一段

52、时间喝的，问他身体中的酒精含量怎样变化？若喝酒是分段脉冲式的，问他身体中的酒精含量怎样变化？瞬间喝下瞬间喝下2 2瓶啤酒瓶啤酒记体内酒精的含量按指数率下降为t时刻体内酒精的含量，由于在一个晚上一个人的身体没有太大的变化，认为酒精的浓度和含量成比例，我们仅考虑酒精含量的变化。 持续一段时间喝下2瓶啤酒 体内酒精含量先上升，再按指数率下降分段脉冲式饮酒类似于服药问题问题的推广问题的推广 血液中的酒精浓度的变化的建模和分析方法与服药问题完全相同，在这里我们将它用在酒后驾车的分析中。饮酒的历史与现状 酒后驾车所带来的交通事故概述不同类型的酒及不同饮酒方式的建模与分析对交通管理中酒精的测定与标准的建议迄

53、今为止，全世界约有共7亿辆乘用车、货车和大客车，3亿辆摩托车，每年大约有50万人死于交通事故，3500万人受伤。美国有专家预测，交通事故将很快成为世界第三大死亡原因，仅次于心脏病和中风（脑血管意外）。 交通事故的惨痛教训交通事故的惨痛教训世界交通安全报告列出的十大“影响事故发生的危险因素”中第二大因素就是使用酒精酒精、药品和毒品。瑞典、芬兰2529 的严重损伤事故与酒精有关，英国2549 的严重损伤事故与酒精有关，加拿大、美国4363 严重损伤事故与酒精有关。法国6个事故死亡者中有1个与酒精有关，德国6个事故死亡者中有1个与酒精有关，英国7个事故死亡者中有1个与酒精有关，荷兰12个事故死亡者中

54、有1个与酒精有关，日本20个事故死亡者中有1个与酒精有关；加拿大43.5%的死亡事故与酒精有关；美国41% 的死亡事故与酒精有关。酒精的危害有多大？酒精的危害有多大？我国交通部副部长冯正霖指出，根据事故统计分析，驾驶员驾驶员的原因造成的道路交通事故死亡人数占总死亡人数的比例高达71.6%71.6%。驾驶员素质不高、操作技能差、安全意识不强，是导致道路交通事故多发的最主要因素。 中国疾病预防控制中心的专家李丹也着重指出了人为的原因中首要的原因就是酒后驾车酒后驾车。 我国的情况又如何呢？我国的情况又如何呢？酒精的影响酒精的影响酒精是一种神经中枢抑制剂。随着酒精浓度的加大，它由皮层开始向下逐渐扩展，

55、至皮质下中枢、小脑，最后逐渐抑制延髓血管运动中枢和呼吸中枢。不少研究结果显示，醉酒驾驶往往与贪婪、冒险、有野心、精力充沛、性格外向、不负责任、反社会倾向等人格特点有关。酒精对驾驶能力的影响酒精对驾驶能力的影响酒精会使驾驶员的感觉能力降低酒精会使驾驶员的视觉敏锐度下降酒精会使驾驶员的视觉神经系统短暂麻痹酒精会使驾驶员的触觉能力降低酒精会使驾驶员的运动反射神经迟钝血液酒精浓度血液酒精浓度BACBAC饮酒会影响驾驶员的正常的驾驶能力以及良好的驾驶心态，从而增加事故发生的可能性。饮酒后驾驶能力的损害程度虽然存在种族、性别、个体等的差异。但有一点是一致的：总体上血液酒精浓总体上血液酒精浓度（度（bloo

56、d alcohol concentrationblood alcohol concentration，BACBAC）越高，驾驶行为损害越严重）越高，驾驶行为损害越严重。 呼气中酒精呼气中酒精浓浓度度（血液中酒精（血液中酒精浓浓度）度）行行为为表表现现或状或状态态肇事率肇事率0.10mg/L(0.02%;20mg/dL)注意力分散、注意力分散、视觉视觉功能功能损损害害0.15mg/L(0.03%;30mg/dL)驾驶驾驶能力开始有能力开始有损损害害0.25mg/L(0.05%；50mg/dL)复复杂杂技巧障碍、技巧障碍、驾驶驾驶能力明能力明显损显损害害2倍倍0.40mg/L(0.08%；80mg

57、/dL)多多话话、感、感觉觉障碍障碍6倍倍0.50mg/L(0.10%；100mg/dL)说话说话含糊、脚步不含糊、脚步不稳稳7倍倍0.55mg/L(0.11%；110mg/dL)平衡感与判断力障碍度升高平衡感与判断力障碍度升高10倍倍0.75mg/L(0.15%；150mg/dL)明明显显酒醉、步履酒醉、步履蹒蹒跚跚25倍倍0.85mg/L(0.17%；170mg/dL)恶恶心、步履心、步履蹒蹒跚跚50倍倍1.50mg/L(0.30%；300mg/dL)呆滞木僵、可能昏迷呆滞木僵、可能昏迷迷醉迷醉2.00mg/L(0.40%；400mg/dL)呼吸中枢麻痹、呼吸中枢麻痹、渐渐近死亡近死亡无法

58、开无法开车车2.50mg/L(0.50%；500mg/dL)致死致死身体酒精浓度与肇事率（行为表现）之关系身体酒精浓度与肇事率（行为表现）之关系 酒精的吸收与代谢酒精的吸收与代谢1口腔口腔: 酒精被摄取。2胃胃:酒精到达胃部。少量的酒精通过胃壁进入血 液，但是大部分的酒精会下到达小肠。3小肠小肠: 酒精从胃到达小肠。然后绝大部分的酒精 通过小肠壁被吸收进血液。 4血液血液:然后血液会携带酒精到身体的各个部位， 比如说脑部，心脏和肝脏。 5肝脏肝脏:当血液携带着酒精在体内循环时，也会携 带酒精通过肝脏。肝脏会将酒精转化为水，二氧 化碳和能量。这个过程被称为氧化。肝脏氧化酒 精的速率很低，这意味着

59、直到肝脏有时间可以氧 化所有的酒精，酒精会一直持续通过身体的各个 部位，包括脑部。6脑部脑部:酒精几乎一经吸收就会马上到达脑部。血 液携带酒精到达。酒精会一直持续的通过脑部直 到肝脏有时间可以排除所有的酒精。 影响吸收代谢的因素影响吸收代谢的因素酒精的吸收速度与吸收部位、酒的浓度、种类、摄食状态以及性别有关。酒精代谢受各种因素的影响，每个人的代谢速度都不尽相同，男女差别、种族差别使得个体酒精代谢速度有一定的差别。 新英格兰医学杂志上的一份研究报告指出，人喝进体内的酒，一部分在胃里被胃酶分解，不会全进入血液。不同饮酒方式下血液中酒精浓度变化过程建模与分析血液酒精浓度随时间变化的模型时间的起点就为

60、酒精进入口腔的瞬时饮酒方式： 单次饮酒单次饮酒、连续饮酒连续饮酒、间歇饮酒间歇饮酒模型的假设模型的假设饮酒后，不考虑摄入的酒量对人体体液的改变；假设整个过程中人没有摄入任何影响代谢的药类物质和作剧烈性运动；酒精在体内消化、吸收过程同时进行；酒精在血液中的含量与在体液中的含量是相等的；假设人体体液对酒精的吸收速率与胃肠中的酒精含量成正比；假设胃肠中酒精的消除速率与胃肠中的酒精含量成正比；假设肝脏消除酒精的速率与血液中的酒精含量成正比；酒精一经吸收后在人体内即达到平衡，分布均衡；假设单次饮酒方式下酒是在瞬时喝下去并立即到达胃肠中；假设连续饮酒方式下酒是匀速进入人体内的；假设间歇饮酒方式下每隔T0时

61、间饮酒一次，饮酒时酒是瞬时喝下去的； 一个人血液中的酒精含量取决于他的饮酒量、体内原有的酒精含量和饮酒的方式。酒精在人体中的吸收代谢过程可用下图简单表示：模模型型建建立立血液酒精浓度u连续饮酒方式连续饮酒方式u单次饮酒方式单次饮酒方式u间歇饮酒方式间歇饮酒方式酒精浓度酒精浓度参考数据第一组数据：（摘自第一组数据：（摘自20042004年全国大学生数学建模竞赛年全国大学生数学建模竞赛C C题）题）人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是一样的。体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含

62、量（毫克百毫升），得到数据如下表：时间(小时)0.250.50.7511.522.533.544.55酒精含量306875828277686858515041时间(小时)678910111213141516酒精含量3835282518151210774参数估计1)单次饮酒方式单次饮酒方式模型求解1)单次饮酒方式单次饮酒方式体重为70kg的体内不含酒精的人单次摄入N瓶啤酒的血液酒精浓度变化函数为： 而体重为Bkg的人单次摄入Mmg的酒精后血液酒精浓度变化的函数为： 若将体液的比例不固定，则单次摄入Mmg的酒精后血液酒精浓度变化的函数为： 2)连续饮酒方式连续饮酒方式3)间歇饮酒方式间歇饮酒方式参

63、考数据第二组数据：（从一饮酒实验中所得到的一组数据）第二组数据：（从一饮酒实验中所得到的一组数据） 该实验的情况是：实验对象实验对象：16名健康男性自愿者（平均年龄24.5岁，平均体重 59.7kg），无习惯性饮酒史，受试前三天内禁酒。饮酒剂量与方式饮酒剂量与方式：然后将16人随机分为两组，分别按含量酒精 0.65g/kg和0.80g/kg体重给予酒精饮酒，模拟普通社会饮 酒习惯，辅以较高蛋白和脂肪饮食，但限其在30分钟内结 束饮酒。所得0.65g/kg组的实验数据（以开始喝酒起计时，取平均值）如下表所示：时间（小时）0.751.01.522.54.56.58.510.512.5BAC(mg/

64、dL)53.1780.3585.2079.7374.0451.6229.7915.259.295.05模型预测与对照2)连续饮酒方式连续饮酒方式三种饮酒方式的比较假设人体内没有酒精，然后以饮一瓶啤酒为例，比较体重为70千克的人在三种饮酒方式下饮酒后血液酒精浓度变化的情况 :三种饮酒方式下血液酒精浓度变化的曲线图三种饮酒方式下血液酒精浓度变化的曲线图 单次连续间歇考虑最大代谢速率的情况考虑最大代谢速率的情况代谢速率酒精含量代谢速率与酒精含量的关系最大代谢速率函数的图形 系数的调整BAC低于何值时，酒精在体内的消除方式开始转化，文献报道不一。有的认为BAC4mg/dL，有的认为10mg/dL，有的

65、认为20mg/dL，有的实验显示为16.77mg/dL。初始模型（考虑最大代谢速率）初始模型（考虑最大代谢速率）调整后的模型此模型的进一步说明u优点：优点： 这个模型很好的体现了在浓度高时血液中酒这个模型很好的体现了在浓度高时血液中酒精浓度的减精浓度的减 少呈直线下降，而降到一定浓度少呈直线下降，而降到一定浓度（约（约20mg/100ml20mg/100ml）时，下降速率越来越小，跟实）时，下降速率越来越小，跟实际中酒精的代谢机理十分吻合。际中酒精的代谢机理十分吻合。 u缺点缺点: : 在用在用MathematicaMathematica软件求解时只能给出数值软件求解时只能给出数值解，而不能像

66、上一组模型可以求出明确的解析解，解，而不能像上一组模型可以求出明确的解析解，因此在模型的运用上还存在一些困难，有待进一因此在模型的运用上还存在一些困难，有待进一步的研究和推敲。步的研究和推敲。大多数BAC预测的原理 大多数BAC 预测模型是基于Widmark在 1932年提出的一个数学方程式。通过假定消除速率（为平均消除速率）和多长时间完全吸收之后，来估计血液酒精浓度BAC的值 。体液的计算 为了根据摄入的酒精量给出对BAC的合理准确的估计就必须计算总的体液。男性：TBW (升) = 2.447 - 0.09516 年龄 (年) + 0.1074 身高(厘米) + 0.3362 体重 (千克)

67、女性： TBW (升) = -2.097 + 0.1069 身高(厘米) + 0.2466 体重 (千克)国外关于BAC的估计所建立的模型同国外所建立的模型同国外BACBAC估计模型比较说明估计模型比较说明二者存在两大区别 ：u在体液的计算上在体液的计算上，国外存在两种方法，一个是给定Widmark常数（即体液所占比率），另一种是用体重、身高、年龄和性别来确定体液的方法，而我所建立的模型是以一固定的体液所占体重的比率来求体液的；u在酒精消除的问题上在酒精消除的问题上，国外许多方法为了简单起见，都是基于固定的消除速率，即假定一固定的线性消除的速率，我所建立的模型是依据药物动力学的原理，根据消除速

68、率与酒精含量成正比的关系所得到的。 模型的改进模型的改进如果知道更详细的饮酒者的个人数据，可以在我们的模型中引入国外计算体液的方程式从而得到更准确的BAC估计值。测定方法测定方法 可由血液、尿液、唾液与呼气进行测量，但较可行的方法是以呼气酒精浓度（BrAC ）换算血液中酒精浓度（BAC）来检测。两者的换算比率两者的换算比率 我国的标准：呼气和血液酒精含量换算比例为1：2000。交通管理中血液酒精浓度的规定交通管理中血液酒精浓度的规定 国家标准规定车辆驾驶人员每百毫升血液中的酒精含量大于或者等于20毫克、小于80毫克为饮酒驾车，每百毫升血液中的酒精含量大于或者等于80毫克为醉酒驾车。交通管理中酒

69、精的规定根据规定应用模型根据规定应用模型根据前面得到的饮酒模型，结合我国对车辆驾驶人员的血液酒精浓度的最低规定来给出一系列针对司机朋友的饮酒数量的建议针对我国对酒精标准的规定从模型得到针对我国对酒精标准的规定从模型得到的一系列建议的一系列建议 体重为70kg的司机不同方式饮入不同数量的啤酒后多长时间后开车不会违规的建议 :啤酒瓶数啤酒瓶数（瓶）（瓶）0.511.522.533.544.55单单次次1.535.788.099.7311.0012.0412.9213.6814.3514.95连续连续1.746.038.359.9811.2612.3013.1713.9314.6115.21间间歇歇

70、1.746.038.359.9911.2612.3013.1813.9414.6115.21所所需需时时间间（小小时时）（其中连续方式时间为半小时，间歇方式时间为半小时且次数为三次）选取40度的白酒为代表 ,体重为70kg的司机不同方式饮入不同数量的白酒后多长时间后开车不会违规的建议 :白酒数量白酒数量（毫升）（毫升）50100125150175200225250275300单单次次4.378.339.6010.6411.5212.2812.9513.5514.0914.59连续连续4.638.589.8510.8911.7712.5313.2013.8014.3414.84间间歇歇4.638.589.5810.8911.7712.5313.2013.8014.3414.84所所需需时时间间（小小时时）（其中连续方式时间为半小时，间歇方式时间为半小时且次数为三次）本周作业本周作业 利用微分方程解决一个实际问题,写一个完整的报告.问题和背景知识分析、假设、模型模型的求解参数的确定和模型应用讨论和推广