数学模型的应用.docxVIP

数学建模 数模作业（第一章） P21 第一章 6、利用 1.5 节药物中毒施救模型确定对于孩子(血液容量为 2000ml)以及成人(血液容量为 4000ml)服用氨茶碱能引起严重中毒和致命的最小剂量。 解：设孩子服用氨茶碱能引起严重中毒的最小剂量为A ，则由 1.5 节中的药物中毒施救模 1 型可知: 在 胃 肠 道 中 药 物 的 量 为 x(t) ? Ae?0 . 1 t3 ，8 1  而 在 血 液 系 统 中 药 物 的 量 为 y(t) ? 6A (e?0.1155t ?e?0.1386t ) 1 y(t)的图像如下： ， 再 令 y(t )? y t ( )A 1 /? 6e?0 . 1t (?1 e5? 5  t 再0 做)出 0.06 0.06 0.05 0.04 0.03 0.02 0.01 5 10 15 20 由图可知 y(t) 具有最大值，设在这个最大值 y(t) 在孩子血液中容量的比例为严重中 max 毒的比例100? g / ml 以及致命的比例200? g / ml 即为孩子服用氨茶碱的最小剂量。于是可 以去求这个最小剂量。由上图可知最大值位于t ? 8h 左右, 利用 Mathematics 去找出这个最大值。求得 y(t) =0.0669 ，而t ? 7.892h 。于是孩子服用氨茶碱引起严重中毒的最小剂 max 量 A 有式子6 A y(t) / 2000ml ? 100? g / ml ，从而得此时 A 1 1 max 1 ? 498256.1? g 同理可以求 的孩子服用氨茶碱致命的最小剂量为 996512.2? g 。而成人服用氨茶碱严重中毒与致命的 最小剂量分别为996512.2?g、1993024.4?g 。 7、对于 1.5 节的模型，如果采用的是体外血液透析的办法，求解药物中毒施救模型的血液中药量的变化并作图。 解：由题可算得： t=0:2:20 y=275*exp(-0.1386*t)+112.3*exp(-0.6930*t) plot(t,y, b:) 400 350 300 250 200 150 100 50 0 0 2 4 6 8 10 12 14 16 18 20 第二章 3、根据 2.5 节中的流量数据(表 2)和(2)式作插值的数值积分，按照连续模型考虑均流池的容量(用到微积分的极值条件)。 解：可以将表 2 中的数据建立散点图以及平均值，如下： h=0:1:23 y=[150.12,115.56,84.96,66.60,68.04,71.64,82.08,132.84,185.04,226.80,2 46.60,250.92,261.00,271.44,273.96,279,291.60,302.04,310.68,290.52,281 .16,248.40,210.24,186.84] x1=0:0.01:23; t=sum(y)/24; plot(h,y,-,x1,t) hold on plot(h,y,x1,b.) 350 300 250 200 150 100 50 0 5 10 15 20 25 另一方面由c(t ?1) ? c(t) ? f (t) ? g ，经过转化c?(t) ? lim c(t ?1) ? c(t) ? f (t) ? g ，从而 1 即可转为 c(t) ? ?t t f (x)dx ? g(t ? t 0 ), f (x)是f (t)的插值函数，t 是某个初始时刻。又 0 0 因为要求出均流池的最大容量 c(t)  max ，就要令 c?(t)=0，即 f (t) ? g. 从中求出时间t 的值， 再去求 c(t)  max 。从书中可知 g ? 1 ?23 f (t) ? 203.67m3 / h ，又有散点图中可知存在两个 t24 t 时间点t 1 ?(8,9), t 2 ?(22，23) 使得 f (t?1) ? g. 接下来我们来求出这两个时间t , t 1 2 ，不妨在时 间段(8，9)、(22，23)做插值并求出 t , t 1 2 即可求得 t 1 ? 8.45h, t 2 ? 22.208h. 于是 c(t) 在t 时刻 1 或者t 2 时刻达到最大值，显然不可能在t 1 时刻。事实上，在t 1 之前 f (t) 均小于 g 所以c(t) 不 可能达到最大值，故只能在t 2 ? 22.208h 达到最大值。利用插值后的数值以及以直代曲的方 法来求积分 ?t f (x)dx ，从而可以利用数学软件 MATLAB 求得最大值 (代码见附录 4)为 t 0 917.08m3.若要考虑25% 的裕量，可按照1146.4m3 来设