数学建模冰山运输

1、拖运冰山的费用估算摘要据资料显示：地球上的水资源有97%被盐化，仅有3%是淡水资源，且大部分是主要分布在南北两极地区的固体冰川，众所周知，水是生命的源泉，是人类赖以生存和发展的重要物质资源之一。而在一些淡水资源非常匮乏的地方，各种产水方式应运而生，本文主要针对以盛产石油著称的波斯湾地区，在采用高成本的海水淡化下，专家建议用拖船从相距9600km外的南极把冰山运到波斯湾地区，取代海水淡化的方法，在拖运过程中采用三种不同的船型、船速会产生不同的费用，为了使获得淡水资源费用合理化。需要对总费用进行合理的估算，为了解决这个问题，我们建立初等函数模型，计算船只以最大运量运一次的总费用和剩余的冰山体积；总

2、费用=日租金天数+总燃料消耗费用，考虑到实际运输过程中，冰山的融化会引起燃料消耗费用的变化，我们应该求出体积变化的临界值（即所运的冰山何时会缩小到原来体积的工），再根据表2可求出对应的消耗原料。剩10余冰山的体积=船的最大运量-运输过程中融化的冰。运输过程中冰的融化速率与冰山到达处与南极的距离有关。根据表3,在特定的船速条件下，将冰山融化速率与南极距离看成线性关系yaxb,根据曲线积分可以求出融化的冰，然后求剩余的冰转换成水的体积，计算每立方米水的成本，经过计算选择大型船，船速为5km/h,每立方米水的成本0.054英镑/m3。关键字：燃料消耗费用冰山融化的体积融化速率临界值;.1、问题重述在

3、以盛产石油著称的波斯湾地区，浩瀚的沙漠覆盖着大地，水资源十分缺乏，不得不采用淡化海水的办法为国民提供用水，成本大约是每立方米0.1英镑，有些专家提出从相距9600km外的南极用拖船运送冰山到波斯湾，以取代淡化海水的办法。在运送冰山的过程中，拖船的租金、运量、燃料消耗及冰山运送过程中融化速率等方面的数据如下：(1)三种拖船的日租金和最大运量如表1所示。表1船型小中大日租金/(英镑)4.06.28.0取大运里/m5105106107(2)燃料消耗(英镑/km),主要依赖于船速和所运冰山的体积，船型的影响可以忽略，如表2所小表2冰山体积/m3、船速/km/h10510610718.410.512.6

4、210.813.516.2313.216.519.8(3)冰山运输过程中的融化速率(m/d),指在冰山与海水接触处每天融化的深度，融化速率除与船速有关外，还与运输过程中冰山到达处与南极的距离有关,这是由于冰山要从南极运往赤道附近的缘故，如表3所示。表3极距离/km船速/km/h01000>40001230000.10.150.20.30.450.6根据所给数据，建立数学模型，解决问题：试选择拖船的船型与船速，使冰山到达目的地后，可以得到的每立方米水所花的费用最低，并与海水淡化的费用相比较2、问题分析此题研究的是每立方米水成本的数学建模问题，题目给出的数据有三种船型，三种船速，那么可分9种

5、情况去求解问题。求解每立方米水的成本就要求解出运输一次的费用和剩余的冰山体积。由于在运输过程中，单位路程燃料消耗与冰山的体积成正相关关系。但是冰山的体积是动态变化的，从而单位路程燃料消耗也随之改变。为了计算方便，利用均值求一段路程的燃料的消耗，而冰山的融化速率与冰山到达处与南极的距离有关，根据图1可以得出当船速确定时，融化速率与距南极距离成线性相关。同样地可采用均值求解，可以确定到达目的地后冰山的剩余体积。将冰转化为水，从而估算出每立方米水的成本。船速E/knih船速1,km/h1000200030004000500060007000S000900010000与南极的距离乐m融化速率与南极距震

6、的关系图船速京km/h图23、模型的假设及符号说明拖船在托运冰山的过程中，有以下假设：假设1、假设拖船航行过程中的船速不变，航行不考虑天气等任何因素的影响,总航行距离9600km假设2、假设冰山形状为球形，球面各点的融化速率相同；假设3、假设冰山到达目的地后，1m3的冰可以融化成0.85m3的水。假设4、假设冰山在运输过程中融化的速率与距离为线性关系。假设5、三种型号的船是按最大运输量拖运冰山。假设6、冰山到达目的地后融化过程中不考虑损耗。假设7、假设拖船所在地就在南极，不用考虑返程费用。符号说明：符合符号说明Vi拖船的速度i1,2,3ZiZ拖船的日租金Wj燃料消耗速度i1,2,3j1,3,5

7、li品巨离南极的品目离i1,2,3uij融化速度i1,2,3R冰山的半径R耗运动过程中冰山融化的深度j1,3,5Rj0,Vj初始时刻冰山的半径j1,3,5最后剩余的冰山体积j1,3,5R末最后时刻冰山的半径j1,3,5ZiZi拖船的总租金LL为总长度mjP全程燃料费用i1,2,3j1,3,5运冰的总费用i1,2,34、模型建立为了求解每立方米水的成本，就要计算出拖船运输一次的总费用和剩余的冰山体积，首先，根据题目所给的冰山融化速率与船速、距南极距离的相关数据，我们可以求解出当船速确定的时候，冰山融化速率与距南极距离成线性关系，利用均值和曲线积分可以求解出船只行驶一段距离与融化冰山的体积的关系。

8、RUL2L3Ui1-Ui2-Ui3ViViR。后（1）R末R0R耗VjRj末0.85其次，由上述计算可以得出燃料消耗与冰山体积的关系，因为冰山体积是动态变化的，所以燃料消耗也是动态的，为了计算的简单化，我们把计算一段路程的燃料消耗简化为固定两个端点的单位路程燃料消耗与行驶路程的乘积和的均值。W1jIW2jIW3jI24ViZi(2)最后：由上述条件利用总费用和冰融化成水后的单位成本:PijpVjVjR末0.855、模型求解首先，计算出不同船型以不同速度拖运冰山所需的总费用。其次，计算出不同船型以不同速度拖运冰山到达目的地，融化成的水。最后，由上述条件求出冰山拖运到目的地融化为水后每立方米的成本

9、，与淡化每立方米水的成本相比较得出成本最低的方案。第一步：计算总费用mijwijliW2jI2W3jI3(4)不同船型以不同船速拖冰山，所需的燃料费：表4,33小5105/m中106/m3大107/m31km/h无意义无意义118780英镑3km/h100800英镑121070英镑155520英镑5km/h120960英镑158202英镑190080英镑乙拖船的总租金：(5)L24vi用不同船型拖运所用的总租金:表5型r5.3小510/m中106/m3大107/m31km/h无意义无意义3200英镑3km/h533.3332英镑826.6665英镑:1066.7英镑5km/h320英镑496英

10、镑640英镑总费用:«Zimj(6)表6r_.53小510/m中106/m3大107/m31km/h无意义无意义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)计算

11、可知用小船拖运时以1km/h到达目的地时冰块已经完全融化，中船拖运时以1km/h到达目的地时冰块已经完全融化所以不再考虑这两种情况。不同船型以不同船速拖冰山，到达目的地剩余体积：R末R。R耗,3小5*10A5/m,3中10A6/m,3大10A7/m1km/h无意义无意义3494840m3km/h31403.5m332236m33194000m5km/h315310m393954m34165900m冰山到达后无其他损耗全部都融化成水的体积：VjRj末0.85表9小5*10人5/m33中10A6/m3大10A7/m1km/h无意义无意义3420614m3km/h31193m327401m32714

12、900m5km/h313014m379861m33541015m每立方米冰的成本piLVj从南极运冰到目的地的每立方米成本为:表10,3小5*10人5/m»3中10A6/m3大10A7/m1km/h无意义无意义30.283英镑/m3km/h384.940英镑/m一,一34.449英镑/m30.058英镑/m5km/h39.320英镑/m31.987英镑/m30.054英镑/m根据计算结果与已知条件（成本每立方米0.1英镑）对比得出0.054英镑/m3的成本价格是最低的，所以我们选择大型船，船速为5km/h6、模型评价优点：1、构造的模型比较简单2、该模型的建立将实际问题模型化，复杂问

13、题简单化；3、我们采用了均值与曲线积分的计算方法，使计算变得相对简单；4、此模型的全局规划比较合理。缺点：1、利用均值求解融化冰山的体积与燃料消耗不准确；2、运输过程融化速率与距南极距离并不成线性关系。7、模型的改进与推广模型改进：虽然得到的最低费用0.054英镑/m3小于海水淡化需要的成本，但是模型中未考虑天气等其它因素会导致拖船过程中成本变高。我们在计算过程中应该以冰山拖运过程中最大损耗、燃料的最大费用，来避免由自然因素带来的船只到达目的地的时间延期增加的费用。方法改进冰山融化速率r(m/d)与船速u(km/h)与南极距离d(km)的关系r是u的线性函数；0d4000时u与d成正比，d40

14、00时u与d无关,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又因为t大冰山的半径为tRRortt0(13)故剩余冰山的体积为Vt3VRt3(14)燃料消耗费用q(英镑/km)z对u线性，对logov线性,zC!uc2igvC3代入数据得出Ci=0.3,c2=6C3="1z3V航行第t

15、大的燃料消耗t-，4，z7.2uu6Ig-(43廉)0而日租金与船的型号有关故f(Vo)4.06.28.0_5V05105105V0106106V0107所以船型运输一次的总费用为400f(v0)T4P-027.2u(u6)Ig一ut13Trt1t1此模型的建立会减少误差，由第一次模型计算的结果,5km/h0的数据代入上述公式，计算结果。参考文献(1)、姜启源、谢金星，叶俊，数学模型(第三版)北京,(15)(16)(17)(18)(19)将大型船，船速为高等教育出版社，2003;(2)、韩忠庚，数学建模方法及其应用，北京高等教育出版社，2005;(3)、吴建国，数学建模案例精编北京，中国水利水

16、电出版社2005;附录模型求解在Matlab中的实现：画图1:船速固定时，冰山体积与燃料的关系图>>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:船速固定时，燃料消耗与冰山体积的关系图> >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')计算过程>>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