运筹学案例分析报告示例

食油生产问题（案例一）分析报告一、模型构造1.1 变量设置设两种硬质油代号分别为HD1、HD2（HD代表Hard）,三种软质油代号分别为SF1、SF2、SF3（SF代表Soft）。每种油的采购（Buy）、耗用（Use）和储存（Store）量分别在油品的代号前加B、U和S表示。16月份5种油品的采购、耗用和储存量分别在油品代号后面加16表示。总产量用PROD（Product）表示。第一种硬质油六个月的采购量、耗用量、月末储存量共有17变量，其中，六月末的存储量为500吨。BHD11，BHD12，BHD13，BHD14，BHD15，BHD16；UHD11，UHD12，UHD13，UHD14，UHD15，UHD16；SHD11，SHD12，SHD13，SHD14，SHD15；第二种硬质油六个月的采购量、耗用量、月末储存量共有17变量，其中，六月末的存储量为500吨。BHD21，BHD22，BHD23，BHD24，BHD25，BHD26；UHD21，UHD22，UHD23，UHD24，UHD25，UHD26；SHD21，SHD22，SHD23，SHD24，SHD25；第一种软质油六个月的采购量、耗用量、月末储存量共有17变量，其中，六月末的存储量为500吨。BSF11，BSF12，BSF13，BSF14，BSF15，BSF16；USF11，USF12，USF13，USF14，USF15，USF16；SSF11，SSF12，SSF13，SSF14，SSF15；第二种软质油六个月的采购量、耗用量、月末储存量共有17变量，其中，六月末的存储量为500吨。BSF21，BSF22，BSF23，BSF24，BSF25，BSF26；USF21，USF22，USF23，USF24，USF25，USF26；SSF21，SSF22，SSF23，SSF24，SSF25；第三种软质油六个月的采购量、耗用量、月末储存量共有17变量，其中，六月末的存储量为500吨。BSF31，BSF32，BSF33，BSF34，BSF35，BSF36；USF31，USF32，USF33，USF34，USF35，USF36；SSF31，SSF32，SSF33，SSF34，SSF35；六个月的油品产量共有六个变量PROD1，PROD2，PROD3，PROD4，PROD5，PROD6；这样,一共有517+6=91个变量。1.2 模型构建这个问题的建模可以分为两步，首先建立单期问题模型，然后建立多期问题模型。1.2.1单期模型以一月份的模型为例，如果原料油没有库存，一月份的采购量和耗用量相等，生产计划线性规划模型如下，其中有5个约束条件：产品油硬质油1硬质油2软质油1软质油2软质油3利润max200 PROD1-110 BHD11-120 BHD21-130 BSF11-110 BSF21-115 BSF31硬质生产能力s.t.UHD11+UHD21200软质生产能力USF11+USF21+USF31250硬度上限约束-6.0 PROD1+8.8 UHD11+6.1 UHD21+2.0 USF11+4.2 USF21+5.0 USF310硬度下限约束-3.0 PROD1+8.8 UHD11+6.1 UHD21+2.0 USF11+4.2 USF21+5.0 USF310物料平衡约束-PROD1+UHD11+UHD21+USF11+USF21+USF31=0变量非负约束PROD1,UHD11,UHD21,USF11,USF21,USF310二月份的生产计划模型与以上模型的结构完全相同，不同的是模型中的变量全部换成5种油二月份的耗用量，目标函数中5种油品的价格换成二月份的价格。二月份的生产计划模型如下，其中有5个约束条件：产品油硬质油1硬质油2软质油1软质油2软质油3利润max200 PROD2-130 BHD12-130 BHD22-110 BSF12-90 BSF22-115 BSF32硬质生产能力s.t.UHD12+UHD22200软质生产能力USF12+USF22+USF32250硬度上限约束-6.0 PROD2+8.8 UHD12+6.1 UHD22+2.0 USF12+4.2 USF22+5.0 USF320硬度下限约束-3.0 PROD2+8.8 UHD12+6.1 UHD22+2.0 USF12+4.2 USF22+5.0 USF320物料平衡约束-PROD2+UHD12+UHD22+USF12+USF22+USF32=0变量非负约束PROD2,UHD12,UHD22,USF12,USF22,USF320类似地，可以分别建立三月份、四月份、五月份和六月份的生产计划线性规划模型。六月份的生产计划模型如下：产品油硬质油1硬质油2软质油1软质油2软质油3利润max200 PROD6-90 BHD16-110 BHD26-140 BSF16-80 BSF26-135 BSF36硬质生产能力s.t.UHD16+UHD26200软质生产能力USF16+USF26+USF36250硬度上限约束-6.0 PROD6+8.8 UHD16+6.1 UHD26+2.0 USF16+4.2 USF26+5.0 USF360硬度下限约束-3.0 PROD6+8.8 UHD16+6.1 UHD26+2.0 USF16+4.2 USF26+5.0 USF360物料平衡约束-PROD6+UHD16+UHD26+USF16+USF26+USF36=0变量非负约束PROD6,UHD16,UHD26,USF16,USF26,USF360以上6个单期模型，共有30约束条件。1.2.2 多期模型以上6个单期模型中，目标函数中的变量是当月各种原料油的采购量和成品油的产量，而约束条件中的变量是当月原料油的耗用量和成品油的产量。但根据题意，每一种原料油上月末的存储量，本月的采购量，本月耗用量，本月末存储量有以下关系：（上月末存储量）+（本月采购量）-（本月耗用量）-（本月末存储量）=0对于一月份，有：（一月份采购量）-（一月份耗用量）-（一月末存储量）=500对于六月份，有：（五月末存储量）+（六月份采购量）-（六月份耗用量）=500用以上的变量代号表示，第一种原料油HD1一月份到六月份的采购量、耗用量和存储量有以下约束： BHD11-UHD11-SHD11=-500SHD11+BHD12-UHD12-SHD12= 0SHD12+BHD13-UHD13-SHD13= 0SHD13+BHD14-UHD14-SHD14= 0SHD14+BHD15-UHD15-SHD15= 0SHD15+BHD16-UHD16=500第二种原料油HD2一月份到六月份的采购量、耗用量和存储量有以下约束： BHD21-UHD21-SHD21=-500SHD21+BHD22-UHD22-SHD22= 0SHD22+BHD23-UHD23-SHD23= 0SHD23+BHD24-UHD24-SHD24= 0SHD24+BHD25-UHD25-SHD25= 0SHD25+BHD26-UHD26=500第三种原料油SF1一月份到六月份的采购量、耗用量和存储量有以下约束： BSF11-USF11-SSF11=-500SSF11+BSF12-USF12-SSF12= 0SSF12+BSF13-USF13-SSF13= 0SSF13+BSF14-USF14-SSF14= 0SSF14+BSF15-USF15-SSF15= 0SSF15+BSF16-USF16=500第四种原料油SF2一月份到六月份的采购量、耗用量和存储量有以下约束： BSF21-USF21-SSF21=-500SSF21+BSF22-USF22-SSF22= 0SSF22+BSF23-USF23-SSF23= 0SSF23+BSF24-USF24-SSF24= 0SSF24+BSF25-USF25-SSF25= 0SSF25+BSF26-USF26=500第五种原料油SF3一月份到六月份的采购量、耗用量和存储量有以下约束： BSF31-USF31-SSF31=-500SSF31+BSF32-USF32-SSF32= 0SSF32+BSF33-USF33-SSF33= 0SSF33+BSF34-USF34-SSF34= 0SSF34+BSF35-USF35-SSF35= 0SSF35+BSF36-USF36=500以上5组反映5种不同的原料油一月到六月的采购量、耗用量、存储量之间平衡关系的约束称为耦合约束。以上5组耦合约束条件共有56=30个，连同6个月单期模型的30个约束条件，总体模型共有60个约束条件。由于油罐容量的限制，每个月5种原料的存储量不能超过1000吨，25个存储变量的上界为1000。变量的上界限制在实际问题中一般不用约束条件表示，在LINDO中有为变量设定上界的语句。形式为：SUB变量名上界值如SUBSHD111000由于考虑了每个月的存储量，在目标函数中还必须引进存储变量，系数为5元/吨。这样，91个变量，60个约束条件的总体线性规划模型的结构如下：二、LINDO模型和求解结果2.1、LINDO模型TITLE 案例1 食油生产问题（1）!-!目标函数=六个月成品油（PROD1-PROD6）的销售收入减去六个月五种原料油的成本!减去每个月五种原料油的库存费用max 200 PROD1-110 BHD11-120 BHD21-130 BSF11-110 BSF21-115 BSF31+200 PROD2-130 BHD12-130 BHD22-110 BSF12- 90 BSF22-115 BSF32+200 PROD3-110 BHD13-140 BHD23-130 BSF13-100 BSF23- 95 BSF33+200 PROD4-120 BHD14-110 BHD24-120 BSF14-120 BSF24-125 BSF34+200 PROD5-100 BHD15-120 BHD25-150 BSF15-110 BSF25-105 BSF35+200 PROD6- 90 BHD16-110 BHD26-140 BSF16- 80 BSF26-135 BSF36-5 SHD11-5 SHD21-5 SSF11-5 SSF21-5 SSF31-5 SHD12-5 SHD22-5 SSF12-5 SSF22-5 SSF32-5 SHD13-5 SHD23-5 SSF13-5 SSF23-5 SSF33-5 SHD14-5 SHD24-5 SSF14-5 SSF24-5 SSF34-5 SHD15-5 SHD25-5 SSF15-5 SSF25-5 SSF35st!-!一月份硬质油生产线的生产能力约束（JAN-HPC）!一月份软质油生产线的生产能力约束（JAN-SPC）!一月份成品油的硬度上限约束（JAN-HUB）!一月份成品油的硬度下限约束（JAN-HLB）!一月份原料和成品的平衡约束（JAN-BLS）JAN-HPC) UHD11+UHD21200JAN-SPC) USF11+USF21+USF31250JAN-HUB)-6.0PROD1+8.8UHD11+6.1UHD21+2.0USF11+4.2USF21+5.0 USF310JAN-BLS) -PROD1+UHD11+UHD21+USF11+USF21+USF31=0!-!二月份硬质油生产线的生产能力约束（FEB-HPC）!二月份软质油生产线的生产能力约束（FEB-SPC）!二月份成品油的硬度上限约束（FEB-HUB）!二月份成品油的硬度下限约束（FEB-HLB）!二月份原料和成品的平衡约束（FEB-BLS）FEB-HPC) UHD12+UHD22200FEB-SPC) USF12+USF22+USF32250FEB-HUB)-6.0PROD2+8.8UHD12+6.1UHD22+2.0USF12+4.2USF22+5.0 USF320FEB-BLS) -PROD2+UHD12+UHD22+USF12+USF22+USF31=0!-!三月份硬质油生产线的生产能力约束（MAR-HPC）!三月份软质油生产线的生产能力约束（MAR-SPC）!三月份成品油的硬度上限约束（MAR-HUB）!三月份成品油的硬度下限约束（MAR-HLB）!三月份原料和成品的平衡约束（MAR-BLS）MAR-HPC) UHD13+UHD23200MAR-SPC) USF13+USF23+USF33250MAR-HUB)-6.0PROD3+8.8UHD13+6.1UHD23+2.0USF13+4.2USF23+5.0USF330MAR-BLS) -PROD3+UHD13+UHD23+USF13+USF23+USF33=0!-!四月份硬质油生产线的生产能力约束（APR-HPC）!四月份软质油生产线的生产能力约束（APR-SPC）!四月份成品油的硬度上限约束（APR-HUB）!四月份成品油的硬度下限约束（APR-HLB）!四月份原料和成品的平衡约束（APR-BLS）APR-HPC) UHD14+UHD24200APR-SPC) USF14+USF24+USF34250APR-HUB)-6.0PROD4+8.8UHD14+6.1UHD24+2.0USF14+4.2USF24+5.0USF340APR-BLS) -PROD4+UHD14+UHD24+USF14+USF24+USF34=0!-!五月份硬质油生产线的生产能力约束（MAY-HPC）!五月份软质油生产线的生产能力约束（MAY-SPC）!五月份成品油的硬度上限约束（MAY-HUB）!五月份成品油的硬度下限约束（MAY-HLB）!五月份原料和成品的平衡约束（MAY-BLS）MAY-HPC) UHD15+UHD25200MAY-SPC) USF15+USF25+USF35250MAY-HUB)-6.0PROD5+8.8UHD15+6.1UHD25+2.0USF15+4.2USF25+5.0USF350MAY-BLS) -PROD5+UHD15+UHD25+USF15+USF25+USF35=0!-!六月份硬质油生产线的生产能力约束（JUN-HPC）!六月份软质油生产线的生产能力约束（JUN-SPC）!六月份成品油的硬度上限约束（JUN-HUB）!六月份成品油的硬度下限约束（JUN-HLB）!六月份原料和成品的平衡约束（JUN-BLS）JUN-HPC) UHD16+UHD26200JUN-SPC) USF16+USF26+USF36250JUN-HUB)-6.0PROD6+8.8UHD16+6.1UHD26+2.0USF16+4.2USF26+5.0USF360JUN-BLS) -PROD6+UHD16+UHD26+USF16+USF26+USF36=0!-!第一种硬质油（HD1）相邻两个月的采购量，耗用量和存储量的平衡HD1-J ) BHD11-UHD11-SHD11=-500HD1-JF) SHD11+BHD12-UHD12-SHD12=0HD1-FM) SHD12+BHD13-UHD13-SHD13=0HD1-MA) SHD13+BHD14-UHD14-SHD14=0HD1-AM) SHD14+BHD15-UHD15-SHD15=0HD1-MJ) SHD15+BHD16-UHD16-SHD16=0!-!第二种硬质油（HD2）相邻两个月的采购量，耗用量和存储量的平衡HD2-J ) BHD21-UHD21-SHD21=-500HD2-JF ) SHD21+BHD22-UHD22-SHD22=0HD2-FM ) SHD22+BHD23-UHD23-SHD23=0HD2-MA ) SHD23+BHD24-UHD24-SHD24=0HD2-AM ) SHD24+BHD25-UHD25-SHD25=0HD2-MJ ) SHD25+BHD26-UHD26-SHD26=0!-!第一种软质油（SF1）相邻两个月的采购量，耗用量和存储量的平衡 SF1-J ) BSF11-USF11-SSF11=-500SF1-JF )SSF11+BSF12-USF12-SSF12=0SF1-FM )SSF12+BSF13-USF13-SSF13=0SF1-MA )SSF13+BSF14-USF14-SSF14=0SF1-AM )SSF14+BSF15-USF15-SSF15=0SF1-MJ )SSF15+BSF16-USF16-SSF16=0!-!第二种软质油（SF2）相邻两个月的采购量，耗用量和存储量的平衡 SF2-J ) BSF21-USF21-SSF21=-500SF2-JF)SSF21+BSF22-USF22-SSF22=0SF2-FM)SSF22+BSF23-USF23-SSF23=0SF2-MA)SSF23+BSF24-USF24-SSF24=0SF2-AM)SSF24+BSF25-USF25-SSF25=0SF2-MJ)SSF25+BSF26-USF26-SSF26=0!-!第三种软质油（SF3）相邻两个月的采购量，耗用量和存储量的平衡SF3-J) BSF31-USF31-SSF31=-500SF3-JF)SSF31+BSF32-USF32-SSF32=0SF3-FM)SSF32+BSF33-USF33-SSF33=0SF3-MA)SSF33+BSF34-USF34-SSF34=0SF3-AM)SSF34+BSF35-USF35-SSF35=0SF3-MJ)SSF35+BSF36-USF36-SSF36=0!-!期末存储量约束SHD16=500SHD26=500SSF16=500SSF26=500SSF36=500!-end!存储量的上限约束SUB SHD11 1000SUB SHD12 1000SUB SHD13 1000SUB SHD14 1000SUB SHD15 1000SUB SHD16 1000SUB SHD21 1000SUB SHD22 1000SUB SHD23 1000SUB SHD24 1000SUB SHD25 1000SUB SHD26 1000SUB SSF11 1000SUB SSF12 1000SUB SSF13 1000SUB SSF14 1000SUB SSF15 1000SUB SSF16 1000SUB SSF21 1000SUB SSF22 1000SUB SSF23 1000SUB SSF24 1000SUB SSF25 1000SUB SSF26 1000SUB SSF31 1000SUB SSF32 1000SUB SSF33 1000SUB SSF34 1000SUB SSF35 1000SUB SSF36 10002.2、LINDO求解结果 LP OPTIMUM FOUND AT STEP 9 OBJECTIVE FUNCTION VALUE 1) 250731.5VARIABLEVALUEREDUCEDCOSTPROD1450.0000000.000000BHD110.00000025.000000BHD210.00000035.000000BSF110.00000025.000000BSF210.00000025.000000BSF310.00000030.000000PROD2450.0000000.000000BHD120.00000040.000000BHD220.00000040.000000BSF120.0000000.000000BSF22750.0000000.000000BSF320.00000025.000000PROD3450.0000000.000000BHD130.00000015.000000BHD230.00000045.000000BSF130.00000015.000000BSF230.0000005.000000BSF330.0000000.000000PROD4450.0000000.000000BHD140.00000020.000000BHD240.00000010.000000BSF140.0000000.000000BSF240.00000020.000000BSF340.00000025.000000PROD5450.0000000.000000BHD15159.2592620.000000BHD250.00000015.000000BSF150.00000025.000000BSF250.0000005.000000BSF350.0000000.000000PROD6450.0000000.000000BHD16659.2592770.000000BHD26381.4814760.000000BSF160.00000010.000000BSF26750.0000000.000000BSF360.00000025.000000SHD11340.7407530.000000SHD21459.2592470.000000SSF11500.0000000.000000SSF21250.0000000.000000SSF31500.0000000.000000SHD12181.4814760.000000SHD22418.5185240.000000SSF12500.0000000.000000SSF22750.0000000.000000SSF32500.0000000.000000SHD1322.2222210.000000SHD23377.7777710.000000SSF13500.0000000.000000SSF23500.0000000.000000SSF33500.0000000.000000SHD140.0000005.000000SHD24200.0000000.000000SSF14500.0000000.000000SSF24250.0000000.000000SSF34500.0000000.000000SHD150.00000015.000000SHD25159.2592620.000000SSF15500.0000000.000000SSF250.00000030.000000SSF35500.0000000.000000UHD11159.2592620.000000UHD2140.7407420.000000USF110.00000020.000000USF21250.0000000.000000USF310.0000000.000000UHD12159.2592620.000000UHD2240.7407420.000000USF120.00000020.000000USF22250.0000000.000000USF320.0000000.000000UHD13159.2592620.000000UHD2340.7407420.000000USF130.00000020.000000USF23250.0000000.000000USF330.0000000.000000UHD1422.2222210.000000UHD24177.7777710.000000USF140.00000020.000000USF24250.0000000.000000USF340.0000000.000000UHD15159.2592620.000000UHD2540.7407420.000000USF150.00000015.925926USF25250.0000000.000000USF350.0000001.481481UHD16159.2592620.000000UHD2640.7407420.000000USF160.00000033.703705USF26250.0000000.000000USF360.00000035.925926SHD16500.0000000.000000SHD26500.0000000.000000SSF16500.0000000.000000SSF26500.0000000.000000SSF36500.0000000.000000ROWSLACK OR SURPLUSDUAL PRICESJAN-HPC)0.000000115.000000JAN-SPC)0.000000115.000000JAN-HUB)0.0000000.000000JAN-HLB)1350.0000000.000000JAN-BLS)0.000000-200.000000FEB-HPC)0.000000110.000000FEB-SPC)0.000000110.000000FEB-HUB)0.0000000.000000FEB-HLB)1350.0000000.000000FEB-BLS)0.000000-200.000000MAR-HPC)0.000000105.000000MAR-SPC)0.000000105.000000MAR-HUB)0.0000000.000000MAR-HLB)1350.0000000.000000MAR-BLS)0.000000-200.000000APR-HPC)0.000000100.000000APR-SPC)0.000000100.000000APR-HUB)370.0000000.000000APR-HLB)980.0000000.000000APR-BLS)0.000000-200.000000MAY-HPC)0.00000094.814812MAY-SPC)0.00000098.333336MAY-HUB)0.0000001.851852MAY-HLB)1350.0000000.000000MAY-BLS)0.000000-211.111115JUN-HPC)0.00000089.259262JUN-SPC)0.000000133.333328JUN-HUB)0.0000007.407407JUN-HLB)1350.0000000.000000JUN-BLS)0.000000-244.444443HD1-J)0.000000-85.000000HD1-JF)0.000000-90.000000HD1-FM)0.000000-95.000000HD1-MA)0.000000-100.000000HD1-AM)0.000000-100.000000HD1-MJ)0.000000-90.000000HD2-J)0.000000-85.000000HD2-JF)0.000000-90.000000HD2-FM)0.000000-95.000000HD2-MA)0.000000-100.000000HD2-AM)0.000000-105.000000HD2-MJ)0.000000-110.000000SF1-J)0.000000-105.000000SF1-JF)0.000000-110.000000SF1-FM)0.000000-115.000000SF1-MA)0.000000-120.000000SF1-AM)0.000000-125.000000SF1-MJ)0.000000-130.000000SF2-J)0.000000-85.000000SF2-JF)0.000000-90.000000SF2-FM)0.000000-95.000000SF2-MA)0.000000-100.000000SF2-AM)0.000000-105.000000SF2-MJ)0.000000-80.000000SF3-J)0.000000-85.000000SF3-JF)0.000000-90.000000SF3-FM)0.000000-95.000000SF3-MA)0.000000-100.000000SF3-AM)0.000000-105.000000SF3-MJ)0.000000-110.00000062)0.000000-90.00000063)0.000000-110.00000064)0.000000-130.00000065)0.000000-80.00000066)0.000000-110.000000 NO. ITERATIONS= 9 RANGES IN WHICH THE BASIS IS UNCHANGED: OBJ COEFFICIENT RANGESVARIABLECURRENTALLOWABLEALLOWABLECOEFINCREASEDECREASEPROD1200.000000INFINITY115.000000BHD11-110.00000025.000000INFINITYBHD21-120.00000035.000000INFINITYBSF11-130.00000025.000000INFINITYBSF21-110.00000025.000000INFINITYBSF31-115.00000030.000000INFINITYPROD2200.000000INFINITY110.000000BHD12-130.00000040.000000INFINITYBHD22-130.00000040.000000INFINITYBSF12-110.00000015.9259250.000000BSF22-90.00000030.0000000.000000BSF32-115.00000025.000000INFINITYPROD3200.000000INFINITY105.000000BHD13-110.00000015.000000INFINITYBHD23-140.00000045.000000INFINITYBSF13-130.00000015.000000INFINITYBSF23-100.0000005.000000INFINITYBSF33-95.0000000.000000INFINITYPROD4200.000000INFINITY100.000000BHD14-120.00000020.000000INFINITYBHD24-110.00000010.000000INFINITYBSF14-120.0000000.000000INFINITYBSF24-120.00000020.000000INFINITYBSF34-125.00000025.000000INFINITYPROD5200.000000INFINITY94.814812BHD15-100.00000015.0000004.999999BHD25-120.00000015.000000INFINITYBSF15-150.00000025.000000INFINITYBSF25-110.0000005.000000INFINITYBSF35-105.0000000.000000INFINITYPROD6200.000000INFINITY89.259262BHD16-90.00000041.36362815.000000BHD26-110.0000004.99999910.000000BSF16-140.00000010.000000INFINITYBSF26-80.000000INFINITY30.000000BSF36-135.00000025.000000INFINITYSHD11-5.0000000.00000024.545452SHD21-5.00000024.5454520.000000SSF11-5.00000025.00000020.000000SSF21-5.0000000.000000INFINITYSSF31-5.00000030.0000000.000000SHD12-5.0000000.00000024.545452SHD22-5.00000024.5454520.000000SSF12-5.00000015.9259250.000000SSF22-5.0000000.0000000.000000SSF32-5.00000025.0000000.000000SHD13-5.0000000.0000000.000000SHD23-5.0000000.0000000.000000SSF13-5.00000015.9259250.000000SSF23-5.00000030.0000000.000000SSF33-5.0000000.0000000.000000SHD14-5.0000005.000000INFINITYSHD24-5.00000010.0000005.000000SSF14-5.00000015.92592510.000000SSF24-5.00000030.0000001.481480SSF34-5.0000001.4814800.000000SHD15-5.00000015.000000INFINITYSHD25-5.00000010.0000004.999999SSF15-5.00000033.70370110.000000SSF25-5.00000030.000000INFINITYSSF3