數(shù)學(xué)模型微分方程_第1頁(yè)
數(shù)學(xué)模型微分方程_第2頁(yè)
數(shù)學(xué)模型微分方程_第3頁(yè)
數(shù)學(xué)模型微分方程_第4頁(yè)
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

PAGEPAGE2數(shù)學(xué)模型13.人體注射葡萄糖溶液時(shí),血液中葡萄糖濃度g(t)的增長(zhǎng)率與注射速率r成正比,與人體血液容積v成反比,而由于人體組織的吸收作用,g(t)的減少率與g(t)本身成正比。分別在以下幾種假設(shè)下建立模型,并討論穩(wěn)定情況。(1)人體血液容積v不變。(2)v隨著注入溶液而增加。(3)由于排泄等因素v的增加有極限值解:模型假設(shè):本模型中主要符號(hào)說(shuō)明為:葡萄糖濃度g(t)注射速率r人體血液容積v基本模型為:(,,常數(shù))⑴(1)V為常數(shù)時(shí),平衡點(diǎn)穩(wěn)定。如果以g為橫軸、為縱軸作出方程的圖形(圖1),可以分析葡萄糖濃度增長(zhǎng)速度隨著g的增加而變化的情況,從而大概地看出g(t)的變化規(guī)律。令,,利用Mathematica在操作窗口中輸入以下代碼命令:Plot[0.2/100-0.5g,{g,0,100},PlotStyle->{RGBColor[1,0,0]}]得到:圖1~g曲線再利用matlab在操作窗口中輸入以下代碼命令:g=dsolve('Dg=k1*r/v-k2*g','g(0)=g0','t')其解為g=k1*r/v/k2+exp(-k2*t)*(-k1*r+g0*v*k2)/v/k2整理得到:⑵令,,利用Mathematica在操作窗口中輸入以下代碼命令:Plot[0.2/100+Exp[-0.5t],{t,0,100},PlotStyle->{RGBColor[1,0,0]}]得到:圖2g~t曲線由圖可以知道它在平衡點(diǎn)穩(wěn)定。(2)不妨設(shè)(,常數(shù))⑶方程⑴,⑵不存在平衡點(diǎn)。若由⑵解出代入⑴,得到⑷則⑷不能是自治方程。因?yàn)槠胶恻c(diǎn)及穩(wěn)定性的概念只是對(duì)自治方程而言才有意義,而⑷不能是自治方程,所以不能考慮它的穩(wěn)定性。(3)不妨設(shè)(,常數(shù))⑸如果以V為橫軸、為縱軸作出方程的圖形(圖3),可以分析人體血液容積V增長(zhǎng)速度隨著V的增加而變化的情況,從而大概地看出V(t)的變化規(guī)律。令,,利用Mathematica在操作窗口中輸入以下代碼命令:Plot[0.5(2000-v),{v,0,2000},PlotStyle->{RGBColor[0,1,0]}]得到:圖3~V曲線再利用matlab在操作窗口中輸入以下代碼命令:v=dsolve('Dv=u*(v1-v)','v(0)=v0','t')得到:v=v1+exp(-u*t)*(-v1+v0)整理得到:⑹令,,,利用Mathematica在操作窗口中輸入以下代碼命令:Plot[2000-Exp[-0.5t](2000-1000),{t,0,100},PlotStyle->{RGBColor[0,1,0]}]得到:圖4V~t曲線由圖可以知道它在平衡點(diǎn)穩(wěn)定.方程

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論