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

《药物中毒数模作业matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、课堂程序提供给大家参考：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=(110

2、0*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('

3、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*

4、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(

5、231/1250))/exp((231*t)/1000)+(1650*exp((231*t)/2500))/exp((231*t)/1000)解得:（1）根据模型（8）计算出施救后血液达到最大值的时间。（即算结果=5.26）;symsxx>>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('74358

6、9/(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（1）要使孩子在施救后立即下降，算出排除率多大？（即算结果=0.4886）（2）如果使用体外血液透析的方法，药物排除率可增加到=0.1155*6=0.693，用这个重新求解模型（7）并作图。z=112.2738*exp(-0.693*t)-28

7、5*exp(-0.1386*t);