药物中毒数模作业matlab程序

课堂程序提供给大家参考：clear% 求方程1 2syms lamda t x mu yx=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)grid on% 计算药量达到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.89syms xx y=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)解得:（2） 根据模型（8）计算出施救后血液达到最大值的时间。（即算结果=5.26）;syms xx 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）（4） 如果使用体外血液透析的方法，药物排除率可增加到=0.1155*6=0.693，用这个重新求解模型（7）并作图。z=112.2738*exp(-0.693*t)-285*exp(-0.1386*t);