1、汽车租赁调度问题摘要国内汽车租赁市场兴起于 1900 年北京亚运会，随后在北京、上海、广州及深圳等国际化程度较高的城市率先发展直至 2000 年左右，汽车租赁市场开始在其他城市发展。为了对某市的一家租赁公司获利情况进行分析并确定汽车调度方案，本文我们以非线性规划为基础，通过 matlab，excel 等软件对数据进行处理，最小二乘法对缺失数据进行预测，最终使用 lingo 软件进行编程求解得到最终的优化方案。在问题一中，我们基于对题目中尽量满足需求的理解，考虑到总的车辆数和总的需求量之间的关系，用最小偏差法和分段考虑法进行了计算，分别建立多目标规划模型和非线性规划模型，通过对转运后各代理点最终
2、的车辆数进行分析，比较两种结果得到更优的转运方案。在问题二中，我们一方面要对其短缺损失进行理解，另一方面要考虑，是否应该考虑在尽量满足需求的条件下求其最低的转运费用和短缺损失，此问题中我们同样分两种情况对其进行考虑，通过比较两者最低费用并且结合实际情况，得到更合理的转运方案。在问题三中，首先我们分析数据，剔除了其中一场的部分，并用最小二乘法对缺失数据进行预测，得到完整的单位租赁费用与短缺损失费用，然后综合考虑各种因素后，我们将公司获利最大作为最终目标函数通过非线性规划的模型求得最佳方案。在问题四中，我们没有直接对是否购买新车作出判断，而是直接以其八年获利最大为目标进行非线性规划，购买的车辆数成
3、为其目标函数中的一个未知数，用 lingo 可直接求得在获利最大时的购车数量，将其与不购车时的利润进行比较可得到最佳的购买方案。关键词：非线性规划 全局最优 短缺损失 最小二乘法 一 问题重述国内汽车租赁市场兴起于 1990 年北京亚运会，随后在北京、上海、广州及深圳等国际化程度较高的城市率先发展，直至 2000 年左右，汽车租赁市场开始在其他城市发展。某城市有一家汽车租赁公司，此公司年初在全市范围内有 379 辆可供租赁的汽车，分布于 20 个代理点中。每个代理点的位置都以地理坐标 X 和 Y 的形式给出，单位为千米。假定两个代理点之间的距离约为他们之间欧氏距离（即直线距离）的 1.2 倍。
4、要求根据附件所给数据计算如下问题：1给出未来四周内每天的汽车调度方案，在尽量满足需求的前提下，使总的转运费用最低；2考虑到由于汽车数量不足而带来的经济损失，给出使未来四周总的转运费用及短缺损失最低的汽车调度方案；3综合考虑公司获利、转运费用以及短缺损失等因素，确定未来四周的汽车调度方案；4为了使年度总获利最大，从长期考虑是否需要购买新车？如果购买的话，确定购买计划（考虑到购买数量与价格优惠幅度之间的关系，在此假设如果购买新车，只购买一款车型） 。二问题分析汽车租赁调度问题是一个典型的数学规划问题，需要综合考虑转运费用，短缺损失，公司获利等多方面因素，在掌握了各代理点实际需求下，根据一定要求，寻
5、找到使目标函数满意的优化解。问题一中，要求在尽量满足需求的前提下，使未来四周的总转运费用最低。对数据进行处理后，对尽量满足需求这一约束条件，认为其在需求量大于供应量时应保证每辆车都能够被利用，在需求量小于供应量时应保证每个代理点的需求都能被满足。然后据此约束建立多目标规划模型求全局最优解，使得未来四周总的转运费用最小。针对问题二，我们需要考虑在汽车数量不足的情况下所带来的短缺损失，所谓短缺损失是指，在某代理点某天经过转运后最终的车辆数比需求量少时，少的车辆数与单位短缺损失的乘积。在此基础上建立两种模型，第一种是尽量满足需求条件下的模型，第二种是不考虑尽量满足需求这一条件下的模型。然后分别建立非
6、线性规划模型求全局最优，使得未来四周的转运费与短缺损失之和最小。针对问题三，综合考虑公司获利、转运费用以及短缺损失等因素，以公司获利最多作为目标函数，考虑到前期尽量满足需求对公司后续的租赁需求影响，在此仅分析在尽量满足需求条件下获利最多。对于附录中丢失的数据，我们将平均需求量与租赁收入之间的关系曲线采用最小二乘法进行拟合，预测出缺失的数据以及异常数据。最后将其考虑为非线性规划问题对其进行规划求全局最优，得到最佳的调度方案。针对问题四，由于一年中最大需求量要比实际供应量多 66 辆车，故我们将购买车的数量 m 取小于 66 的值，然后分别计算每增加一辆能够获得的最大的利润，然后求得最优的 m 值
7、，该 m 的取值区间会有一个值使得获利最大。由于车型不影响租赁收入，所以在考虑车型时，选择是 8 年成本和维修费用之和最低的一款。三符号说明kix第 k 天第 i 个代理点转运之后最终的车辆数kia第 k 天第 i 个代理点的需求车量数kijn第 k 天第 i 个代理点转到第 j 个代理点的车辆数ijp第 i 个代理点到第 j 个代理点转运一辆的运费ib第 i 个代理点的单位短缺损失（万元/天*辆）A该公司拥有的总的车辆数kB每天所有代理点总的需求量e 选中的车型每辆总的花费m 需要购买的车辆数kid第 k 天第 i 个代理点短缺的车辆数四模型假设1.假设租赁车辆不会损坏，且不会产生维修保养费
8、用。2.假设当天租出去的车会当天归还，不影响第二天租赁。3.假设每次车辆转运发生在一天的结束后，第二天之前。4.假设附件 2 所给一年各代理点的汽车需求量代表未来八年的汽车需求量。5.假设购买新车的周期为 8 年。6.假设价格不考虑涨价等情况。7.假设前期的不满足需求不会影响到后续的需求量。五模型建立与求解5.1 问题一5.1.1 对问题一的理解问题一要求在尽量满足需求的前提下，使总的转运费用最低。对于尽量满足需求，我们对其有两种理解。一是使每天每个代理点转运后最终的车辆数与其需求量的偏差最小。二是认为其在需求量大于供应量时应保证每辆车都能够被利用，在需求量小于供应量时应保证每个代理点的需求都
9、能被满足。5.1.2 基于偏差最小的多目标规划模型的建立与求解首先用 matlab 对附件 1 和附件 6 中数据进行处理，得到两两代理点之间每转运一辆车的转运费用。具体结果见附件 1。用 表示其偏差，建立多目标规划如下：kiixamin 2901kiikixas.t. 1379nkiiX上式可求得当其偏差和最小时每天每个代理点经过转运后的最终车辆数。在此基础上以其转运费用最低为目标函数建立如下模型：min 2901kijkijnps.t. 2020()11kijkjikikijj x利用 lingo 软件编程解得最小的转运费用为 70.4987 万元，以下是前 11 天各代理点转运之后最终的
10、车辆数。由下表数据可知，在该模型下，虽然大部分代理点几乎完全满足需求，但是一些代理点经过转运之后一辆车也没有，这违背了尽量满足需求这一条件，也不符合实际情况，同时求解得到其运输费用最小为 70.4987 万元，远高于第二个模型的最小运输费用，所以该模型被舍弃。表 5.1.2 前 11 天各代理点转运之后最终的车辆数日期代理点1 2 3 4 5 6 7 8 9 10 11A 22 15 9 8 0 22 19 12 13 28 26B 18 22 28 17 23 19 39 0 29 20 16C 19 22 25 21 30 54 27 25 0 25 14D 18 27 15 20 20
11、20 20 12 15 11 0E 24 15 28 0 17 15 12 28 18 13 21F 16 20 24 18 19 13 15 29 14 17 20G 19 15 25 22 21 12 22 16 24 22 14H 17 4 19 14 18 26 17 24 23 13 28I 22 19 0 17 25 25 18 20 27 14 11J 15 16 0 18 0 23 19 30 14 29 11K 18 27 17 11 23 16 16 30 27 27 12L 23 24 18 29 13 17 14 12 16 22 23M 14 30 26 23 20
12、11 20 15 29 19 30N 18 13 19 27 15 16 12 15 11 12 18O 18 17 15 15 20 14 15 13 14 13 19P 17 24 30 28 29 11 30 17 18 23 16Q 21 16 28 18 30 11 30 12 15 14 25R 23 13 12 15 25 14 11 18 22 15 29S 18 12 28 28 30 14 12 21 22 16 25T 19 28 13 30 1 26 11 30 28 26 215.1.3 基于分段考虑的非线性规划模型的建立与求解对该公司拥有的总的车辆数和总的需求量进行
13、比较，通过对两者大小的判断，以此述判断为分段约束条件，直接以转运费用最低位目标函数建立非线性规划模型如下：min 2901kijkijnps.t. 2020()11kijkjikikijj x当 时 kABkiia当 时 kkikix利用 lingo 软件编程对该模型进行求解得，最小的转运费用为 40.4916 万元，下表给出前 11 天各代理点转运之后的最终车辆数，完整表格见附件 3：表 5.1.3 前 11 天各代理点转运之后最终的车辆数1 2 3 4 5 6 7 8 9 10 11A 22 15 15 15 14 22 21 15 22 28 28B 18 22 22 17 21 19
14、12 11 25 20 11C 19 22 22 21 27 27 27 25 11 25 14D 18 26 15 20 20 20 20 12 12 11 12E 24 15 15 12 17 15 15 18 18 13 21F 16 20 24 18 18 14 15 29 14 17 20G 19 15 21 21 21 13 22 16 24 22 14H 17 12 15 14 18 26 24 24 23 13 22I 22 19 18 17 25 25 20 20 24 14 11J 15 16 13 18 13 23 19 26 14 29 11K 18 21 17 11
15、12 16 16 30 27 27 12L 23 23 18 18 13 17 14 12 16 22 22M 1 3 2 2 2 1 2 1 2 1 34 0 1 3 0 2 0 5 0 9 0N 18 13 19 20 15 16 14 15 11 12 18O 18 17 15 15 15 15 15 13 13 13 17P 17 24 28 28 29 25 30 17 18 23 16Q 21 16 28 18 24 14 30 12 15 14 25R 23 13 12 15 15 14 11 18 22 15 29S 18 12 28 28 28 20 20 21 22 16
16、25T 19 28 13 30 14 26 14 30 28 26 21由上表数据可以看出，该模型在尽量满足需求的条件下分配的也比较合理，远远优于第一种模型，而且转运费用也远远低于第一种模型。转运方案：1-20 分别代表 A-T20 个租赁代理点第一天：12(7) 213(3) 510(9) 74(5) 84(1) 820(4) 911(3) 103(3) 106(4) 10 7(1) 1413(5) 154(1) 1613(2) 1720(5)184(1) 1816(9) 1 913(6)第二天：47(2) 414(9) 910(1) 107(4) 116(4) 1214(5) 1316(4
17、) 1319(5) 1419(8) 119(2) 18 19(1) 208(3) 2017(12)第十四天：411(12) 611(2) 820(3) 910(16) 1413(2) 1416(4) 152(6) 153(2) 1 5 4(5) 1514(1) 173(3) 1918(1) 203(2)完整转运方案可见附件 2分析第十四天的转运方案在代理点 4 既有转入又有转出，表面上看如此周折会产生多余费用，实际上这样是节省了转运费用，由附件 1 可知，1 511 转运费用0.04 余万元每千米，而 1 54 再从 411 转运费用是 0.03 余万元每千米。因此如此周转节省了转运费用。5.
18、1.4 模型比较通过对以上两种方案最小转运费用的的比较，发现第二种基于全局优化的非线性规划模型得到的最小转运费用远远小于第一种模型，并且由第一种方案得到的转运后的最终车辆数在一部分代理点中出现了零，这是不符合实际情况的，也偏离了尽量满足需求这一要求。所以我们选择了第二种方案作为最终的调度方案。5.2 问题二5.2.1 对问题二的理解问题二要求在考虑短缺损失的情况下，求得使四周总的转运费用及短缺损失最低的最佳汽车调度方案。本题我们同样分两种情况考虑，第一种是考虑在尽量满足需求的前提进行求解，即在问题一的基础上保留使得尽量满足需求的约束条件。第二种是不考虑尽量满足需求即在问题一的基础上去掉使得尽量
19、满足需求的约束条件。并且应该明确的是：在尽量满足需求的前提下，只有在该公司拥有的总的车辆数小于总的需求量时才存在短缺损失，并且是转运费用与短缺损失之和最低。由此我们可以以和费用最低为目标函数建立如下模型。5.2.2 考虑尽量满足需求时模型的建立与求解考虑尽量满足需求，即要对该公司拥有的总的车辆数和总的需求量之间进行比较，在需求量大于供应量时应保证每辆车都能够被利用，在需求量小于供应量时应保证每个代理点的需求都能被满足。建立非线性规划模型如下：29029011minbkij kikij kinpds.t. 2020()11kijkjikikijjnx当 时 kiiaxkikikida当 时 ki
20、i0ki当 时 kABkiixa当 时 kkiki利用 lingo 软件编程对该模型进行求解，最小的转运费用与短缺损失共为 70.1639万元5.2.3 不考虑尽量满足需求时的模型建立与求解不考虑尽量满足需求，即使得损失与运输费之和最小而不用考虑尽量满足需求这一约束条件，建立非线性规划模型如下：29029011minbkij kikij kinpds.t. 2020()11kijkjikikijj x当 时 kiiaxkikikida当 时 kii0ki运用 lingo 软件编程对该模型进行求解，最小的转运费用为 64.2085 万元。汽车调度方案见附件 4通过分析各代理点的转运费用和短缺损失
21、费用，可以知道造成两种模型的差价原因，由数据可知，部分代理点中从一个代理点到另一个代理点的转运费用是大于后者代理点的短缺损失费用的。因此在尽量满足需求的前提下，调动进行的多一些，忽略了差价，使得总的费用高于第二种模型。5.2.4 模型比较通过比较以上两种情况下的最小转运费用，我们会发现考虑尽量满足需求与不考虑尽量满足需求的最小转运费用相差并不大，我们综合考虑能让每个代理点尽可能的正常运行，该公司的信誉等多方面因素，考虑尽量满足需求这种情况更加合理。5.3 问题三5.3.1 对问题三的理解问题三要在综合考虑公司获利、转运费用以及短缺损失等因素的情况下，确定未来四周的汽车调度方案，使总的获利最大化
22、，同样需要考虑尽量满足需求，在此与问题 2 相似，仅考虑尽量满足需求这一情况，其余不再赘述。由于附件 5所给的单位租赁收入数据不完整，我们需要对其缺少的数据进行预测，从而计算其总的租赁收入。5.3.2 基于最小二乘法的预测模型对于该题中的数据缺失，需要对其数据进行预测，我们认为其租赁收入预期地理位置有关，用 MATLAB 做出代理点分布图如下，对地理位置邻近的代理点的租赁收入分析后发现并无必然联系。图 5.3.2.1再次经过对数据的比对分析后我们认为其与未来四周各代理点平均每天的需求量有一定关系，计算出每个代理点前四周平均需求量与租赁收入表如下表 5.3.2.1代 理 点 A B C D E
23、F G H I J租 赁 收 入 （ 万 元 /天 辆 ） 0.43 0.39 0.35 0.49 0.3 0.41 0.395 0.36 0.37 0.4各 代 理 点 平 均 需 求 20.5517218.6896620.1379318.82759 17.206920.0689719.1724118.6206918.2413819.51724代 理 点 K L M N O P Q R S T租 赁 收 入 （ 万 元 /天 辆 ） 0.33 0.379 0.365 0.4 0.42各 代 理 点 平 均 需 求 22.3793119.4137921.41379 16.7931 16.206
24、922.0344818.1724118.0344818.6551723.24138其数据基本满足线性关系，在 MATLAB 下用最小二乘法分别进行一次函数拟合，二次函数拟合，三次函数拟合， 可得其拟合函数均为一次线性关系，如下图图 5.3.2.2所以采用一次拟合函数补全缺失数据如下表表 5.3.2.2代 理 点 A B C D E F G H I J租 赁 收 入 （ 万 元 /天 辆 ） 0.43 0.39 0.35 0.49 0.3 0.41 0.395 0.36 0.37 0.4各 代 理 点 平 均 需 求 20.5517218.6896620.1379318.8275917.2069
25、20.0689719.1724118.6206918.2413819.51724代 理 点 K L M N O P Q R S T租 赁 收 入 （ 万 元 /天 辆 ） 0.33 0.379 0.365 0.4 0.42 0.3748 0.3897 0.3902 0.3878 0.3702各 代 理 点 平 均 需 求 22.3793119.4137921.4137916.7931 16.206922.0344818.1724118.0344818.6551723.241385.3.3 基于公司获利最多的模型建立与求解考虑到公司获利要多，转运费用及短缺损失要尽量少，我们建立非线性规划模型如下
26、：290290290111max bki kij kiki kij kiXCnpds.t. 2020()11kijkjikikijjnx当 时 kiiaxkikikida当 时 kiiax0kid当 时 kABkiixa当 时 kkiki利用 lingo 软件编程对该模型进行求解，可以得到在尽量满足需求的条件下公司在这四周内获利最多为 4103.29 万元。5.4 问题四5.4.1 对问题四的理解首先我们要考虑购买新车是否能够使获利增加。其次，如果购买，则要对其购买哪种车型及购买多少辆进行确定。在此我们假设购买了 m 辆车，使其总车辆数达到 379+m 辆，然后对其最大获利进行计算，从而可得到
27、最大获利时 m 的值。如果 m=0 则表示购买车并不能够使获利增加。5.4.2 问题四模型的建立与求解a.对车型的选择通过 excel 对附件 4 中数据进行处理，得到购买每种车型八年所需总费用。结果如下表：表 5.4.2.1车型 1 2 3 4 5 6 7 8 9 10八年总花费（万元）43.1946.6648.7544.2758.353.354.7141.2160.2456.68由上表可得第 8 种车型八年总的花费最低，因为不管车好车坏我们都只考虑其只能使用 8 年，且不考虑八年后车辆的残存费用，较贵的车型并不能使得租赁费用增加或使需求量有所改变，故在对车型的选择时我们只考虑车费和 8 年
28、总的维修费之和最低，根据上表，如需购买车辆，则选择第 8 种型号的车。b.模型建立假设购买 m 辆第 8 种型号的车，由附件 2 中一年的数据分析可得，需求量最多时缺少 66 辆车，所以购买车辆数 m 的上限为 66 辆，为使其八年总的获利最多，即建立目标函数使得平均每年的的获利最多，以减少计算的数据量，将每辆车的总费用平均分摊到每一年，建立非线性规划模型如下： 365203652036520111ax( b)8ki kij kiki kij kiXCnpdems.t. 2020()11kijkjikikijj x0m当 时 kiiaxkikikidax当 时 kii0kic.模型求解用 li
29、ngo 软件编程对该模型进行求解，得到最大获利与够买车辆数量关系图如下图 5.4.2.1可得出该公司八年中平均每年获利最大值为 52855.52 万元，获利最大时对应 m值为 46 即购买 46 辆第 8 种型号的车能够使每年获利达到最大值 52855.52 万元，所以应该购买新车，购买 46 辆第八种型号的车可以使总获利最大。六.模型评价与改进问题一通过比较对尽量满足需求的两种理解，得出分段考虑总的车辆数和总的需求量大小的结果更优，但我们并不知道该公司具体对尽量满足需求的要求。在本题中分别建立了多目标规划模型和非线性规划模型，用 lingo 求解运行时间过长，如果可以建立线性规划模型那么可以
30、大大缩短运行时间。问题二同样建立非线性规划模型，对转运费用与短缺费用之和求最小值，模型简单明了，易于执行，并且考虑了尽量满足需求与不尽量满足需求两种情况，通过比较分析，得到的结果更加符合实际。问题三巧妙地使用最小二乘法求得缺失的数据，虽然不是非常精确，但方法简单，执行方便，可靠程度较高。通过对问题的分析，将获利最大作为最终目标进行规划求解，简化了求解过程。考虑到预测数据的精确性，本题还可用 BP 神经网络，综合多个因素影响，结果会更加可靠。问题四引入变量购买车辆数将两问结合成一个非线性规划问题，一次性可以解决两个问题，模型简单，可行性高，便于推广。七参考文献1李庆杨等，数值分析（第四版） ，华
31、中科技大学出版社， 2006.72胡运权，运筹学教程（第三版）,清华大学出版社，2007.43吴建国等，数学建模案例精编，中国水利水电出版社，20054姜启源等，数学模型（第三版） ，高等教育出版社， 20035赵东方，数学模型与计算，科学出版社，2007.2八附录程序问题 1：ti1.lg4model:sets:i/120/:shu;k/129/:shuzi;lin(k):c;link(k,i):a,b;links(i,i,k):p;linkp(i,i):yun;endsetsmin=sum(k(z):sum(i(x):sum(i(y):p(x,y,z)*yun(x,y);for(k(z)|
32、z#le#28:for(i(y):sum(i(x):p(y,x,z+1)-p(x,y,z+1)=b(z,y)-b(z+1,y); for(k(z):c(z)=sum(i(x):a(z,x);for(link(z,x):if(c(z)#gt#379,b(z,x),-b(z,x)=if(c(z)#gt#379,a(z,x),-a(z,x); for(k(t):sum(i(x):b(t,x)=379); for(link:gin(b); for(links:gin(p);data:enddataendti1.m%计算两两代理点之间的距离clear all;clc;%读取坐标coordinate=xl
33、sread(C:UserslifengDesktop2014模拟练习2A题汽车租赁调度问题附件1：代理点的位置及年初拥有车辆数.xls, Sheet1,B3:U4);distance=zeros(20);length=size(distance);% 为距离矩阵赋权值for i=1:lengthfor j=1:idistance(i,j)=sqrt(coordinate(1,i)-coordinate(1,j)2+(coordinate(2,i)-coordinate(2,j)2);endend% 计算每辆车在代理点间的转运价格%单位转运成本cost=xlsread(C:Userslifeng
34、Desktop2014模 拟练习2A 题汽车租赁调度问题附件6：不同代理点之间的转运成本.xls ,Sheet1,C3:V22);price=cost.*distance*1.2;for i=1:lengthfor j=i:lengthprice(i,j)=price(j,i);endendxlswrite(C:UserslifengDesktop2014模拟练习2代码每辆车在代理点间的转运价格.xls ,price,C3:V22);问题 2：ti2.lg4!第二题;model:!i代表站点，k代表日期，lost代表站点的单位损失，c代表每天的需求总量，a，b，d分别为每个站点需求量,实际量和
35、缺口量，p为转移数量，yun代表运费;sets:i/120/:shu;k/129/:shuzi;li(i):lost;lin(k):c;link(k,i):a,b,d;links(i,i,k):p;linkp(i,i):yun;endsetsmin=sum(k(z):sum(i(x):sum(i(y):p(x,y,z)*yun(x,y)+sum(link(z,x):d(z,x)*lost(x);for(k(z)|z#le#28:for(i(y):sum(i(x):p(y,x,z+1)-p(x,y,z+1)=b(z,y)-b(z+1,y); for(link(z,x):d(z,x)=if(a(z
36、,x)#gt#b(z,x),a(z,x)-b(z,x),0);for(k(z):c(z)=sum(i(x):a(z,x);for(link(z,x):if(c(z)#gt#379,b(z,x),-b(z,x)=if(c(z)#gt#379,a(z,x),-a(z,x); for(k(t):sum(i(x):b(t,x)=379); for(link:gin(b); for(links:gin(p);data:enddataend问题 3：ti3.lg4!第三题;model:!i代表站点，k代表日期，lost代表站点的单位损失，c代表每天的需求总量，a，b，d分别为每个站点需求量,实际量和缺口量
37、，p为转移数量，yun代表运费;sets:i/120/:shu;k/129/:shuzi;li(i):lost,win;lin(k):c;link(k,i):a,b,d,e;links(i,i,k):p;linkp(i,i):yun;endsetsmax=sum(link(z,x):e(z,x)*win(x)-sum(k(z):sum(i(x):sum(i(y):p(x,y,z)*yun(x,y)-sum(link(z,x):d(z,x)*lost(x);for(k(z)|z#le#28:for(i(y):sum(i(x):p(y,x,z+1)-p(x,y,z+1)=b(z,y)-b(z+1,
38、y); for(link(z,x):e(z,x)=if(a(z,x)#gt#b(z,x),b(z,x),a(z,x);for(link(z,x):d(z,x)=if(a(z,x)#gt#b(z,x),a(z,x)-b(z,x),0);for(k(z):c(z)=sum(i(x):a(z,x);for(link(z,x):if(c(z)#gt#379,b(z,x),-b(z,x)=if(c(z)#gt#379,a(z,x),-a(z,x); for(k(t):sum(i(x):b(t,x)=379); for(link:gin(b); for(links:gin(p);data:enddatae
39、nd问题 4：ti4.lg4!第四题;model:!i代表站点，k代表日期，lost代表站点的单位损失，c代表每天的需求总量，a，b，d分别为每个站点需求量,实际量和缺口量，p为转移数量，yun代表运费;sets:i/120/:shu;k/1365/:shuzi;li(i):lost,win;lin(k):c;link(k,i):a,b,d,e;links(i,i,k):p;linkp(i,i):yun;endsetsmax=sum(link(z,x):e(z,x)*win(x)-sum(k(z):sum(i(x):sum(i(y):p(x,y,z)*yun(x,y)-sum(link(z,x
40、):d(z,x)*lost(x)-q*5.15125;for(k(z)|z#le#28:for(i(y):sum(i(x):p(y,x,z+1)-p(x,y,z+1)=b(z,y)-b(z+1,y); for(link(z,x):e(z,x)=if(a(z,x)#gt#b(z,x),b(z,x),a(z,x);for(link(z,x):d(z,x)=if(a(z,x)#gt#b(z,x),a(z,x)-b(z,x),0);for(k(z):c(z)=sum(i(x):a(z,x);for(link(z,x):if(c(z)#gt#379+q,b(z,x),-b(z,x)=if(c(z)#gt#379+q,a(z,x),-a(z,x); for(k(t):sum(i(x):b(t,x)=379+q); for(link:gin(b); for(links:gin(p);gin(q);data:enddataend数据附件 1：每辆车在代理点间的转运价格.xls附件 2：问题 1 转运方案.xlsx附件 3：问题 1 转运后的.xlsx附件 4：问题 2 转运方案.xlsx附件 5：问题 3 转运方案.xlsx
