數(shù)學(xué)建模冰山運(yùn)輸

1、拖運(yùn)冰山的費(fèi)用估算摘要據(jù)資料顯示：地球上的水資源有97%被鹽化，僅有3%是淡水資源，且大部分是主要分布在南北兩極地區(qū)的固體冰川，眾所周知，水是生命的源泉，是人類賴以生存和發(fā)展的重要物質(zhì)資源之一。而在一些淡水資源非常匱乏的地方，各種產(chǎn)水方式應(yīng)運(yùn)而生，本文主要針對(duì)以盛產(chǎn)石油著稱的波斯灣地區(qū)，在采用高成本的海水淡化下，專家建議用拖船從相距9600km外的南極把冰山運(yùn)到波斯灣地區(qū)，取代海水淡化的方法，在拖運(yùn)過(guò)程中采用三種不同的船型、船速會(huì)產(chǎn)生不同的費(fèi)用，為了使獲得淡水資源費(fèi)用合理化。需要對(duì)總費(fèi)用進(jìn)行合理的估算，為了解決這個(gè)問(wèn)題，我們建立初等函數(shù)模型，計(jì)算船只以最大運(yùn)量運(yùn)一次的總費(fèi)用和剩余的冰山體積；總

2、費(fèi)用=日租金天數(shù)+總?cè)剂舷馁M(fèi)用，考慮到實(shí)際運(yùn)輸過(guò)程中，冰山的融化會(huì)引起燃料消耗費(fèi)用的變化，我們應(yīng)該求出體積變化的臨界值（即所運(yùn)的冰山何時(shí)會(huì)縮小到原來(lái)體積的工），再根據(jù)表2可求出對(duì)應(yīng)的消耗原料。剩10余冰山的體積=船的最大運(yùn)量-運(yùn)輸過(guò)程中融化的冰。運(yùn)輸過(guò)程中冰的融化速率與冰山到達(dá)處與南極的距離有關(guān)。根據(jù)表3,在特定的船速條件下，將冰山融化速率與南極距離看成線性關(guān)系yaxb,根據(jù)曲線積分可以求出融化的冰，然后求剩余的冰轉(zhuǎn)換成水的體積，計(jì)算每立方米水的成本，經(jīng)過(guò)計(jì)算選擇大型船，船速為5km/h,每立方米水的成本0.054英鎊/m3。關(guān)鍵字：燃料消耗費(fèi)用冰山融化的體積融化速率臨界值;.1、問(wèn)題重述在

3、以盛產(chǎn)石油著稱的波斯灣地區(qū)，浩瀚的沙漠覆蓋著大地，水資源十分缺乏，不得不采用淡化海水的辦法為國(guó)民提供用水，成本大約是每立方米0.1英鎊，有些專家提出從相距9600km外的南極用拖船運(yùn)送冰山到波斯灣，以取代淡化海水的辦法。在運(yùn)送冰山的過(guò)程中，拖船的租金、運(yùn)量、燃料消耗及冰山運(yùn)送過(guò)程中融化速率等方面的數(shù)據(jù)如下：(1)三種拖船的日租金和最大運(yùn)量如表1所示。表1船型小中大日租金/(英鎊)4.06.28.0取大運(yùn)里/m5105106107(2)燃料消耗(英鎊/km),主要依賴于船速和所運(yùn)冰山的體積，船型的影響可以忽略，如表2所小表2冰山體積/m3、船速/km/h10510610718.410.512.6

4、210.813.516.2313.216.519.8(3)冰山運(yùn)輸過(guò)程中的融化速率(m/d),指在冰山與海水接觸處每天融化的深度，融化速率除與船速有關(guān)外，還與運(yùn)輸過(guò)程中冰山到達(dá)處與南極的距離有關(guān),這是由于冰山要從南極運(yùn)往赤道附近的緣故，如表3所示。表3極距離/km船速/km/h01000>40001230000.10.150.20.30.450.6根據(jù)所給數(shù)據(jù)，建立數(shù)學(xué)模型，解決問(wèn)題：試選擇拖船的船型與船速，使冰山到達(dá)目的地后，可以得到的每立方米水所花的費(fèi)用最低，并與海水淡化的費(fèi)用相比較2、問(wèn)題分析此題研究的是每立方米水成本的數(shù)學(xué)建模問(wèn)題，題目給出的數(shù)據(jù)有三種船型，三種船速，那么可分9種

5、情況去求解問(wèn)題。求解每立方米水的成本就要求解出運(yùn)輸一次的費(fèi)用和剩余的冰山體積。由于在運(yùn)輸過(guò)程中，單位路程燃料消耗與冰山的體積成正相關(guān)關(guān)系。但是冰山的體積是動(dòng)態(tài)變化的，從而單位路程燃料消耗也隨之改變。為了計(jì)算方便，利用均值求一段路程的燃料的消耗，而冰山的融化速率與冰山到達(dá)處與南極的距離有關(guān)，根據(jù)圖1可以得出當(dāng)船速確定時(shí)，融化速率與距南極距離成線性相關(guān)。同樣地可采用均值求解，可以確定到達(dá)目的地后冰山的剩余體積。將冰轉(zhuǎn)化為水，從而估算出每立方米水的成本。船速E/knih船速1,km/h1000200030004000500060007000S000900010000與南極的距離樂(lè)m融化速率與南極距震

6、的關(guān)系圖船速京km/h圖23、模型的假設(shè)及符號(hào)說(shuō)明拖船在托運(yùn)冰山的過(guò)程中，有以下假設(shè)：假設(shè)1、假設(shè)拖船航行過(guò)程中的船速不變，航行不考慮天氣等任何因素的影響,總航行距離9600km假設(shè)2、假設(shè)冰山形狀為球形，球面各點(diǎn)的融化速率相同；假設(shè)3、假設(shè)冰山到達(dá)目的地后，1m3的冰可以融化成0.85m3的水。假設(shè)4、假設(shè)冰山在運(yùn)輸過(guò)程中融化的速率與距離為線性關(guān)系。假設(shè)5、三種型號(hào)的船是按最大運(yùn)輸量拖運(yùn)冰山。假設(shè)6、冰山到達(dá)目的地后融化過(guò)程中不考慮損耗。假設(shè)7、假設(shè)拖船所在地就在南極，不用考慮返程費(fèi)用。符號(hào)說(shuō)明：符合符號(hào)說(shuō)明Vi拖船的速度i1,2,3ZiZ拖船的日租金Wj燃料消耗速度i1,2,3j1,3,5

7、li品巨離南極的品目離i1,2,3uij融化速度i1,2,3R冰山的半徑R耗運(yùn)動(dòng)過(guò)程中冰山融化的深度j1,3,5Rj0,Vj初始時(shí)刻冰山的半徑j(luò)1,3,5最后剩余的冰山體積j1,3,5R末最后時(shí)刻冰山的半徑j(luò)1,3,5ZiZi拖船的總租金LL為總長(zhǎng)度mjP全程燃料費(fèi)用i1,2,3j1,3,5運(yùn)冰的總費(fèi)用i1,2,34、模型建立為了求解每立方米水的成本，就要計(jì)算出拖船運(yùn)輸一次的總費(fèi)用和剩余的冰山體積，首先，根據(jù)題目所給的冰山融化速率與船速、距南極距離的相關(guān)數(shù)據(jù)，我們可以求解出當(dāng)船速確定的時(shí)候，冰山融化速率與距南極距離成線性關(guān)系，利用均值和曲線積分可以求解出船只行駛一段距離與融化冰山的體積的關(guān)系。

8、RUL2L3Ui1-Ui2-Ui3ViViR。后（1）R末R0R耗VjRj末0.85其次，由上述計(jì)算可以得出燃料消耗與冰山體積的關(guān)系，因?yàn)楸襟w積是動(dòng)態(tài)變化的，所以燃料消耗也是動(dòng)態(tài)的，為了計(jì)算的簡(jiǎn)單化，我們把計(jì)算一段路程的燃料消耗簡(jiǎn)化為固定兩個(gè)端點(diǎn)的單位路程燃料消耗與行駛路程的乘積和的均值。W1jIW2jIW3jI24ViZi(2)最后：由上述條件利用總費(fèi)用和冰融化成水后的單位成本:PijpVjVjR末0.855、模型求解首先，計(jì)算出不同船型以不同速度拖運(yùn)冰山所需的總費(fèi)用。其次，計(jì)算出不同船型以不同速度拖運(yùn)冰山到達(dá)目的地，融化成的水。最后，由上述條件求出冰山拖運(yùn)到目的地融化為水后每立方米的成本

9、，與淡化每立方米水的成本相比較得出成本最低的方案。第一步：計(jì)算總費(fèi)用mijwijliW2jI2W3jI3(4)不同船型以不同船速拖冰山，所需的燃料費(fèi)：表4,33小5105/m中106/m3大107/m31km/h無(wú)意義無(wú)意義118780英鎊3km/h100800英鎊121070英鎊155520英鎊5km/h120960英鎊158202英鎊190080英鎊乙拖船的總租金：(5)L24vi用不同船型拖運(yùn)所用的總租金:表5型r5.3小510/m中106/m3大107/m31km/h無(wú)意義無(wú)意義3200英鎊3km/h533.3332英鎊826.6665英鎊:1066.7英鎊5km/h320英鎊496英

10、鎊640英鎊總費(fèi)用:«Zimj(6)表6r_.53小510/m中106/m3大107/m31km/h無(wú)意義無(wú)意義119100英鎊3km/h101333.33英鎊121896.66英鎊156586.7英鎊5km/h121280英鎊158698英鎊190720英鎊第二步:R耗£U1£52£幾不同船型以不同船速拖冰山，冰山融化的深度:表7，一一53小510/m中106/m3大107/m31km/h84.583384.583384.58333km/h42.291742.291742.29175km/h33.833333.833333.8333R末%R耗(9)計(jì)算

11、可知用小船拖運(yùn)時(shí)以1km/h到達(dá)目的地時(shí)冰塊已經(jīng)完全融化，中船拖運(yùn)時(shí)以1km/h到達(dá)目的地時(shí)冰塊已經(jīng)完全融化所以不再考慮這兩種情況。不同船型以不同船速拖冰山，到達(dá)目的地剩余體積：R末R。R耗,3小5*10A5/m,3中10A6/m,3大10A7/m1km/h無(wú)意義無(wú)意義3494840m3km/h31403.5m332236m33194000m5km/h315310m393954m34165900m冰山到達(dá)后無(wú)其他損耗全部都融化成水的體積：VjRj末0.85表9小5*10人5/m33中10A6/m3大10A7/m1km/h無(wú)意義無(wú)意義3420614m3km/h31193m327401m32714

12、900m5km/h313014m379861m33541015m每立方米冰的成本piLVj從南極運(yùn)冰到目的地的每立方米成本為:表10,3小5*10人5/m»3中10A6/m3大10A7/m1km/h無(wú)意義無(wú)意義30.283英鎊/m3km/h384.940英鎊/m一,一34.449英鎊/m30.058英鎊/m5km/h39.320英鎊/m31.987英鎊/m30.054英鎊/m根據(jù)計(jì)算結(jié)果與已知條件（成本每立方米0.1英鎊）對(duì)比得出0.054英鎊/m3的成本價(jià)格是最低的，所以我們選擇大型船，船速為5km/h6、模型評(píng)價(jià)優(yōu)點(diǎn)：1、構(gòu)造的模型比較簡(jiǎn)單2、該模型的建立將實(shí)際問(wèn)題模型化，復(fù)雜問(wèn)

13、題簡(jiǎn)單化；3、我們采用了均值與曲線積分的計(jì)算方法，使計(jì)算變得相對(duì)簡(jiǎn)單；4、此模型的全局規(guī)劃比較合理。缺點(diǎn)：1、利用均值求解融化冰山的體積與燃料消耗不準(zhǔn)確；2、運(yùn)輸過(guò)程融化速率與距南極距離并不成線性關(guān)系。7、模型的改進(jìn)與推廣模型改進(jìn)：雖然得到的最低費(fèi)用0.054英鎊/m3小于海水淡化需要的成本，但是模型中未考慮天氣等其它因素會(huì)導(dǎo)致拖船過(guò)程中成本變高。我們?cè)谟?jì)算過(guò)程中應(yīng)該以冰山拖運(yùn)過(guò)程中最大損耗、燃料的最大費(fèi)用，來(lái)避免由自然因素帶來(lái)的船只到達(dá)目的地的時(shí)間延期增加的費(fèi)用。方法改進(jìn)冰山融化速率r(m/d)與船速u(km/h)與南極距離d(km)的關(guān)系r是u的線性函數(shù)；0d4000時(shí)u與d成正比，d40

14、00時(shí)u與d無(wú)關(guān),aid(1bu)a2(1bu)0d4000d4000(10)代入(u,d)=(1,3);r(u,d)=(1,5)得出a1=7.5105,a2=0.225b=0.33航行t天后，融化速率為a1ud(1bu)trta2(1bu)40000tu4000tu(11)0.225(10.33u)(12)540007.5105u(10.33u)t0tu4000又因?yàn)閠大冰山的半徑為tRRortt0(13)故剩余冰山的體積為Vt3VRt3(14)燃料消耗費(fèi)用q(英鎊/km)z對(duì)u線性，對(duì)logov線性,zC!uc2igvC3代入數(shù)據(jù)得出Ci=0.3,c2=6C3="1z3V航行第t

15、大的燃料消耗t-，4，z7.2uu6Ig-(43廉)0而日租金與船的型號(hào)有關(guān)故f(Vo)4.06.28.0_5V05105105V0106106V0107所以船型運(yùn)輸一次的總費(fèi)用為400f(v0)T4P-027.2u(u6)Ig一ut13Trt1t1此模型的建立會(huì)減少誤差，由第一次模型計(jì)算的結(jié)果,5km/h0的數(shù)據(jù)代入上述公式，計(jì)算結(jié)果。參考文獻(xiàn)(1)、姜啟源、謝金星，葉俊，數(shù)學(xué)模型(第三版)北京,(15)(16)(17)(18)(19)將大型船，船速為高等教育出版社，2003;(2)、韓忠庚，數(shù)學(xué)建模方法及其應(yīng)用，北京高等教育出版社，2005;(3)、吳建國(guó)，數(shù)學(xué)建模案例精編北京，中國(guó)水利水

16、電出版社2005;附錄模型求解在Matlab中的實(shí)現(xiàn)：畫圖1:船速固定時(shí)，冰山體積與燃料的關(guān)系圖>>clear>>a=0100040009600;>>b1=00.10.30.3;>>b3=00.150.450.45;>>c=0:100:9600;>>interp1(a,b1,G'cubic);>>interp1(a,b3,c,'cubic);>>interp1(a,b5,G'cubic);>>plot(a,b1,'o',a,b1,'g

17、9;,a,b3,'o',a,b3,'r',a,b5,'o',a,b5,'b')>>gtext(船速1/km/h')>>gtext(船速3/km/h')>>gtext(船速5/km/h')畫圖2:船速固定時(shí)，燃料消耗與冰山體積的關(guān)系圖> >a=10A510A610A7;> >b1=8.410.512.6;> >b3=10.813.516.2;> >b5=13.216.519.8;> >c=10A5:100:10A7

18、;> >interp1(a,b1,c,'cubic');> >interp1(a,b3,c,'cubic');> >interp1(a,b5,c,'cubic');>>plot(a,b1,'o',a,b1,'g',a,b3,'o',a,b3,'r',a,b5,'o',a,b5,'b')> >gtext('船速1/km/h')> >gtext('船速3/km/

19、h')>>gtext('船速5/km/h')計(jì)算過(guò)程>>x=(0.75*10A5/(pi)A(1/3)x1=(0.75*10A5/(pi)A(1/3)x=(0.75*10A6/(pi)A(1/3)x=(0.75*10A7/(pi)A(1/3)1000/(5*24)8.3333*0.13000/(5*24)*0.25600/(5*24)*0.6x=(0.75*10A5/(pi)A(1/3)x=(0.75*5*10A5/(pi)A(1/3)>>133.6505-2.0833-12.5ans=119.0672>>119.0672

20、-62.0350>>84.5833-2.0833-12.5>>60-2.0833-12.5ans=45.4167>>12.9678/0.3*24ans=1.0374e+003>>9600-1.0374e+003ans=8.5626e+0036441.1*13.5+3158.9*10.8ans=1.2107e+005>>9540*16.5+60*13.2ans=158202>>8562.6*12.6+1037.4*10.5ans=1.1878e+005>>9600*16.2ans=155520>>9600*19.8ans=190080>>1000/(5*24)ans=8.3333>>8.3333*0.1ans=0.8333