常微分方程建模方法及案例分析

沈延鋒 姜永慧 沈延琦

摘? 要：常微分方程建模是數學建模中一類十分重要的方法，使用它通常需要建立含多個變量及導數信息的常系數微分方程。本文首先給出了此類建模問題的基本思路、步驟和建模方法，然后通過最速降線、懸鏈線及藥物擴散衰減三個問題對該建模方法進行了分析。分析過程中強調了變量及其變量間關系的確定在常微分方程建立過程中的重要作用。

關鍵詞：微分方程? 微元分析法? 最速降線? 懸鏈線? 藥物擴散衰減

中圖分類號：O175? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A? ? ? ? ? ? ? ? ? ? 文章編號：1674-098X（2020）08（c）-0199-04

Abstract： Ordinary differential equation modeling is a very important method in mathematical modeling. Using it， it is usually necessary to establish constant coefficient differential equations with multiple variables and derivative information. In this paper， we shall discuss the basic idea， steps and several methods about thus modeling problems firstly. Then three practical problems will be studied， as the brachistochrone problems， catenary problems and medicament diffusion problems. The important role of the determination of variables and their relations in the establishment of ordinary differential equations is emphasized.

Key Words： Ordinary differential equations; Microelement analysis method; Brachistochrone; Catenary; Medicament diffusion

函數的本質是兩個變量之間的依賴關系，而對事物變化相互影響的關系研究是工程計算、醫療衛生和金融經濟等眾多領域中的核心問題。通常情況下這些領域中所用到的函數比較復雜，很難用解析表達式直接給出因此用包含多變量（因變量、自變量）的常系數方程是可取的方式之一。有時還需要考慮含未知函數導數的方程，這類方程就是常微分方程。建立適當的常微分方程是數學建模中解決連續性問題的一類十分重要的方法。常微分方程建模的應用領域也非常廣泛，如生物學模型、經濟模型、航空航天模型和物理學模型等等。

常微分方程建模[1-5]是利用常微分方程來近似模擬某些事物或現象隨時間而連續發生變化的方法。所遇問題中提到的“變化”、“增加”、“減少”和“改變”等類似詞語，通常都與導數有關，這也是常微分方程建模時需要特別注意的地方。總體來說，常微分方程建模的思路是：首先尋找問題中的條件得到自變量和因變量及其導數之間的聯系，建立合適的常微分方程;然后利用分析或數值方法求解對應常微分方程;再結合問題本身對結果進行分析;最后對所研究的事物或現象的推演過程給出分析或預測等工作。

一般的常微分方程模型建立的步驟及方法[1-4]總結如下：

（1）分清變量的類型，明確自變量和因變量，理清多個因變量之間的相互聯系。從問題出發，利用“奧卡姆”剃刀，篩選出直接反映所需結論的核心因變量。這里可能會有多個因變量反映結論的不同方面。結合題目所創設的環境，挖掘影響因變量的因素，注意區分確定性因素和隨機性因素。對于隨機性的因素還需要概率統計等基礎知識，甚至還應考慮建立隨機微分方程模型。另外某些因變量之間可能會有明確的（導數關系）關聯，比如位移和速度以及加速度之間的聯系。對于這類關系應該盡早明確，以備后續建立方程時候直接使用。

（2）探索各變量之間的相互關系。通常問題的因變量與自變量之間不是線性關系。對于具有確定性關系的情況，可以參考已有模型進行處理。例如，指數（對數）關系模型、多項式型關系模型和引力模型（在經濟領域中應用也非常多）等等。而對于關系不是很明確的情況，需要建立含兩類變量的函數方程進行處理。部分可以參考人口增長模型、傳染病模型等等，而另外一些還需要分析題設中的不變量。

（3）分析并建立合適的常微分方程。通過上面兩個步驟，已經可以建立起粗糙的常微分方程，而這步通常需要添加一些影響因素對模型進行改進，并得到合理的模型。建立常微分方程的方法多種多樣，大致可以歸結為以下幾種：

①按規律列方程。簡單地說，含自變量和因變量的常系數等式就是常微分方程，而部分等式可以通過數學、化學和物理等學科中的規律獲取。牛頓第二定律就描述了力和加速度的關系。另外放射規律、折射定律等等都可以幫助我們建立較為合適的常微分方程。

②微元分析法。這里所謂的微元分析法不同于定積分應用的微元法，詳見參考文獻[8]。求解微分方程所用的積分幾乎都是變上限的積分，其本質上也是一個函數。因此這里的微元可以理解為微分。我們知道，微分是對應改變量的一種近似。固定自變量的兩點和，分析對應因變量，及其改變量所滿足的等式條件。這里的等式通常由確定的規律或者題設中的不變量得到。然后近似轉化為各變量的微分，就構造出了所需要的常微分方程。例如分析勻速圓周運動的向心加速度，對圓周運動的速度變化進行微元分析，便可以推導出向心加速度的表達式。

③模擬近似法。工程應用領域中的問題通常比較復雜，且受到許多不可控制的因素影響。全面考慮這些因素往往使得微分方程十分復雜甚至無法建立。因此有必要刪繁就簡，添加一些假設條件，用簡單些的微分模型去近似。在Logistic人口阻滯增長模型中，人口的增長率被簡化為只由當前人口總量和其他一些阻滯增長率的因素（如資源環境等）所決定。

1? 常微分方程案例分析

利用常微分方程模型解決實際問題，第一步就是要利用數學理論、初等模型以及其他學科中的定律，建立起較為合適的常微分方程。下面筆者通過幾個典型實例對上述常微分方程建立方法進行闡述。

1.1 最速降線問題

意大利科學家伽利略曾提出這樣的問題：一個質點受重力作用，從點A滑到不在它垂直下方的點B，若不計摩擦，沿什么曲線所需時間最短？取時間t為自變量，并建立坐標系如圖1所示，位置、速度和P點切線傾角都是關于t的因變量。另設質點的質量為m，初速度為。因為不計摩擦，速度滿足? ? ? ? ? ?。

由費馬原理可知，光沿著遵守折射定律的路線傳播，所用的時間最短。將AB曲線看作光在非均勻介質中的傳播路徑，要使得時間最短，每點都應滿足折射定律，即

1.2 懸鏈線問題

懸鏈線是指受重力作用，兩端水平固定的柔軟的繩子自由下垂所形成的曲線[6]。如圖2所示，O點為繩子的最低點，此點受到一個水平向左的拉力。取橫坐標為自變量，懸鏈線上點的高度、拉力、拉力與水平方向夾角和累計重力均為因變量，分別記為、、和。任取繩子上極小的一段，且設和 的橫坐標分別為和。另設繩子的線密度為，下面對繩子上的一小段進行受力分析。

1.3 藥物擴散衰減問題

在使用藥物治療疾病時，通常有三種給藥方式，分別為快速靜脈注射（方式a）、恒速靜脈注射（方式b）和口服或肌肉注射（方式c）[7]。下面通過模擬近似分析血藥濃度隨時間的變化趨勢。為此，我們將身體看成一個房室，給藥方式不同，藥物的吸收和消除速率不同。

取時間為自變量，設為表觀分布容積，為? 時刻體內的藥物量，為時刻體內的藥物濃度，則滿足關系體內的藥物濃度，則滿足關系。令 為一級消除速率。

下面我們將建立描述三種給藥方式的常微分方程。利用微元分析法，任取小區間，設藥物的改變量為，其值可以利用不同的吸收和消除速率計算：

（1） 方式a。因吸收和消除速率分別為0和，故改變量，從而借助函數微分即可得到常微分方程為，且有初值條件。為注射總藥量。解得血藥濃度變化曲線為，近似為指數分布的密度函數。

（2）方式b。吸收速率為定值a，得改變量，可得模型為，且。解得濃度變化曲線為。近似為Logistic人口阻滯函數曲線。

（3）方式c。此時改變量計算稍微復雜，因為吸收速率也是的函數。先考慮吸收子模型。設藥物被吸收的速率與存量藥物的數量成正比，記比例系數為，若記時刻殘留藥物量為，類似地得到常微分方程為，，其中為方式c的藥物總量。解得。進而建立方式c藥量的常微分方程，。解得血藥濃度變化曲線為

近似為卡方分布的密度函數。

（1） 方式a。因吸收和消除速率分別為0和，故改變量，從而借助函數微分即可得到常微分方程為，且有初值條件。為注射總藥量。解得血藥濃度變化曲線為，近似為指數分布的密度函數。

（2）方式b。吸收速率為定值a，得改變量，可得模型為，且。解得濃度變化曲線為。近似為Logistic人口阻滯函數曲線。

近似為卡方分布的密度函數。

圖3給出了三種給藥方式的血藥濃度變化曲線。進一步分析可以根據血藥濃度要求控制再次給藥的時間和劑量，以便達到維持藥效的作用。

2? 結語

常微分方程建模在利用數學模型解決醫療、經濟、工程等領域問題中占據著很大的優勢。要建立適合的微分方程模型，就要從題設、目標中仔細挖掘各變量及其可能存在的關系，還要考慮相應領域中的定律、因變量微分的近似計算特性和微分方程教材中的經典模型，綜合分析構建出較為適合的模型。通常還需要把求得的解與實際情況進行比較，不斷改進和優化模型，盡可能地讓最終模型既能滿足實際問題的精確性要求，同時也具有一定的穩定性要求。

參考文獻

[1] 廖莉.常微分方程在數學建模中的應用[J].課程教育研究，2017（38）：139-140.

[2] 陳尹剛.數學建模在常微分方程中的應用[J].湖北開放職業學院學報，2019（6）：103-104.

[3] 鄭玉敏.對常微分方程數學建模案例的研究[J].黑龍江科學，2017（13）：66-67.

[4] 李明偉.數學建模思想融入常微分方程教學的探討[J].高教學刊，2018（1）：93-95.

[5] 武海輝.一類競爭擴散系統的穩定性分析[J].河南科學，2019，37（6）：869-873.

[6] 周義倉.懸鏈線模型在系泊系統設計中的應用——2016年全國大學生數學建模競賽A題解答評述[J].數學建模及其應用，2016，5（4）：26-33.

[7] 王蓮招.三種給藥方式的微分方程及函數性態比較 探討[J].高教學刊，2016（5）：111-112.