全文預覽已結束
下載本文檔
版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
收稿日期 2005 08 29 作者簡介 楊錦偉 1983 男 河南省葉縣人 平頂山學院數學與信息科學學院助教 飲酒駕車的數學模型 楊錦偉 刁 群 平頂山學院 河南 平頂山467002 摘 要 定義了血醇濃度 借鑒藥物動力學中的房室模型 建立了血醇濃度的微分方程模型 準確地模 擬出血醇濃度變化趨勢的曲線 該模型可以預測喝酒后任一時刻的血醇濃度 同時對于血醇濃度達到最大值的 時間也給出了一個精確的解析解 最后還討論了周期性連續(xù)喝酒條件下血醇濃度的變化情況 關 鍵 詞 飲酒 微分方程模型 數據擬合 中圖分類號 O141 4 文獻標識碼 A 文章編號 1673 1670 2006 02 0056 03 1 問題提出 針對交通事故高發(fā)的嚴重情況 新的酒精含量標準規(guī) 定 車輛駕駛人員血液中的酒精含量大于或等于20mg 100mL 小于80mg 100mL為飲酒駕車 血液中的酒精含量 大于或等于80mg 100mL為醉酒駕車 大李在中午12點喝 了一瓶啤酒 下午6點檢查時符合標準 緊接著他在吃晚飯 時又喝了一瓶啤酒 凌晨2點檢查時卻被定為飲酒駕車 兩 次喝同樣多的酒 檢查結果卻不一樣 筆者建立飲酒后血液 中酒精含量的數學模型 并討論以下問題 1 對大李碰到的情況做出解釋 2 在喝了3瓶啤酒或者半斤低度白酒后多長時間內 駕車就會違反上述標準 a 酒是在很短時間內喝的 b 酒是在較長一段時間 比如 2h 內喝的 3 怎樣估計血液中的酒精含量在什么時間最高 4 根據你的模型論證 如果天天喝酒 是否還能開車 2 問題分析 當人們喝酒時 酒精 即乙醇 就會隨著血液擴散到全 身的器官和組織中 然后不斷地被吸收 代謝和排出 在這 里我們把乙醇在血液中的濃度 即每100mL血液中乙醇含 量 mg 稱血醇濃度 它隨時間和所處范圍的變化而變化 為了建立飲酒后血液中酒精含量的數學模型 我們需 要討論乙醇在體內吸收 分布和排出的動態(tài)過程以及此過 程中的定量關系 為了達到這一目的 我們借鑒藥物動力學 中的房室模型 所謂房室 是指由具有相近的藥物轉運速率 的器官 組織組合而成 同一房室內各部分的藥物處于動態(tài) 平衡 可以認為乙醇在體內各部位間均有較高及相近的轉 運速率 可在體內迅速達到分布平衡 即血醇濃度具有均衡 性 但是在不同的房室則可按一定的規(guī)律進行乙醇的轉移 在這里我們簡單地認為機體是一個血液均勻分布的中心 室 喝酒后 乙醇的擴散與分布可看作在進入中心室之前有 一個將乙醇由腸胃吸收入血液的過程 將之簡化為一個吸 收室 這種簡化在醫(yī)學界和藥理學界已被接受和承認 我 們稱之為一室模型 1 下面在吸收室和中心室之間建立一 個血醇濃度的微分方程來描述其動態(tài)特征 3 模型的假設與符號表示 在不失真的前提下 我們可以建立線性常系數方程 為 此 我們作如下假設 1 機體可看作一個中心室 中心室的容積 即血液容 積 在過程中保持不變 乙醇由腸胃吸收入血液的過程簡化 為一個吸收室 2 乙醇從吸收室向中心室的轉移速率與吸收室的血醇 濃度成正比 中心室向體外的排出速率與中心室的血醇濃 度成正比 3 乙醇全部由吸收室進入中心室 又從中心室排出 在 此過程中 吸收室中轉換的乙醇可忽略不計 4 假設大李在下午6點接受檢查后由于現實原因等8 點開始吃晚飯 5 假設酒在很短時間內喝完即將酒瞬時導入胃腸 沒 有耽擱 沒有擴散 6 假設酒在較長時間內喝是勻速的 7 不考慮人體的個體差異 酒精擴散過程示意如下 吸收室c1 t k1 中心室 c t k 符號表示 c1 t 吸收室中血醇濃度 c t 中心室 中血醇濃度 k1 轉移速率系數 k 排除速率系數 v1 吸 收室體積 v 中心室體積 第21卷第2期 2006年4月 平頂山學院學報 Journal of Pingdingshan University Vol 21 No 2 Apr 2006 4 模型的建立與求解 4 1 短時間內喝酒血醇濃度的確定 4 1 1 微分方程的建立與求解 設t 0時刻 一次喝的酒中乙醇含量為D 則c1 t 滿 足方程 2 dc1 dt k1c1 c1 0 D v1 1 中心室酒精濃度 c t 的變化率由兩部分組成 與c成 正比的排出 比例系數為k 與c1成正比的吸收 比例系數 為k1 則可得 c t 滿足方程 dc dt kc v1 v k1c1 c 0 c0 2 由 1 可解出 c1 t D v1 e k1t 3 將 3 式代入 2 中可解得 c t D v k1 k1 k e kt e k1t c0e kt 4 可以看出 t c1 t 0 c t 0 4 1 2 c t 的最大值問題 在 4 中 記b D v 取初始濃度c0 0 則 c t b g k1 k1 k e kt e k1t 容易求出 當時間t取tm 1 k k1ln k k1 時 c t 取得最大 值 cm c tm b k1 k1 k e k tm e k1tm b k1 k1 k k k1 k k k1 k k1 k1 k k1 5 4 2 較長時間內連續(xù)飲酒血醇濃度的確定 在較長一段時間內連續(xù)飲酒 持續(xù)時間為 我們可以 認為酒精從腸胃吸收入血液的過程是一個勻速變化的過 程 即可假設從吸收室到中心室的轉移速率為常數k0 則 此時 dc1 dt k0 c1 0 0 6 聯合 2 式可解得 c t k0 vk 1 e kt 0 t k0 vk 1 e kt e k t t 7 4 3 每天喝酒情況下血醇濃度的分析 我們假設以時間T為周期喝酒 每周期開始時連續(xù)飲 酒 持續(xù)時間為 則在此情況下 我們根據 7 中的 c t 表達式可以得出當t 時 c t 的穩(wěn)態(tài)濃度為 3 c nT n k0 vk 1 e k 1 e kT c nT n k0 vk ek 1 ekT 1 同時可以求得當 t0 T時 c nT t0 k0 vk 1 ek 1 ekTe k T t0 8 5 模型的檢驗與結果分析 1 短時喝酒模型中 4 式中的3個參數k1 k和D V 即 b 可依據參考資料中喝2瓶酒的試驗數據擬合 我們 用MATLAB對其進行非線性最小二乘法擬合 在確定初始 值時考慮到排出速率k遠小于吸收速率k1 那么k的初值 可設為0 b的初值可認為與實驗數據中c的最大值相當 在此不妨設為85 k1的初值不易確定 可取1或10進行試 驗 在此之前體內不含酒精 c0取值為0 用MATLAB中的優(yōu)化工具箱中的leastsq函數進行擬 合運算 4 最后的擬合結果為 吸收速率k1 2 007 9 排出速率k 0 185 5 b D v 103 860 7 又瓶裝啤酒可設其體積為650mL 酒精濃度為4 乙 醇密度0 8g cm3 則每瓶中乙醇的質量d 20 8g 兩瓶 D 41 6g 于是v D b 0 400 5 則短時喝2瓶酒的血 醇濃度函數 圖1中曲線 2 為 c t 114 432 6 e 0 185 5t e 2 007 9 t 短時喝一瓶酒的情況下可參考以上參數 得出血醇濃度函 數 圖1中曲線1 為 c t 57 216 3 e 0 185 5 t e 2 007 9 t 9 對大李來說 t 6時 c 6 18 799 3 20 因此 凌晨2時被檢查會被定為飲酒駕車 大李出錯的 原因在于忽視了第2次喝酒前 體內仍有部分酒精殘余 會 影響血醇濃度 2 a 短時喝3瓶啤酒 其血醇濃度函數可由 9 得到 為 c t 3 57 216 3 e 0 185 5t e 2 007 9t 75 第2期 楊錦偉 刁 群 飲酒駕車的數學模型 易求得 t 11 588 5時 c t 20 所以 在喝酒后的11 588 5h內駕車會違規(guī) b 若在2h內喝了3瓶啤酒 則依 7 式的參數可取值 如下 2 k0 20 8 3 2 31 2 V 0 400 5 對于參數k仍取0 185 5 則有 c t 419 960 2 1 e 0 185 5t 0 t 2 130 187 7e 0 185 5 t 2 t 2 該函數可在MATLAB中編程實現 當t 12 098 35時 c 20 故在此模型下喝酒后的 12 098 35h內駕車會違規(guī) 3 要知道血醇濃度何時達到最大值 即是求 c t 的最 大值 該問題在4 1 2中已有論述 有公式可知 無論瞬時喝多少酒 體內血醇濃度達到最 高值所用的時間是不變的 如圖1中的曲線1與曲線2 下面引入參數驗證 例如 在短時喝2瓶酒的過程中 tm 1 k k11n k k1 1 0 185 5 2 007 91n 0 185 5 2 007 9 1 307 0 此時cm 81 500 7 對比試驗數據 其最大值在t 1 和t 1 5處取到 為82 可以看出 該模型擬合度較高 比較精確 當參數改變 時 只需代入即可 4 天天喝酒時 根據上面分析 即使只喝一瓶啤酒也需 要5h血醇濃度才能降到20以下 因此對于一個上白班的 司機早上和中午是不宜喝酒的 要想天天喝酒 只能在晚上 下班之后喝適量的酒 并且在第二天早上上班以前血醇濃 度必須降到20以下 這樣 我們來取參數 假設每天晚上8點開始 連續(xù)喝2h的酒 第二天早8 點上班 則為使早8點時的血醇濃度為20 取t0 12 c 20 T 24 2 v 0 400 5 k 0 185 5 利用式 8 可求得k0 30 283 1 即相當于2 911 8瓶 啤酒 也就是說 司機師傅要想天天喝酒 那么 他每天最多 可以在晚上10點以前喝2 911 8瓶啤酒或者含等量酒精的 白酒 6 模型的評價與改進 6 1 模型評價 1 筆者建立了描述酒精在血液中含量變化規(guī)律的微 分方程模型 討論了短時間內喝酒和較長時間內連續(xù)喝酒 2種情況 分別得出了血醇濃度的變化函數 通過此模型可 以較準確的說明血液中酒精濃度的變化規(guī)律 2 所建立的模型比較簡明 可以很方便的使用MAT2 LAB來編程實現其結果 并且采用了準確度較高的非線性 最小二乘進行擬合 曲線吻合度非常高 3 模型具有較好的穩(wěn)定性 可以根據每天飲酒量來計 算可以安全駕車的時間 6 2 模型改進 1 筆者對酒精在血液中的擴散 吸收 排出機理并未 作過多探討 2 只是近似的把全身看作一個分布均勻的機體 事實 上 全身血液分布是不均的 可以分為血液較豐富的中心室 心 肺 腎等 和血液較貧乏的周邊室 四肢 肌肉組織等 可以建立一個兩室模型進一步探討 3 本模型存在近似誤差 系數是通過擬合產生的 而 原始數據組有限 7 經常開汽車 喝酒須適量 自新交通法規(guī)頒布以來 對飲酒駕車的懲罰更加嚴厲 為了減少危害 避免無謂的損失 司機不得不控制自己的飲 酒習慣 但是仍不免有不少困惑 在此提出幾點建議 1 不要以為第一次喝酒沒事就認為每次喝同樣的酒間 隔相同的時間是安全的 要知道經常飲酒體內酒精代謝不 完全 是有殘余 2 經常出車的司機 飲酒頻率應該降低 并且每次的飲 酒量也不應該太多 3 不經常開車的司機 要注意出車前的較長時間內不 要飲酒 至少不能飲太多的酒 參考文獻 1 蕭樹鐵 數學試驗 M 北京 高等教育出版社 1999 2 姜啟源 謝金星 葉 俊 數學模型 M 第3版 北京 高等教育出版社 2003 3 姜啟源 謝金星 葉 俊 數學模型 第3版 習題參考解 答 M 北京 高等教育出版社 2003 4 張智星 MATLAB程序設計與應用 M 北京 清華大 學出版社 2002 Mathematical Model of Drinking and Driving YANGJin wei DIAO Qun Pingdingshan University Pingdingshan Henan467002 China Abstract In this paper we give two differential equations concerning situations in which a certain amount of alcohol drinks are given in a short time vs in a long time by analyzing the
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025版智能電網建設與運營入股合同范本3篇
- 2025年度個人委托代繳社保代理合同樣本3篇
- 二零二五年度地下管線探測與測繪分包合同精準實施范本3篇
- 2025年水泥編織袋市場拓展與品牌戰(zhàn)略合作框架協議3篇
- 2025年度制片人知識產權聘用合同規(guī)范
- 二零二五年度倉儲用地租賃合同簡易范本3篇
- 二零二五年度農行電子商務平臺技術支持與維護合同
- 2025年離婚協議簽訂時效與婚姻解除后續(xù)子女監(jiān)護權協議合同3篇
- 二零二五版廢輪胎膠粉回收及橡膠制品生產合同3篇
- 二零二五年度品牌酒店用品采購合同
- JTG∕T E61-2014 公路路面技術狀況自動化檢測規(guī)程
- 高中英語短語大全(打印版)
- 2024年資格考試-對外漢語教師資格證筆試參考題庫含答案
- 軟件研發(fā)安全管理制度
- 三位數除以兩位數-豎式運算300題
- 寺院消防安全培訓課件
- 比摩阻-管徑-流量計算公式
- GB/T 42430-2023血液、尿液中乙醇、甲醇、正丙醇、丙酮、異丙醇和正丁醇檢驗
- 五年級數學應用題100道
- 西方經濟學(第二版)完整整套課件(馬工程)
- GB/T 33688-2017選煤磁選設備工藝效果評定方法
評論
0/150
提交評論