數(shù)學(xué)建模-微分方程模型課件

1、微分方程模型常微分方程的基本方法微分方程基礎(chǔ)微分方程是含有函數(shù)及其導(dǎo)數(shù)的方程。如果方程(組)只含有一個(gè)自變量(通常是時(shí)間t)，則稱為常微分方程。否則稱為偏微分方程。例：下面的方程都是微分方程：微分方程的解是函數(shù)，對(duì)應(yīng)一個(gè)變化過程。常微分方程的解是隨時(shí)間t變化的函數(shù)，比如一輛汽車在公路上飛馳，一個(gè)球從空中落下等。偏微分方程不但描述物體隨時(shí)間變化發(fā)生位置的改變，而且物體各部分之間的位置的相對(duì)變化。如水的流動(dòng)，煙霧的擴(kuò)散，公路上車流的涌動(dòng)等。微分方程解決的主要問題： (1)描述對(duì)象特征隨時(shí)間(空間)的演變過程 (2)分析對(duì)象特征的變化規(guī)律 (3)預(yù)報(bào)對(duì)象特征的未來性態(tài) (4)研究控制對(duì)象特征的手段導(dǎo)

2、數(shù)的意義：瞬時(shí)變化率在實(shí)際上我們遇到的描述變化的詞有 速率(物理) 增長率(經(jīng)濟(jì)，生物，人口等) 衰變(原子反應(yīng)) 邊際的(經(jīng)濟(jì))瞬時(shí)變化率的描述： 絕對(duì)增加率：單位時(shí)間增加的量。 相對(duì)增加率：單位時(shí)間增加的百分比。 變化率= 增加率-減少率由于是瞬時(shí)的，其量的關(guān)系只有在很短的時(shí)間間隔中才能夠利用靜態(tài)的方法分析。(微元法)微分方程的建模方法： (1)利用導(dǎo)數(shù)的意義，建立含有導(dǎo)數(shù)的方程(微分方程)。 (2)微元法。微分方程的穩(wěn)定性理論：對(duì)微分方程組若f(x0)=0，則稱x0是方程組的平衡點(diǎn)。穩(wěn)定的平衡點(diǎn)的實(shí)際意義：如果微分方程存在穩(wěn)定的平衡點(diǎn)，設(shè)x(t)是微分方程的解，則當(dāng)t時(shí)， x(t)趨向于

3、某個(gè)穩(wěn)定的平衡點(diǎn)。例：對(duì)Logistic方程，它有兩個(gè)平衡點(diǎn) x=0和x=N。其中x=0是不穩(wěn)定的平衡點(diǎn)，x=N是穩(wěn)定的平衡點(diǎn)。 例1：某人的食量是2500卡/天。其中1200卡用于基本的新陳代謝。在健身訓(xùn)練中，他每公斤體重所消耗的熱量大約是16卡/天。設(shè)以脂肪形式貯存的熱量100%有效，且1公斤脂肪含熱量10000卡，分析這個(gè)人體重的變化。分析：問題研究人體重量隨時(shí)間的變化w(t)。條件給出的是熱量單位時(shí)間的變化 2500-1200-16w(t)轉(zhuǎn)換成體重為 (2500-1200-16w(t)/10000因此得到變化關(guān)系 常微分方程建模的物理方法熱傳導(dǎo)：牛頓冷卻定律(加熱定律)：例：將一只讀

4、數(shù)為25度的溫度計(jì)放在室外，10分鐘后度數(shù)為30度，又過了10分鐘，讀數(shù)變?yōu)?3度，問室外溫度是多少？如果遇到我們不熟悉的問題時(shí)，應(yīng)該怎么辦？答案：不要回避，到網(wǎng)上查一下相關(guān)的概念你就會(huì)發(fā)現(xiàn)：這個(gè)不熟悉的問題可能是比較簡單的!分析：上網(wǎng)查一下熱傳導(dǎo)，我們可以了解到：熱的傳導(dǎo)從溫度高的地方向溫度低的地方傳導(dǎo)，單位時(shí)間傳送的熱量與溫差T成正比，與兩個(gè)熱源的距離成反比。即對(duì)于兩個(gè)固定熱源，距離d是常數(shù)，則 在我們的問題中，室外溫度可以看做常數(shù)T0，大于室內(nèi)溫度，而熱量正比于溫差，從而變化規(guī)律為問題：現(xiàn)有4000毫升溫度為10度的化學(xué)溶液，將一個(gè)體積40毫升溫度為90度的玻璃球放在溶液中。求溶液溫度的

5、變化規(guī)律。(平均溫度)模型的解為這里有三個(gè)參數(shù)，其中T0=25。還剩兩個(gè)參數(shù)，利用剩下的兩個(gè)條件可以確定。動(dòng)力學(xué)：牛頓第二定律能量守恒定律歐拉-拉格朗日方程空氣和水的阻力例1：求單擺的運(yùn)動(dòng)：擺長L，擺錘質(zhì)量m的單擺的運(yùn)動(dòng)方程例2：一只裝滿水的圓柱形桶，底半徑3m，高6m。底部有一個(gè)直徑0.02米的孔。(1)水多長時(shí)間可以流光？(2)如果孔在側(cè)面，而桶放在距地面3m 的高度。求水流噴出距離的變化規(guī)律。解：直接利用Newton第二定律建模比較困難，我們利用能量的轉(zhuǎn)換。在流水的過程中，桶的頂部減少的勢能化為水的動(dòng)能。(如圖)hhdsdhhds設(shè)桶的水平面積為A，孔的面積為B，則由于質(zhì)量守恒，則 Ad

6、h=-Bds符號(hào)反映了此消彼長。設(shè)水的流速是v則 根據(jù)能量轉(zhuǎn)換關(guān)系，水失去的勢能轉(zhuǎn)化為動(dòng)能，即或 練習(xí)題：1、在一所大學(xué)，某個(gè)教師每天從圖書館借出一本書，而圖書館每周收回所借圖書的10%。幾年后，這個(gè)教師手中有大約多少本圖書館的書？2、某學(xué)院的教育基金，最初投資P元，以后按利率r的連續(xù)復(fù)利增長。另外，每年在基金開算的時(shí)間，都要投入新的資本A/年求7年的累計(jì)資金數(shù)量。另外，如果每年在基金開算的時(shí)間，把其中20%用于獎(jiǎng)學(xué)金的發(fā)放，求7年后累計(jì)資金數(shù)量。3、一場降雪開始于中午前的某個(gè)時(shí)刻，降雪量穩(wěn)定。某人從正午12點(diǎn)開始清掃人行道，他的鏟雪速度(m3/小時(shí))和路面寬度都不變，到下午2點(diǎn)他掃了1000

7、米，到下午4點(diǎn)又清掃了500米。雪是什么時(shí)間開始下的？另外，如果他在下午4點(diǎn)開始回頭清掃，什么時(shí)間回到開始清掃的地點(diǎn)？2004C題 飲酒駕車據(jù)報(bào)載，2003年全國道路交通事故死亡人數(shù)為10.4372萬，其中因飲酒駕車造成的占有相當(dāng)?shù)谋壤?。針?duì)這種嚴(yán)重的道路交通情況，國家質(zhì)量監(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

8、點(diǎn)喝了一瓶啤酒，下午6點(diǎn)檢查時(shí)符合新的駕車標(biāo)準(zhǔn)，緊接著他在吃晚飯時(shí)又喝了一瓶啤酒，為了保險(xiǎn)起見他呆到凌晨2點(diǎn)才駕車回家，又一次遭遇檢查時(shí)卻被定為飲酒駕車，這讓他既懊惱又困惑，為什么喝同樣多的酒，兩次檢查結(jié)果會(huì)不一樣呢？請(qǐng)你參考下面給出的數(shù)據(jù)（或自己收集資料）建立飲酒后血液中酒精含量的數(shù)學(xué)模型，并討論以下問題：1. 對(duì)大李碰到的情況做出解釋；2. 在喝了3瓶啤酒或者半斤低度白酒后多長時(shí)間內(nèi)駕車就會(huì)違反上述標(biāo)準(zhǔn)，在以下情況下回答：酒是在很短時(shí)間內(nèi)喝的；酒是在較長一段時(shí)間（比如2小時(shí)）內(nèi)喝的。3. 怎樣估計(jì)血液中的酒精含量在什么時(shí)間最高。4. 根據(jù)你的模型論證：如果天天喝酒，是否還能開車？5. 根據(jù)

9、你做的模型并結(jié)合新的國家標(biāo)準(zhǔn)寫一篇短文，給想喝一點(diǎn)酒的司機(jī)如何駕車提出忠告。時(shí)間(小時(shí))0.250.50.7511.522.533.544.55酒精含量306875828277686858515041時(shí)間(小時(shí))678910111213141516酒精含量3835282518151210774參考數(shù)據(jù)1. 人的體液占人的體重的65%至70%，其中血液只占體重的7%左右；而藥物（包括酒精）在血液中的含量與在體液中的含量大體是一樣的。2. 體重約70kg的某人在短時(shí)間內(nèi)喝下2瓶啤酒后，隔一定時(shí)間測量他的血液中酒精含量（毫克百毫升），得到數(shù)據(jù)如下：在這個(gè)問題中，有兩個(gè)問題：(1)酒精在身體中吸收和排

10、除過程的描述。(2)利用所給出的數(shù)據(jù)，分析大李的具體例子。這是某人喝下兩瓶啤酒后，測得的血液中酒精含量數(shù)據(jù)我們可以把喝酒后酒精的變化過程描述為喝酒酒精進(jìn)入腸胃消化后進(jìn)入血液排出。這里，血液循環(huán)系統(tǒng)可以看作中心室，腸胃可以看作吸收室。M1克酒精在很短時(shí)間進(jìn)入吸收室，從吸收室逐漸進(jìn)入中心室，最后逐漸排出。吸收室 x中心室 y排出k1k2 m1吸收室 x中心室 y排出k1k2 m1設(shè)從吸收室到中心室的酒精轉(zhuǎn)移速率為k1，中心室中的酒精排出速率為k2，則2011 C題 企業(yè)退休職工養(yǎng)老金制度的改革養(yǎng)老金也稱退休金，是一種根據(jù)勞動(dòng)者對(duì)社會(huì)所作貢獻(xiàn)及其所具備享受養(yǎng)老保險(xiǎn)的資格，以貨幣形式支付的保險(xiǎn)待遇，用

11、于保障職工退休后的基本生活需要。我國企業(yè)職工基本養(yǎng)老保險(xiǎn)實(shí)行“社會(huì)統(tǒng)籌”與“個(gè)人賬戶”相結(jié)合的模式，即企業(yè)把職工工資總額按一定比例（20%）繳納到社會(huì)統(tǒng)籌基金賬戶，再把職工個(gè)人工資按一定比例（8%）繳納到個(gè)人賬戶。這兩個(gè)賬戶我們合稱為養(yǎng)老保險(xiǎn)基金。退休后，按職工在職期間每月（或年）的繳費(fèi)工資與社會(huì)平均工資之比（繳費(fèi)指數(shù)），再考慮到退休前一年的社會(huì)平均工資等因素，從社會(huì)統(tǒng)籌賬戶中撥出資金（基礎(chǔ)養(yǎng)老金），加上個(gè)人工資賬戶中一定比例的資金（個(gè)人賬戶養(yǎng)老金），作為退休后每個(gè)月的養(yǎng)老金。養(yǎng)老金會(huì)隨著社會(huì)平均工資的調(diào)整而調(diào)整。如果職工死亡，社會(huì)統(tǒng)籌賬戶中的資金不退給職工，個(gè)人賬戶中的余額可繼承。個(gè)人賬戶儲(chǔ)

12、存額以銀行當(dāng)時(shí)公布的一年期存款利率計(jì)息，為簡單起見，利率統(tǒng)一設(shè)定為3%。養(yǎng)老金的發(fā)放與職工在職時(shí)的工資及社會(huì)平均工資有著密切關(guān)系；工資的增長又與經(jīng)濟(jì)增長相關(guān)。近30年來我國經(jīng)濟(jì)發(fā)展迅速，工資增長率也較高；而發(fā)達(dá)國家的經(jīng)濟(jì)和工資增長率都較低。我國經(jīng)濟(jì)發(fā)展的戰(zhàn)略目標(biāo)，是要在21世紀(jì)中葉使我國人均國民生產(chǎn)總值達(dá)到中等發(fā)達(dá)國家水平?，F(xiàn)在我國養(yǎng)老保險(xiǎn)改革正處于過渡期。養(yǎng)老保險(xiǎn)管理的一個(gè)重要的目標(biāo)是養(yǎng)老保險(xiǎn)基金的收支平衡，它關(guān)系到社會(huì)穩(wěn)定和老齡化社會(huì)的順利過渡。影響?zhàn)B老保險(xiǎn)基金收支平衡的一個(gè)重要因素是替代率。替代率是指職工剛退休時(shí)的養(yǎng)老金占退休前工資的比例。按照國家對(duì)基本養(yǎng)老保險(xiǎn)制度的總體思路，未來基本養(yǎng)

13、老保險(xiǎn)的目標(biāo)替代率確定為58.5%. 替代率較低，退休職工的生活水準(zhǔn)低，養(yǎng)老保險(xiǎn)基金收支平衡容易維持；替代率較高，退休職工的生活水準(zhǔn)就高，養(yǎng)老保險(xiǎn)基金收支平衡較難維持，可能出現(xiàn)缺口。所謂缺口，是指當(dāng)養(yǎng)老保險(xiǎn)基金入不敷出時(shí)出現(xiàn)的收支之差。附件1是山東省職工歷年平均工資數(shù)據(jù)；附件2是2009年山東省某企業(yè)各年齡段職工的工資分布情況，附件3是養(yǎng)老金的計(jì)算辦法。請(qǐng)建立數(shù)學(xué)模型，解決如下問題：問題一：對(duì)未來中國經(jīng)濟(jì)發(fā)展和工資增長的形勢做出你認(rèn)為是簡化、合理的假設(shè)，并參考附件1，預(yù)測從2011年至2035年的山東省職工的年平均工資。問題二：根據(jù)附件2計(jì)算2009年該企業(yè)各年齡段職工工資與該企業(yè)平均工資之比

14、。如果把這些比值看作職工繳費(fèi)指數(shù)的參考值，考慮該企業(yè)職工自2000年起分別從30歲、40歲開始繳養(yǎng)老保險(xiǎn)，一直繳費(fèi)到退休（55歲，60歲，65歲），計(jì)算各種情況下的養(yǎng)老金替代率。問題三：假設(shè)該企業(yè)某職工自 2000年起從30歲開始繳養(yǎng)老保險(xiǎn)，一直繳費(fèi)到退休（55歲，60歲，65歲），并從退休后一直領(lǐng)取養(yǎng)老金，至75歲死亡。計(jì)算養(yǎng)老保險(xiǎn)基金的缺口情況，并計(jì)算該職工領(lǐng)取養(yǎng)老金到多少歲時(shí)，其繳存的養(yǎng)老保險(xiǎn)基金與其領(lǐng)取的養(yǎng)老金之間達(dá)到收支平衡。問題一：對(duì)未來中國經(jīng)濟(jì)發(fā)展和工資增長的形勢做出你認(rèn)為是簡化、合理的假設(shè)，并參考附件1，預(yù)測從2011年至2035年的山東省職工的年平均工資。下面的圖是題目所給出

15、的職工工資隨時(shí)間增長數(shù)據(jù)。許多同學(xué)分析這個(gè)問題的一個(gè)誤區(qū)是把這組數(shù)據(jù)看作無法分析機(jī)理的問題，即只能見黑箱模型。因此，一眼看上去像一個(gè)3次多項(xiàng)式，馬上利用三次多項(xiàng)式擬合，并利用這個(gè)多項(xiàng)式做預(yù)測得到下圖人均達(dá)到五十五萬！先不談?lì)A(yù)測的正確與否，我們來分析問題是否沒有機(jī)理可以在建模中起作用。首先讀一下題目中下面段落：養(yǎng)老金的發(fā)放與職工在職時(shí)的工資及社會(huì)平均工資有著密切關(guān)系；工資的增長又與經(jīng)濟(jì)增長相關(guān)。近30年來我國經(jīng)濟(jì)發(fā)展迅速，工資增長率也較高；而發(fā)達(dá)國家的經(jīng)濟(jì)和工資增長率都較低。我國經(jīng)濟(jì)發(fā)展的戰(zhàn)略目標(biāo)，是要在21世紀(jì)中葉使我國人均國民生產(chǎn)總值達(dá)到中等發(fā)達(dá)國家水平。這一段告訴我們：(1)職工工資逐年增

16、長；(2)發(fā)達(dá)國家的經(jīng)濟(jì)和工資增長率都較低，而要在21世紀(jì)中葉使我國人均國民生產(chǎn)總值達(dá)到中等發(fā)達(dá)國家水平，因此工資增長率要逐年降低，到2050年到較低的水平。利用數(shù)學(xué)模型描述，就是：其中r(x)是單調(diào)減函數(shù)，且r(2050)很小。由于x(1978)=566，是一個(gè)很小的數(shù)，我們可以作變換y=x-1978. 一個(gè)最直接的模型是Logistic模型其精確解為 利用最小二乘法可以求出參數(shù)r,y0和N。但計(jì)算得到的y(50)=200萬擬合結(jié)果(好像不錯(cuò))預(yù)測結(jié)果：到2050年人均650萬！(合理嗎？)模型的修改：修改后的模型與數(shù)據(jù)產(chǎn)生偏移。所以，模型和數(shù)據(jù)擬合的不太好。是不是可以改進(jìn)模型呢？可以考慮更

17、一般的模型：問題：這個(gè)方程的解比較復(fù)雜，難以寫成解析形式。那么，如何對(duì)解作最小二乘擬合呢？ 也可以預(yù)先確定N的值。由于我國在2050年達(dá)到中等發(fā)達(dá)國家的工資水平。中等發(fā)達(dá)國家工資水平在2050年應(yīng)在15-18萬元。(可以利用網(wǎng)上查出的數(shù)據(jù)加上通貨膨脹因素計(jì)算得到)取N=150000,即可得到合理的數(shù)據(jù)。練習(xí)題：2003D題 搶渡長江“渡江”是武漢城市的一張名片。1934年9月9日，武漢警備旅官兵與體育界人士聯(lián)手，在武漢第一次舉辦橫渡長江游泳競賽活動(dòng)，起點(diǎn)為武昌漢陽門碼頭，終點(diǎn)設(shè)在漢口三北碼頭，全程約5000米。有44人參加橫渡，40人達(dá)到終點(diǎn)，張學(xué)良將軍特意向冠軍獲得者贈(zèng)送了一塊銀盾，上書“力

18、挽狂瀾”。2001年，“武漢搶渡長江挑戰(zhàn)賽”重現(xiàn)江城。2002年，正式命名為“武漢國際搶渡長江挑戰(zhàn)賽”，于每年的5月1日進(jìn)行。由于水情、水性的不可預(yù)測性，這種競賽更富有挑戰(zhàn)性和觀賞性。2002年5月1日，搶渡的起點(diǎn)設(shè)在武昌漢陽門碼頭，終點(diǎn)設(shè)在漢陽南岸咀，江面寬約1160米。據(jù)報(bào)載，當(dāng)日的平均水溫16.8, 江水的平均流速為1.89米/秒。參賽的國內(nèi)外選手共186人（其中專業(yè)人員將近一半），僅34人到達(dá)終點(diǎn)，第一名的成績?yōu)?4分8秒。除了氣象條件外，大部分選手由于路線選擇錯(cuò)誤，被滾滾的江水沖到下游，而未能準(zhǔn)確到達(dá)終點(diǎn)。假設(shè)在競渡區(qū)域兩岸為平行直線, 它們之間的垂直距離為 1160 米, 從武昌漢

19、陽門的正對(duì)岸到漢陽南岸咀的距離為 1000米，見示意圖。請(qǐng)你們通過數(shù)學(xué)建模來分析上述情況, 并回答以下問題: 1. 假定在競渡過程中游泳者的速度大小和方向不變，且競渡區(qū)域每點(diǎn)的流速均為 1.89 米/秒。試說明2002年第一名是沿著怎樣的路線前進(jìn)的，求她游泳速度的大小和方向。如何根據(jù)游泳者自己的速度選擇游泳方向，試為一個(gè)速度能保持在1.5米/秒的人選擇游泳方向，并估計(jì)他的成績。(2)在（1）的假設(shè)下，如果游泳者始終以和岸邊垂直的方向游, 他(她)們能否到達(dá)終點(diǎn)？根據(jù)你們的數(shù)學(xué)模型說明為什么 1934年 和2002年能游到終點(diǎn)的人數(shù)的百分比有如此大的差別；給出能夠成功到達(dá)終點(diǎn)的選手的條件。(3)

20、若流速沿離岸邊距離的分布為 (設(shè)從武昌漢陽門垂直向上為 y軸正向) ： 游泳者的速度大?。?.5米/秒）仍全程保持不變，試為他選擇游泳方向和路線，估計(jì)他的成績。4. 若流速沿離岸邊距離為連續(xù)分布, 例如 或你們認(rèn)為合適的連續(xù)分布，如何處理這個(gè)問題。5.用普通人能懂的語言，給有意參加競渡的游泳愛好者寫一份競渡策略的短文。6.你們的模型還可能有什么其他的應(yīng)用？ 偏微分方程的建模含有多個(gè)自變量的微分方程稱為偏微分方程偏微分方程的自變量一般有時(shí)間變量t 和位形變量x。如果運(yùn)動(dòng)過程只與時(shí)間有關(guān)，物體中各點(diǎn)的相對(duì)位置不變，通?？梢岳贸Ｎ⒎址匠探M描述。如質(zhì)點(diǎn)組或剛體的運(yùn)動(dòng)等。如果運(yùn)動(dòng)的物質(zhì)的相對(duì)空間還發(fā)生

21、變化，則需要利用偏微分方程(組)描述。如水的流動(dòng)，氣體的擴(kuò)散和彈性體的振動(dòng)等。偏微分方程的建模有以下特點(diǎn):(1) 由于自變量是高維空間，邊界條件比較復(fù)雜；(2) 高維空間的變化描述的微元不是簡單的區(qū)間，而是任意的小區(qū)域，處理需要曲面積分和多重積分的技巧。流動(dòng)現(xiàn)象的建模流動(dòng)是一類普遍的自然現(xiàn)象。江河里的水流，公路上的車流，大廈中的人流，企業(yè)中的資金流和物流等，都是流動(dòng)現(xiàn)象。流動(dòng)是物質(zhì)的轉(zhuǎn)移的基本形式之一。流動(dòng)是由于力的作用造成的。如水流從高處向低處流動(dòng)，桶中的水由于壓力的作用而發(fā)生流速的變化等。流動(dòng)的一個(gè)基本特點(diǎn)是物質(zhì)不滅，即在一個(gè)封閉的區(qū)域中，流入的量減去流出的量等于區(qū)域中物質(zhì)的增量。流動(dòng)的描

22、述：利用什么樣的量描述流動(dòng)呢？(我們通過什么樣的量的變化能夠刻畫出流動(dòng)的過程呢？)我們可以從河流中的污染物的隨波逐流可以看出，刻畫流動(dòng)有以下基本的量：(1)物質(zhì)的分布。(由于是連續(xù)分布，可以利用分布函數(shù)或密度函數(shù)刻畫分布)流動(dòng)的刻畫實(shí)際上是由物質(zhì)分布的改變而顯示出來的。(2)流速。在一點(diǎn)附近的流動(dòng)在時(shí)間很短時(shí)，可以近似看作小的固體沿某方向流動(dòng)，因此流速是利用向量描述的。(3)流量。(單位時(shí)間流過的物質(zhì)數(shù)量)流量是標(biāo)量。例2：公路上的車流考慮一段公路上的汽車流動(dòng)。設(shè)這段公路的長度與車輛之間的距離相比長得多。設(shè)在t 時(shí)刻，x位置 車流密度函數(shù)為p(x,t) 車的流量為q(x,t)。(單位時(shí)間流過的

23、物質(zhì))在x,x+x段，t,t+t時(shí)段車輛總數(shù)的改變等于進(jìn)入這一區(qū)域的車輛數(shù)減去離開這一區(qū)域的車輛數(shù)，即 可以將上述方程化簡：取極限得到：問題：一個(gè)方程，兩個(gè)未知函數(shù)，解無法確定，如何處理？解決的方法有兩個(gè)： (1)尋找另一個(gè)獨(dú)立的方程； (2)尋找兩個(gè)變量之間的關(guān)系。我們從第2個(gè)問題入手。直觀上可以看出，當(dāng)密度很小的時(shí)候，速度較大，而密度達(dá)到極限時(shí)，速度很小(可看作速度等于零，堵車了).因此，設(shè)pm為極限密度，則密度和速度的關(guān)系模型為 最簡單的是線性化模型 從而得到流量和密度的關(guān)系式把這個(gè)方程代入原方程得到模型的結(jié)構(gòu)就確定了。但模型還沒有最后確定，還必須給出初邊值條件。 問題：如果已知v=60

24、m/h，pm=280輛/mile,并設(shè)開始時(shí)刻公路上均勻分布著150輛車，車輛以30m/h進(jìn)入公路。計(jì)算t=3時(shí)的車輛分布。首先考慮煙草內(nèi)毒素的流動(dòng)。要分析物質(zhì)的流動(dòng)，首先要分析兩個(gè)量在t時(shí)刻x位置 毒素的密度 w(x,t) 毒素的流量 q(x,t)考慮在t,t+t時(shí)段，x,x+x上毒素?cái)?shù)量的變化。由質(zhì)量守恒得到例3：煙卷里毒素的流動(dòng)和沉積。類似上一例子得到描述流動(dòng)的方程下面分析流量和密度的關(guān)系。煙草中的毒素是由于本來含有的毒素加上由于煙草中的毒素隨煙霧流動(dòng)，經(jīng)過x點(diǎn)被煙草吸附的毒素。而煙霧中毒素的含量與燃燒處毒素的密度有關(guān)?？梢约僭O(shè)： 毒素的流量與燃燒處燃燒掉的煙草所含有的毒素成比例。比例系

25、數(shù)。 煙霧經(jīng)過x點(diǎn)被煙草吸附的毒素與經(jīng)過x點(diǎn)的煙霧所含有的毒素量成比例，比例系數(shù)。 進(jìn)一步分析。設(shè)在x0處在燃燒的煙草毒物沉積在xx0處的數(shù)量。由于沉積率為，則這些毒素隨著煙霧的流動(dòng)而減少，變化規(guī)律為其中w1是煙霧中的毒物密度。設(shè)煙霧的運(yùn)行速度為u，則從x0到x的時(shí)間為 得到下面分析在t 時(shí)刻x0的位置。注意到x0的位置與煙草的燃燒速度有關(guān)，如果燃燒速度是常數(shù)v，則x0=vt。綜合得到 上面的式子建立了密度w和流量q的關(guān)系。問題在于兩者不同步。我們可以設(shè)計(jì)算法求解。例4：人口分布的動(dòng)力學(xué)方程從廣義上看，人口的發(fā)展也是一個(gè)流動(dòng)過程，即隨著時(shí)間的發(fā)展，從低齡向高齡變化。差別在于：由于死亡的發(fā)生，人

26、數(shù)是不守恒的。設(shè)t時(shí)刻r齡人口的密度為f(t,r)，經(jīng)過時(shí)間dt，r,r+dr區(qū)段的人數(shù)減少f(t,r)(r)dr而年齡變化為r+dt,r+dr+dt。其中(r)是r齡人口的死亡率。從而得到令dt,dr0得到高維流動(dòng)一維流動(dòng)的特點(diǎn)是流動(dòng)方向唯一，容易考慮。下面考慮三維流動(dòng)。首先引入 流量 q(x,y,z,t) (由于流的方向可以變化，流量成為向量)。 密度 w(x,y,z,t)流動(dòng)的特點(diǎn)是：從區(qū)域邊界流出的物質(zhì)的總數(shù)等于內(nèi)部物質(zhì)數(shù)量的變化。考慮一個(gè)區(qū)域上物質(zhì)在t,t+t時(shí)段的流動(dòng)。流出的流量總和為這一段時(shí)間內(nèi)物質(zhì)總量的變化兩者是相等的。令t0得到什么？練習(xí)：對(duì)二維問題，導(dǎo)出流動(dòng)的質(zhì)量守恒方程。

27、2011 A題 城市表層土壤重金屬污染分析隨著城市經(jīng)濟(jì)的快速發(fā)展和城市人口的不斷增加，人類活動(dòng)對(duì)城市環(huán)境質(zhì)量的影響日顯突出。對(duì)城市土壤地質(zhì)環(huán)境異常的查證，以及如何應(yīng)用查證獲得的海量數(shù)據(jù)資料開展城市環(huán)境質(zhì)量評(píng)價(jià)，研究人類活動(dòng)影響下城市地質(zhì)環(huán)境的演變模式，日益成為人們關(guān)注的焦點(diǎn)。按照功能劃分，城區(qū)一般可分為生活區(qū)、工業(yè)區(qū)、山區(qū)、主干道路區(qū)及公園綠地區(qū)等，分別記為1類區(qū)、2類區(qū)、5類區(qū)，不同的區(qū)域環(huán)境受人類活動(dòng)影響的程度不同?，F(xiàn)對(duì)某城市城區(qū)土壤地質(zhì)環(huán)境進(jìn)行調(diào)查。為此，將所考察的城區(qū)劃分為間距1公里左右的網(wǎng)格子區(qū)域，按照每平方公里1個(gè)采樣點(diǎn)對(duì)表層土（010 厘米深度）進(jìn)行取樣、編號(hào)，并用GPS記錄采樣

28、點(diǎn)的位置。應(yīng)用專門儀器測試分析，獲得了每個(gè)樣本所含的多種化學(xué)元素的濃度數(shù)據(jù)。另一方面，按照2公里的間距在那些遠(yuǎn)離人群及工業(yè)活動(dòng)的自然區(qū)取樣，將其作為該城區(qū)表層土壤中元素的背景值。附件1列出了采樣點(diǎn)的位置、海拔高度及其所屬功能區(qū)等信息，附件2列出了8種主要重金屬元素在采樣點(diǎn)處的濃度，附件3列出了8種主要重金屬元素的背景值?，F(xiàn)要求你們通過數(shù)學(xué)建模來完成以下任務(wù)：(1) 給出8種主要重金屬元素在該城區(qū)的空間分布，并分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度。(2) 通過數(shù)據(jù)分析，說明重金屬污染的主要原因。(3) 分析重金屬污染物的傳播特征，由此建立模型，確定污染源的位置。(4) 分析你所建立模型的優(yōu)缺點(diǎn)，

29、為更好地研究城市地質(zhì)環(huán)境的演變模式，還應(yīng)收集什么信息？有了這些信息，如何建立模型解決問題？這里只考慮第3問。由于城市表層土壤中，地表徑流沖刷為主要傳播途徑，而擴(kuò)散的作用較小，因此只考慮對(duì)流的影響。設(shè)u(x,y,t)是t時(shí)刻在(x,y)處的重金屬元素的濃度，對(duì)平面上任意區(qū)域，在t,t+t時(shí)段，通過的邊界流出的總和為其中q是流量向量，n是區(qū)域的外法線方向。在同一時(shí)段內(nèi)，中物質(zhì)數(shù)量的變化為 從而有令t0得到等式兩邊分別是曲面積分和三重積分。不容易比較。為了方便比較，我們利用Gauss公式把它換一下：代入上一個(gè)式子得到由于區(qū)域的任意性得到這是流動(dòng)的基本方程，反映了流動(dòng)的共性。流量沒有數(shù)據(jù)，但由于考慮重

30、金屬的傳播的成因主要是由于雨水的沖刷。因此傳播的快慢應(yīng)該與重金屬的濃度及地面的陡峭程度有關(guān)。設(shè)地面高度為h(t,x,y)，則可以假設(shè) 代入上面的式子即得到最終模型。練習(xí)題：( 2005)長江水質(zhì)的評(píng)價(jià)和預(yù)測水是人類賴以生存的資源，保護(hù)水資源就是保護(hù)我們自己，對(duì)于我國大江大河水資源的保護(hù)和治理應(yīng)是重中之重。專家們呼吁：“以人為本，建設(shè)文明和諧社會(huì)，改善人與自然的環(huán)境，減少污染?！?長江是我國第一、世界第三大河流，長江水質(zhì)的污染程度日趨嚴(yán)重，已引起了相關(guān)政府部門和專家們的高度重視。2004年10月，由全國政協(xié)與中國發(fā)展研究院聯(lián)合組成“保護(hù)長江萬里行”考察團(tuán)，從長江上游宜賓到下游上海，對(duì)沿線21個(gè)重

31、點(diǎn)城市做了實(shí)地考察，揭示了一幅長江污染的真實(shí)畫面，其污染程度讓人觸目驚心。為此，專家們提出“若不及時(shí)拯救，長江生態(tài)10年內(nèi)將瀕臨崩潰”（附件），并發(fā)出了“拿什么拯救癌變長江”的呼喚（附件2）。附件3給出了長江沿線17個(gè)觀測站（地區(qū)）近兩年多主要水質(zhì)指標(biāo)的檢測數(shù)據(jù)，以及干流上個(gè)觀測站近一年多的基本數(shù)據(jù)（站點(diǎn)距離、水流量和水流速）。通常認(rèn)為一個(gè)觀測站（地區(qū)）的水質(zhì)污染主要來自于本地區(qū)的排污和上游的污水。一般說來，江河自身對(duì)污染物都有一定的自然凈化能力，即污染物在水環(huán)境中通過物理降解、化學(xué)降解和生物降解等使水中污染物的濃度降低。反映江河自然凈化能力的指標(biāo)稱為降解系數(shù)。事實(shí)上，長江干流的自然凈化能力可

32、以認(rèn)為是近似均勻的，根據(jù)檢測可知，主要污染物高錳酸鹽指數(shù)和氨氮的降解系數(shù)通常介于0.10.5之間，比如可以考慮取0.2(單位：1/天)。附件4是“19952004年長江流域水質(zhì)報(bào)告”給出的主要統(tǒng)計(jì)數(shù)據(jù)。下面的附表是國標(biāo)(GB3838-2002) 給出的地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)中4個(gè)主要項(xiàng)目標(biāo)準(zhǔn)限值，其中、類為可飲用水。請(qǐng)你們研究下列問題：(1)對(duì)長江近兩年多的水質(zhì)情況做出定量的綜合評(píng)價(jià)，并分析各地區(qū)水質(zhì)的污染狀況。(2)研究、分析長江干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的污染源主要在哪些地區(qū)?(3)假如不采取更有效的治理措施，依照過去10年的主要統(tǒng)計(jì)數(shù)據(jù)，對(duì)長江未來水質(zhì)污染的發(fā)展趨勢做出預(yù)測分析

33、，比如研究未來10年的情況。(4)根據(jù)你的預(yù)測分析，如果未來10年內(nèi)每年都要求長江干流的類和類水的比例控制在20%以內(nèi)，且沒有劣類水,那么每年需要處理多少污水？ (5)你對(duì)解決長江水質(zhì)污染問題有什么切實(shí)可行的建議和意見。擴(kuò)散模型擴(kuò)散是自然界的另一個(gè)基本現(xiàn)象。如煙霧向空中擴(kuò)散，污染物在水中擴(kuò)散等。在社會(huì)活動(dòng)中，疾病的傳播、謠言的泛濫也有類似的動(dòng)態(tài)特征。熱量的傳播也是一種擴(kuò)散過程。擴(kuò)散的特點(diǎn)：物質(zhì)從濃度高的地方向濃度低的地方轉(zhuǎn)移。即物質(zhì)沿濃度的負(fù)梯度方向流動(dòng)。在曲面一點(diǎn)(x,y,z)的流量為其中C(x,y,z,t)是t時(shí)刻的濃度，k稱為擴(kuò)散系數(shù)。 在同一時(shí)段內(nèi)，中物質(zhì)數(shù)量的變化為 從而有令t0得到再利用Gauss積分公式得到從而 由的任意性得到熱的傳導(dǎo)和擴(kuò)散過程類似，熱的傳導(dǎo)是熱量由溫度高的地方向溫度低的地方流動(dòng)的結(jié)果。(可以稱作熱的擴(kuò)散！)引入變量 u(x,t)：x位置t時(shí)刻的溫度；J(x,t): x位置t時(shí)刻的熱流速度。熱力學(xué)定律：(1)熱流從溫度高的地方向溫度低的地方流動(dòng)，熱流速度正比于溫度函數(shù)的梯度；(2)體積元素的溫度的升高正比于流入它的熱量。第一個(gè)定律的數(shù)學(xué)描述為第二個(gè)定律的數(shù)學(xué)描述為把第一個(gè)式子代入第二個(gè)式子得到利用Gauss公式得到當(dāng)溫度穩(wěn)定在某個(gè)分布時(shí)，熱量的變化為0，方程化為這個(gè)方