(上机实验WINQSB)运筹学上机指导手册.doc

《(上机实验WINQSB)运筹学上机指导手册.doc》由会员分享，可在线阅读，更多相关《(上机实验WINQSB)运筹学上机指导手册.doc（17页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、一、软件下载、安装1、下载地址：ftp:/2、将文件夹WinQSB拷贝到硬盘打开硬盘中的文件夹WinQSB运行Set.up文件安装程序二、线性规划、整数规划、0-1规划上机程序1、运行“Linear and integer programming”，出现图1所示界面图 12、运行file菜单下的new problem 命令，出现图2所示界面。图 2问题名称决策变量个数约束条件个数（不含变量约束）目标函数类型数据类型输入数据格式：选择Spreadsheet Matrix From非负连续变量非负整数变量0-1整数变量不定义如：秋解下面线性规划问题图2输入为：图3所示图 33、按图2所示输入完成确

2、定后出现图4所示界面。图 4目标函数系数约束条件系数变量类型：双击改变约束形式：双击改变右端项图 5运行键4、输入完成后，按图5所示运行键。5、运行结果如图6所示图 6图6中各列的含义为：Decision Variable：决策变量Solution Value：决策方案取值Solution Value：决策变量对目标的单位贡献/目标函数系数Total Contribution：总贡献=（Solution Value）（Solution Value）Reduced Cost：检验数Allowable Min c(j) / Allowable Man c(j)：目标系数的灵敏度范围Objectiv

3、e Function：目标函数Constraint：约束条件（C1，C2，C3分别表示约束条件1、2、3）Left Hand Side：左端项，将决策变量取值代入约束方程左端计算的结果Right Hand Side：右端项，表示目前资源的拥有量Slack or Surplus：左端项与右端项的差额：资源的不足/slack或剩余/surplusShadow Price：资源的影子价格Allowable Min. RHS/ Allowable Max. RHS：右端项的灵敏度范围P65 2.9 已知线性规划问题：已知用单纯形法求得最优解的单纯形表如下，试分析在下列各种条件单独变化时，最优解如何变化

4、，看看与你的分析是否一致？x1x2x3x4x5x62x24/3012/3-1/3003x110/310-1/32/3000x5300-11100x62/300-2/31/301CjZj00-1/3-3/400(a) 第1个和第2个约束条件的右端项分别由6变成7，由8变成4；(b) 目标函数变为maxZ=2x1+5x2(c) 增加一个变量x3，其在目标函数中系数C34，在约束系数矩阵中列P3（1，2，3，2）T;3)整数规划P100 习题4.6P101 习题4.8（1），分别直接求解和用分枝定界法求解，比较结果。4)目标规划书P108 例3 (Goal Programming)四、图与网络分析上

5、机程序1、运行“Network Modeling”，出现图7所示界面图132、运行file菜单下的new problem 命令，出现图14所示界面。 图14图14中各项目含义：Problem Type(问题类型)如下：Transportation Problem:运输问题 Assignment Problem： 分配问题Shortest Path Problem：最短路问题 Maximal Flow Problem：最大流问题Minimal Spanning Tree：最小分布树 Traveling Salesman Problem：旅行商问题Objective Criterion选择目标为求

6、最大或最小Data Entry Format选择数据输入是以矩阵或图形输入Number of Nodes共有几个节点数3、如求一个最短路 P137 习题6.8 （b） 则选择Shortest Path Problem，Objective Criterion选最小化，Number of Nodes后写上节点数：7，Problem Title后给文件命名，OK确定则得矩阵形式输入如图15所示：图15根据P137 习题6.8 （b）图，按照从左到右，从上到下顺序给节点标上号，共有17个节点，填上节点间的直接距离，两点间如果没有直接连线，则空着不填，自己到自己的节点也可空着不填，如果是无向图（即i节点

7、到j节点距离与j到i相等）可只填上三角矩阵，填好后则得图16图16运行，出现图17如下，选择求最短路的起始节点和最后节点，点solve运行，也可点solve and Display Steps看求解过程，最后得结果，如图18。图17图18图18中表示1节点到7节点的最短路径为从Node1Node2Node5Node7，总距离13。4、求最大流、最小树输入与此类似，最大流中权数为最大流量。5、上机任务最小树 P137 习题6.5 （d） Minimal Spanning Tree最短路 P137 习题6.8 （a） Shortest Path Problem最大流 P139 习题6.13 （c）

8、 Maximal Flow Problem分配问题 P100 习题4.6 （a）（b） Assignment Problem五、计划评审方法和关键路径法上机程序1、运行“PERT_CPM”，出现图19所示界面图192、运行file菜单下的new problem 命令，出现图20所示界面。 图20图20中各项目含义：Problem Type(问题类型)如下：Deterministic CPM： 确定型关键路线法 Probabilistic PERT： 概率型网络计划技术Data Entry Format选择数据输入是以矩阵或图形输入Select CPM Data FieldNormal Time

9、 正常时间 Crash Time 赶工时间 Normal Cost 正常费用 Crash Cost 赶工费用3、如求一个PERT问题， P143 例1 则Problem Title后给文件命名，Number of Activities后给出作业数9，Time Unit后给出时间单位hour，Problem Type选择Deterministic CPM， Select CPM Data Field选Normal Time， 输入界面如图21所示，OK确定后出现输入矩阵如图22所示，图21图22图22中表格各项含义如下：Activity Number：作业编号，按1、2、3等依次对各项作业编号A

10、ctivity Name：作业名称，可自行取名填入Immediate Predecessor：紧前工序，填入该项作业的紧前作业，可以填紧前作业的编号或名称，若有多项紧前作业，每项之间用西文状态下的逗号,格开Normal Time：作业时间按P143例1填入后如图23所示图23按运行图标运行，得结果如图24所示图24图24中从左到右各列含义依次如下：1、作业编号 2、作业名称 3、该作业是否是关键路径上的关键作业，若是则为Yes，若不是则no4、作业时间 5、作业最早可能开始时间 6、作业最早可能完成时间7、作业最迟必须开始时间 8、作业最迟必须完成时间9、作业总时差图24中最后两行给出了项目完

11、成时间，本题为17.50，和关键路径的数量，本题为1。4、求一个概率型PERT，书P150，例3（1）按顺序给书P140表7-4中的9项作业编号为19，（2）选择File菜单下New Problem，（3）在出现的对话框中：Number of Activities填入9， Problem Type选择Problistic PERT, 点击OK, (4) 在出现的表格中填入书P150表7-4的数据后如图25所示，点击运行键运行，（5）结果意义与上例同，其中Activity Mean Time和Standard Deviation给出了完成各项作业的期望时间和标准差。图255、上机任务P152 习

12、题7.4 P153 习题7.6 六、动态规划1、运行“Dynamic Programming”，进入动态规划程序，菜单栏上选择“File”“New”，进入图26所示对话框；对话框中列出了本程序可以求解的三种动态规划问题，分别为“最短路径问题”，“背包问题”，“生产与存储计划问题”，下面将分别举例说明 图26 2、最短路径问题，以P155例1为例。（1）在图26的对话框中，点选第一项，并在“Number of Nodes”中填入节点个数，例1中共有节点数“10”， 完成后点击“OK”。（2）在出现的界面中，按照各个节点间连接关系，完成距离矩阵，填完后的矩阵如图27所示，按运行键运行，出现一个对话

13、框如图28所示，在左边的列表中选择路径的开始节点，在右边的列表中选择路径的结束节点，点“solve”求解问题，也可点solve and Display Steps看求解过程，最后得结果，如图29。（3）图29中给出了用“顺序法”分四个阶段求解的每个阶段的决策，及当阶段指标值“Distance”，及该阶段的最优指标值”Cumulative Distance”，并在求得最后一个阶段即第四阶段的最优指标值时即得到整个过程上的最短路径“11” 图27图28图293、背包问题，以P169例8为例。（1）运行file菜单下的new problem 命令，在出现的对话框中，“Problem Type”点选“

14、Knapsack Problen”，“Number of Items”根据例8填入“3”，点击“ok”。（2）在出现的界面中，按照例8数据填入后如图30所示，其中从左到右各项含义为：Item背包中装载货物数目，也即为动态规划的阶段数；Item Identification货物名称，可以改变其命名；Units Avaliable可提供的该货物的最大数量，因本例中未提及货物的数量限制，因此均为M；Unit Capacity Required单位货物的体积或重量，本题为重量；Return Function单位货物的价值，用函数式表示，如单位Item1货物价值为50，则写为“50Item1”，以此类推

15、；在底行的“Knapsack Capacity”背包最大容量；填完后运行得结果如图31。（3）图31结果中各项含义可参照上例最短路径的结果来理解。 图30 图314、生产与存储计划问题，以P171习题8.3为例。（1）运行file菜单下的new problem 命令，在出现的对话框中，“Problem Type”点选“Production and Inventory Scheduling Problem”，“Number of Periods”根据例8填入“4”，点击“ok”。（2）在出现的界面中，按照例8数据填入后如图32所示，其中从左到右各项含义为：Period生产与存储划分的时期，即动态

16、规划的阶段数，本题为4阶段生产与存储问题；Period Identification各时期的名称，可以改变其命名；Demand各时期需求量，本题4个时期需求量分别为2、3、2、4；Production Capacity各时期生产能力，本题由C函数（x大于6时，C为无穷）分析为6；Storage Capacity各时期存储能力，本题无限制，取默认值M；Production Setup Cost生产准备费用，本题由C函数（0x6时，C为3x）分析为3；Variable Cost Function费用函数，用变量P、H、B三个变量的函数式表示，P、H、B分别表示本阶段末的产量、存储量、缺货量，本题中

17、单位生产成本3，单位存储费用0.5，无缺货损失，因此费用函数为P+0.5H,填完后运行得结果如图33。（3）图33结果中给出了最优生产计划，每个时期初始存储量（“Starting Inventory”） 相当于动态规划的状态变量分别为0、3、0、4，每个时期的产量(“Production Quantity”)分别为5、0、6、0，该策略对应的总变量成本（不包括生产准备费用）为14.50，总成本（加上生产准备费用）为20.50。图32图335、上机任务上述三个例子的练习七、存储论1、运行“Inventory Theory and System”，菜单栏上选择“File”“New”，进入图34所示

18、对话框；对话框中列出了本程序可以求解的问题类型，根据所学内容，将应用两种类型，即“EOQ(经济订货批量)模型”，“单时期随机存储模型”。下面将分别举例说明。图342、EOQ(经济订货批量)模型，以P190习题1为例。在图34中点选第一个，即EOQ模型，并填入Problem Title, Time Unit后，点击OK, 进入图35界面，并依据习题1填入参数，需求量=10 000, 订货费用=2000，单位存储费用1000.2=20, 具体如图35所示，填入后运行，得到图36所示结果，结果显示最佳订货批量为1414，总的订货费用为14142.13，总存储费用14142.14，总的费用为28284

19、.27。图35图363、单时期随机存储模型，以P187例8为例。在图34中点选第3个，即模型Single-period Stochastic Demand(Newsboy) Problem，点击OK，依据例8填入参数：需求分布Normal（正态分布，可以双击Normal，在弹出的对话框中改变随机分布类型）, 需求期望300，需求方差50，单位成本Unit acquistion cost0.3, 单位售价Unit Selling Price0.5，单位处理价Unit salvage value0.2，具体如图37所示，填入后运行。图37得到图38所示结果，结果显示最佳订货批量为321.5373，

20、期望的利润为54.55。图384、上机任务 P191 习题9.2，9.3，9.12。 八、排队论1、运行“Qeueing Analysis”，菜单栏上选择“File”“New”，进入图39所示对话框；根据问题类型选择“Simple M/M System” 或“General Queuing System”。图392、M/M/S排队模型，以P209例5为例，点选“Simple M/M System” “OK”。在图40中，依据例5填入参数，具体如图40所示，填入后运行，得到图41所示结果。图40图41图41结果显示：系统有效到达率eff19.3135人/h，系统忙率为92.5506，平均队长Ls9.5563，平均排队长Lq5.8543，平均等待时间Wq=0.3031h等。3、上机任务 P224 习题10.7，10.9。