药物中毒数模作业matlab程序

1、药物中毒数模作业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)>

2、;> 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=fze

3、ro('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 xxy=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*

4、exp(-231/2000*xx)+22869/25*exp(-693/5000*xx)=0e。1386t 1609.5* e ,t _ 2,= 回View Insert Tools Desktop Window Help盘m、工踮Q房二 DE T otmax=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)作业

5、：对于该案例，结合教材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.1386t2.310tz(t) =1650* e-1609.5* e ,t _ 2EJ Figure 1； I 回File Edit View Insert Jools Desktop Window Help七ON。的瀚年/总|日国|口|120010

6、00BOO600I 4002000 0510152025(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/(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)/231 ans =5.2595T3=5.2595(3) 要R使孩子在施救后z(t)立即下降，算出排除率多大？(即算结果 R =0.4886)(4) 如