飲酒駕車問(wèn)題的數(shù)學(xué)模型(初步修正版)

1、飲酒駕車問(wèn)題的數(shù)學(xué)模型摘要本問(wèn)題是生活中的飲酒駕車問(wèn)題，酒精對(duì)人體的作用過(guò)程實(shí)際上類似于生物醫(yī)學(xué)中的藥用過(guò)程，針對(duì)飲酒方式的不同，本文將飲酒過(guò)程分成快速飲酒、某時(shí)間段內(nèi)勻速飲酒和周期飲酒三種形式來(lái)討論。并分別建立了一室快速飲酒、二室勻速飲酒以及周期飲酒三種系統(tǒng)動(dòng)力學(xué)模型，并運(yùn)用非線性最小二乘法進(jìn)行數(shù)據(jù)擬合得到相關(guān)參數(shù),從而得到了血液中酒精含量與時(shí)間的函數(shù)關(guān)系。結(jié)合模型，運(yùn)用MATLAB工具得到了快速飲用三瓶啤酒時(shí)的違規(guī)時(shí)間分布，t:0.0650.24小時(shí)內(nèi)飲酒駕車；t:0.244.5小時(shí)內(nèi)醉酒駕車；t:4.512小時(shí)內(nèi)飲酒駕車。結(jié)合模型，得到了在2個(gè)小時(shí)內(nèi)均勻飲用三瓶啤酒的違規(guī)時(shí)間分布，t:2

2、4.5小時(shí)內(nèi)為醉酒駕車；當(dāng)t為4.5-12小時(shí)為飲酒駕車。模型的建立,使問(wèn)題一以及問(wèn)題三得到了較為確切的解釋。關(guān)鍵詞: 動(dòng)力學(xué)；吸收速率；消除速率；模型一、問(wèn)題重述在2003年全國(guó)道路交通事故死亡人數(shù)中，飲酒駕車造成的占有相當(dāng)比例，為此，國(guó)家發(fā)布了新的車輛駕駛?cè)藛T血液、呼氣酒精含量閾值與檢驗(yàn)國(guó)家標(biāo)準(zhǔn)。在新標(biāo)準(zhǔn)下，大李在中午12點(diǎn)喝了一瓶啤酒，下午6點(diǎn)檢查時(shí)符合標(biāo)準(zhǔn)，接著晚上又喝了一瓶，但凌晨2點(diǎn)檢查時(shí)卻被定為飲酒駕車，為什么喝同樣多的酒，兩次檢查結(jié)果不一樣？建立飲酒時(shí)人體內(nèi)酒精含量與時(shí)間關(guān)系的數(shù)學(xué)模型，并討論快速或慢速飲3瓶啤酒在多長(zhǎng)時(shí)間內(nèi)駕車就會(huì)違反新標(biāo)準(zhǔn)，估計(jì)血液中的酒精含量在什么時(shí)間最高

3、，如果某人天天喝酒，是否還能開(kāi)車，并根據(jù)你所做的結(jié)合新國(guó)家標(biāo)準(zhǔn)寫一篇短文，給想喝一點(diǎn)酒的司機(jī)如何駕車提出忠告。二、符號(hào)說(shuō)明及模型假設(shè)2.1符號(hào)說(shuō)明-人體飲入酒精總量t-飲用酒的時(shí)間-t時(shí)刻血液中的酒精量-t時(shí)刻人體吸收的酒精量M-人的體重-人的體液占人的體重的百分含量-人的血液占人體重的百分含量-酒精在人體中的吸收速率常數(shù)1-酒精在人體中的消除速率常數(shù)1-t時(shí)刻血液中的酒精濃度F-酒精在人體中的吸收度V-人體的血液體積V酒-喝酒的體積-酒中的酒精含量-飲酒持續(xù)時(shí)間2.2基本假設(shè)1. 酒精在血液中的含量與在體液中的含量大至相同；2. 每瓶啤酒的酒精含量、體積基本相同；3. 酒精進(jìn)入人體后，不考慮

4、其他因素對(duì)酒精的分解作用；4. 如果在很短時(shí)間內(nèi)飲酒，認(rèn)為是一次性飲入，中間的時(shí)間差不計(jì)；5. 確定是否飲酒駕車或醉酒駕車以新的國(guó)家標(biāo)準(zhǔn)為界；6. 不管喝的是什么酒，只以涉入的酒精總量納入計(jì)算；7. 酒精按一級(jí)吸收過(guò)程進(jìn)入體內(nèi)；8. 正常情況下，酒精在各人體中的吸收和消除速率基本相同；9. 將慢速飲酒看作是一個(gè)勻速過(guò)程。三、問(wèn)題分析與模型建立3.1模型（快速飲酒模型）同藥物一樣，酒精進(jìn)入機(jī)體后，作用于機(jī)體而影響某些器官組織的功能；另一方面酒精在機(jī)體的影響下，可以發(fā)生一系列的運(yùn)動(dòng)和體內(nèi)過(guò)程：自用藥部位被吸收進(jìn)入血液循環(huán)；然后分布于各器官組織、組織間隙或細(xì)胞內(nèi)；有部分酒精則在血漿、組織中與蛋白質(zhì)結(jié)

5、合；或在各組織（主要是肝臟）發(fā)生化學(xué)反應(yīng)而被代謝；最后，酒精可通過(guò)各種途徑離開(kāi)機(jī)體（排泄）；即吸收、分布、代謝和排泄過(guò)程。它們可歸納為兩大方面：一是酒精在體內(nèi)位置的變化，即酒精的轉(zhuǎn)運(yùn)，如吸收、分布、排泄；二是酒精的化學(xué)結(jié)構(gòu)的改變，即酒精的轉(zhuǎn)化亦即狹義的代謝。由于轉(zhuǎn)運(yùn)和轉(zhuǎn)化以致形成酒精在體內(nèi)的量或濃度（血漿內(nèi)、組織內(nèi)）的變化，而且這一變化可隨時(shí)間推移而發(fā)生動(dòng)態(tài)變化。又因?yàn)榫凭写龠M(jìn)血液循環(huán)的作用2。而藥物動(dòng)力學(xué)模型中的一室模型3是指給藥后，藥物一經(jīng)進(jìn)入血液循環(huán)，即均勻分布至全身,故快速飲酒情況可通過(guò)建立一室模型求解。雖然酒精在體內(nèi)的分布狀況復(fù)雜，但酒精的吸收、分解等則都在系統(tǒng)內(nèi)部進(jìn)行，酒精進(jìn)入人

6、體后，經(jīng)一段時(shí)間進(jìn)入血液，進(jìn)入血液后，當(dāng)在血液中達(dá)最高濃度時(shí)，隨后開(kāi)始消除3，把酒精在體內(nèi)的代謝過(guò)程看為進(jìn)與出的過(guò)程，這樣便會(huì)使問(wèn)題得到簡(jiǎn)化。用和分別表示酒精輸入速率和輸出速率。由于單位時(shí)間內(nèi)血液中酒精的改變即變化率就等于輸入與輸出速率之差，所以其動(dòng)力學(xué)模型為：=- （1）又因?yàn)榫凭谘褐械南俾逝c當(dāng)時(shí)血液內(nèi)的藥量成正比，所以=kx，代入（1）式得： =-kx （2）則由（2）式可知x(t)的變化規(guī)律由飲酒速率而定。而酒精在人體內(nèi)的代謝可簡(jiǎn)單的由圖一表示：Fx0吸收室x1(t)kx(t)V (圖一)則t時(shí)刻吸收室的藥量為x1(t)，又藥物是按一級(jí)吸收過(guò)程進(jìn)入體內(nèi)的，對(duì)于吸收室有： =-k1

7、x1 (3)對(duì)于房室，=，于是（2）式變?yōu)椋?（4）（3）、（4）兩式構(gòu)成一階線性方程組，當(dāng)t=0時(shí)，x(0)=0,解（3）式得：，將其代入（4）式得一階線性非齊次方程：解之得：從而，人體內(nèi)酒精含量為：在這種情況下，酒精含量最大值出現(xiàn)的時(shí)間：使時(shí)t的值。一般情況下，又因?yàn)榫凭谘褐械暮颗c在體液中的含量大至相同。則有： (F為常數(shù)且0<F<1)X0=V酒則人體內(nèi)酒精含量與時(shí)間函數(shù)關(guān)系為: （一般情況）3.1.1模型的求解根據(jù)圖一中酒精含量實(shí)測(cè)數(shù)據(jù)擬合,顯然它無(wú)法化為線性最小二乘,我們直接作非線性最小二乘擬合4。用MATLAB優(yōu)化工具箱的Leastsq5計(jì)算,擬合參數(shù)程序見(jiàn):JM2

8、004C1.m。擬合得：=2.0079mg·ml-1·h-1,k=0.1855mg·ml-1·h-1,=11.2423（毫克/毫升），又由，得到問(wèn)題中隱含的一瓶啤酒的酒精量約為：27543.635毫克。3.1.2問(wèn)題一的解答雖然大李喝等量的酒，并且相隔的時(shí)間也相同的情況下,兩次檢查的結(jié)果不一樣是因?yàn)榈谝淮魏认氯サ木?，?小時(shí)內(nèi)并沒(méi)有完全分解，還殘留有相當(dāng)一部分在血液中，并且這一部分在較長(zhǎng)時(shí)間內(nèi)不能完全分解。由圖（二）喝一瓶啤酒的酒精含量隨時(shí)間變化的函數(shù)圖像可知，6小時(shí)后，第一次喝的酒的酒精含量x(6)約19.5毫克/百毫升。也就是說(shuō)此時(shí)血液中已有一定的酒

9、精量，這樣雖然第二次喝的是同樣多的酒，由于第一次殘留部分的存在，相當(dāng)于涉入的酒精量已增大了，使其再過(guò)8小時(shí)，酒精含將會(huì)大于20毫克/百毫升。這樣大李碰到的情況也就很自然的解釋了（圖二）現(xiàn)通過(guò)實(shí)際計(jì)算證明：設(shè)A為第一次喝酒在六個(gè)小時(shí)后，殘留在血液中的酒精量。則第二次喝酒時(shí)，x(0)=A=19.5毫克/百毫升,酒精含量C（t）與時(shí)間t的函數(shù)關(guān)系為： 代值用MATLAB計(jì)算（程序見(jiàn)JM2004C5.LOG），此時(shí)再過(guò)8小時(shí)酒精含量為：27.4毫克/百毫升。這就大于了第一次的值。3.1.3問(wèn)題二第一問(wèn)的求解當(dāng)一個(gè)人在很短時(shí)間內(nèi)喝3瓶啤酒時(shí)，相當(dāng)于在瞬間使其吸收室的酒精濃度達(dá)最大值，運(yùn)用模型并用計(jì)算機(jī)擬

10、合得其函數(shù)圖象如(圖三)（程序見(jiàn)JM20004C2.LOG）， (圖三)由圖可得： t:0.0650.24小時(shí)內(nèi)飲酒駕車；t:0.344.5小時(shí)內(nèi)醉酒駕車；t:4.512小時(shí)內(nèi)飲酒駕車。3.2模型 （勻速飲酒二室模型）3.2.1模型的建立及求解由于在兩小時(shí)內(nèi)慢速喝酒有多種方式，為了便于計(jì)算我們考慮為勻速飲入等量酒精，即在時(shí)間內(nèi)勻速吸收酒精。若體內(nèi)酒精含量不超過(guò)一級(jí)消除動(dòng)力學(xué)范圍，假設(shè)人的酒精含量未達(dá)到平衡狀態(tài)，隨著人體吸收次數(shù)增多，血液中酒精濃度逐漸升高，當(dāng)在時(shí)間內(nèi)飲完后，由于時(shí)刻前一小段時(shí)間內(nèi)飲入的酒精在短時(shí)刻內(nèi)沒(méi)有被吸收，故在時(shí)刻后一段時(shí)間內(nèi)，酒精濃度將繼續(xù)升高，某一時(shí)刻將達(dá)到最大值。這一

11、時(shí)刻后，酒精含量就會(huì)逐漸減小?？紤]到一室模型在求慢速喝酒情況下的局限性，我們采用建立二室模型6。對(duì)二室模型我們將建立兩個(gè)關(guān)于酒精濃度的動(dòng)力系統(tǒng)模型來(lái)描述其動(dòng)態(tài)特性。3.2.1.1模型二假設(shè)：1、 機(jī)體分為中心室和周邊室，兩個(gè)室的容積在過(guò)程中保持不變。2、 藥物從一室向另一室的轉(zhuǎn)移速率，及向體外的排除速率，與該室的酒精濃度成正比。3、 只在中心室一體外有酒精交換，即酒精從體外進(jìn)入中心室，最后又從中心室排出體外，與轉(zhuǎn)移和排除的數(shù)量相比，酒精的吸收可以忽略。3.2.1.2模型建立K12二室模型的示意圖如(圖四)所示： K21周邊室c2(t),x2(t) V2中心室c1（t）,x1(t)V1 飲酒 (

12、圖四)兩個(gè)房室中酒精量滿足的微分方程。的變化率由一室向二室的轉(zhuǎn)移，一室向體外排除，二室向一室的轉(zhuǎn)移及酒精組成；的變化率由一室向二室的轉(zhuǎn)移及二室向一室的轉(zhuǎn)移組成，于是有： （1）與血液中酒精含量、房室容積顯然有關(guān)系式（2）將（2）式代入（1）式可得： （3）喝酒相當(dāng)于在酒精進(jìn)入中心室之前先有一個(gè)將酒精吸收入血液的過(guò)程，可以簡(jiǎn)化為有一個(gè)吸收室，如下圖，為吸收室的酒精，酒精由吸收室進(jìn)入中心室的轉(zhuǎn)移速率系數(shù)為，于是滿足： （4） 當(dāng)時(shí)，（3）可以化為： （5）中心室吸收室 （圖四） 酒精經(jīng)吸收室進(jìn)入中心室是飲入的酒精量，而酒精進(jìn)入中心室的速率為： （6）將方程（5）式代入（6）式得： （7）在這種情況

13、下方程（3）的解的一般形式為： （8）此時(shí)，這種情況下，酒精含量最大值出現(xiàn)的時(shí)間為：使時(shí)的時(shí)間。其中 （9） 3.2.1.3參數(shù)估計(jì) 不妨設(shè)，于是當(dāng)t充分大時(shí)可近似為： 對(duì)于適當(dāng)大的和相應(yīng)的，用最小二乘法估計(jì)出的值。然后計(jì)算（10）再利用（5）式得：對(duì)于較小的和（10）式算出的，仍用最小二乘法可得到。由上可得參數(shù)值：3.2.1.4模型的求解將題中所給參考數(shù)據(jù)代入，運(yùn)用MATLAB工具，對(duì)前參數(shù)估計(jì)進(jìn)行代值計(jì)算得:程序見(jiàn)JM2004C6-C8A=69.4908;a=0.3234;B=-48.0251;b= 2.77;將所求得的參數(shù)代入（8）當(dāng)中可得，并用MATABT畫出勻速喝三瓶啤酒的酒精含量-

14、時(shí)間曲線圖如圖（五）程序見(jiàn)JM2004C9.M （圖五）由上函數(shù)圖象可知：當(dāng)t為24.5小時(shí)內(nèi)為醉酒駕車；當(dāng)t為4.5-12小時(shí)為飲酒駕車。3.3模型（周期飲酒模型）3.3.1模型假設(shè)1、每次喝酒相隔時(shí)間相同，且為常量。2、每次喝酒量相同3.3.2模型建立設(shè)表式t時(shí)刻酒精在人體內(nèi)的濃度，表示t=0時(shí)飲入酒精量。我們知道，飲酒后，隨時(shí)間推移，酒精在體內(nèi)逐漸被吸收，也就是體內(nèi)酒精的濃度逐漸降低，根據(jù)藥物學(xué)理論，酒精濃度的變化率與飲酒量成線性比，則有： （1）其中k>0為吸收速率常數(shù)，其解為： 當(dāng)時(shí)，由于經(jīng)過(guò)時(shí)間間隔T，又第二次飲酒，飲入量為，所以時(shí)類似有，當(dāng)時(shí)，酒精含量為：并且當(dāng)時(shí)，又第三次

15、飲酒，飲酒量仍為，所以， 那么當(dāng)時(shí)，體內(nèi)酒精濃度應(yīng)為：按此類推，則當(dāng)時(shí)，體內(nèi)酒精濃度達(dá)到： （2）（2）式的右邊為一等比數(shù)列之和，利用等比數(shù)列求和得： 當(dāng)時(shí)，根據(jù)自己要求，如果我們需要酒精含量接近時(shí)，我們近似地有如果間隔時(shí)間T確定，那么飲入酒精量可由： （3）體內(nèi)酒精的分布可由圖六看出多飲幾次酒后，體內(nèi)酒精濃度緩慢趨于極限值。（圖六）由上模型，代入?yún)?shù)計(jì)算，也可以解決問(wèn)題一，取T=6時(shí)，酒精含量為19.98毫克/毫升，取=12時(shí)，酒精含量為27.6毫克/毫升。此結(jié)果與前面用一次模型求解的結(jié)果吻合，這樣更全而的解釋了相等飲量，間隔相等時(shí)間，前次達(dá)標(biāo)，而后一次超標(biāo)的原因。運(yùn)用此模型，更方便解決天天

16、飲酒是否能駕車的問(wèn)題，令T=24，毫克/百毫升，代入（3）式得，毫克，也就是說(shuō)只要每天的飲酒量小于20000毫克就不會(huì)出現(xiàn)違規(guī)情況，可以駕車；令毫克/百毫升，毫克，即飲酒量大于80000毫克時(shí)就會(huì)出現(xiàn)醉酒駕車情況，就不能駕車。當(dāng)每飲酒量為20000<<80000毫克時(shí)，就會(huì)出現(xiàn)飲酒駕車情況。類似，本模型還可推廣到更長(zhǎng)周期情況的求解。綜上所述，現(xiàn)對(duì)第三問(wèn)解答：三種情況下，因考慮到擬合函數(shù)只存在一個(gè)駐點(diǎn)，即一個(gè)極值，則可對(duì)擬合函數(shù)求導(dǎo)，,然后計(jì)算駐點(diǎn)值即可解得各種情況下酒精含量達(dá)最高時(shí)的時(shí)間，在快速飲酒情況下，飲三瓶啤酒為例，解得Tmax=1.3145小時(shí)。四、結(jié)果分析與檢驗(yàn)由于酒精不

17、同于一般的藥物，人飲用之后，酒精能迅速的進(jìn)入人體，所以在快速飲酒情況下用一室模型來(lái)研究酒精的作用過(guò)程及代謝過(guò)程，與幾室模型來(lái)研究酒精的作用過(guò)程及代謝過(guò)程并無(wú)太大的區(qū)別，（圖七）酒精含量C（t）的擬合曲線與實(shí)測(cè)數(shù)據(jù)的比較,可以看出,我們的結(jié)果是較為準(zhǔn)確的，并且擬合的曲線是一條較為理想的曲線。（圖七）對(duì)于勻速飲酒情況，由于勻速飲酒要考慮飲酒過(guò)程中的代謝情況，與快速飲酒有一定差異，所以采用二室模型來(lái)建立勻速飲酒模型，達(dá)到減小誤差的目的。 五、模型的評(píng)價(jià)與改進(jìn)5.1模型的優(yōu)點(diǎn)1.綜合運(yùn)用MATLAB和LINGO兩個(gè)軟件，準(zhǔn)確求解，在運(yùn)用MATLAB進(jìn)行數(shù)據(jù)擬合時(shí)，得到了較理想化的曲線。在表示喝三瓶啤酒

18、的人什么時(shí)候是飲酒駕車，什么時(shí)候是醉酒駕車時(shí)，運(yùn)用MATLAB準(zhǔn)確的做出了函數(shù)據(jù)圖像，使結(jié)果一目了然。2.本模型從三種情況分別建立模型，模型穩(wěn)定性高，適用性強(qiáng)。3.本模型計(jì)算步驟清晰，引用了醫(yī)藥動(dòng)力學(xué)的二室模型進(jìn)行計(jì)算，可靠性較高4.從問(wèn)題出發(fā)，分析了應(yīng)該考慮的各種情況，建立了一般的數(shù)學(xué)模型，并進(jìn)行實(shí)例驗(yàn)證，從而證明我們建立的數(shù)學(xué)模型可以較好的解決實(shí)際問(wèn)題。5.此模型具有極為廣泛的應(yīng)用性，對(duì)每一個(gè)具體的情況，都可以通過(guò)模型求解。5.2模型的缺點(diǎn)1.本文的模型參數(shù)僅是依靠題中給出的一組數(shù)據(jù)擬合求解得出，可能有一點(diǎn)偏差。2.模型為使計(jì)算簡(jiǎn)便，使所得的結(jié)果更理想化，忽略了一些次要的因素。如：酒進(jìn)入身

19、體后隨著血液流動(dòng)，人體對(duì)酒精的吸收率是隨時(shí)間變化的，而本模型是在吸收率恒定的情況下，進(jìn)行求解的。對(duì)于這些問(wèn)題，由于時(shí)間關(guān)系本模型還未能更好的研究，有待以后的改進(jìn)和完善。5.3模型的改進(jìn)由于人體內(nèi)部的復(fù)雜性，及各器官對(duì)酒精轉(zhuǎn)化的多樣性，用一室或二室都較為初級(jí)，三室或多室的情況模型更準(zhǔn)確，但考慮起來(lái)會(huì)很復(fù)雜，又由于短期收集資料的局限，實(shí)行起來(lái)較為困難，可留著時(shí)間充裕時(shí)考慮。六、模型的推廣我們建立模型的方法和思想可以推廣到其它類似的問(wèn)題。本文所建立的模型不僅估算出了t時(shí)刻人體內(nèi)的酒精含量，而且還能給飲酒駕車的司機(jī)們?nèi)绾物嬀铺峁┮恍├碚搮⒖?。又如：現(xiàn)在有些感冒藥對(duì)人的大腦有刺激作用，當(dāng)其血藥濃度高于某

20、標(biāo)準(zhǔn)的時(shí)候就不能進(jìn)行駕車等一系列安全操作，于是我們也可以用此模型的研究結(jié)論來(lái)對(duì)其相關(guān)問(wèn)題進(jìn)行分析。下面是給想喝點(diǎn)酒的司機(jī)如何駕車的忠告：酒后駕車危害多給想喝點(diǎn)酒駕車的司機(jī)們的忠告俗語(yǔ)說(shuō)：“美酒佳肴”美酒自古以來(lái)對(duì)人的誘惑從未衰減。多少人因貪杯而命喪黃泉。據(jù)統(tǒng)計(jì)，酒后駕車發(fā)生事故的比率為沒(méi)有飲酒情況下的16倍，幾率高達(dá)27，為了你的安全，請(qǐng)你注意以下信息。一.人的健康飲酒量肝臟處理酒精的能力，按體重每公斤每小時(shí)計(jì)算可處理0.125ml。體重為70公斤的人1小時(shí)能處理8.75ml，即相當(dāng)于能處理清酒約60毫升，啤酒約200ml，威士忌酒約20ml?，F(xiàn)在綜合對(duì)酒的處理能力與免疫學(xué)調(diào)查，可以得出以下結(jié)

21、論：健康的安全性飲純酒量每日為50ml以內(nèi)，有害量是每日100ml，危險(xiǎn)量是每日150ml以上。二.過(guò)量飲酒對(duì)人體的危害飲酒駕車，是造成交通事通行證的重要原因之一，酒精被胃、腸吸收后深于血液當(dāng)中，當(dāng)血液中酒精濃度達(dá)到一定程度時(shí)，中樞神經(jīng)系統(tǒng)活動(dòng)逐漸遲鈍，致使大腦判斷發(fā)生障礙，手腳遲鈍不靈活，甚至喪失操作能力。1.在血液中，酒精含量在0.5-2mg/毫升時(shí)，造成微醉。表現(xiàn)為臉紅、話多、反應(yīng)遲鈍、做事不顧后果，但尚未忘記自我。2.酒精含量在2-3mg/毫升時(shí)，造成輕醉。表現(xiàn)為言語(yǔ)不清、哭笑失常。3.酒精含量在3-4mg/毫升時(shí)，造成深醉。表現(xiàn)為腿腳發(fā)軟，動(dòng)作失調(diào)，陷入麻痹狀態(tài)。4.酒精含量在4-5

22、mg/毫升時(shí)，造成泥醉。表現(xiàn)為陷入昏睡狀態(tài)，四肢無(wú)力，甚至造成大小便失禁，呼吸困難，最終可能導(dǎo)致死亡。雖然飲酒駕車危害甚多，但并不是說(shuō)一點(diǎn)都不能喝酒。甚至還可以天天喝，但一定要注意控制自己的飲酒量和出車時(shí)間，結(jié)合上面的信息，注意以下幾點(diǎn)，想喝一點(diǎn)酒的司機(jī)們也能過(guò)一把酒癮。1.如果你想每天即飲酒又駕車，而又不違規(guī)，請(qǐng)你一定記住你每天涉入的酒精量不要超過(guò)20000毫克。2.一次性飲酒的酒精量越大，到達(dá)標(biāo)時(shí)的時(shí)間會(huì)越長(zhǎng)，所以你等待時(shí)間的長(zhǎng)短應(yīng)根據(jù)你飲酒量的多少而定。比如說(shuō)一次飲一瓶啤酒，大約6個(gè)小時(shí)后酒精含量就可達(dá)標(biāo)；一次性喝2瓶啤酒，大概要等9.5小時(shí)后才能達(dá)標(biāo)；而一次性喝3瓶啤酒，則大概要等12

23、小時(shí)后才能達(dá)標(biāo)。3.連續(xù)飲酒次數(shù)越多，每次間隔時(shí)間應(yīng)越長(zhǎng)。以第一題為例，第一次飲啤酒一瓶，過(guò)六個(gè)小時(shí)達(dá)標(biāo)，但第二次飲同樣多的酒，同樣再過(guò)六個(gè)，酒精含量增加到27毫克/百毫升，要使第二次飲酒后，不超標(biāo)，則至少應(yīng)在7.5小時(shí)后再駕車。參考文獻(xiàn)：1 卓先義.血中酒精消除速度與濃度推算關(guān)系的研究2中華醫(yī)學(xué)網(wǎng).健康常識(shí). 3 張雙德.醫(yī)用高等數(shù)學(xué).天津：天津科學(xué)技術(shù)出版社，1999-84 蕭樹(shù)鐵.大學(xué)數(shù)學(xué)實(shí)驗(yàn).北京：高等教育出版社，1999-7，199-20015王沫然.MATLAB6.0與科學(xué)計(jì)算.北京:電子工業(yè)出版社,20006姜啟源.數(shù)學(xué)模型，第三版.北京：高等教育出版社程序附頁(yè):1. JM200

24、4C1.mfunction f=fun(x)t=0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16;c=0.3 0.68 0.75 0.82 0.82 0.77 0.68 0.68 0.58 0.51 0.50 0.41 0.38 0.35 0.28 0.25 0.18 0.15 0.12 0.10 0.07 0.07 0.04;f=c-x(1)*x(3)*(exp(-x(2)*t)-exp(-x(1)*t)/(x(1)-x(2);x0=1:0.25:20;x=leastsq('fun',x0

25、);y=sum(fun(x).*fun(x)2、程序JM2004C2.LOGmodel:sets:endsetsc=k1*(exp(-k*t)-exp(-k1*t)*d/(k1-k);DATA:k1=2.0079;k=0.1855;d=168.65;t=12;ENDDATAEND3、JM2004C3.LOGmodel:sets:endsetsc=k1*(exp(-k*t)-exp(-k1*t)*d/(k1-k);DATA:k1=2.0079;k=0.1855;d=56.217;t=12;ENDDATAEND4、JM2004C4.LOGmodel:sets:endsets82630.905/(k

26、*v)=k0/k*log(v)-k0/k*log(v)*exp(-k*t);data:k=0.185;v=6938.082065;enddata5、JM2004C5.LOGmodel:sets:endsetsc=k1*(27543.635+9555)*(exp(-k*t)-exp(-k1*t)/(V*(k1-k);DATA:k1=2.0079;k=0.1855;V=490;t=6;ENDDATAEND6、Jm2004C6.mfunction f=fun(x)t=0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16;c=