工业大学课程设计报告课程设计名称运筹课程设计专业班级学生姓名指导教师2013年6月28日课程设计任务书运筹学课程设计报告组别：第十二组设计人员：设计时间：2013年6月17日至2013年6月28日1.设计进度本课程设计时间分为两周：1.1 第一周（2013年6月17日----2013年6月21日）：建模阶段。
此阶段各小组根据给出的题目完成模型的建立。
主要环节包括：1.1.16月17日上午：发指导书；按组布置设计题目；说明进度安排。
1.1.2 6月18日下午至6月20日：各小组审题，查阅资料，进行建模前的必要准备(包括求解程序的编写与查找)。
1.1.3 6月21日至6月22日：各个小组进行建模，并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。
1.2 第二周（2013年6月24日---6月28日）：上机求解，结果分析及答辩。
主要环节包括:1.2.1 6月24日至6月25日：上机调试程序1.2.2 6月26日：完成计算机求解与结果分析。
1.2.3 6月27日：撰写设计报告。
1.2.4 6月28日：设计答辩及成绩评定。
2.设计题目第三十题某糖果厂用原料A、B、C加工成三种不同品牌的糖果甲、乙、丙。
已知各种品牌糖果中A、B、C的含量，原料成本，各种原料的每月限制用量，三种牌号糖果的单位加工费及售价如表1所示。
问该厂每月生产这三种糖果各多少公斤，使该厂获利最大。
并按要求分别完成下列分析：（1）乙产品的售价在何范围内变化时最优生产方案不变？（2）B原料的成本在何范围内变化时最优生产方案不变？（3）C原料的每月限制量在何范围内变化时最优基不变？（4）甲产品的加工费在何范围内变化时最优生产方案不变？表13.建模过程3.1分析过程用i=1,2,3分别代表原材料A、B、C，用j=1,2,3分别表示甲、乙、丙三种糖果。
设Xij为生产第j种糖果使用的第i种原料的公斤数。
甲糖果的质量为Y1乙糖果的质量为Y2丙糖果的质量为Y3生成甲糖果使用的A原料的公斤数为：X11生成甲糖果使用的B原料的公斤数为：X12生成甲糖果使用的C原料的公斤数为：X13生成乙糖果使用的A原料的公斤数为：X21生成乙糖果使用的B原料的公斤数为：X22生成乙糖果使用的C原料的公斤数为：X23生成丙糖果使用的A原料的公斤数为：X31生成丙糖果使用的B原料的公斤数为：X32生成丙糖果使用的C原料的公斤数为：X33X11+X21+X31 =Y1表示A、B、C三种原料质量之和为甲糖果的质量。
X11>=0.6Y1，化简为0.6X21+0.6X31-0.4X11<=0表示甲糖果中A原料含量要求大于或等于甲糖果总重量的60%。
X31<=0.2Y1化简为0.8X31-0.2X11-0.2X21<=0表示甲糖果中C原料含量要求小于或等于甲糖果总重量的20%。
X12+X22+X32=Y2表示A、B、C三种原料质量之和为乙糖果的质量。
X12>=0.3Y2化简为0.3X22+0.3X32-0.7X12<=0表示乙糖果中A原料含量要求大于或等于乙糖果总重量的30%。
X32<=0.5Y2化简为0.5X32-0.5X12-0.5X22<=0表示乙糖果中C原料含量要求小于或等于乙糖果总重量的50%。
X13+X23+X33=Y3表示A、B、C三种原料质量之和为丙糖果的质量。
X33<=0.6Y3化简为0.4X33-0.6X13-0.6X23<=0表示丙糖果中C原料含量要求小于或等于丙糖果总重量的60%。
三种原料的限制用量第一种原料的资源限量X11+X12+X13<=2000第二种原料的资源限量X21+X22+X23<=2500第三种原料的资源限量X31+X32+X33<=1200而此时的总利润为Z=（3.4-0.5）Y1+（2.85-0.4）Y2+（2.25-0.3）Y3-2（X11+X12+X13）-1.5（X21+X22+X23））-1（X31+X32+X33）3.2模型该问题的LP模型：Max Z =0.9X11+1.4 X21+1.9 X31+0.45X12+0.95 X22+1.45 X32 -0.05X13 +0.45X23-0.95X33 X11+X12+X13<=2000X21+X22+X23<=2500X31+X32+X33<=12000.6X21+0.6X31-0.4X11<=00.8X31-0.2X11-0.2X21<=00.3X22+0.3X32-0.7X12<=00.5X32-0.5X12-0.5X22<=00.4X33-0.6X13-0.6X23<=0X ij>=0 (i=1,2,3, j=1,2,3,)4.求解程序功能介绍4．1 程序流程图4.2 C语言程序介绍C 语言具有各种各样的数据类型,可使程序效率更高。
并且计算功能、逻辑判断功能也比较强大, 可以实现决策目的。
这种结构化方式可使程序层次清晰, 便于使用、维护以及调试。
C语言是以函数形式提供给用户的,这些函数可方便的调用,并具有多种循环、条件语句控制程序流向,从而使程序完全结构化。
4.2.1 数据录入（1）输入方程组的系数矩阵A（8行17列）：（2）输入初始基变量的数字代码num矩阵：（3）输入方程组右边的值矩阵b：（4）输入目标函数各个变量的系数所构成的系数阵C：图 14.3 利用LINDO软件进行灵敏度分析4.3.1 LINDO软件介绍LINDO是一种专门用于求解数学规划问题的软件包。
由于LINDO执行速度快，易于方便地输入、求解和分析数学规划问题，因此在教学、科研和工业界得到广泛应用。
LINDO主要用于求解线性规划、非线性规划、二次规划和整数规划等问题，也可以用于一些线性和非线性方程组的求解以及代数方程求根等。
LINDO中包含了一种建模语言和许多常用的数学函数(包括大量概率函数)，可供使用者建立数学规划问题模型时调用。
4.3.2数据录入Max 0.9X11+1.4X21+1.9X31+0.45X12+0.95X22+1.45X32 -0.05X13+0.45X23-0.95X33 stX11+X12+X13<=2000X21+X22+X23<=2500X31+X32+X33<=12000.6X21+0.6X31-0.4X11<=00.8X31-0.2X11-0.2X21<=00.3X22+0.3X32-0.7X12<=00.5X32-0.5X12-0.5X22<=00.4X33-0.6X13-0.6X23<=0end4.3.3 注意事项：. 目标函数及各约束条件之间一定要有“Subject to (ST) ”分开。
. 变量名不能超过8个字符。
. 变量与其系数间可以有空格，单不能有任何运算符（如乘号“＊”等）。
. 要输入<=或>=约束，相应以<或>代替即可。
. 一般LINDO中不能接受括号“（）”和逗号“，”，例：400(X1+X2) 需写成400X1+400X2；10,000需写成10000。
5．结果分析5.1 C语言程序输出结果图 2根据所求最优解得知：甲糖果的质量为Y1 =X11+X21+X31=966.667乙糖果的质量为Y2=X12+X22+X32=4733丙糖果的质量为Y3=X13+X23+X33=05.2 LINDO分析的结果5.2.1 OBJECTIVE FUNCTION VALUE （给出目标函数的最优值）5450.00 （目标函数的最优值为 5450.00）VARIABLE（变量）VALUE （变量值）REDUCED COST（检验数所在行变量系数） X11 580.000000 0.000000X21 386.666656 0.000000X31 0.000000 0.000000X12 1420.000000 0.000000X22 2113.333252 0.000000X32 1200.000000 0.000000X13 0.000000 1.550000X23 0.000000 0.050000X33 0.000000 1.9500005.2.2 SLACK OR SURPLUS：给出松弛变量的值。
列出最优单纯形表中判别数所在行的松弛变量的系数ROW SLACK OR SURPLUS DUAL PRICES2) 0.000000 1.5000003) 0.000000 0.5000004) 0.000000 1.0000005) 0.000000 1.5000006) 193.333328 0.0000007) 0.000000 1.5000008) 1166.666626 0.0000009) 0.000000 0.0000005.2.3 RANGES IN WHICH THE BASIS IS UNCHANGED: 当目标函数的变量系数什么变化范围内时，最优基不变。
CURRENT COEF：初始目标函数的系数。
ALLOWABLE INCREASE：允许变量系数增加的范围。
ALLOWABLE DECREASE：允许变量系数减少的范围。
当目标函数的系数C在 [初始目标函数的系数-允许变量系数减少的范围，初始目函数的系数+允许变量系数增加的范围] 内变化时，最优基不变，最优解也不变，由于目标函数的系数发上改变了，所以最优值有可能改变。
如：目标函数中X1的变量系数为0.9，当它在 [0.15，0.98]内变化时，最优基保持不变。
Objective Coefficient Ranges:VARIABLE CURRENT ALLOWABLE ALLOWABLECOEF INCREASE DECREASEX11 0.900000 0.083333 0.750000X21 1.400000 0.125000 0.000000X31 1.900000 0.000000 INFINITYX12 0.450000 0.750000 0.083333X22 0.950000 0.000000 0.035714X32 1.450000 INFINITY 0.000000X13 -0.050000 1.550000 INFINITYX23 0.450000 0.050000 INFINITYX33 -0.950000 1.950000 INFINITY5.2.4 约束条件右端项在什么范围内变化时，最优基不变：CURRENT RHS：初始约束条件右端项的值；ALLOWABLE INCREASE：允许b值增加的范围；ALLOWABLE DECREASE：允许b值减少的范围。