数学模型的应用.docVIP

数学建模 数模作业（第一章） P21 第一章 6、利用1.5节药物中毒施救模型确定对于孩子(血液容量为2000ml)以及成人(血液容量为4000ml)服用氨茶碱能引起严重中毒和致命的最小剂量。 解：设孩子服用氨茶碱能引起严重中毒的最小剂量为，则由1.5节中的药物中毒施救模型可知: 在胃肠道中药物的量为，而在血液系统中药物的量为，再令再做出 由图可知具有最大值，设在这个最大值在孩子血液中容量的比例为严重中毒的比例以及致命的比例即为孩子服用氨茶碱的最小剂量。于是可以去求这个最小剂量。由上图可知最大值位于左右, 利用Mathematics去找出这个最大值。求得，而。于是孩子服用氨茶碱引起严重中毒的最小剂量有式子，从而得此时同理可以求的孩子服用氨茶碱致命的最小剂量为。而成人服用氨茶碱严重中毒与致命的最小剂量分别为。 7、对于1.5节的模型，如果采用的是体外血液透析的办法，求解药物中毒施救模型的血液中药量的变化并作图。 解：由题可算得： t=0:2:20 y=275*exp(-0.1386*t)+112.3*exp(-0.6930*t) plot(t,y,b:) 第二章 3、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,246.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.) 另一方面由，经过转化，从而即可转为。又因为要求出均流池的最大容量，就要令即从中求出时间的值，再去求。从书中可知，又有散点图中可知存在两个时间点使得接下来我们来求出这两个时间，不妨在时间段做插值并求出即可求得于是在时刻或者时刻达到最大值，显然不可能在时刻。事实上，在之前均小于所以不可能达到最大值，故只能在达到最大值。利用插值后的数值以及以直代曲的方法来求积分，从而可以利用数学软件MATLAB求得最大值(代码见附录4)为若要考虑的裕量，可按照来设计均流池。 数模作业（第二章插值法） P56 题目：根据2.5节中的流量数据（表2）和（2）式作插值和数值积分，按照连续模型考虑考虑均流池的容量（用到微积分的极值条件）。 时间/h 0 1 2 3 4 5 6 7 流量/(m3.h-1) 150.12 115.56 84.96 66.60 68.04 71.64 82.08 132.84 时间/h 8 9 10 11 12 13 14 15 流量/(m3.h-1) 185.04 226.80 246.60 250.92 261.00 271.44 273.96 279.00 时间/h 16 17 18 19 20 21 22 23 流量/(m3.h-1) 291.60 302.04 320.68 290.52 281.16 248.40 210.24 186.84 分析：我们已知的只有数据的散点。通过已学知识，用matlab画图，画出散点所形成的曲线。 建立matlab文件e。m文件，输入的代码为： 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,246.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,r+) 在matlab工作区间运行结果为： 现用插值进行运算： 在matlab中建立M文件 x = 0:23; y = [150.12 115.56 84.96 66.6 68.04 71.64 82.08 32.84 185.04 226.8 246.6 250.92 261 271.44 273.96 279 291.6 302.04 310.68 290.52 281.16 248.4 210.24 186.84]; h = 0:0.001:23; t = interp1(x, y, h, spline) %一维插值 利用插值后的数值来求积分，从而利用如下MATLAB代码求得最大值为 若要考虑的裕量，可按照来设计均流池。 在matl