1、數學建模數模作業（第一章）P21第一章6、利用1.5節藥物中毒施救模型確定對于孩子(血液容量為2000ml)以及成人(血液容量為4000ml)服用氨茶堿能引起嚴重中毒和致命的最小劑量。解：設孩子服用氨茶堿能引起嚴重中毒的最小劑量為，則由1.5節中的藥物中毒施救模型可知:在胃腸道中藥物的量為，而在血液系統中藥物的量為，再令再做出由圖可知具有最大值，設在這個最大值在孩子血液中容量的比例為嚴重中毒的比例以及致命的比例即為孩子服用氨茶堿的最小劑量。于是可以去求這個最小劑量。由上圖可知最大值位于左右, 利用Mathematics去找出這個最大值。求得，而。于是孩子服用氨茶堿引起嚴重中毒的最小劑量有式子，

2、從而得此時同理可以求的孩子服用氨茶堿致命的最小劑量為。而成人服用氨茶堿嚴重中毒與致命的最小劑量分別為。7、對于1.5節的模型，如果采用的是體外血液透析的辦法，求解藥物中毒施救模型的血液中藥量的變化并作圖。解：由題可算得：t=0:2:20y=275*exp(-0.1386*t)+112.3*exp(-0.6930*t) plot(t,y,b:)第二章3、根據2.5節中的流量數據(表2)和(2)式作插值的數值積分，按照連續模型考慮均流池的容量(用到微積分的極值條件)。解：可以將表2中的數據建立散點圖以及平均值，如下： h=0:1:23y=150.12,115.56,84.96,66.60,68.0

3、4,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.84x1=0:0.01:23;t=sum(y)/24;plot(h,y,-,x1,t)hold onplot(h,y,x1,b.) 另一方面由，經過轉化，從而即可轉為。又因為要求出均流池的最大容量，就要令即從中求出時間的值，再去求。從書中可知，又有散點圖中可知存在兩個時間點使得接下來我們來求出這兩個時間，不妨在時間段做插值并求出即可求得于是在時

4、刻或者時刻達到最大值，顯然不可能在時刻。事實上，在之前均小于所以不可能達到最大值，故只能在達到最大值。利用插值后的數值以及以直代曲的方法來求積分，從而可以利用數學軟件MATLAB求得最大值(代碼見附錄4)為若要考慮的裕量，可按照來設計均流池。數模作業（第二章插值法）P563、 題目：根據2.5節中的流量數據（表2）和（2）式作插值和數值積分，按照連續模型考慮考慮均流池的容量（用到微積分的極值條件）。時間/h01234567流量/(m3.h-1)150.12115.5684.9666.6068.0471.6482.08132.84時間/h89101112131415流量/(m3.h-1)185.

5、04226.80246.60250.92261.00271.44273.96279.00時間/h1617181920212223流量/(m3.h-1)291.60302.04320.68290.52281.16248.40210.24186.84分析：我們已知的只有數據的散點。通過已學知識，用matlab畫圖，畫出散點所形成的曲線。建立matlab文件e。m文件，輸入的代碼為：h=0:1:23y=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,

6、279,291.60,302.04,310.68,290.52,281.16,248.40,210.24,186.84x1=0:0.01:23;t=sum(y)/24;plot(h,y,*-,x1,t)hold onplot(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 2

7、81.16 248.4 210.24 186.84;h = 0:0.001:23;t = interp1(x, y, h, spline) %一維插值利用插值后的數值來求積分，從而利用如下MATLAB代碼求得最大值為 若要考慮的裕量，可按照來設計均流池。在matlab中建立M文件clear; a = 876.15;x = 0:23;y = 150.12 115.56 84.96 66.6 68.04 71.64 82.08 132.84 185.04 226.8 246.6 250.92 261 271.44 273.96 279 291.6 302.04 310.68 290.52 281.

8、16 248.4 210.24 186.84;h = 0:0.001:23;t = interp1(x, y, h, spline); %一維插值t1 = t(2:22209);m = 0.001 * (sum(t1) - 203.67 * 22.208;Max = m + a %容量數模作業（第四章）1、 （1）解：根據題意及表格信息，可列出下列關系試。設投資證券A，B,C,D,E的證券的金額分別為，則： （1）整理后得： （2）在LINGO中輸入如下命令：model:max=0.043*x1+0.027*x2+0.025*x3+0.022*x4+0.045*x5;x2+x3+x4=4;6*

9、x1+6*x2-4*x3-4*x4+36*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;x1+x2+x3+x4+x5=4;6*x1+6*x2-4*x3-4*x4+36*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;x1+x2+x3+x4+x5=4;6*x1+6*x2-4*x3-4*x4+36*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;x1+x2+x3+x4+x5=4;6*x1+6*x2-4*x3-4*x4+36*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;x1+x2+x3+x4+x5=4;x1+x2+y1+y2=3;x1+x

10、2+y1+y2+y3=4;x2+y1+y2+y3+y4=6;x1+y2+y3+y4+y5=5;x1+x2+y3+y4+y5=6;x1+x2+y4+y5=8;x1+x2+y5=8;y1+y2+y3+y4+y5=6000+5*x1;65*s2=7500+5*x2;65*s3=5500+5*x3;65*s4=9000+5*x4;s1=120+x1;s2=0.85*s1+x2;s3=0.85*s2+x3;s4=0.85*s3+x4;end運行結果如下：Global optimal solution found. Objective value: 478.5107 Total solver iterat

11、ions: 1 Variable Value Reduced Cost S1 120.0000 0.000000 S2 116.5000 0.000000 S3 99.02500 0.000000 S4 142.9857 0.000000 X1 0.000000 0.8732231 X2 14.50000 0.000000 X3 0.000000 0.9291667 X4 58.81448 0.000000 Row Slack or Surplus Dual Price 1 478.5107 -1.000000 2 1800.000 0.000000 3 0.000000 -0.2982986

12、E-01 4 936.6250 0.000000 5 0.000000 -0.1666667E-01 6 0.000000 -0.8732231 7 0.000000 0.1491493 8 0.000000 -0.9291667 9 0.000000 0.8333333E-01因此，4個季度開始時，公司新招聘的保姆數量分別為0,15,0,59人。（2） 允許解雇的條件下，設4個季度開始時的保姆數量分別為人，4個季度公司新招的保姆數分別為人，相對應的4個季度結束時解雇的保姆數分別為人，重新建立模型如下：model:min=s1+s2+s3+s4;65*s1=6000+5*x1;65*s2=75

13、00+5*x2;65*s3=5500+5*x3;65*s4=9000+5*x4;s1=120+x1;s2=0.85*s1+x2-y1;s3=0.85*s2+x3-y2;s4=0.85*s3+x4-y3;end運行結果如下：Global optimal solution found. Objective value: 465.1218 Total solver iterations: 0 Variable Value Reduced Cost S1 120.0000 0.000000 S2 116.5000 0.000000 S3 84.61538 0.000000 S4 144.0064 0.

14、000000 X1 0.000000 0.9291667 X2 14.50000 0.000000 X3 0.000000 0.7147436E-01 X4 72.08333 0.000000 Y1 0.000000 0.8333333E-01 Y2 14.40962 0.000000 Y3 0.000000 0.8333333E-01 Row Slack or Surplus Dual Price 1 465.1218 -1.000000 2 1800.000 0.000000 3 0.000000 -0.1666667E-01 4 0.000000 -0.1429487E-01 5 0.0

15、00000 -0.1666667E-01 6 0.000000 -0.9291667 7 0.000000 0.8333333E-01 8 0.000000 0.000000 9 0.000000 0.8333333E-01因此，四個季度新招聘保姆人數分別為：0、15、0、72人。四個季度結束后解雇保姆人數分別為：0、15、0.6、 假設混合池中甲、乙、丁所占的比例分別為，A類產品中來自混合池和原料丙的噸數分別為噸，B類產品來自混合池和原料丙的噸數分別為噸。優化目標總利潤最大。建立模型如下：model:max=(9-6*x1-16*x2-15*x4)*y1+(15-16*x1-16*12-15*x4)*y2+(9-10)*z1+(15-10)*z2;x4*(y1+y2)=50;y1+y2=100;y2+z2=200;(3*x1+x2+x4-2.5)*y1-0.5*z1=0;(3*x1+x2+x4-1.5)*y2+0.5*z2 L=0 4 3;1/2 0 0;0 1/4 0; A=L3A = 0.37500000000000 8.00000000000000 6.00000000000000 1.00000000000000 0.37500000000000 0 0 0.50000000000000 0.37500000000000 B=1000 1000 1