运筹学课程设计报告书

1、西安建筑科技大学 课程设计（论文）任务书一、本次课程设计应达到的目的1. 掌握运筹学知识在管理问题中应用的基本方法与步骤；2. 巩固和加深对所学运筹学理论知识及方法的理解与掌握；3. 培养与锻炼学生从管理实践中提炼问题、分析问题、构建模型求解问题的综合应用能力；4. 上机练习，了解与掌握几种常用的运筹学计算软件及其使用与操作方法；5. 锻炼并初步掌握运筹学模型求解程序的编写方法与技术。6. 初步了解学术研究的基本方法与步骤，并通过设计报告的撰写，了解学术报告的写作方法。 二、本次课程设计任务的主要内容和要求1. 结合专业知识，对某一实际管理问题进行分析，调查收集相关数据，并整理出符合问题特征的

2、数据，包括目标因素、约束因素以及必须的参数与系数等等；2. 在上一步分析基础上，按照运筹学建模的基本方法与要求，通过抽象处理，建立所研究问题的运筹学模型，判断模型的类型并选择求解方法；3. 上机练习，学习常用运筹学计算软件的使用与基本操作方法，并选择其中一种对所建运筹学模型进行求解，得出最优解、灵敏度计算等相关计算结果；4. 结合理论课以及计算机程序设计课程所学的基本知识，编写线性规划单纯形法的计算程序，别用所编写程序和已学习的某种运筹学计算软件，并分求解相关课后习题，对所编写的算程序进行验证；5. 总结设计过程，整理与记录设计中的关键工作与成果，撰写设计报告。三、应收集的资料及主要参考文献：

3、 1. 应收集的资料：1研究对象的现状数据材料2与所建模型的参数、系数、约束条件等因素相关的数据材料2. 主要参考文献：1杨茂盛.运筹学(第三版).陕西科学技术出版社，20062运筹学编写组. 运筹学(第三版).清华大学出版社，20053徐玖平, 胡知能, 王緌. 运筹学(第二版). 北京: 科学出版社, 20044胡运权. 运筹学基础及应用. 哈尔滨: 哈尔滨工业大学出版社, 19985陈汝栋，于延荣. 数学模型与数学建模(第2版).国防工业出版社，20096刘建永.运筹学算法与编程实践：delphi实现.清华大学出版社，20047谢金星，薛毅.建优化建模lindo/lingo软件.清华大学

4、出版社，2005设计总说明运筹学在实际生活中的应用主要解决两类问题，而本次设计是研究烟草种植的成本最小化问题。在一个系统中经常遇到两类决策问题：一类是如何利用确定的资源去完成最大的任务，另一类是如何利用最少的资源来完成一个确定的任务。问题的背景是某烟草基地的烟草种植对肥料的合理分配问题，可以有多种方案，但要使成本资金达到最小，节约成本，增加经济效益。从表面现象难以做出准确的决策。这就需要用运筹学中的有关知识，通过对一个问题条件分析，建立模型，求解问题，提高对解进行分析与评价的综合运用能力。并利用lindo计算软件求解问题的操作掌握计算软件的基本操作方法，了解计算机在运筹学中的运用，提高自身分析

5、问题、解决问题的能力。到最后求解出最优的研究路线，来得出问题的最优方案。关键词：运筹学，烟草肥料，最优搭配，成本最小化，lindo 目 录1绪论- 4 -1.1研究背景- 4 -1.2研究意义与目的- 4 -1.3研究的内容- 5 -1.4研究主要方法与思路- 6 -2理论方法的选择- 7 -2.1所研究问题特点- 7 -2.2拟采用的运筹学理论方法特点- 7 -2.3理论方法的实用性及有效性论证- 7 -3模型的建立- 8 -3.1基础数据的确定- 8 -3.2变量的设定- 8 -3.3目标函数的建立- 8 -3.4限制条件的确定- 9 -3.5模型的建立- 10 -4模型的求解及解的分析-

6、 11 -4.1模型的求解- 11 -4.2解的分析及评价- 15 -5总结 - 16 -代码：- 16 -参考文献：- 22 -1绪论1.1研究背景近两年来，在玛纳斯县境内、塔西河冲积扇的中部进行了烤烟的试验和试种。试种结果表明，玛纳斯县种植的g140，nc89，中烟14、15等新品种，均获得了可喜的成果。19851986年全县试种烤烟700亩，收获优质烟叶71490.5公斤，总产值达184157.84元。1987年扩种面积1000亩。预计比往年有较大幅度的增长。经河南省烟草研究所化验，烟叶指标基本上达到了国家规定的质量标准。 而肥料是作物增产的物质基础，合理施用肥料是提高烟叶产量和质量的一

7、项重要技术措施。烤烟施肥的目的不但要提高单位面积上的产量和经济效益，更重要的是要有利于烟叶品质的形成与提高。本课程设计通过对烟草施肥进行研究，运用运筹学的有关知识，对研究数据进行规划、建模实现施肥搭配优化。1.2研究意义与目的研究的意义： 通过本次研究，可以得到此题的最佳解决方案，并可以建立数学模型，把其推广至同类问题，为其他类似问题提供快捷、高效的方法。让烟草种植者再以后遇到类似问题时可以一定的参考文献，可以从一定程度上提高种植者的收入，作到合理的资金分配和最有效的资源的配制。研究的目的:1)巩固和学习曾经学习过的运筹学课程，学会用理论解决实际问题。2)种植花费最小，获得最大收益。 3)通过

8、运筹设计，一定程度上解决种植成本分配，助决策者或决策机构对其所控制的活动进行实现优化决策。加快烟草种植的标准化步伐，增强产品的市场竞争力，提高种植的效益，刺激并提高烟草种植者的种植积极性。4）通过利用运筹学软件求解模型的操作,掌握运筹学计算机软件的基本操作方法。1.3研究的内容通过对资金的计算和生产的需要的分析可以制定出有效的肥料搭配，使农民既可以提高烟草产量又可以合理改善其品质。本次研究针对以下存在问题：（一）、改善烟草的肥料搭配，并且不影响烟草种植新技术的推广工作，提高烟草生产效益。（二）、改善标准化生产水平低的状况。 （三）、不同时期的烟草所用的肥料搭配不同。本次研究课题是：某烟草种植基

9、地所用的化肥是由各种肥料配混而成。包括氮肥、磷肥、钾肥、复合肥，且各种烟草的肥料配比成分及比例是不一样的，这里以每亩烟草为单位。具体需要如下表：表1： 各种肥料及配比关系表种类比例氮肥磷肥钾肥复合肥经济烟草2/92/93/92/9高档烟草2/91/91/31/3烟草苗1/31/301/3虽然各种烟草的肥料比例各不一样，但各种原料每周的最大供应量是确定的。根据近期调查和统计，具体数据如下：表2： 各种原料的单价及供应量表原料种类氮肥磷肥钾肥复合肥单价（元/斤）0.80.032.21.5周供应量（斤）42323050再根据烟草的种类及数目确定每周需要的肥料数量，数据如下表：表3： 各种肥料需求表：

10、烟草种类种类经济烟草高档烟草烟草苗周需求量7447331.4研究主要方法与思路 在现实生活中，如果若干项活动都要使用数量有限的某种资源，或者一项或几项活动生产出同一种资源要供另外一些活动来使用，那么这些活动之间就存在着相互抵赖的关系。这种有限的资源可以是资金、材料、劳力、或生产设备或经营设备。在任何实际活动中，这些资源中总有一些或几种必然要限制利润的最大化或成本的最小化。在这中情况下找出最好的活最优的分配方案就是线性规划的任务。本次研究用运筹学线性规划解决的实际问题是如何使种植费用最小，首先根据自己所需要结局的问题建立线性规划的数学模型，其次对已得模型利用计算机求解这个线性规划，得出最优解，确

11、定符合实际再施于实践。本次研究将采用运筹学中线性规划的有关思想方法，从而取得问题的最优解决方案。主要思路是：从要求和条件入手，分析数据结构，建立恰当的数学模型，用lindo软件在计算机上求解。2理论方法的选择2.1所研究问题特点既然是合理搭配必定会遇到各种线性规划的条件。利用并不是盲目的，而是在一定的条件下进行的。要求成本最小，肥料选用不能为负数。求解过程在lindo软件上进行计算。2.2拟采用的运筹学理论方法特点将采用线性规划的思想方法对此题求解。线性规划是运筹学中发展最完善，并且应用最广泛的一个分支，其研究的主要对象有：一类是给定了人力、物力资源，研究如何用这些资源完成任务，另一类是研究如

12、何统筹安排，尽量以最少的人力、物力资源完成该项任务。在本设计中研究的是资金一定的条件下，合理规划资金，使烟农种植费用花费最少。在本次设计中，我利用到了处理线性规划所需的电脑软件：lindo，lindo一种专门用于求解数学规划问题的软件包。由于lindo执行速度很快、易于方便输入、求解和分析数学规划问题。因此在数学、科研和工业界得到广泛的应用。lindo主要用于解线性规划、非线性规划、二次规划和整数规划等问题。也可以用于一些非线性和线性方程组得球界以及代数方程求根等。lindo中包括了一种建模语言和许多常用的数学函数（包括大量概论函数），可供使用者建立规划问题是调用。2.3理论方法的实用性及有效

13、性论证此题的特点与运筹学中线性规划的第二类问题的特点完全一致，故采用线性规划求解此题是可行的。3模型的建立3.1基础数据的确定 各种肥料配比关系及供应和需求关系3.2变量的设定为了方便直观起见，我们将表1、表2和表3放在一起，来确定变量。表4： 各种肥料配比关系及供应、需求关系表原料单价（元/斤）周供应量（斤）肥料配比关系经济烟草a1高档烟草a2烟草苗 a3氮肥 b10.8422/92/91/3磷肥b20.03322/91/91/3钾肥 b32.2301/31/30化肥 b41.5502/91/31/3各种肥料周需求量744733设配制m1斤a1，m2斤a2和m3斤a3，所用的b1、b2、b3

14、、b4四种原料数量分别为xij（i=1，2，3；j=1，2，3，4）斤。再设z，z配料的总成本。其中，z，xij的取值范围都是大于或等于0。3.3目标函数的建立 根据问题的特点，我们要求的目标函数是使z=0.8（x11+x21+x31）+0.03（x12+x22+x32）+2.2（x13+x23+x33）+1.5（x14+x24+x34）的最小的解。这里的是单位（一千克）最优目标，因为种植肥料的供应量总量一定，需要确定是配比的比例，通过确定最优配比比例达到种植肥料的最低花费。3.4限制条件的确定 假设混合调配的过程中没有重量损失，可以得出以下式子：x11+x12+x13+x14=m1 （1-1

15、）x21+x22+x23+x24=m2 （1-2）x31+x32+x33+x34=m3 （1-3）由配比关系可得：x112/9m1，x122/9m1，x131/3m1，x142/9m1 （1-4）x212/9m2，x221/9m2，x231/3m2，x241/3m2 （1-5）x31m3，x322/3m3，x33=0m3，x341/3m3 （1-6）在式（1-1）两边同乘以2/9得：2/9x11+2/9x12+2/9x13+2/9x14=2/9m1又因（1-4）的第一项为：x112/9m1所以有：2/9x11+2/9x12+2/9x13+2/9x14x11移项整理得：-7/9x11+2/9x1

16、2+2/9x13+2/9x140同理将式（1-1）（1-2）（1-3）分别与式（1-4）（1-5）（1-6）相结合，得其它约束条件：5/9x11-4/9x12+5/9x13+5/9x140-2/9x11-2/9x12-2/9x13+7/9x140-7/9x21+2/9x22+2/9x23+2/9x2401/9x21-8/9x22+1/9x23+1/9x240-2/3x21-2/3x22+1/3x23-2/3x240-2/3x31+1/3x32-2/3x33-3/3x3401/3x31+1/3x32+1/3x33-2/3x340因为供应量也是有限的，所以还有：x11+x21+x3142x12+x

17、22+x3232x13+x23+x3330x14+x24+x3450又因为需求量也是有底线的，所以：x11+x12+x13+x1474x21+x22+x23+x2447x31+x32+x33+x3433以及：x31=0，x13=0，x33=0，x24=0和xij03.5模型的建立综上所述，得该问题的线性规划模型min z=0.8（x11+x21+x31）+0.03（x12+x22+x32）+2.2（x13+x23+x33）+1.5（x14+x24+x34）s.t. -7/9x11+2/9x12+2/9x13+2/9x1405/9x11-4/9x12+5/9x13+5/9x140-2/9x11-

18、2/9x12-2/9x13+7/9x140-7/9x21+2/9x22+2/9x23+2/9x2401/9x21-8/9x22+1/9x23+1/9x240-2/3x21-2/3x22+1/3x23-2/3x240-2/3x31+1/3x32-2/3x33-3/3x3401/3x31+1/3x32+1/3x33-2/3x340x11+x21+x3142x12+x22+x3232x13+x23+x3330x14+x24+x3450x11+x12+x13+x1474x21+x22+x23+x2447x31+x32+x33+x3433x31=0，x13=0，x33=0，x24=0xij04模型的求解

19、及解的分析4.1模型的求解在此由于约束条件及变量比较多，我们运用lindo软件来进行计算。lindo 是一个专门对线性规划问题求解的软件，它的运用简化了人工计算时的复杂性。在原来的约束条件经过整理后在lindo输入：min 0.53x11+0.53x21+0.53x31+0.76x12+0.76x22+0.76x32+0.3x13+0.3x23+0.3x33+1.5x14+1.5x24+1.5x34st-0.778x11+0.222x12+0.222x13+0.222x14=00.556x11-0.444x12+0.556x13+0.556x14=0-0.222x11-0.222x12-0.2

20、22x13+0.778x14=0-0.778x21+0.222x22+0.222x23+0.222x24=00.111x21-0.889x22+0.111x23+0.111x24=0-0.667x21-0.667x22+0.333x23-0.667x24=0-0.667x31+0.333x32-0.667x33-0.667x34=00.333x31+0.333x32+0.333x33-0.667x34=0x11+x21+x31=120x12+x22+x32=320x13+x23+x33=60x14+x24+x34=420x21+x22+x23+x24=94.5x31+x32+x33+x34=1

21、07.1x11=0x12=0x13=0x14=0x21=0x22=0x23=0x24=0x31=0x32=0x33=0x34=0end运行后的结果为： lp optimum found at step 10 objective function value 1) 507.2000variable value reduced cost x11 99.020996 0.000000 x21 20.979000 0.000000 x31 0.000000 0.000000 x12 306.479004 0.000000 x22 13.521000 0.000000 x32 0.000000 0.00

22、0000 x13 0.000000 0.000000 x23 60.000000 0.000000 x33 0.000000 0.000000 x14 14.500000 0.000000 x24 0.000000 0.000000 x34 107.099998 0.000000 row slack or surplus dual prices 2) 5.780999 0.000000 3) 72.959000 0.000000 4) 78.739998 0.000000 5) 0.000000 0.000000 6) 3.031500 0.000000 7) 3.031500 0.00000

23、0 8) 71.435699 0.000000 9) 71.435699 0.000000 10) 0.000000 0.970000 11) 0.000000 0.740000 12) 0.000000 1.200000 13) 8.400000 0.000000 14) 0.000000 -1.500000 15) 0.000000 -1.500000 16) 0.000000 -1.500000 17) 99.020996 0.000000 18) 306.479004 0.000000 19) 0.000000 0.000000 20) 14.500000 0.000000 21) 2

24、0.979000 0.000000 22) 13.521000 0.000000 23) 60.000000 0.000000 24) 0.000000 0.000000 25) 0.000000 0.000000 26) 0.000000 0.000000 27) 0.000000 0.000000 28) 107.099998 0.000000 no. iterations= 10 ranges in which the basis is unchanged: obj coefficient ranges variable current allowable allowable coef

25、increase decrease x11 0.530000 0.000000 infinity x21 0.530000 infinity 0.000000 x31 0.530000 infinity infinity x12 0.760000 0.000000 0.000000 x22 0.760000 0.000000 1.200000 x32 0.760000 infinity 0.000000 x13 0.300000 infinity infinity x23 0.300000 1.200000 infinity x33 0.300000 infinity infinity x14

26、 1.500000 0.000000 0.000000 x24 1.500000 infinity 0.000000 x34 1.500000 0.000000 1.500000 righthand side ranges row current allowable allowable rhs increase decrease 2 0.000000 infinity 5.780999 3 0.000000 infinity 72.959000 4 0.000000 infinity 78.739998 5 0.000000 20.979000 3.031500 6 0.000000 infi

27、nity 3.031500 7 0.000000 infinity 3.031500 8 0.000000 infinity 71.435699 9 0.000000 infinity 71.435699 10 120.000000 14.500000 5.780999 11 320.000000 14.500000 8.400000 12 60.000000 3.031500 8.400000 13 130.000000 infinity 8.400000 14 420.000000 8.400000 14.500000 15 94.500000 8.400000 4.544977 16 1

28、07.099998 8.400000 107.099998 17 0.000000 99.020996 infinity 18 0.000000 306.479004 infinity 19 0.000000 60.000000 0.000000 20 0.000000 14.500000 infinity 21 0.000000 20.979000 infinity 22 0.000000 13.521000 infinity 23 0.000000 60.000000 infinity 24 0.000000 0.000000 0.000000 25 0.000000 5.780999 0

29、.000000 26 0.000000 0.000000 infinity 27 0.000000 60.000000 0.000000 28 0.000000 107.099998 infinity4.2解的分析及评价 由此可见，计算机在7次迭代后算出每个月的最少花费为507.2元，这时：x11=99.021， x21=20.979， x31=0，x12=235.043， x22=13.521， x32=71.436，x13=0，x23=60， x33=0，x14=85.936， x24=0， x34=35.664通过以上计算结果可以看出，最优值符合实际情况，如果种植户按这个结果来配料，将在

30、满足肥料需求及其供给的情况下，使自己的利益得到最大化的满足。5 总结 通过完成本次课程设计，将所学的运筹学知识进行了简单实现。从题目的选定、背景的调研、实际数据的查证，到问题分析、问题建模和最后的解决问题，巩固了运筹知识。在使用计算机软件解决问题时，能将实际问题联系到计算机上，提高了解决问题的效率。通过lindo的使用，加强了运筹问题的计算能力。在使用程序编译并将问题实现后，给以后遇到类似问题时可以直接通过程序将其解决。在调研烟草数据过程中，更加意识到实践的重要性。本次设计不仅是巩固了运筹知识，更是开阔了视野，初步了解烟草行业。在实现设计过程中遇到过复杂问题，比如问题如何分析、模型如何建立，都

31、是本次设计的难点。在查阅资料后虽解决的部分问题，可模型还有不严整的地方，在以后建模中将会加强。代码：#include#include#define m 3 /*定义约束条件方程组的个数*/#define n 5 /*定义未知量的个数*/float m=1000000.0;float amn; /*用于记录方程组的数目和系数;*/float cn; /*用于存储目标函数中各个变量的系数*/float bm; /*用于存储常约束条件中的常数*/float cbm; /*用于存储基变量的系数*/float setam; /*存放出基与入基的变化情况*/float deltan; /*存储检验数矩阵*

32、/float xn;int numm; /*用于存放出基与进基变量的情况*/float zb=0; /*记录目标函数值*/void input();void print();int danchunxing1();int danchunxing2(int a);void danchunxing3(int a,int b);int danchunxing1()int i,k=0;int flag=0;float min=0;for(i=0;i=0)flag=1;else flag=0;break;if(flag=1)return -1;for(i=0;ideltai) min=deltai;k=i

33、;return k;int danchunxing2(int a)int i,k,j;int flag=0;float min;k=a;for(i=0;im;i+)if(aik=0)flag=1;else flag=0;break;if(flag=1)printf(n该线性规划无最优解!n); return -1;for(i=0;i0) setai=bi/aik; else setai=m;min=m;for(i=0;i=setai) min=setai;j=i;numj=k+1;cbj=ck;return j;void danchunxing3(int p,int q)int i,j,c,l

34、;float temp1,temp2,temp3;c=p;/*行号*/l=q;/*列号*/temp1=acl;bc=bc/temp1;for(j=0;jn;j+)acj=acj/temp1; for(i=0;im;i+) if(i!=c) if(ail!=0) temp2=ail; bi=bi-bc*temp2; for(j=0;jn;j+) aij=aij-acj*temp2; temp3=deltal;for(i=0;in;i+)deltai=deltai-aci*temp3;void print()int i,j=0;printf(n-n);for(i=0;im;i+) printf(%8.2ftx(%d) %8.2f ,cbi,numi,bi); for(j=0;jn;j+) printf(%8.2f ,aij); printf(n);printf(n-n);printf(ttt);for(i=0;in;i+)printf( %8.2f,deltai);printf(n-n);void input()int i,j