数学建模敏感性分析幻灯片.ppt

《数学建模敏感性分析幻灯片.ppt》由会员分享，可在线阅读，更多相关《数学建模敏感性分析幻灯片.ppt（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数学建模敏感性分析数学建模敏感性分析课件课件第1页，共20页，编辑于2022年，星期六内容提要内容提要:例例2.4 例例2.5 2.2 敏感性分析敏感性分析第2页，共20页，编辑于2022年，星期六例例2.4 某某家家具具公公司司制制造造书书桌桌、餐餐桌桌和和椅椅子子，所所用用的的资资源源有有三三种种：木料、木工和漆工。生产数据如下表所示。木料、木工和漆工。生产数据如下表所示。每个书桌每个餐桌每个椅子现有资源总数木料8单位6单位1单位48单位漆工4单位2单位1.5单位20单位木工2单位1.5单位0.5单位8单位成品单价60单位30单位20单位 若若要要求求桌桌子子的的生生产产量量不不超超过过5

2、件件，如如何何安安排排三三种种产产品品的的生生产可使利润最大？产可使利润最大？第3页，共20页，编辑于2022年，星期六解：解：用用DESKS、TABLES和和CHAIRS分别表示三种产品的分别表示三种产品的生产量（决策变量），容易建立生产量（决策变量），容易建立LP模型。模型。在在LINDO模型窗口中输入模型：模型窗口中输入模型：MAX 60 DESKS+30 TABLES+20 CHAIRSSUBJECT TO 2)8 DESKS+6 TABLES+CHAIRS=48 3)4 DESKS+2 TABLES+1.5 CHAIRS=20 4)DESKS+1 5 TABLES+O 5 CHAIR

3、S=8 5)TABLES=5END解这个模型，并对弹出的对话框解这个模型，并对弹出的对话框 “DO RANGE(SENSITIVITY)ANALYSIS?”选择选择“是（是（Y）”按钮，这表示你需要做灵敏性分析。然按钮，这表示你需要做灵敏性分析。然后，查看输出结果。后，查看输出结果。第4页，共20页，编辑于2022年，星期六LP OPTIMUM FOUND AT STEP 1 OBJECTIVE FUNCTION VALUE 1)280.0000 VARIABLE VALUE REDUCED COST DESKS 2.000000 0.000000 TABLES 0.000000 5.0000

4、00 CHAIRS 8.000000 0.000000 ROW SLACK OR SURPLUS DUAL PRICES 2)24.000000 0.000000 3)0.000000 10.000000 4)0.000000 10.000000 5)5.000000 0.000000 NO.ITERATIONS=1输出结果的前半部分：输出结果的前半部分：第5页，共20页，编辑于2022年，星期六前半部分的输出结果的解释与前一节例前半部分的输出结果的解释与前一节例2.1的结果类似：的结果类似：“LP OPTIMUM FOUND AT STEP2”表示两次迭代（旋转变换）后得到最优解。表示两次迭

5、代（旋转变换）后得到最优解。“OBJECTIVE FUNCTION VALUE 1）280.000000”表示最优目标值为表示最优目标值为280。“VALUE”给出最优解中各变量的值：造给出最优解中各变量的值：造2个书桌（个书桌（desks）,0个餐桌（个餐桌（tables）,8个椅子个椅子（chairs）。所以）。所以desks、chairs是基变量（取值非是基变量（取值非0），），tables是非基变量（取是非基变量（取值为值为0）。）。“SLACK OR SURPLUS”给出松驰变量的值：给出松驰变量的值：第第2行松驰变量行松驰变量=24 （第（第1行表示目标函数，第行表示目标函数，第2

6、行对应第行对应第1个约束）个约束）第第3行松驰变量行松驰变量=0 第第4行松驰变量行松驰变量=0 第第5行松驰变量行松驰变量=5第6页，共20页，编辑于2022年，星期六“REDUCED COST”列出最优单纯形表中判别数所在行的变量的系数，表示当变量有微小列出最优单纯形表中判别数所在行的变量的系数，表示当变量有微小变动时变动时,目标函数的变化率目标函数的变化率.其中基变量的其中基变量的reduced cost值应为值应为0，对于非基变量对于非基变量 Xj（请注意，非基变量的取值一定是（请注意，非基变量的取值一定是0），），相应的相应的 reduced cost值表示当某个变量值表示当某个变量

7、Xj 增加一个单位时目标函数减少的增加一个单位时目标函数减少的量量(max型问题型问题)。本例中：变量本例中：变量TABLES对应的对应的reduced cost值为值为5，表示当非基变量，表示当非基变量TABLES 的值从的值从0变为变为 1时（此时假定其他非基变量保持不变时（此时假定其他非基变量保持不变,但为了但为了满足约束条件，基变量显然会发生变化），最优的目标函数值满足约束条件，基变量显然会发生变化），最优的目标函数值=280-5=275。第7页，共20页，编辑于2022年，星期六“DUAL PRICE”（对偶价格）表示当对应约束有微小变动时（对偶价格）表示当对应约束有微小变动时,目标

8、函数目标函数的变化率的变化率.输出结果中对应于每一个约束有一个对偶价格输出结果中对应于每一个约束有一个对偶价格.若其数值若其数值为为p,表示对应约束中不等式右端项若增加表示对应约束中不等式右端项若增加1 个单位个单位,目标函数将增目标函数将增加加p个单位（个单位（max型问题型问题)。显然，如果在最优解处约束正好取等号（也就是显然，如果在最优解处约束正好取等号（也就是“紧约束紧约束”，即起作用，即起作用约束），对偶价格值才可能不是约束），对偶价格值才可能不是0。本例中：第。本例中：第3、4行是紧约束，对应的行是紧约束，对应的对偶价格值为对偶价格值为10，表示当紧约束，表示当紧约束 3)4 DE

9、SKS+2 TABLES+1.5 CHAIRS=20变为变为 3)4 DESKS+2 TABLES+1.5 CHAIRS=21时，目标函数值时，目标函数值=280+10=290。对第。对第4行也可类似解释。行也可类似解释。对于非紧约束（如本例中第对于非紧约束（如本例中第2、5行是非紧约束），行是非紧约束），DUAL PRICE 的值为的值为0,表示对应约束中不等式右端项的微小扰动不影响目标函表示对应约束中不等式右端项的微小扰动不影响目标函数。有时数。有时,通过分析通过分析DUAL PRICE,也可对产生不可行问题的原因也可对产生不可行问题的原因有所了解。有所了解。第8页，共20页，编辑于202

10、2年，星期六输出结果的后半部分：输出结果的后半部分：RANGES IN WHICH THE BASIS IS UNCHANGED:OBJ COEFFICIENT RANGES VARIABLE CURRENT ALLOWABLE ALLOWABLE COEF INCREASE DECREASE DESKS 60.000000 20.000000 4.000000 TABLES 30.000000 5.000000 INFINITY CHAIRS 20.000000 2.500000 5.000000 RIGHTHAND SIDE RANGES ROW CURRENT ALLOWABLE ALL

11、OWABLE RHS INCREASE DECREASE 2 48.000000 INFINITY 24.000000 3 20.000000 4.000000 4.000000 4 8.000000 2.000000 1.333333 5 5.000000 INFINITY 5.000000（报告中（报告中INFINITY表示正无穷表示正无穷）第9页，共20页，编辑于2022年，星期六 敏感性分析的作用是给出敏感性分析的作用是给出“RANGES IN WHICH THE BASIS IS UNCHANGED”，即研究当目标函数的系数和约束右端项在什么，即研究当目标函数的系数和约束右端项在什么

12、范围变化（此时假定其他系数保持不变）时，最优基（矩阵）范围变化（此时假定其他系数保持不变）时，最优基（矩阵）保持不变。保持不变。1.目标函数中系数变化的范围目标函数中系数变化的范围 （OBJ COEFFICIENT RANGES）如本例中：目标函数中如本例中：目标函数中DESKS变量当前的系数（变量当前的系数（CURRENT COEF）=60，允许增加（，允许增加（Allowable Increase）=4、允许减少、允许减少（Allowable Decrease）=2，说明当这个系数在，说明当这个系数在60-4，60+20=56，80范围变化时，最优基保持不变。对范围变化时，最优基保持不变。

13、对TABLES、CHAIRS变量，可以类似解释。由于此时约束没有变化（只是目标变量，可以类似解释。由于此时约束没有变化（只是目标函数中某个系数发生变化），所以最优基保持不变的意思也就是函数中某个系数发生变化），所以最优基保持不变的意思也就是最优解不变（当然，由于目标函数中系数发生了变化，所以最优最优解不变（当然，由于目标函数中系数发生了变化，所以最优值会变化）。值会变化）。这个部分包括两方面的敏感性分析内容：这个部分包括两方面的敏感性分析内容：第10页，共20页，编辑于2022年，星期六2.约束右端项变化的范围（约束右端项变化的范围（Right Hand Side RANGES）如本例中：第如

14、本例中：第2行约束中当前右端项（行约束中当前右端项（CURRENT RHS）=48，允许增加（允许增加（Allowable Increase）=INFINITY（无穷）、允许（无穷）、允许减少（减少（Allowable Decrease）=24，说明当它在，说明当它在 48-24，48+)=24，）范围变化时，最优基保持不变。第范围变化时，最优基保持不变。第3、4、5行可以类似解释。行可以类似解释。不过由于此时约束发生变化，最优基即使不变，最优解、最优不过由于此时约束发生变化，最优基即使不变，最优解、最优值也会发生变化。如何变化呢？我们将在本节后面结合第值也会发生变化。如何变化呢？我们将在本节

15、后面结合第1章章1.2.1节例节例1.1给出的实际问题来进行说明。给出的实际问题来进行说明。第11页，共20页，编辑于2022年，星期六最后，如果你对单纯形法比较熟悉，你可以直接查看最优解时的单纯形最后，如果你对单纯形法比较熟悉，你可以直接查看最优解时的单纯形表，这只要选择菜单命令表，这只要选择菜单命令Reports|Tableau(Alt+7)执行即可执行即可,输出输出结果如下：结果如下：ROW (BASIS)DESKS TABLES CHAIRS 1 ART 0.000 5.000 0.000 2 SLK 2 0.000 -2.000 0.000 3 DESKS 1.000 1.250 0

16、.000 4 CHAIRS 0.000 -2.000 1.000 5 SLK 5 0.000 1.000 0.000 ROW SLK 2 SLK 3 SLK 4 SLK 5 1 0.000 10.000 10.000 0.000 280.000 2 1.000 2.000 -8.000 0.000 24.000 3 0.000 1.500 -0.500 0.000 2.000 4 0.000 -4.000 2.000 0.000 8.000 5 0.000 0.000 0.000 1.000 5.000第12页，共20页，编辑于2022年，星期六在输出结果中，在输出结果中，基变量为基变量为BV

17、=SLK2,Chairs,DESKS,SLK5，ART是人工变量（是人工变量（artificial variable）,即相应的目标值即相应的目标值z;这样，你就可以知道这样，你就可以知道 z=5 TABLES+10 SLK3+10 SLK4=280。敏感性分析结果表示的是最优基保持不变的系数范围。由此，敏感性分析结果表示的是最优基保持不变的系数范围。由此，也可以进一步确定当目标函数的系数和约束右端项发生小的变化也可以进一步确定当目标函数的系数和约束右端项发生小的变化时，最优解、最优值如何变化。下面我们通过求解第时，最优解、最优值如何变化。下面我们通过求解第1章章1.2.1节节例例1.1的实际

18、问题来进行说明。的实际问题来进行说明。第13页，共20页，编辑于2022年，星期六例例2.5 继续讨论例继续讨论例1.11.1 在在LINDO模型窗口中输入模型模型窗口中输入模型 问题的数学模型：问题的数学模型：MAX 72 x1+64 x2SUBJECT TO 2)x1+x2=50 3)12 x1+8 x2=480 4)3 x1 =100END第14页，共20页，编辑于2022年，星期六求解这个模型并做灵敏性分析，查看报告窗口（求解这个模型并做灵敏性分析，查看报告窗口（Reports Window）。输出结果告诉我们：这个线性规划的最优解为）。输出结果告诉我们：这个线性规划的最优解为x1=2

19、0，x2=30，最优值为，最优值为z=3360，即用，即用20桶牛奶生产桶牛奶生产A1,30桶牛奶生产桶牛奶生产A2，可获最大利润，可获最大利润3360元。元。输出中除了告诉我们问题的最优解和最优值以外，还有许多对输出中除了告诉我们问题的最优解和最优值以外，还有许多对分析结果有用的信息，下面结合题目中提出的分析结果有用的信息，下面结合题目中提出的3个附加问题给予个附加问题给予说明。说明。3个约束条件的右端不妨看作个约束条件的右端不妨看作3种种“资源资源”：原料、劳动时间、车间甲：原料、劳动时间、车间甲的加工能力。输出中的加工能力。输出中SLACK OR SURPLUS（松弛或剩余）给出（松弛或

20、剩余）给出这这3种资源在最优解下是否有剩余：原料、劳动时间的剩余均为种资源在最优解下是否有剩余：原料、劳动时间的剩余均为零（即约束为紧约束），车间甲尚余零（即约束为紧约束），车间甲尚余40公斤加工能力（不是紧约公斤加工能力（不是紧约束）。束）。第15页，共20页，编辑于2022年，星期六目标函数可以看作目标函数可以看作“效益效益”，成为紧约束的，成为紧约束的“资源资源”一旦增加，一旦增加，“效益效益”必然跟着增长。必然跟着增长。输出中输出中DUAL PRICES（对偶价格）（对偶价格）给出这给出这3种资源在最优解下种资源在最优解下“资资源源”增加增加1个单位时个单位时“效益效益”的增量：原料增

21、加的增量：原料增加1个单位（个单位（1桶牛奶）时桶牛奶）时利润增长利润增长48（元），劳动时间增加（元），劳动时间增加1个单位（个单位（1小时）时利润增长小时）时利润增长2（元），而增加非紧约束车间甲的能力显然不会使利润增长。（元），而增加非紧约束车间甲的能力显然不会使利润增长。这里，这里，“效益效益”的增量可以看作的增量可以看作“资源资源”的潜在价值，经济学上称为影子价格的潜在价值，经济学上称为影子价格（shadow price），即），即1桶牛奶的影子价格为桶牛奶的影子价格为48元，元，1小时劳动的影子价格为小时劳动的影子价格为2元，元，车间甲生产能力的影子价格为零。车间甲生产能力的影子价

22、格为零。可以用直接求解的办法验证上面的结论，即将输入文件中原料约束可以用直接求解的办法验证上面的结论，即将输入文件中原料约束milk）右端的）右端的50改为改为51，看看得到的最优值（利润）是否恰好增长，看看得到的最优值（利润）是否恰好增长48（元）。用影子价格的概念（元）。用影子价格的概念很容易回答附加问题很容易回答附加问题1）：用）：用35元可以买到元可以买到1桶牛奶，低于桶牛奶，低于1桶牛奶的影子价格桶牛奶的影子价格48，当，当然应该作这项投资。回答附加问题然应该作这项投资。回答附加问题2）：聘用临时工人以增加劳动时间，付给的工资）：聘用临时工人以增加劳动时间，付给的工资低于劳动时间的影

23、子价格才可以增加利润，所以工资最多是每小时低于劳动时间的影子价格才可以增加利润，所以工资最多是每小时2元。元。第16页，共20页，编辑于2022年，星期六目标函数的系数发生变化时（假定约束条件不变），最优解和最目标函数的系数发生变化时（假定约束条件不变），最优解和最优值会改变吗？优值会改变吗？这个问题不能简单地回答。这个问题不能简单地回答。上面的输出结果给出了最优基不变条件下目标函数系数的允许变化上面的输出结果给出了最优基不变条件下目标函数系数的允许变化范围：范围：x1的系数范围为（的系数范围为（72-8，72+24）=（64，96）；）；x2的系数的系数范围为（范围为（64-16，64+8）

24、=（48，72）。）。注意：注意：x1系数的允许范围需要系数的允许范围需要x2的系数的系数64不变，反之亦然。不变，反之亦然。由于目标函数的系数变化并不影响约束条件，因此此时最优基由于目标函数的系数变化并不影响约束条件，因此此时最优基不变可以保证最优解也不变，但最优值变化。用这个结果很容不变可以保证最优解也不变，但最优值变化。用这个结果很容易回答附加问题易回答附加问题3）：若每公斤）：若每公斤A1的获利增加到的获利增加到30元，则元，则x1系数系数变为变为303=90，在允许范围内，所以不应改变生产计划，但最，在允许范围内，所以不应改变生产计划，但最优值变为优值变为9020+6430=3720

25、。第17页，共20页，编辑于2022年，星期六下面对下面对“资源资源”的影子价格作进一步的分析。的影子价格作进一步的分析。影子价格的作用影子价格的作用（即在最优解下（即在最优解下“资源资源”增加增加1个单位时个单位时“效益效益”的增的增量）量）是有限制的是有限制的。每增加每增加1桶牛奶利润增长桶牛奶利润增长48元（影子价格），但是，约束的右端项元（影子价格），但是，约束的右端项（CURRENT RHS）的的“允许增加允许增加”（ALLOWABLE INCREASE）和和“允许减少允许减少”（ALLOWABLE DECREASE）给出了影子价格给出了影子价格有意义条件下约束右端的限制范围（因为此

26、时最优基不变，所以有意义条件下约束右端的限制范围（因为此时最优基不变，所以影子价格才有意义；如果最优基已经变了，那么结果中给出的影影子价格才有意义；如果最优基已经变了，那么结果中给出的影子价格也就不正确了）。子价格也就不正确了）。具体对本例来说：具体对本例来说：milk）原料最多增加）原料最多增加10（桶牛奶），（桶牛奶），time）劳动）劳动时间最多增加时间最多增加53（小时）。（小时）。现在可以回答附加问题现在可以回答附加问题1）的第）的第2问：虽然应该批准用问：虽然应该批准用35元买元买1桶牛奶的投桶牛奶的投资，但每天最多购买资，但每天最多购买10桶牛奶。顺便地说，可以用低于每小时桶牛奶

27、。顺便地说，可以用低于每小时2元的工资元的工资聘用临时工人以增加劳动时间，但最多增加聘用临时工人以增加劳动时间，但最多增加53.3333小时。小时。第18页，共20页，编辑于2022年，星期六需要注意的是：需要注意的是：灵敏性分析给出的只是最优基保持不变的充灵敏性分析给出的只是最优基保持不变的充分条件，而不一定是必要条件分条件，而不一定是必要条件。比如对于上面的问题，。比如对于上面的问题，“原料最原料最多增加多增加10（桶牛奶）（桶牛奶）”的含义只能是的含义只能是“原料增加原料增加10（桶牛奶）（桶牛奶）”时最优基保持不变，所以影子价格有意义，即利润的增加时最优基保持不变，所以影子价格有意义，

28、即利润的增加大于牛奶的投资。反过来，原料增加超过大于牛奶的投资。反过来，原料增加超过10（桶牛奶），最优（桶牛奶），最优基是否一定改变？影子价格是否一定没有意义？一般来说，这是基是否一定改变？影子价格是否一定没有意义？一般来说，这是不能从灵敏性分析报告中直接得到的。此时，应该重新用新数据不能从灵敏性分析报告中直接得到的。此时，应该重新用新数据求解规划模型，才能做出判断。所以严格来说，我们上面回答求解规划模型，才能做出判断。所以严格来说，我们上面回答“原料最多增加原料最多增加10（桶牛奶）（桶牛奶）”并不是完全科学的。并不是完全科学的。第19页，共20页，编辑于2022年，星期六THE ENDTHE ENDThank you very much!Thank you very much!第20页，共20页，编辑于2022年，星期六