冰山运输数学模型

1、(1)三种拖船的日租金和最大运量如表1.所示冰山运输数学模型摘要当今社会，水资源短缺已成为世界性问题，水资源紧张地区正不断扩大，除淡 化海水的方法外，专家提出从相距9600千米以外的南极托运冰山到波斯湾，将其化成冰水从而取代淡化海水作为国民用水。本文所要解决的是选择合适的拖船与船 速使得冰山到达目的地后得到每立方米水所花的费用最低的问题，由此建立了一个 关于费用y的数学模型。首先，根据表3中的拖船速率V和拖船与南极的距离可知冰山融化速率，从而确定剩余的冰山体积。然后，根据表2中的船速V和运输过程中剩余冰山的体积N可知每千米燃料消耗量q0，从而可以求出所需燃料总消耗量Q，再分别选取小、中、大三种

2、船型确定拖船的租金总费用M,则运输总费用Y Q M，运输每立方米水所花费用即为yOYN0.0626。根据运输每立方米水所花的费用最低，将该问题归结为优化问题，运用积分方法，通过Matlab计算，得到最优解确定船型和船速，再与海水淡化的费用相比较，确定其可行性。关键字：冰山体积融化速率燃料消耗量最优化1问题重述在以石油着称的波斯湾地区，浩瀚的沙漠覆盖着大地，水资源十分缺乏，不得 不采用淡化海水的办法为国民提供用水。成本大约是每立方米英镑。有些专家提出 从相距9600km外的南极用拖船运送冰山到波斯湾，以取代淡化海水的办法。在运送冰山的过程中，拖船的租金、运量、燃料消耗以及冰山运送过程中融化 速率

3、等方面的数据如下：表1.船型小中大日租金/（英镑）最大运里/m（2）燃料消耗（英镑/km）,主要依赖于船速和所运冰山的体积，船型的影响 可以忽略，如表2.所示。深度。融化速率除与船速有关，还与运输过程中冰山到达与南极的距离有关，这是 由于冰山要从南极运往赤道附近的缘故。如表3.所示。表3.与南极的距离km船速/km/h010004000103050本文所要解决的问题是： 选择拖船的船型与船速，使冰山到达目的地后，可以得到的每立方米水所花的费用最低，并与海水淡化的费用相比较。拖船在拖运冰山 的过程中，有以下假设：（1）拖船航行过程中船速不变，航行不考虑天气等任何因素的影响，总航行 距离9600k

4、m；（2）冰山形状为球形，球面各点的融化速率相同；（3）冰山到达目的地后，1m3的冰可以融化成 m3的水。2.问题分析 为更好地计算冰山运输的费用，我们对问题进行了分析。 根据题目已给的资料和数据， 我们发现： 冰山的运输主要和拖船的租金、 运量、 燃料消耗及冰山运输过程中融化速率有关，因此，我们可以把问题分成以下五步来 分析解决：1、 冰山的融化规律。 冰山融化是冰山体积的变化，而冰山体积的变化N实质上是冰山半径的变化。由表3中数据可发现，冰山运输过程中的融化速率S与船速v和拖船距离南极d呈 线性关系，因此，根据数据拟合公式可以得出三者的函数关系式，然后联系拖船航 行的天数t，就得到冰山的融

5、化规律。2、 燃料消耗费用Q。由表二数据可发现，燃料消耗q与船速v和冰山体积N呈线性关系，但考虑到N的值比较大，我们可以取体积N的对数lgN，再根据数据拟合公式可以得出三 者的函数关系式，然后联系拖船航行的天数t，就得到燃料消耗费用。3、冰山运输的总费用。 冰山运输的总费用由拖船的租金和燃料消耗两部分组成，燃料消耗费用上面已 经分析解决，拖船的租金则由表1的数据可确定，进而得到冰山运输的总费用。4、冰山到达目的地可获得的水体积。最终获得的水体积是冰山到达目的地所需天数T与冰山融化规律得到。5、每立方米水的费用。上面得到的冰山运输费用的函数与得到的冰山到达目的地后可获得水的体积 的函数相除， 即

6、可得到每立方米水的费用。3.模型假设与符号说明模型假设：假设一：拖船航行过程中船速不变，航行不考虑天气等任何因素的影响，总航行距离9600km;假设二：冰山形状为球形，球面各点的融化速率相同；假设三：冰山到达目的地后，1m3的冰可以融化成 m3的水。假设四：冰山运输距离南极4000m之后忽略温度对冰山融化的影响。符号说明：符号符号说明拖船的日租金（小、中、大）i 1,2,3拖船的速度拖船到达波斯湾所用的总时间距离南极为d时所需的时间冰山融化速率拖船与南极的距离:南极到波斯湾的总距离每千米燃料消耗费用每天燃料消耗费用燃料消耗的总费用拖船到达波斯湾所用租金冰山原来的半径冰山融化后的半径冰山幵始的体

7、积变化后冰山体积每立方水所需费用拖船到达波斯湾时所需的总费用4.模型建立本文的目的是选择出拖船的最适船型与船速，使冰山到达目的地后，可以 得到的每立方米水所花的费用最低，通过问题分析可以分别从冰山的融化规律、燃 料消耗费用、冰山运输的总费用、冰山到达目的地可获得的水体积和每立方米水的 费用建立数学模型。冰山的融化规律假设运输过程中距离南极的距离为dkm，拖船的船速为vkm h，融化速率为S， 由表三可以得到如图1。图1通过图1,很明显看出：当距离0 d 4000km时，融化速率S不仅与船速v成 线性关系，而且也与距离d呈线性关系，也就是说，船速越快，距离越远，冰山接 触的温度越高，导致它融化越

8、快；当d 4000km时，距离对融化速率的影响忽略， 没有影响，所以这时融化速率S只与船速v呈线性关系。那么通过数据拟合公式可设融化速率S的函数关系式为:(1)(6)S %d k2v b 0 d 4000 k2v b d 4000然后将表2相关数据代入可得到:当拖船从南极出发行走t天时，与南极距离为d 24vt由（1）、（2）及k1,k2,b的值得:燃料消耗费用由表2中的数据，通过Matlab拟合可得图形如图2所示:(2)糅vt 0.075V0.225S1250.075v 0.2255005003v400v(3)接下来求变化后冰山体积设第t天冰山半径为 &，体积为N，则tRtRO0Sdt(4)

9、N牛R3，N屯戍(5)其中R）,V0为从南极启运时冰山的幵始半径和体积。由（4）、（5）得冰山体积为:N 却34NOt0Sdt)33N冰山体积由图可知每千米燃料消耗量q与船速v和冰山体积N的对数呈线性关系，则可 设其函数关系式为：q。C v c2lg N三（7）其中c1,c2，c3为待定参数，根据表2中数据可求得所以q00.3 v 6 lg N 1（8）则每天燃料消耗量为q 24v q0（9）把（2）式代入（3）式，再根据（6）式得q 7.2v v 6 lg台（3丁。；Sdt）31（10）所以，总燃料消耗量为QTqdt；7.2v v 6 lg笛（3?0；Sdt）31dt(11)冰山运输总费用冰

10、山运输的总费用是由租金的总费用和燃料消耗总费用相加。由表1知船的日租金取决于船型，船型又由冰山的初始体积N0决定，记日租金为m（N0）贝V有4.0 No5 105m(N0)6.2 5 105No1068.0 106N0107又冰山运输总时间为：T _D 9600 40024v24vV所以租金总费用为：由（11）、（8）式可得冰山运输总费用为Y M Q400m(N0)；7.2vv 6 lg等(書0；Sdt)3冰山运到目的地后化成水的体积为：W 0.85N 0.85 （。亍。t）Sdt）3（16）每立方米水所需费用为：(12)(13)m(N)T40m(N0)(14)1dt(15)冰山运输的总费用是

11、由租金的总费用和燃料消耗总费用相加。由表1知船的日首先，对（3）式进行积分运算，有5.模型求解t500400Sdt3v-vt 0.075v 0.225 dt爲0.075v 0.225 dt（18）进而求o o1255300出不同的船速下冰山到达目的地之后总的融化深度，如表4所示：表4船速v12345其次，由（17）式可以发现，这个模型最终由船速V和冰山体积N0决定。由（12）、（13）、（17）式以及表2、表4,通过Matlab计算可以得到每立方米水所需费用, 结果如表5所示：表512345负值负值负值负值负值负值负值负值负值由表5数据可知，当船速v 3km h,N0106时，不符合实际，所得

12、结果为负值，所以，应选择大型拖船N0107，速度v 5km h，则每立方米水所需费用为y 0.0626英镑。因为淡化海水每立方米约为英镑，所以，冰山运输比淡化海水的 成本更低。6.模型评价优点：模型中运用图表分析数据，使各因素间的关系明朗化，运输费用只考虑租费和燃料消耗费，并将冰山体积看成球形，简化了计算，思路清晰简单缺点：模型中只考虑了拖船的租费和燃料消耗费用，没有考虑到影响航行的天气因素和社会因素（如人工费和税收等） ，所以冰山运输方法是否可行还得取决于 综合考虑因素。7.参考文献1姜启源.数学模型（第三版）M.北京：高等教育出版社，2003.2王文波.数学建模及其基础知识详解（第一版）M

13、.武汉大学出版社，2006.3刘玉琏.数学分析（上、下）M.北京：高等教育出版社，1988.8.附 录附录一： a=1 3 5;b= ; %融化速率与船速、与南极距离的关系图 x0=1:5; y0=interp1（a,b,x0,linear）; c=0 1 4 5;d=0 ; x1=0:5; y1=interp1（c,d,x1,linear）; gtext（船速） gtext（与南极距离x1000） a=1 3 5;b= ;%燃料消耗与船速、冰山体积的关系图 x0=1:5; y0=interp1（a,b,x0,linear）; c=5 6 7;d= ;s =%当 v 为4时冰山的融 x1=10

14、A5:1:10A7; y1=interp1(c,d,x1,linear); plot(a,b,r,c,d,-) gtext(船速) gtext(冰山体积 10i) syms k2 b;k2,b=solve(=k2*1+b,=k2*3+b,k2,b) %达式冰山融化速率的表k2 = syms k1;k1=solve(=k1*1000*1+,k1)k1 = syms t,s=int(*t,t,0,500/3)+int(,500/3,400) %融化速率当 v 为1时冰山的s = syms t,s=int(*2/125*t,t,0,250/3)+int(,250/3,200)s = %化速率 sym

15、s t,s=int(*3/125*t,t,0,500/9)+int(,500/9,400/3)s = %化速率当 v 为2时冰山的融当 v 为3时冰山的融 syms t,s=int(*4/125*t,t,0,125/3)+int(,125/3,100)化速率 syms t,s=int(*5/125*t,t,0,100/3)+int(,100/3,80) %山的融化速率s =(400*4+*i nt(log(4*pi/3*(3*10A5/4*pi)A(1/3)F3)-1,0,400)/*4*pi/3*(4*pi/3*(3*10A5/4*piF(1/3)F3)ans =(200*4+* in t(

16、log(4*pi/3*(3*10A5/4*pi)A(1/3)A3)-1,0,200)/*4*pi/3*(4*pi/3*(3*10A5/4*pi)A(1/3)A3)ans =+002 +002i(400/3*4+*int(log(4*pi/3*(3*10A5/4*pi)A(1/3)A3)-1,0,400/3)/*4*pi/3*(4*pi/3*(3*10A5/4*pi)A(1/3)A3)(400/3*4+*int(log(4*pi/3*(3*10A5/4*pi)A(1/3)A3)-1,0,400/3)/*4*pi/3*(4*pi/3*(3*10A5/4*pi)A(1/3)A3)(100*4+288

17、*int(log(4*pi/3*(3*10A5/4*pi)A(1/3)A3)-1,0,100)/*4*pi/3*(4*pi/3*(3*10A5/4*pi)A(1/3)A3)(80*4+396*int(log(4*pi/3*(3*10A5/4*pi)A(1/3)A3)-1,0,80)/*4*pi/3*(4*pi/3*(3*10A5/4*pi)A(1/3)A3)(80*+396*int(log(4*pi/3*(3*10A6/4*pi)A(1/3)A3)-1,0,80)/*4*pi/3*(4*pi/3*(3*10A6/4*pi)A(1/3)A3)ans =当 v 为5时冰(100*+288*i nt

18、(log(4*pi/3*(3*10A6/4*pi)A(1/3)F3)-1,0,100)/*4*pi/3*(4*pi/3*(3*10W/4*pi)A(1/3)A3)ans =(400/3*+*i nt(log(4*pi/3*(3*10A6/4*piF(1/3)F3)-1,0,400/3)/*4*pi/3*(4*pi/3*(3*10A6/4*pi)A(1/3)A3)ans =(400*8+*int(log(4*pi/3*(3*10A7/4*pi)A(1/3)A3)-1,0,400)/*4*pi/3*(4*pi/3*(3*10A7/4*pi)A(1/3)A3)ans =(200*8+*int(log(