数模论文的撰写方法1. 题目2.摘要3. 问题重述4. 问题分析5. 模型假设与约定6. 符号说明及名词定义7. 模型建立与求解①补充假设条件，明确概念，引进参数; ②模型形式(可有多个形式的模型);8. 进一步讨论(参数的变化、假设改变对模型的影响)9. 模型检验(使用数据计算结果，进行分析与检验)10. 模型优缺点(改进方向，推广新思想)11. 参考文献及参考书籍和网站12.附录(计算程序，框图;各种求解演算过程，计算中间结果;各种图形、表格。

)下面是范例：1 问题的提出位于我国西南地区的某个偏远贫困村，年平均降水量不足20mm ，是典型的缺水地区。

过去村民的日常生活和农业生产用水一方面靠的是每家每户自行建造的小蓄水池，用来屯积每逢下雨时获得的雨水，另一方面是利用村里现有的四口水井。

由于近年来环境破坏，经常是一连数月滴雨不下，这些小蓄水池的功能完全丧失。

而现有的四口水井经过多年使用后，年产水量也在逐渐减少，在表1中给出它们在近9年来的产水量粗略统计数字。

2009年以来，由于水井的水远远不能满足需要，不仅各种农业生产全部停止，而且大量的村民每天要被迫翻山越岭到相隔十几里外去背水来维持日常生活。

为此，今年政府打算着手帮助该村解决用水难的问题。

从两方面考虑，一是地质专家经过勘察，在该村附近又找到了8个可供打井的位置，它们的地质构造不同，因而每个位置打井的费用和预计的年产水量也不同，详见表2，而且预计每口水井的年产水量还会以平均每年10%左右的速率减少。

二是从长远考虑，可以通过铺设管道的办法从相隔20公里外的地方把河水引入该村。

铺设管道的费用为L 66Q .0P 0.51（万元），其中Q 表示每年的可供水量（万吨/年），L 表示管道长度（公里）。

铺设管道从开工到完成需要三年时间，且每年投资铺设管道的费用为万元的整数倍。

要求完成之后，每年能够通过管道至少提供100万吨水。

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

表1 现有各水井在近几年的产水量（万吨）表2 8个位置打井费用（万元）和当年产水量（万吨）2 问题的分析题中要求制定一个总费用（决策目标）最小的抗旱（打井，铺设管道）方案，属于优化问题，并且使得该村从2010至2014年这五年间每年分别能至少获得150、160、170、180、190万吨水，每年费用不超过60万元。

（此两点为主要约束条件）其他的约束条件有：a.每口井只能在2010年开始，连续三年中的其中一年施工b.铺设管道费用为万元整数倍c.由于河位于与该村相隔20公里外的地方,所以管道总长度不小于20公里d.铺设管道需要3年时间,故前3年管道供水量为0,而第4,5年供水量不小于100万吨。

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

3 模型的假设a．忽略小蓄水池的作用和利息因素b．不考虑意外情况导致所需经费增加c．假设井在年初修建且时间很短，修完之后即可利用，管道铺好后即可用于供水d．假设这五年之内村民需水量基本稳定e．假设井供水量呈稳定规律变化，不考虑其他因素对产水量的影响f．从长远利益考虑，打井和铺设管道两个方案应同时协调进行4 符号说明X ij 0—1变量，表示第i号井在第j年的施工情况，X ij=1第i号井在第j年施工，X ij=0表示不施工Z j 第j年的总费用P j 第j年的铺管道费用L j 第j 年铺管道公里数W j 第j 年的水量Q 管道供水量N j 所有新建的水井在第j 年的产水量5 模型建立决策变量为三年间铺设管道和打井的总费用。

0—1变量X ij 表示i 号井j 年是否施工，为1则施工,产生费用，P j 表示第j 年的铺路费用。

所以第j 年的总费用Z j =5*X 1j +7*X 2j +5*X 3j +4*X 4j +6*X 5j +5*X 6j +5*X 7j +3*X 8j +P j三年费用min Z=Z 1+Z 2+Z 3=5*X 11+7*X 21+5*X 31+4*X 41+6*X 51+5*X 61+5*X 71+3*X 81+P 1+ 5*X 12+7*X 22+5*X 32+4*X 42+6*X 52+5*X 62+5*X 72+3*X 82+P 2+ 5*X 13+7*X 23+5*X 33+4*X 43+6*X 53+5*X 63+5*X 73+3*X 83+P 3约束条件：1） 由于第i 号井只能在三年中的某一年打造或者不打造，故应有∑=31j Xij <=1;2) 每年的费用不能超过计划即Z 1=5*X 11+7*X 21+5*X 31+4*X 41+6*X 51+5*X 61+5*X 71+3*X 81+P 1; Z 2=5*X 12+7*X 22+5*X 32+4*X 42+6*X 52+5*X 62+5*X 72+3*X 82+P 2; Z 3=5*X 13+7*X 23+5*X 33+4*X 43+6*X 53+5*X 63+5*X 73+3*X 83+P 3; Z 1〈=60，Z 2〈=60 ， Z 3〈=603) 每年的水量应满足要求，水量有三部分构成：现有水井的产水量，新建水井的产水量， 管 道铺好后的管道水量。

现有水井产水量可根据2001——2009数据拟合出2010——2014年的，程序编码及拟合图 见附录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*+36*X21*+32*X31*+15 *X41*+31*X51*+28*X61*+22*X71*+12*X81*第三年：N3=25*X13+36*X23+32*X33+15*X43+31*X53+28*X63+22*X73+12*X83+25*X12*+36*X22*+32*X32*+15 *X42*+31*X52*+28*X62*+22*X72*+12*X82*+25*X11*+36*X21*+32*X31*+15*X41*+31*X51*+28*X61*+22*X71*+12*X81*;第四年：N4=N3*第五年：N5=N3*管道水量：前三年为0，后两年为Q故每年的总水量W1=+N1W2=+N2W3=+N3W4=+N4+QW5=+N5+Q满足，W1>=150, W2>=160, W3>=170, W4>=180, W5>=1904)每年的铺管道费取整且总管道不小20公里即P j=^*L iP j取整L1+L2+L3 >= 206 模型求解将上述模型输入LINGO可得到【2】Local optimal solution found.Objective value:Extended solver steps: 308Total solver iterations: 10226Variable Value Reduced CostZ1Z2Z3X11X21X31X41X51X61X71X81P1X12 X22 X32 X42 X52 X62 X72 X82 P2 X13 X23 X33 X43 X53 X63 X73 X83 P3 Q L1 L2 L3 W1 W2 W3 W4 W5Row Slack or Surplus Dual Price 123456789101112131415161718192021222324252627282930313233即最小总费用为172万元第一年花费20万元打造1,3,6,7号井；花费35万元铺管道公里，共计55万元；第二年花费7万元打造2号井，花费53万元铺管道公里，共计60万元；第三年花费6万元打造5号井，花费51万元铺管道公里，共计577结果分析由结果可知第一年打井1，3，6，7号。

产生水量万吨。

由各井的产水量可知无论是减少井量，或是替换成其他的井，在保证费用不增加的情况下都会使产水量减小，所以第一年只能打井1，3，6，7号。

第二年新增水井2号，总水量，可供替换的井为4，5，7号，与2号水量之差分别为21，5，24皆大于4万吨，故也无法满足水量只能打2号井。

同理第三年也只能打5号井。

这样方案费用是最小的。

8方案评价1）本文把所解决的问题归结为优化问题，建立的数学模型清晰合理。

2）运用MATLAB和LINGO软件处理数据和进行运算，降低运算量，简单易行，有很大的可操作性。

且所得数据较为合理可靠。

3）运用0—1模型解题，全面可靠4）但在实际运用本方案中还应考虑自然因素对产水量的影响，还有需水量的变化，根据实际情况进行灵活改变。

9参考资料1 姜启源谢金星叶俊《数学模型》，20032 穆国旺MATLAB课件LINGO课件3 陈綖《决策分析》19874杨启帆《数学建模中的优化问题》199010 附录附录一：一号井：x=1:1:9y=[,,,,,,,, ]plot(x,y,'k.','markersize',25)a=polyfit(x,y,1)t=1:1:14s=polyval(a,t)hold onplot(t,s,'r-','linewidth',2)grid二号井：x=1:1:9y=[,,,,,,,, ]plot(x,y,'k.','markersize',25) a=polyfit(x,y,3)t=1:1:14s=polyval(a,t)hold onplot(t,s,'r-','linewidth',2) grid 一号井水量模拟图线（万吨）（年份减去2000）三号井：x=1:1:9y=[ ,,,,,,,,]plot(x,y,'k.','markersize',25) a=polyfit(x,y,1)t=1:1:14s=polyval(a,t)hold onplot(t,s,'r-','linewidth',2) grid 二号井水量模拟图线（万吨）（年份减去2000）四号井：x=2:1:9y=[ ,,,,,,,]plot(x,y,'k.','markersize',25) a=polyfit(x,y,3)t=2:1:14s=polyval(a,t)hold onplot(t,s,'r-','linewidth',2) grid 三号井水量模拟图线（万吨）（年份减去2000）（万吨）（年份减去2000）四号井水量模拟图线附录二：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=*Q^*L1;P2=*Q^*L2;P3=*Q^*L3;L1+L2+L3>=20;@gin(P1);@gin(P2);@gin(P3);Z1<=60;Z1>=0;Z2>=0;Z2<=60;Z3<=60;Z3>=0;W1=+25*X11+36*X21+32*X31+15*X41+31*X51+28*X61+22*X71+12*X81;W2=+25*X12+36*X22+32*X32+15*X42+31*X52+28*X62+22*X72+12*X82+25*X11*+36*X21*+32*X3 1*+15*X41*+31*X51*+28*X61*+22*X71*+12*X81*;W3=+25*X13+36*X23+32*X33+15*X43+31*X53+28*X63+22*X73+12*X83+25*X12*+36*X22*+32*X3 2*+15*X42*+31*X52*+28*X62*+22*X72*+12*X82*+25*X11*+36*X21*+32*X31*+15*X41*+31*X 51*+28*X61*+22*X71*+12*X81*;W4=+（）*+Q;W5=+（）*+Q;Q>=100;W1>=150;W2>=160;W3>=170;W4>=180;W5>=190;end附录三：。