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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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