冰山运输数学模型

冰山运输数学模型摘 要当今社会,水资源短缺已成为世界性问题,水资源紧张地区正不断扩大,除淡化海水的方法外,专家提出从相距 9600 千米以外的南极托运冰山到波斯湾,将其化成冰水从而取代淡化海水作为国民用水本文所要解决的是选择合适的拖船与船速使得冰山到达目的地后得到每立方米水所花的费用最低的问题,由此建立了一个关于费用 y 的数学模型首先,根据表 3 中的拖船速率 v 和拖船与南极的距离可知冰山融化速率,从而确定剩余的冰山体积然后,根据表 2 中的船速 v 和运输过程中剩余冰山的体积 N 可知每千米燃料消耗量 q0 ,从而可以求出所需燃料总消耗量 Q ,再分别选取小、中、大三种船型确定拖船的租金总费用 M , 则运输总费用Y Q M,运输每立方米水所花费用即为Yy 0.0626 根据运输每立方米水所花的费用最低,将该问题归结为优化问0.85N题,运用积分方法,通过 Matlab 计算,得到最优解确定船型和船速 , 再与海水淡化的费用相比较,确定其可行性关键字:冰山体积 融化速率 燃料消耗量 最优化1. 问题重述在以石油着称的波斯湾地区,浩瀚的沙漠覆盖着大地,水资源十分缺乏,不得不采用淡化海水的办法为国民提供用水。
成本大约是每立方米 0.1 英镑有些专家提出从相距9600km外的南极用拖船运送冰山到波斯湾,以取代淡化海水的办法在运送冰山的过程中,拖船的租金、运量、燃料消耗以及冰山运送过程中融化速率等方面的数据如下:(1)三种拖船的日租金和最大运量如表 1. 所示表 1.船型 小 中 大日租金 /( 英镑 ) 4.0 6.2 8.0最大运量 / m3(2)燃料消耗(英镑 /km),主要依赖于船速和所运冰山的体积,船型的影响可以忽略,如表 2. 所示表 2.冰山体积 / m3船速/km/h18.410.512.6310.813.516.2513.216.519.8(3)冰山运输过程中的融化速率( m/d), 指在冰山与海水接触处每天融化的深度融化速率除与船速有关,还与运输过程中冰山到达与南极的距离有关,这是由于冰山要从南极运往赤道附近的缘故如表 3. 所示表 3.与南极的距离01000>4000/km船速 /km/h100.10.3300.150.45500.20.6本文所要解决的问题是: 选择拖船的船型与船速,使冰山到达目的地后,可以得到的每立方米水所花的费用最低,并与海水淡化的费用相比较拖船在拖运冰山的过程中,有以下假设:(1) 拖船航行过程中船速不变,航行不考虑天气等任何因素的影响,总航行距离9600km;(2) 冰山形状为球形,球面各点的融化速率相同;(3) 冰山到达目的地后, 1 m3 的冰可以融化成 0.85 m3 的水。
2. 问题分析为更好地计算冰山运输的费用,我们对问题进行了分析根据题目已给的资料和数据,我们发现:冰山的运输主要和拖船的租金、运量、燃料消耗及冰山运输过程中融化速率有关,因此,我们可以把问题分成以下五步来分析解决:1、冰山的融化规律冰山融化是冰山体积的变化,而冰山体积的变化 N 实质上是冰山半径的变化由表 3中数据可发现,冰山运输过程中的融化速率 S 与船速 v 和拖船距离南极 d 呈线性关系,因此,根据数据拟合公式可以得出三者的函数关系式,然后联系拖船航行的天数 t ,就得到冰山的融化规律2、燃料消耗费用 Q 由表二数据可发现,燃料消耗 q 与船速 v 和冰山体积 N 呈线性关系,但考虑到 N 的值比较大,我们可以取体积 N 的对数 lg N ,再根据数据拟合公式可以得出三者的函数关系式,然后联系拖船航行的天数 t ,就得到燃料消耗费用3、冰山运输的总费用冰山运输的总费用由拖船的租金和燃料消耗两部分组成,燃料消耗费用上面已经分析解决,拖船的租金则由表 1 的数据可确定,进而得到冰山运输的总费用4、冰山到达目的地可获得的水体积最终获得的水体积是冰山到达目的地所需天数 T 与冰山融化规律得到。
5、每立方米水的费用上面得到的冰山运输费用的函数与得到的冰山到达目的地后可获得水的体积的函数相除,即可得到每立方米水的费用3. 模型假设与符号说明3.1 模型假设:假设一: 拖船航行过程中船速不变,航行不考虑天气等任何因素的影响,总航行距离 9600km;假设二: 冰山形状为球形,球面各点的融化速率相同;假设三: 冰山到达目的地后, 1 m3 的冰可以融化成 0.85 m3 的水假设四: 冰山运输距离南极 4000m之后忽略温度对冰山融化的影响3.2 符号说明:符号 符号说明拖船的日租金(小、中、大)i 1,2,3拖船的速度拖船到达波斯湾所用的总时间距离南极为 d 时所需的时间冰山融化速率拖船与南极的距离南极到波斯湾的总距离每千米燃料消耗费用每天燃料消耗费用燃料消耗的总费用拖船到达波斯湾所用租金冰山原来的半径冰山融化后的半径冰山开始的体积变化后冰山体积每立方水所需费用拖船到达波斯湾时所需的总费用4. 模型建立本文的目的是选择出拖船的最适船型与船速,使冰山到达目的地后,可以得到的每立方米水所花的费用最低,通过问题分析可以分别从冰山的融化规律、燃料消耗费用、冰山运输的总费用、冰山到达目的地可获得的水体积和每立方米水的费用建立数学模型。
4.1 冰山的融化规律假设运输过程中距离南极的距离为 dkm ,拖船的船速为 v km h ,融化速率为 S ,由表三可以得到如图 1图 1通过图 1,很明显看出:当距离 0 d 4000km时,融化速率 S 不仅与船速 v 成线性关系,而且也与距离 d 呈线性关系,也就是说,船速越快,距离越远,冰山接触的温度越高,导致它融化越快;当 d 4000km时,距离对融化速率的影响忽略,没有影响,所以这时融化速率 S 只与船速 v呈线性关系那么通过数据拟合公式可设融化速率S 的函数关系式为:k1d k2vb0 d4000S(1)k2vbd4000然后将表 2 相关数据代入可得到:当拖船从南极出发行走 t 天时,与南极距离为d 24vt(2)由( 1)、(2)及 k1,k2 ,b的值得:1vt 0.075v 0.2250t500S1253v0.075v 0.225500t4003vv接下来求变化后冰山体积设第 t 天冰山半径为 Rt ,体积为 N ,则Rt R0tSdt0N 43 Rt3, N0 43 R03其中 R0 ,V0 为从南极启运时冰山的开始半径和体积由( 4)、(5)得冰山体积为:N4(3 3N0tSdt3340)4.2 燃料消耗费用由表 2 中的数据,通过 Matlab 拟合可得图形如图2 所示:图 2(3)(4)(5)(6)冰山 体积x由图可知每千米燃料消耗量 q0 与船速 v 和冰山体积 N 的对数呈线性关系,则可设其函数关系式为:q c1 v c2 lg Nc3( 7)0其中 c1,c2,c3 为待定参数,根据表2 中数据可求得所以q0 0.3 v 6 lg N 1(8)则每天燃料消耗量为q24v q0(9)把( 2)式代入( 3)式,再根据( 6)式得433N0t3( 10)q7.2v v 6lg3 (40Sdt) 1所以,总燃料消耗量为Qqdt7.2v v6 lg4(3 3N0Sdt)31dtTTt00340(11)4.3 冰山运输总费用冰山运输的总费用是由租金的总费用和燃料消耗总费用相加。
由表 1 知船的日租金取决于船型,船型又由冰山的初始体积 N0 决定,记日租金为 m(N0 ) ,则有4.0N05105m(N0 ) 6.2 5 105N0106(12)8.0106N0107又冰山运输总时间为:TD9600400(13)24v24vv所以租金总费用为:Mm(N0 )T400 m(N0)(14)v由( 11)、(8)式可得冰山运输总费用为400m(N0 )T43N 0t3dt ( 15)Y M Qv07.2v v 6 lg3(340Sdt) 14.4 冰山运到目的地后化成水的体积为:433N0t3(16)W 0.85N 0.853 (40Sdt)4.5 每立方米水所需费用为:Yy W400v m( N0 )T43N0t37.2v v6lg3(340Sdt) 1 dt0(17)43N0t30.85( 3Sdt)3405. 模型求解首先,对( 3)式进行积分运算,有t5001 vt 0.075v 0.225 dt4003v500v0.075v 0.225 dt (18)进而求出Sdt001253v不同的船速下冰山到达目的地之后总的融化深度,如表 4 所示:表 4船速 v1233.544.55103.33364.583351.666747.976245.208343.055641.33333其次,由( 17)式可以发现,这个模型最终由船速v 和冰山体积 N0 决定。
由( 12)、(13)、( 17)式以及表2、表 4,通过 Matlab 计算可以得到每立方米水所需费用,结果如表 5所示:表 51233.544.55负值负值负值负值负值负值负值负值负值16.9470 12.3749 7.80275.86423.59620.85870.07280.07140.06660.06420.06300.0626由表 5 数据可知,当船速 v 3km h, N0 106 时,不符合实际,所得结果为负值, 所以,应选择大型拖船N0107 ,速度v5km h,则每立方米水所需费用为y0.0626英镑因为淡化海水每立方米约为0.1 英镑,所以,冰山运输比淡化海水的成本更低6. 模型评价优点:模型中运用图表分析数据,使各因素间的关系明朗化,运输费用只考虑租费和燃料消耗费,并将冰山体积看成球形,简化了计算,思路清晰简单缺点:模型中只考虑了拖船的租费和燃料消耗费用,没有考虑到影响航行的天气因素和社会因素(如人工费和税收等) ,所以冰山运输方法是否可行还得取决于综合考虑因素7. 参考文献[1] 姜启源 . 数学模型(第三版) [M]. 北京:高等教育出版社, 2003.[2] 王文波 . 数学建模及其基础知识详解(第一版) [M]. 武汉大学出版社, 2006.[3] 刘玉琏 . 数学分析(上、下) [M]. 北京:高等教育出版社, 1988.8.附 录附录一:>> a=[1 3 5];b=[0.1 0.15 0.2]; % 融化速率与船速、与南极距离的关系图>> x0=1:0.1:5;>> y0=interp1(a,b,x0,'linear');>> c=[0 1 4 5];d=[0 0.1 0.3 0.3];>> x1=0:0.1:5;>> y1=interp1(c,d,x1,'linear');>> gtext(' 船速 ')>> gtext(' 与南极距离× 1000')>> a=[1 3 5];b=[8.4 10.8 13.2]; % 燃料消耗与船速、冰山体积的关系图>> x0=1:0.1:5;>> y0=interp1(a,b,x0,'linear');>> c=[5 6 7];d=[8.4 10.5 12.6];>> x1=10^5:1:10^7;>> y1=interp1(c,d,x1,'linear');>> plot(a,b,'r',c,d,'--')>> gtext(' 船速 ')>> gtext(' 冰山体积 10 i ')>> syms k2 b;[k2,b]=solve('0.3=k2*1+b','0.45=k2*3+b',k2,b) % 冰山融化速率的表达式k2 =0.0750b =0.2250>> syms k1;k1=solve('0.1=k1*1000*(0.075*1+0.225)',k1)k1 =0.00033333333333333333333333333333333>> syms t,s=int('0.0024*t',t,0,500/3)+int('0.3',500/3,400) % 当 v 为1 时冰山的融化速率s =103.33333333333333333333333333333>> syms t,s=int('0.375*2/125*t',t,0,250/3)+int('0.375',250/3,200)s = % 当 v 为 2 时冰山的融化速率64.583333333333333333333333333333>> syms t,s=int('0.45*3/125*t',t,0,500/9)+int('0.45',500/9,400/3)s = % 当 v 为 3 时冰山的融化速率51.666666666666666666666666666667>> syms t,s=int('0.525*4/125*t',t,0,125/3)+int('0.525',125/3,100)s = % 当 v 为 4 时冰山的融化速率45.208333333333333333333333333333>> syms t,s=int('0.6*5/125*t',t,0,100/3)+int('0.6',100/3,80) % 当 v 为5 时冰山的融化速率s =41.333333333333333333333333333333>>(400*4+50.4*int('log(4*pi/3*(((3*10^5/4*pi)^(1/3))-103.3333)^3)-1',0,400))/(0.85*4*pi/3*(4*pi/3*(((3*10^5/4*pi)^(1/3))-103.3333)^3))ans =-0.2201 - 0.0591i>>(200*4+115.2*int('log(4*pi/3*(((3*10^5/4*pi)^(1/3))-64.5833)^3)-1',0,200))/(0.85*4*pi/3*(4*pi/3*(((3*10^5/4*pi)^(1/3))-64.5833)^3))ans =-2.4661e+002 -2.1663e+002i>>(400/3*4+115.2*int('log(4*pi/3*(((3*10^5/4*pi)^(1/3))-51.6667)^3)-1',0,400/3))/(0.85*4*pi/3*(4*pi/3*(((3*10^5/4*pi)^(1/3))-51.6667)^3))>>(400/3*4+194.4*int('log(4*pi/3*(((3*10^5/4*pi)^(1/3))-51.6667)^3)-1',0,400/3))/(0.85*4*pi/3*(4*pi/3*(((3*10^5/4*pi)^(1/3))-51.6667)^3))>>(100*4+288*int('log(4*pi/3*(((3*10^5/4*pi)^(1/3))-45.2083)^3)-1',0,100))/(0.85*4*pi/3*(4*pi/3*(((3*10^5/4*pi)^(1/3))-45.2083)^3))>>(80*4+396*int('log(4*pi/3*(((3*10^5/4*pi)^(1/3))-41.3333)^3)-1',0,80))/(0.85*4*pi/3*(4*pi/3*(((3*10^5/4*pi)^(1/3))-41.3333)^3))>>(80*6.2+396*int('log(4*pi/3*(((3*10^6/4*pi)^(1/3))-41.3333)^3)-1',0,80))/(0.85*4*pi/3*(4*pi/3*(((3*10^6/4*pi)^(1/3))-41.3333)^3))ans =3.596>>(100*6.2+288*int('log(4*pi/3*(((3*10^6/4*pi)^(1/3))-45.2083)^3)-1',0,100))/(0.85*4*pi/3*(4*pi/3*(((3*10^6/4*pi)^(1/3))-45.2083)^3))ans =7.8027>>(400/3*6.2+194.4*int('log(4*pi/3*(((3*10^6/4*pi)^(1/3))-51.6667)^3)-1',0,400/3))/(0.85*4*pi/3*(4*pi/3*(((3*10^6/4*pi)^(1/3))-51.6667)^3))ans =16.9470>>(400*8+50.4*int('log(4*pi/3*(((3*10^7/4*pi)^(1/3))-103.3333)^3)-1',0,400))/(0.85*4*pi/3*(4*pi/3*(((3*10^7/4*pi)^(1/3))-103.3333)^3))ans =0.8587>>(200*8+115.2*int('log(4*pi/3*(((3*10^7/4*pi)^(1/3))-64.5833)^3)-1',0,200))/(0.85*4*pi/3*(4*pi/3*(((3*10^7/4*pi)^(1/3))-64.5833)^3))ans =0.0728>>(400/3*8+194.4*int('log(4*pi/3*(((3*10^7/4*pi)^(1/3))-51.6667)^3)-1',0,400/3))/(0.85*4*pi/3*(4*pi/3*(((3*10^7/4*pi)^(1/3))-51.6667)^3))ans =0.0714>>(400/4.5*8+340.2*int('log(4*pi/3*(((3*10^7/4*pi)^(1/3))-43.0556)^3)-1',0,400/4.5))/(0.85*4*pi/3*(4*pi/3*(((3*10^7/4*pi)^(1/3))-43.0556)^3))ans =0.0626。