药物中毒数模作业matlab程序VIP专享

m^#a^^^ikmaf-ab^^ 课堂程序提供给大家参考：clear %%%%%% 求方程 1 2 syms lamda t x mu y x=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)) %作图 1 t=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.87 tt=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.89 syms xx y=6600*(exp(-0.1155*xx)-exp(-0.1386*xx))-400; z=diff(y,xx) z = 22869/(25*exp((693*xx)/5000)) - 7623心 0*exp((231*xx)/2000)) %%% 算出来式子 -7623/10*exp(-231/2000*xx)+22869/25*exp(-693/5000*xx)=0 tmax=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.8927 t=-10000/231*log(5/6) 作业：对于该案例，结合教材P12 ，给岀具体的施救方案，解决如下问题：(1)作图，即求岀口服活性炭药物后模型( 8 ) 后，参考教材，画岀图形2；((473*exp(231/500))/2 - 1650*exp(231/1250))/exp((231*t)/1000) + (1650*exp((231*t)/2500))/exp((231*t)/1000) 解得 : 0.1386t 2.310t z(t) =1650* e 「1609.5* e ,t _ 2 （2）根据模型（ 8）计算出施救后血液达到最大值的时间。（即算结果 t3 =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心 00*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))/231 ans = 5.2595 T3=5.2595 （3） 要」使孩子在施救后Z（t）立即下降，算出排除率多大？（即算结果」=0.4886 ）（4）如果使用体外血液透析的方法，药物排除率可增加到」=0.1155*6=0.693，用这个」重新求解模型（ 7）并作图。z=112.2738*exp(-0.693*t)-285*exp(-0.1386*t);