飲酒與駕車的關(guān)系設(shè)計論文

1、 攀枝花學院本科數(shù)學模型設(shè)計飲酒駕車的數(shù)學模型 學生姓名： 王 哲 辰 同組學生姓名： 王祝鵬，譙志成 學生學號： 200910801070 200910801071 200910801057 院 （系）： 計算機學院 年級專業(yè)：2009級計算機科學與技術(shù)2班 指導教師： 李 思 霖 助理指導教師： 無 二一年十二月三十日飲酒與駕車的關(guān)系摘 要針對酒后駕車問題，本文從實際情況出發(fā)并結(jié)合題設(shè)作出了合理的假設(shè)，并考慮了影響體液中酒精含量變化的因素吸收過程和代謝過程的同時進行，建立了一般微分方程模型。進而考慮到飲酒的快慢的不同，于是在一般模型的基礎(chǔ)上又得出了兩個具體的模型。接下來用常數(shù)變易法對各個模

2、型進行求解，而且由題得知題中所給的數(shù)據(jù)符合具體模型，于是我們用最小二乘法并借助于Matlab軟件，對模型所得的體液中酒精濃度關(guān)于時間的表達式進行擬合，得到了模型中吸收能力系數(shù)和代謝能力系數(shù)的值分別為2.0261和0.1842，到此時我們得到了可以用來解答題中問題的模型的解。于是我們就根據(jù)已知條件分別判斷各問題符合的模型，再逐步地采用Mathematica軟件來一一解答了各個問題。解答問題1時，得到下午6點的體液中酒精的含量為18.8198毫克/百毫升，小于國家標準規(guī)定的20毫克/百毫升，而凌晨2點的為20.3968毫克/百毫升，大于國家標準規(guī)定的20毫克/百毫升，如此通過定量的分析合理的解釋了

3、大李遇到的情況；在問題2中我們得到兩種情況下分別在飲酒后11.6341小時、12.7169小時內(nèi)駕車就會違反國家新標準，而通過對題中數(shù)據(jù)的分析比較可知此結(jié)果是很符合實際的；問題3中采用了作圖的方法并結(jié)合了高等數(shù)學中關(guān)于求駐點及求極值的相關(guān)知識，在問題2的前提條件下結(jié)合兩種不同的情況來說明如何估計體液中酒精含量達到最大時的時間，如：在問題2中對兩種飲酒方式求得分別在飲酒后1.35067小時和2.62436小時時體液中酒精含量達到最大值：124.638毫克/百毫升和116.682毫克/百毫升；對于第四個小問題，先根據(jù)模型求解出酒精被吸收的百分比，再由求無窮等比數(shù)列的極限得出體內(nèi)酒精濃度的極限值，進

4、一步算出每天酒精涉入量的極限安全值為8399.7毫克,相當于0.387瓶啤酒所含的酒精量。本文的最大特色在于運用合理的假設(shè)建立了一個針對酒后駕車問題的一般的微分方程模型，然后用常數(shù)變易法求出了其通解，這使得我們一開始就站在一個很高的理論高度，在這個理論的指導下，各種具體問題都變得非常容易求解，這就極大的簡化了計算。此外，我們還對一般的模型進行了誤差和靈敏度分析，利用微分方程的穩(wěn)定性理論嚴格的證明了微分方程對初值與非齊次項都是漸進穩(wěn)定的，這表明我們的模型是完全可行的。一、 問題的提出據(jù)報載，2003年全國道路交通事故死亡人數(shù)為10.4372萬，其中因飲酒駕車造成的占有相當?shù)谋壤?。針對這種嚴重的道

5、路交通情況，國家量監(jiān)督檢驗檢疫局2004年5月31日發(fā)布了新的車輛駕駛?cè)藛T血液、呼氣酒精含量閾值與檢驗國家標準。新標準規(guī)定，車輛駕駛?cè)藛T血液中的酒精含量大于或等于20毫克百毫升，小于80毫克百毫升為飲酒駕車（原標準是小于100毫克百毫升），血液中的酒精含量大于或等于80毫克百毫升為醉酒駕車（原標準是大于或等于100毫克百毫升）。大李在中午12點喝了一瓶啤酒，下午6點檢查時符合新的駕車標準，緊接著他在吃晚飯時又喝了一瓶啤酒，為了保險起見他呆到凌晨2點才駕車回家，又一次遭遇檢查時卻被定為飲酒駕車，這讓他既懊惱又困惑，為什么喝同樣多的酒，兩次檢查結(jié)果會不一樣呢？請你參考下面給出的數(shù)據(jù)（或自己收集資料

6、）建立飲酒后血液中酒精含量的數(shù)學模型，并討論以下問題：1. 對大李碰到的情況做出解釋；2. 在喝了3瓶啤酒或者半斤低度白酒后多長時間內(nèi)駕車就會違反上述標準，在以下情況下回答：1） 酒是在很短時間內(nèi)喝的；2） 酒是在較長一段時間（比如2小時）內(nèi)喝的。3. 怎樣估計血液中的酒精含量在什么時間最高。參考數(shù)據(jù)1. 人的體液占人的體重的65%至70%，其中血液只占體重的7%左右；而藥物（包括酒精）在血液中的含量與在體液中的含量大體是一樣的。2. 體重約70kg的某人在短時間內(nèi)喝下2瓶啤酒后，隔一定時間測量他的血液中酒精含量（毫克百毫升），數(shù)據(jù)如下：時間(小時)0.250.50.7511.522.533.

7、544.55酒精含量306875828277686858515041時間(小時)678910111213141516酒精含量3835282518151210774二、 問題假設(shè)為了更簡便的解決問題，我們在研究這個問題的過程中作出以下假設(shè)：1. 假設(shè)整體過程中人沒有攝入任何影響代謝的藥類物質(zhì)和作劇烈性運動。2. 人的吸收速率和代謝速率是恒定的。3. 血液與體液中酒精的濃度相同。4. 酒精代謝速率與當前血液中酒精濃度成正比。5. 人體體液對酒精吸收速率與當前腸胃中酒精含量成正比。6. 人體自身產(chǎn)生的酒精忽略不計。7. 忽略不同人對酒精代謝能力的差異。三、符號說明變 量含 義單 位備 注吸收能力系數(shù)

8、1/小時代謝能力系數(shù)1/小時開始飲酒時人體體液中的酒精濃度毫克/百毫升人喝下酒精的總量毫克一瓶啤酒中酒精的量約為21700毫克5人體液所占的體積百毫升約為420百毫升4時間小時較長時間飲酒時的持續(xù)時間小時酒精由腸胃進入人體體液的速率毫克/小時酒精由口進入腸胃的速率毫克/小時t時刻體液中酒精的量毫克時刻腸胃中酒精的量毫克在t時刻人體體液（或血液）中的酒精濃度毫克/百毫升四、問題分析及模型的建立1問題分析依常識，我們知道酒精無需經(jīng)過消化系統(tǒng)即可被腸胃直接吸收，進入血管，在幾分鐘后迅速擴散到人體全身，且參與代謝，即吸收過程和代謝過程是同時進行的。由此我們可以知道血液中酒精濃度同時受吸收和代謝影響。而

9、由假設(shè)和題目的條件可知，體液中酒精濃度和血液中的酒精濃度是一樣的，所以我們只考慮體液中酒精濃度的變化。接下來我們將分別對吸收過程和代謝過程作分析和討論：（1）吸收過程在此過程中我們考慮的是腸胃吸收酒精進入體液所引起的酒精量的增量，其中用代表時刻酒精由腸胃進入人體體液的速率（單位：毫克/小時），則時間段內(nèi)體液中酒精量的增量可表示為：。再以表示人喝下酒精的總量，以剛開始飲酒的時刻為記時的初始時刻，即=0，由于喝下的酒精在時間無窮長時總會被全部吸收。（2）代謝過程在此過程中我們考慮的是機體代謝排出體液中酒精所引起的量的變化，這里我們認為酒精的代謝的速率與當前體液中含有的酒精量成正比關(guān)系，設(shè)比例系數(shù)為

10、b，以表示t時刻體液中酒精的量，則時間段內(nèi)體液中酒精量的減少量可表示為：。結(jié)合以上兩個過程，依據(jù)量守恒定律，在時間段內(nèi)有：體液中酒精量的變化量 = 吸收進入體液的酒精量 - 代謝了的酒精的量即得方程：兩邊同除后讓取極限得微分方程：根據(jù)體液中酒精的量與體液中酒精濃度間的關(guān)系，在上式兩邊同時除以得酒精濃度滿足的微分方程：2模型建立根據(jù)以上分析，我們可以得出血液中酒精濃度關(guān)于時間變化的一般模型：一般模型： 由題中可知，酒可以在很短時間內(nèi)喝完，也可以在一段較長的時間內(nèi)喝完，這樣將有兩個不同的吸收速率，即可以得到兩個不同的表達式，于是可以得到兩個具體模型。第一個模型是考慮短時飲酒效應，第二個模型考慮的是

11、長時飲酒效應。具體模型的建立：由于考慮的是短時的飲酒效應，可以認為人飲入的酒量在最初時間就全部存在于機體的腸胃中，以表示時刻腸胃中酒精的量，即有,從生物機體對酒精吸收的規(guī)律，可以知道腸胃中酒精的減少速率與剩余量成正比，設(shè)比例系數(shù)為，由此有微分方程：而又很容易想到酒精進入人體體液的速率與腸胃中酒精的減少速率是一個相當過程，這樣便可以得到：再結(jié)合一般模型，我們可以得出如下模型具體模型： 具體模型的建立：在這里研究的是長時飲酒效應，可近似認為在持續(xù)飲酒的過程中酒精是勻速進入腸胃的，參照模型可有，在此我們引入函數(shù)來表示酒精進入腸胃的速率（單位：毫克/小時），表示飲酒時的持續(xù)總時間，則酒精進入腸胃的速率

12、與整個過程中喝入的酒精量有如下關(guān)系： 而腸胃里的酒精量的變化與機體對酒精的吸收和喝入機體的酒精量都有關(guān)系，而機體對酒精的吸收速率可以與模型一樣用來表示，則時間段內(nèi)有：腸胃里的酒精量的變化量 = 喝入機體的酒精量 - 機體對酒精的吸收量即得：兩邊同除后讓取極限得：而與模型一樣有： 再結(jié)合一般模型我們可以得到模型如下：具體模型： 五、模型求解 對于本題的模型求解，我們分以下幾個步驟進行： Step1:對一般模型進行求解。 Step2:根據(jù)Step1求得的結(jié)果分別對具體模型和具體模型進行求解，得到兩個模型的解的表達式。 Step3:根據(jù)已有數(shù)據(jù)應用最小二乘法對具體模型進行數(shù)據(jù)擬合，分別獲得參數(shù)與的值

13、。從而得出具體模型和具體模型的解的具體表達式。 Step4:利用Step3的結(jié)果分別對題中的每一問進行解答。 1一般模型的求解 觀察一般模型，這是一個一階線性常系數(shù)非齊次微分方程，其特征方程為：，故其特征根為：，于是此微分方程對應的齊次方程有通解：為任意常數(shù) 運用常數(shù)變易法，令代入原方程，得： （k為任意常數(shù)） 于是得原方程通解： 由初始條件得：，故特解為：2具體模型的求解由及得：，又，于是有如下結(jié)論：（此時也是滿足的。）將其代入一般模型的特解中并化簡得：3具體模型的求解根據(jù)題設(shè)，我們?nèi)?。由及可得：又，可得：將其代入特解中并化簡得：4參數(shù)估計及具體解根據(jù)題中所給數(shù)據(jù)以及所給的條件可知所給數(shù)據(jù)

14、符合依據(jù)具體模型求得的酒精濃度關(guān)于時間變化的表達式，這里我們通過資料搜集得到信息人體中含有體液的總體積約為420百毫升4，由于常識我們很容易的知道了一瓶啤酒含有的酒精量為21700毫克。接下來我們運用最小二乘法結(jié)合Matlab軟件擬合出、的值分別為：于是具體模型的解為：圖像為：（其中點為題中已知的數(shù)據(jù)點，曲線為擬合的曲線。）具體模型的解為：圖像為：5問題1的解答在此問題中，大李是在下午6點時接受檢查的， 首先考慮到他喝啤酒的時間很短，則我們將此處理為具體模型的情景。在求解過程中我們以喝啤酒時刻為計算的初始時刻，而根據(jù)假設(shè)5，可以知道，又，把數(shù)據(jù)代入具體模型的解中求出第一次他在喝一瓶后6小時的血

15、液中酒精濃度為：（毫克/百毫升）小于20毫克/百毫升,即通過檢查。7小時后血液中酒濃度為15.654毫克/百毫升。此時再喝下一瓶啤酒，我們以此刻為初始時刻，初始濃度為15.654毫克/百毫升，即：（毫克/百毫升）將此代入模型的解，依然借用Mathematica軟件得出凌晨2點時血液中酒精濃度為：（毫克/百毫升）大于20毫克/百毫升，即未通過檢查。6問題2的解答1）從題可知此小題的情況符合模型，且由題可知：將此代入模型的具體解，并計算出的臨界值，也就是當時的值?？山馊缦路匠蹋簯肕athematica軟件解得（小時）即短時間喝完3瓶啤酒后11.3641小時內(nèi)駕車出行就會違反標準。2）從題可知此小

16、題的情況符合模型，且由題可知：將此代入模型的解，并計算出的臨界值，也就是當時的值。這里我們同樣假設(shè)T=2，由圖像觀察可知人體內(nèi)酒精濃度的降低在大于2的區(qū)域，于是可解如下方程：應用Mathematica軟件解得：（小時）即用2小時喝完3瓶啤酒的方式，在開始飲酒后12.7169小時內(nèi)駕車出行就會違反標準。7問題3的解答我們首先假設(shè)只喝3瓶啤酒，用作圖的方式發(fā)現(xiàn)無論啤酒是在短時間內(nèi)喝的還是在較長一段時間內(nèi)喝的，在時只存在一個極值點并且當不斷增大時曲線逐漸趨近于軸。于是我們可以用求的駐點的方法來求的最大值點即酒精含量最高的點。在模型中，我們運用Mathematica軟件求解方程得：（小時）此時最大濃度

17、（毫克/百毫升）。我們在模型的圖像中發(fā)現(xiàn)最大值出現(xiàn)在的區(qū)域內(nèi)，我們運用Mathematica軟件求解方程，得（小時）此時最大濃度（毫克/百毫升）。8問題4的解答在本問中我們假設(shè)每天喝進啤酒的量是，每天只喝一次，是短時間進酒，兩次間隔24小時。一般地講，如果天天飲酒，他們都喜歡短時間飲酒，因此采用模型的結(jié)果：將的值代入上式得到當?shù)谝惶煅褐芯频臐舛葹闀r24小時之后血液中還有未被代謝。進而第二天血液中酒的濃度為：第三天血液中酒的濃度為：第n天血液中酒的濃度為：當時血液中酒的濃度為：若司機還能開車，就要使，即，又，進而計算出毫克。因為一瓶啤酒中酒精的量約為21700毫克，由此計算出每天最多可喝0.3

18、87瓶。六、誤差與靈敏度分析首先，根據(jù)微分方程解的存在唯一性定理可知，只要在上連續(xù)，則對于任意實數(shù)及，方程在上存在唯一解滿足初始條件。其次，由于用微分方程描述的實際問題，其特解密切依賴于初始值和非齊次項，如本文中開始飲酒時人體體液中的酒精濃度和人飲酒的速率的快慢的度量，而初始值往往不能準確得到，于是我們必須考慮其會不會很嚴重的影響微分方程的解。這就是我們要考慮的方程的穩(wěn)定性問題，實際上我們有下面的兩個定理： 定理1：微分方程的解對于初值是漸進穩(wěn)定的。證明：上述微分方程的解為，給初值一個微小的擾動，即，則解為：于是，。由于，故 ，取，有 。（證畢）定理1表明方程的解對初值并不敏感，因此我們不考慮

19、人體內(nèi)的酒精含量0.003%是完全可行的。定理2：微分方程的解對于非齊次項是漸進穩(wěn)定的。證明：上述微分方程的解為考慮不同的非齊次項和，且，有 ，則：因此只要取， 那么，有。（證畢）由定理2可知非齊次項對解的影響不大，即可以知道由一般模型推得具體模型和具體模型是可行的。九、參考文獻1 費培之，數(shù)學模型實用教程，成都，四川大學出版社，1998年。129頁2 張志讓等，數(shù)學實驗，科學出版社，1999年。88頁。3 云舟工作室，MATLAB 6 數(shù)學建?；A(chǔ)教程，北京，人民郵電出版社，2001年。4 龔茜玲，人體解剖生理學，人民衛(wèi)生出版社，2001年。第124頁。5 千龍新聞網(wǎng)，酒在人體內(nèi)的吸收， :/finance.sina ,2003-01-23 。 十、附錄圖