数学建模论文范文xxx

时间:2024.4.14

 1 问题的提出 

位于我国西南地区的某个偏远贫困村,年平均降水量不足20mm,是典型的缺水地区。过去村民的日常生活和农业生产用水一方面靠的是每家每户自行建造的小蓄水池,用来屯积每逢下雨时获得的雨水,另一方面是利用村里现有的四口水井。由于近年来环境破坏,经常是一连数月滴雨不下,这些小蓄水池的功能完全丧失。而现有的四口水井经过多年使用后,年产水量也在逐渐减少,在表1中给出它们在近9年来的产水量粗略统计数字。20##年以来,由于水井的水远远不能满足需要,不仅各种农业生产全部停止,而且大量的村民每天要被迫翻山越岭到相隔十几里外去背水来维持日常生活。

    为此,今年政府打算着手帮助该村解决用水难的问题。从两方面考虑,一是地质专家经过勘察,在该村附近又找到了8个可供打井的位置,它们的地质构造不同,因而每个位置打井的费用和预计的年产水量也不同,详见表2,而且预计每口水井的年产水量还会以平均每年10%左右的速率减少。二是从长远考虑,可以通过铺设管道的办法从相隔20公里外的地方把河水引入该村。铺设管道的费用为(万元),其中表示每年的可供水量(万吨/年),表示管道长度(公里)。铺设管道从开工到完成需要三年时间,且每年投资铺设管道的费用为万元的整数倍。要求完成之后,每年能够通过管道至少提供100万吨水。

    政府从20##年开始,连续三年,每年最多可提供60万元用于该村打井和铺设管道,为了保证该村从2010至20##年这五年间每年分别能至少获得150、160、170、180、190万吨水,请作出一个从20##年起三年的打井和铺设管道计划,以使整个计划的总开支尽量节省(不考虑小蓄水池的作用和利息的因素在内)。

1   现有各水井在近几年的产水量(万吨)

2   8个位置打井费用(万元)和当年产水量(万吨)

2 问题的分析

 题中要求制定一个总费用(决策目标)最小的抗旱(打井,铺设管道)方案,属于优化问题,并且使得该村从2010至20##年这五年间每年分别能至少获得150、160、170、180、190万吨水,每年费用不超过60万元。(此两点为主要约束条件)

 其他的约束条件有:

a.每口井只能在20##年开始,连续三年中的其中一年施工

b.铺设管道费用为万元整数倍

c.由于河位于与该村相隔20公里外的地方,所以管道总长度不小于20公里

d.铺设管道需要3年时间,故前3年管道供水量为0,而第4,5年供水量不小于100万吨。

故此模型即为基于以上约束条件的整数规划(最优决策目标)问题。

3 模型的假设

a.忽略小蓄水池的作用和利息因素

b.不考虑意外情况导致所需经费增加

c.假设井在年初修建且时间很短,修完之后即可利用,管道铺好后即可用于供水

d.假设这五年之内村民需水量基本稳定

e.假设井供水量呈稳定规律变化,不考虑其他因素对产水量的影响

f.从长远利益考虑,打井和铺设管道两个方案应同时协调进行

4 符号说明

Xij       0—1变量,表示第i号井在第j年的施工情况,Xij=1第i号井在第j年施工,Xij=0表示不施工

   

Zj       第j年的总费用

Pj       第j年的铺管道费用

Lj       第j 年铺管道公里数

Wj       第j 年的水量

Q     管道供水量

Nj      所有新建的水井在第j年的产水量

5 模型建立

决策变量为三年间铺设管道和打井的总费用。0—1变量Xij表示i号井j 年是否施工,为1则施工,产生费用,Pj表示第j年的铺路费用。所以第j年的总费用Zj=5*X1j+7*X2j+5*X3j+4*X4j+6*X5j+5*X6j+5*X7j+3*X8j+Pj

三年费用min Z=Z1+Z2+Z3

                   =5*X11+7*X21+5*X31+4*X41+6*X51+5*X61+5*X71+3*X81+P1+

5*X12+7*X22+5*X32+4*X42+6*X52+5*X62+5*X72+3*X82+P2+

5*X13+7*X23+5*X33+4*X43+6*X53+5*X63+5*X73+3*X83+P3

约束条件:

1)由于第i号井只能在三年中的某一年打造或者不打造,故应有

 <=1;

2) 每年的费用不能超过计划即

  Z1=5*X11+7*X21+5*X31+4*X41+6*X51+5*X61+5*X71+3*X81+P1;

  Z2=5*X12+7*X22+5*X32+4*X42+6*X52+5*X62+5*X72+3*X82+P2;

  Z3=5*X13+7*X23+5*X33+4*X43+6*X53+5*X63+5*X73+3*X83+P3;

  Z1〈=60,Z2〈=60 , Z3〈=60

3) 每年的水量应满足要求,水量有三部分构成:现有水井的产水量,新建水井的产水量, 管  道铺好后的管道水量。

 现有水井产水量可根据20##——2009数据拟合出20##——20##年的,程序编码及拟合图 见附录1,拟合结果如下图所示:

新建水井产水量:

第一年:

N1=25*X11+36*X21+32*X31+15*X41+31*X51+28*X61+22*X71+12*X81;

第二年:

N2=25*X12+36*X22+32*X32+15*X42+31*X52+28*X62+22*X72+12*X82+25*X11*0.9+36*X21*0.9+32*X31*0.9+15*X41*0.9+31*X51*0.9+28*X61*0.9+22*X71*0.9+12*X81*0.9

第三年:

N3=25*X13+36*X23+32*X33+15*X43+31*X53+28*X63+22*X73+12*X83+25*X12*0.9+36*X22*0.9+32*X32*0.9+15*X42*0.9+31*X52*0.9+28*X62*0.9+22*X72*0.9+12*X82*0.9+25*X11*0.81+   36*X21*0.81+32*X31*0.81+15*X41*0.81+31*X51*0.81+28*X61*0.81+22*X71*0.81+12*X81*0.81;

第四年:N4=N3*0.9

第五年:N5=N3*0.81

管道水量:前三年为0,后两年为Q

故每年的总水量

W1=43.4754+N1

W2=31.7983+N2

W3=24.0045+N3

W4=20.7178+N4+Q

W5=17.4312+N5+Q

满足,W1>=150, W2>=160, W3>=170, W4>=180, W5>=190

4)每年的铺管道费取整且总管道不小20公里即

Pj=0.66Q^0.51*Li

Pj取整

L1+L2+L3 >= 20

6 模型求解

将上述模型输入LINGO可得到【2】

Local optimal solution found.

   Objective value:                              172.0000

   Extended solver steps:                             308

   Total solver iterations:                         10226

                       Variable           Value        Reduced Cost

                             Z1        55.00000            0.000000

                             Z2        60.00000            0.000000

                             Z3        57.00000            0.000000

                            X11        1.000000            5.000000

                            X21        0.000000            7.000000

                            X31        1.000000            5.000000

                            X41        0.000000            4.000000

                            X51        0.000000            6.000000

                            X61        1.000000            5.000000

                            X71        1.000000            5.000000

                            X81        0.000000            3.000000

                             P1        35.00000            1.000000

                            X12        0.000000            5.000000

                            X22        1.000000            7.000000

                            X32        0.000000            5.000000

                            X42        0.000000            4.000000

                            X52        0.000000            6.000000

                            X62        0.000000            5.000000

                            X72        0.000000            5.000000

                            X82        0.000000            3.000000

                             P2        53.00000            1.000000

                            X13        0.000000            5.000000

                            X23        0.000000            7.000000

                            X33        0.000000            5.000000

                            X43        0.000000            4.000000

                            X53        1.000000            6.000000

                            X63        0.000000            5.000000

                            X73        0.000000            5.000000

                            X83        0.000000            3.000000

                             P3        51.00000            1.000000

                              Q        100.0000            0.000000

                             L1        5.064355            0.000000

                             L2        7.668880            0.000000

                             L3        7.379488            0.000000

                             W1        150.4754            0.000000

                             W2        164.0983            0.000000

                             W3        174.0745            0.000000

                             W4        255.7808            0.000000

                             W5        238.9878            0.000000

                            Row    Slack or Surplus      Dual Price

                              1        172.0000           -1.000000

                              2        0.000000           -1.000000

                              3        0.000000           -1.000000

                              4        0.000000           -1.000000

                              5        0.000000            0.000000

                              6        0.000000            0.000000

                              7        0.000000            0.000000

                              8        1.000000            0.000000

                              9        0.000000            0.000000

                             10        0.000000            0.000000

                             11        0.000000            0.000000

                             12        1.000000            0.000000

                             13      -0.4462500E-07        0.000000

                             14      -0.6757500E-07        0.000000

                             15      -0.6502500E-07        0.000000

                             16       0.1127226            0.000000

                             17        5.000000            0.000000

                             18        55.00000            0.000000

                             19        60.00000            0.000000

                             20        0.000000            0.000000

                             21        3.000000            0.000000

                             22        57.00000            0.000000

                             23        0.000000            0.000000

                             24        0.000000            0.000000

                             25        0.000000            0.000000

                             26        0.000000            0.000000

                             27        0.000000            0.000000

                             28        0.000000            0.000000

                             29       0.4754000            0.000000

                             30        4.098300            0.000000

                             31        4.074500            0.000000

                             32        75.78080            0.000000

                             33        48.98782            0.000000

即最小总费用为172万元

第一年花费20万元打造1,3,6,7号井;花费35万元铺管道5.06公里,共计55万元;

第二年花费7万元打造2号井,花费53万元铺管道7.669公里,共计60万元;

第三年花费6万元打造5号井,花费51万元铺管道7.379公里,共计57

7结果分析 

 由结果可知第一年打井1,3,6,7号。产生水量150.47万吨。由各井的产水量可知无论是减少井量,或是替换成其他的井,在保证费用不增加的情况下都会使产水量减小,所以第一年只能打井1,3,6,7号。第二年新增水井2号,总水量164.098,可供替换的井为4,5,7号,与2号水量之差分别为21,5,24皆大于4万吨,故也无法满足水量只能打2号井。同理第三年也只能打5号井。这样方案费用是最小的。

 8方案评价

 1)本文把所解决的问题归结为优化问题,建立的数学模型清晰合理。

 2)运用MATLAB和LINGO软件处理数据和进行运算,降低运算量,简单易行,有很大的可操作性。且所得数据较为合理可靠。

 3)运用0—1模型解题,全面可靠

4)但在实际运用本方案中还应考虑自然因素对产水量的影响,还有需水量的变化,根据实际情况进行灵活改变。

9参考资料

1 姜启源 谢金星  叶俊 《数学模型》,2003

2 穆国旺    MATLAB课件   LINGO课件

3 陈綖     《决策分析》    1987

4杨启帆    《数学建模中的优化问题》 1990

10 附录

附录一:

一号井:     

x=1:1:9

y=[32.2,31.3,29.7,28.6,27.5,26.1,25.3,23.7, 22.7]

plot(x,y,'k.','markersize',25)

a=polyfit(x,y,1)

t=1:1:14

s=polyval(a,t)

hold on

plot(t,s,'r-','linewidth',2)

grid

二号井:  

x=1:1:9

y=[21.5,15.9,11.8,8.7,6.5,4.8,3.5,2.6, 2.0]

plot(x,y,'k.','markersize',25)

a=polyfit(x,y,3)

t=1:1:14

s=polyval(a,t)

hold on

plot(t,s,'r-','linewidth',2)

grid

三号井:   

x=1:1:9

y=[ 27.9,25.8,23.8,21.6,19.5,17.4,15.5,13.3,11.2

]

plot(x,y,'k.','markersize',25)

a=polyfit(x,y,1)

t=1:1:14

s=polyval(a,t)

hold on

plot(t,s,'r-','linewidth',2)

grid


四号井:

x=2:1:9

y=[ 46.2,32.6,26.7,23.0,20.0,18.9,17.5,16.3

]

plot(x,y,'k.','markersize',25)

a=polyfit(x,y,3)

t=2:1:14

s=polyval(a,t)

hold on

plot(t,s,'r-','linewidth',2)

grid

附录二:

min=Z1+Z2+Z3;

Z1=5*X11+7*X21+5*X31+4*X41+6*X51+5*X61+5*X71+3*X81+P1;

Z2=5*X12+7*X22+5*X32+4*X42+6*X52+5*X62+5*X72+3*X82+P2;

Z3=5*X13+7*X23+5*X33+4*X43+6*X53+5*X63+5*X73+3*X83+P3;

X11+X12+X13<=1;

X21+X22+X23<=1;

X31+X32+X33<=1;

X41+X42+X43<=1;

X51+X52+X53<=1;

X61+X62+X63<=1;

X71+X72+X73<=1;

X81+X82+X83<=1;

@bin(X11);@bin(X12);@bin(X13);@bin(X21);@bin(X22);@bin(x23);

@bin(X31);@bin(X32);@bin(X33);@bin(X41);@bin(X42);@bin(X43);

@bin(X51);@bin(X52);@bin(X53);@bin(X61);@bin(X62);@bin(X63);

@bin(X71);@bin(X72);@bin(X73);@bin(X81);@bin(X82);@bin(X83);

P1=0.66*Q^0.51*L1;P2=0.66*Q^0.51*L2;P3=0.66*Q^0.51*L3;

L1+L2+L3>=20;

@gin(P1);@gin(P2);@gin(P3);

Z1<=60;Z1>=0;Z2>=0;Z2<=60;Z3<=60;Z3>=0;

W1=43.4754+25*X11+36*X21+32*X31+15*X41+31*X51+28*X61+22*X71+12*X81;

W2=31.7983+25*X12+36*X22+32*X32+15*X42+31*X52+28*X62+22*X72+12*X82+25*X11*0.9+36*X21*0.9+32*X31*0.9+15*X41*0.9+31*X51*0.9+28*X61*0.9+22*X71*0.9+12*X81*0.9;

W3=24.0045+25*X13+36*X23+32*X33+15*X43+31*X53+28*X63+22*X73+12*X83+25*X12*0.9+36*X22*0.9+32*X32*0.9+15*X42*0.9+31*X52*0.9+28*X62*0.9+22*X72*0.9+12*X82*0.9+25*X11*0.81+36*X21*0.81+32*X31*0.81+15*X41*0.81+31*X51*0.81+28*X61*0.81+22*X71*0.81+12*X81*0.81;

W4=20.7178+(W3-24.0045)*0.9+Q;

W5=17.4312+(W3-24.0046)*0.81+Q;

Q>=100;

W1>=150;W2>=160;W3>=170;W4>=180;W5>=190;

end

附录三:

更多相关推荐:
数学建模论文20篇第一问方法总结

数学建模论文20篇第一问方法总结C1论文问题一:运用多元线性回归模型并借助Matlab软件建立了机械强度模型;运用包维尔法和通用全局优化法,借助1stOpt软件进行数据拟合,建立了茎秆鲜重模型;根据客观实际和类…

如何撰写数学建模论文

摘要(200-300字,包括模型的主要特点、建模方法和主要结果。)关键词(求解问题、使用的方法中的重要术语)内容较多时最好有个目录1。问题重述2。问题分析3。模型假设与约定4。符号说明及名词定义5。模型建立与求…

数学建模论文写作技巧 论文自评(美赛一等奖获得者从获奖论文评述中总结的经验)

论文自评Successfulteamswouldhavetocombineexistingmodels,data,andnewideasincreativeandoriginalways.(成功的队伍会把现有的…

数学建模论文范本

数学建模论文范本:20##-20##数学建模课程考核论文A题饮酒驾车模型摘要:本题根据已知数据,建立血液中酒精含量随时间变化的微分方程模型。并用Mathematica软件提供的非线性拟合的方法求出了一般解。得出…

数学建模竞赛论文范文样本

宜宾学院数模竞赛论文模版宜宾学院第三届大学生数学建模竞赛20xx年5月19日5月28日参赛题目在所选题目上打勾AB参赛编号竞赛组委会填写论文题目摘要1摘要本文解决什么问题解决问题的方法结论提请大家注意摘要应该是...

数学建模论文范文6

饮酒驾车模型摘要交通事故是目前危害人类生命的第一杀手而酒后驾车已经成为引发交通事故的重要原因之一并日益凸现为社会问题因此必须加强有效防控以保障交通安全和秩序长期以来我国酒后驾车现象一直处于较快增长的态势由酒后驾...

数学建模论文范文

论文题目三号黑体字摘要摘要标题是以最恰当最简明的词语反映论文中主要内容的逻辑组合要求反映内容准确得体外延内涵恰如其分用语凝练醒目题目是给评委的第一印象建议将论文所有模型或者算法加入题目中例如用遗传算法解决XXX...

数学建模论文范文[1]

利用数学建模解数学应用题数学建模随着人类的进步科技的发展和社会的日趋数字化应用领域越来越广泛人们身边的数学内容越来越丰富强调数学应用及培养应用数学意识对推动素质教育的实施意义十分巨大数学建模在数学教育中的地位被...

全国大学生数学建模竞赛论文格式规范

全国大学生数学建模竞赛论文格式规范本科组参赛队从AB题中任选一题专科组参赛队从CD题中任选一题全国评奖时每个组别一二等奖的总名额按每道题参赛队数的比例分配但全国一等奖名额的一半将平均分配给本组别的每道题另一半按...

数学建模论文范文

数学建模论文范文利用数学建模解数学应用题数学建模随着人类的进步科技的发展和社会的日趋数字化应用领域越来越广泛人们身边的数学内容越来越丰富强调数学应用及培养应用数学意识对推动素质教育的实施意义十分巨大数学建模在数...

数模论文范文

Ⅰ、问题的重述石油是重要的战略资源,进入新世纪以来石油价格一路高涨且波动频繁,油价成为全球关注的焦点。成品油的合理定价对国家经济发展及社会和谐稳定具有重要的意义,还关系到民生,石油储备等多方面的问题。石油价格的…

数学建模论文范例

数学建模论文小学数学情境教学中的问题与策略AAaa黄海学院20xx6班号摘要新概念理念下数学教育提倡在情境中教学创设适合学生年龄特征的数学教学情境并在情境中解决问题本文主要探讨了小学数学情境教学中存在的若干误区...

数学建模论文(46篇)