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) 如