MATLAB药丸的计数与识别,药物中毒数模作业matlab程序.doc

药物中毒数模作业matlab程序.doc

课堂程序提供给大家参考：clear%%%%%%求方程12symslamdatxmuyx=dsolve( Dx=-lamda*x , x(0)=1100 , t )y=dsolve( Dy=lamda*1100*exp(-lamda*t)-mu*y , y(0)=0 , t )>>x=dsolve( Dx=-lamda*x , x(0)=1100 , t )x=1100/exp(lamda*t)>>y=dsolve( Dy=lamda*1100*exp(-lamda*t)-mu*y , y(0)=0 , t )y=(1100*lamda)/(exp(mu*t)*(lamda-mu))-(1100*lamda)/(exp(lamda*t)*(lamda-mu))%作图1t=0:0.1:25;y=6600*(exp(-0.1155*t)-exp(-0.1386*t));x=1100*exp(-0.1386*t);plot(t,x,t,y)gridon%计算药量达到400mg的时间t1,即求出t1=4.87tt=fzero( 6600*(exp(-0.1155*x)-exp(-0.1386*x))-400 ,[0,8])>>tt=fzero( 6600*(exp(-0.1155*x)-exp(-0.1386*x))-400 ,[0,8])tt=4.8659%%%%%%%以下程序是求出P12中t2=7.89symsxxy=6600*(exp(-0.1155*xx)-exp(-0.1386*xx))-400;z=diff(y,xx)z=22869/(25*exp((693*xx)/5000))-7623/(10*exp((231*xx)/2000))%%%算出来式子-7623/10*exp(-231/2000*xx)+22869/25*exp(-693/5000*xx)=0tmax=solve( -7623/10*exp(-231/2000*xx)+22869/25*exp(-693/5000*xx)=0 , xx )tmax=-(10000*log(5/6))/231%%%%求出值t=-10000/231*log(5/6)=7.8927t=-10000/231*log(5/6)作业：对于该案例，结合教材P12，给出具体的施救方案，解决如下问题：(1)作图，即求出口服活性炭药物后模型(8)后，参考教材，画出图形2；z=((473*exp(231/500))/2-1650*exp(231/1250))/exp((231*t)/1000)+(1650*exp((231*t)/2500))/exp((231*t)/1000)解得:0.13862.310()5*95*,ttztee(2)根据模型(8)计算出施救后血液达到最大值的时间。(即算结果=5.26);3tsymsxx>>z=1650*exp(-0.1386*xx)-1609.5*exp(-0.2310*xx);a=diff(z,xx)a=743589/(2000*exp((231*xx)/1000))-22869/(100*exp((693*xx)/5000))>>tmax=solve( 743589/(2000*exp((231*xx)/1000))-22869/(100*exp((693*xx)/5000))=0 , xx )tmax=(2500*log(1073/660))/231>>(2500*log(1073/660))/231ans=5.2595T3=5.2595(3)要使孩子在施救后立即下降，算出排除率多大？(即算结果=0.4886)()zt(4)如果使用体外血液透析的方法，药物排除率可增加到=0.1155*6=0.693，用这个重新求解模型(7)并作图。z=112.2738*exp(-0.693*t)-285*exp(-0.1386*t);