数学模型的应用

数学模型的应用数学模型的应用数学模型的应用数学模型的应用编制仅供参考审核批准生效日期地址：电话：传真:邮编:数学建模数模作业（第一章）P21第一章6、利用节药物中毒施救模型确定对于孩子(血液容量为2000ml)以及成人(血液容量为4000ml)服用氨茶碱能引起严重中毒和致命的最小剂量。解：设孩子服用氨茶碱能引起严重中毒的最小剂量为，则由节中的药物中毒施救模型可知:在胃肠道中药物的量为，而在血液系统中药物的量为，再令再做出由图可知具有最大值，设在这个最大值在孩子血液中容量的比例为严重中毒的比例以及致命的比例即为孩子服用氨茶碱的最小剂量。于是可以去求这个最小剂量。由上图可知最大值位于左右,利用Mathematics去找出这个最大值。求得，而。于是孩子服用氨茶碱引起严重中毒的最小剂量有式子，从而得此时同理可以求的孩子服用氨茶碱致命的最小剂量为。而成人服用氨茶碱严重中毒与致命的最小剂量分别为。7、对于节的模型，如果采用的是体外血液透析的办法，求解药物中毒施救模型的血液中药量的变化并作图。解：由题可算得：t=0:2:20y=275*exp*t)+*exp*t)plot(t,y,'b:')第二章3、根据节中的流量数据(表2)和(2)式作插值的数值积分，按照连续模型考虑均流池的容量(用到微积分的极值条件)。解：可以将表2中的数据建立散点图以及平均值，如下：h=0:1:23y=[,,,,,,,,,,,,,,,279,,,,,,,,]x1=0::23;t=sum(y)/24;plot(h,y,'-',x1,t)holdonplot(h,y,x1,'b.')另一方面由，经过转化，从而即可转为。又因为要求出均流池的最大容量，就要令即从中求出时间的值，再去求。从书中可知，又有散点图中可知存在两个时间点使得接下来我们来求出这两个时间，不妨在时间段做插值并求出即可求得于是在时刻或者时刻达到最大值，显然不可能在时刻。事实上，在之前均小于所以不可能达到最大值，故只能在达到最大值。利用插值后的数值以及以直代曲的方法来求积分，从而可以利用数学软件MATLAB求得最大值(代码见附录4)为若要考虑的裕量，可按照来设计均流池。数模作业（第二章插值法）P56题目：根据节中的流量数据（表2）和（2）式作插值和数值积分，按照连续模型考虑考虑均流池的容量（用到微积分的极值条件）。时间/h01234567流量/时间/h89101112131415流量/时间/h1617181920212223流量/分析：我们已知的只有数据的散点。通过已学知识，用matlab画图，画出散点所形成的曲线。建立matlab文件e。m文件，输入的代码为：h=0:1:23y=[,,,,,,,,,,,,,,,279,,,,,,,,]x1=0::23;t=sum(y)/24;plot(h,y,'*-',x1,t)holdonplot(h,y,x1,'r+')在matlab工作区间运行结果为：现用插值进行运算：在matlab中建立M文件x=0:23;y=[261279];h=0::23;t=interp1(x,y,h,'spline')%一维插值利用插值后的数值来求积分，从而利用如下MATLAB代码求得最大值为若要考虑的裕量，可按照来设计均流池。在matlab中建立M文件clear;a=;x=0:23;y=[261279];h=0::23;t=interp1(x,y,h,'spline');%一维插值t1=t(2:22209);m=*(sum(t1)')-*;Max=m+a%容量数模作业（第四章）（1）解：根据题意及表格信息，可列出下列关系试。设投资证券A，B,C,D,E的证券的金额分别为，则：（1）整理后得：（2）在LINGO中输入如下命令：model:max=*x1+*x2+*x3+*x4+*x5;x2+x3+x4>=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<=10;End运行后所得结果：Globaloptimalsolutionfound.Objectivevalue:Totalsolveriterations:4VariableValueReducedCostX1X2X3X4X5RowSlackorSurplusDualPrice12345可得：A，C，E分别投资百万元、百万元、百万元，最大税后收益为百万元。由于（1）可知，增加1百万元收益增加百万元。以%借到1百万元资金需要税收百万元，故借钱合算，列出下列模型：（3）在lingo中输入如下命令：model:max=*x1+*x2+*x3+*x4+*x5;x2+x3+x4>=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<=11;end运行结果如下：Globaloptimalsolutionfound.Objectivevalue:Totalsolveriterations:4VariableValueReducedCostX1X2X3X4X5RowSlackorSurplusDualPrice12345因此解得A，C，E分别投资百万元、百万元、百万元，最大税后收益为百万元。若A税前增加到%，则有：model:max=*x1+*x2+*x3+*x4+*x5;x2+x3+x4>=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<=10;end运行得到：Globaloptimalsolutionfound.Objectivevalue:Totalsolveriterations:0VariableValueReducedCostX1X2X3X4X5RowSlackorSurplusDualPrice12345所以不用改变投资方案。若C税前收益减少为%，则有：model:max=*x1+*x2+*x3+*x4+*x5;x2+x3+x4>=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<=10;end运行结果为：Globaloptimalsolutionfound.Objectivevalue:Totalsolveriterations:0VariableValueReducedCostX1X2X3X4X5RowSlackorSurplusDualPrice12345所以要改变投资方案。3、假设储蓄所每天雇佣的全时服务员中从12点到1点中，为午餐时间的有名，从1点到2点为午餐时间的有名。半时服务员中从9点、10点、11点、12点、1点开始工作的分别为：名。列出相应的模型：model:min=100*x1+100*x2+40*y1+40*y2+40*y3+40*y4+40*y5;x1+x2+y1>=4;x1+x2+y1+y2>=3;x1+x2+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<=3;@gin(x1);@gin(x2);@gin(y1);@gin(y2);@gin(y3);@gin(y4);@gin(y5);end运行结果为：Globaloptimalsolutionfound.Objectivevalue:Extendedsolversteps:0Totalsolveriterations:23VariableValueReducedCostX1X2Y1Y2Y3Y4Y5RowSlackorSurplusDualPrice12345678910所以最优解为x1=3,x2=4,y1=0,y2=0,y3=2,y4=0,y5=1，最小用费为820元。4、（1）假设四个季度开始时，公司招聘的保姆数量分别为人，4个季度开始时保姆总数量分别为人，建立模型如下：model:min=s1+s2+s3+s4;65*s1>=6000+5*x1;65*s2>=7500+5*x2;65*s3>=5500+5*x3;65*s4>=9000+5*x4;s1=120+x1;s2=*s1+x2;s3=*s2+x3;s4=*s3+x4;end运行结果如下：Globaloptimalsolutionfound.Objectivevalue:Totalsolveriterations:1VariableValueReducedCostS1S2S3S4X1X2X3X4RowSlackorSurplusDualPrice123456789因此，4个季度开始时，公司新招聘的保姆数量分别为0,15,0,59人。允许解雇的条件下，设4个季度开始时的保姆数量分别为人，4个季度公司新招的保姆数分别为人，相对应的4个季度结束时解雇的保姆数分别为人，重新建立模型如下：model:min=s1+s2+s3+s4;65*s1>=6000+5*x1;65*s2>=7500+5*x2;65*s3>=5500+5*x3;65*s4>=9000+5*x4;s1=120+x1;s2=*s1+x2-y1;s3=*s2+x3-y2;s4=*s3+x4-y3;end运行结果如下：Globaloptimalsolutionfound.Objectivevalue:Totalsolveriterations:0VariableValueReducedCostS1S2S3S4X1X2X3X4Y1Y2Y3RowSlackorSurplusDualPrice123456789因此，四个季度新招聘保姆人数分别为：0、15、0、72人。四个季度结束后解雇保姆人数分别为：0、15、0.假设混合池中甲、乙、丁所占的比例分别为，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+**z1<=0;(3*x1+x2+*y2+*z2<=0;x1+x2+x4=1;end运行结果为：Localoptimalsolutionfound.Objectivevalue:Totalsolveriterations:4VariableValueReducedCostX1X2X4Y1Y2Z1Z2RowSlackorSurplusDualPrice123