小标
2019-01-23
来源 :
阅读 1393
评论 0
摘要:本文主要向大家介绍了【云计算】蔬菜配送方案制定与实施,通过具体的内容向大家展现,希望对大家学习云计算有所帮助。
本文主要向大家介绍了【云计算】蔬菜配送方案制定与实施,通过具体的内容向大家展现,希望对大家学习云计算有所帮助。
这是今年三月底,人生中交的第一个数学建模作业……现在看来确实很幼稚,问题三的约束条件可能都没写对。相比之下,毕竟是第一篇,也抄了网上的,自己把所有的地方弄明白也很不容易,所以写(抄)的挺好的。嘿嘿
某连锁公司主营蔬菜加工与销售,拟在某省县级及以上城镇继续设立销售连锁店扩大经营规模。该省县级及以上城镇地理位置及道路连接见数据文件:全省交通网络数据.xlsx
问题:
1、目前公司现有2个生产基地、23家销售连锁店,生产基地设在120号和63号城镇,为23家连锁店提供蔬菜,连锁店的日销售量见附录1。若运输成本为0.8元/吨公里,请你为公司设计生产与配送方案,使运输成本最低。
2、公司收集了近5年全省各城镇的蔬菜月度需求数据(文件:各城镇月度需求数据.xlsx)请你分析各城镇需求特征,并预测未来数年,何时全省蔬菜需求达到峰值,达到峰值时需求达到前5位和后5位的城镇是那些?
3、通过广告宣传等手段,未来几年公司在全省的市场占有率可增至3成左右(各城镇对公司产品每日需求预测数据见文件:公司未来各城镇每日需求预测数据.xlsx),调查还发现,公司产品的需求量与销售量并不完全一致,若在当地(同一城镇)购买,则这一部分需求量与销售量相同,若在不足10公里的其他城镇的销售连锁店购买,则这一部分需求量只能实现一半(成为公司产品销售量,由于距离的原因,另一半需求转向购买其他公司或个体工商户的产品),而在超过10公里的其他城镇的销售连锁店购买,销售量只能达到需求量的三成。于是,公司决定在各城镇增设销售连锁店,基于现有条件、成本等的考虑,原有的23家销售连锁店销售能力可在现有销售量的基础上上浮20%,增设的销售连锁店销售能力控制在每日10吨至30吨内,并且要求增设的销售连锁店的销售量必须达到销售能力的下限。同一城镇可设立多个销售连锁店。
请你为公司设计增设销售连锁店方案,使全省销售量达到最大。
结合本题的实际,为了确定模型求解的准确性和合理性,我们排除了一些因素的干扰,提出以下几点假设:
1、运输过程中不考虑分散运输类型,每个销售店只有一辆配送车配送且载货量为无限大;
2、不考虑运输时间对运输成本的影响;
3、运输过程中不产生其他费用,即生产基地到连锁店的距离即为生产基地所在城镇到连锁店所在城镇的距离;
4、设生产基地生产量与配送销售量相等;
5、假设公司环境未来几年无变化;
6、每个连锁店只有一个生产基地来配送。
为了便于问题的求解,给出以下符号说明:
序号 | 符号 | 符号说明 |
1 | An | 各销售店与生产基地最短距离 |
2 | vi | 城镇序号 |
3 | mi | 各销售店日销售量 |
4 | x | 各城镇蔬菜需求量平均值 |
5 | σ | 各城镇蔬菜需求量标准差 |
本问题要求设计运输成本最低的生产与配送方案,建立Floyd模型画出全省交通网络拓扑图。首先求解出所有城镇之间的最短路径,再建立0-1矩阵,选取距离销售连锁店最近的生产基地,即得到各个生产基地的配送方案,将距离加和,之后求解出运输成本即为最低运输成本。
编写MATLAB程序求解画出全省交通网络拓扑图,如图1。代码见附录1。
图 1全省交通网络拓扑图
4.1.2 计算最短距离
建立Floyd模型对问题一进行求解。建模过程如下:
求解最短路径:全省交通网络共有154个位置点,现需要求从生产基地到销售连锁店的的最小运输成本。
设邻接矩阵A0
(4-1)
来存放各边长度,其中:aii=0,1,2,…,n。无向图Ao是对称矩阵,aij=aji。递推产生一个矩阵序列A0,A1,…,An,Akij表示城镇iv到城镇vj路径上所经过的顶点序号不大于k的最短路径的长度。计算时用迭代公式:
(4-2)
k为迭代次数,i,j,k=1,2,…,n。
最后,当k=n时,An即是各城镇之间最短路径值。
选取最近生产基地:将具有销售店的城镇挑选出来,一共23个城镇。
引入0-1变量xi(1,2,...)
令
则运输成本为
Y=0.8miAn (4-3)
其中mi为各销售店日销售量,An为各销售店与生产基地最短距离。
已假设生产量与销售量相等。通过编写MATLAB程序得到各生产基地配送城镇分配方案。对比63和120号基地到连锁店的距离分成两个生产基地分配方案,并用excel计算费用。程序见附录2。
表 163号基地
城镇编号 | 连锁店编号 | 距离 | 日销售量 | 费用 |
10 | 6 | 108.36 | 8481 | 735.200928 |
11 | 16 | 179.15 | 6103 | 874.68196 |
22 | 20 | 168.95 | 6375 | 861.645 |
24 | 17 | 128.94 | 3251 | 335.347152 |
27 | 12 | 135.1 | 9265 | 1001.3612 |
31 | 4 | 114.66 | 23947 | 2196.610416 |
63 | 3 | 0 | 21733 | 0 |
63 | 18 | 0 | 28295 | 0 |
64 | 23 | 7.31 | 1840 | 10.76032 |
65 | 7 | 19.09 | 15570 | 237.78504 |
79 | 8 | 28.17 | 38759 | 873.472824 |
表 2120号基地
城镇编号 | 连锁店编号 | 距离 | 日销售量 | 费用 |
1 | 9 | 134.31 | 14744 | 1584.213312 |
16 | 21 | 103.64 | 14783 | 1225.688096 |
34 | 13 | 119.54 | 451 | 43.130032 |
36 | 11 | 151.19 | 11503 | 1391.310856 |
42 | 14 | 110.58 | 9489 | 839.434896 |
94 | 15 | 170.17 | 12773 | 1738.865128 |
106 | 2 | 63.7 | 38223 | 1947.84408 |
120 | 1 | 0 | 28733 | 0 |
120 | 10 | 0 | 32517 | 0 |
123 | 22 | 5.11 | 18081 | 73.915128 |
141 | 5 | 61.72 | 9258 | 457.123008 |
145 | 19 | 72.85 | 39653 | 2310.97684 |
(1)63号城镇生产基地的货物配送至11个销售连锁店,配送成本为7126.86484
元/天,生产量为163619公斤/天。
(2)120号城镇生产基地的货物配送至12个销售基地,配送成本为11612.50138
元/天,生产量为230208公斤/天。
(3)两个生产基地总运输成本为18739.36622元/天,即为最低运输成本。
具体路线即为生产基地到连锁店的最短距离路线。两个生产基地的配送路线见表3。
表 3配送路线
生产基地 | 配送路线 | 城镇编号 | 连锁店编号 |
63 | 63→6→51→50→10 | 10 | 6 |
63→68→69→3→21→23→11 | 11 | 16 | |
63→68→69→3→21→22 | 22 | 20 | |
63→68→69→3→24 | 24 | 17 | |
63→6→51→50→10→29→28→27 | 27 | 12 | |
63→6→51→50→10→31 | 31 | 4 | |
63→63 | 63 | 3 | |
63→63 | 63 | 18 | |
63→64 | 64 | 23 | |
63→64→65 | 65 | 7 | |
63→6→66→79 | 79 | 8 | |
120 | 120→125→124→133→132→131→15→142→143→16→1 | 1 | 9 |
120→125→124→133→132→131→15→142→143→16 | 16 | 21 | |
120→119→13→45→44→43→34 | 34 | 13 | |
120→119→13→45→40→41→42→35→36 | 36 | 11 | |
120→119→13→45→40→41→42 | 42 | 14 | |
120→125→131→106→91→90→2→84→93→94 | 94 | 15 | |
120→125→131→106 | 106 | 2 | |
120→120 | 120 | 1 | |
120→120 | 120 | 10 | |
120→123 | 123 | 22 | |
120→123→124→133→132→141 | 141 | 5 | |
120→123→134→139→149→146→145 | 145 | 19 |
本问题要求分析各城镇需求特征并预测全省鲜蔬菜何时达到峰值,以及达到峰值时排序前5位和后5位的城镇,在分析特征时,拟求解各城镇蔬菜需求量的平均值、标准差,由于标准差不能完全反映需求特征,因而再求解变异系数,分析这三个数据体现各城镇需求量特征。拟建立拟合模型,对全省各月总需求量进行拟合,求解出何时达到峰值,再对每一个城镇需求量进行拟合,求解出在全省需求量达到峰值的月份各城镇蔬菜需求量,之后对各城镇需求量排序,得到排名前五和后五的城镇。
通过均值和方差以及变异系数的计算,可以看出每个镇对蔬菜的需求量及其波动情况。
(4-4)
(4-5)
<img alt="c*Z喎�" https:="" www.2cto.com="" kf="" yidong="" wp="" "="" target="_blank" class="keylink">WPVxmcmFje3h9e1xzaWdtYSB9" class="mathcode" src="https://www.2cto.com/uploadfile/Collfiles/20180810/2018081009482677.gif" /> (4-6)
4.2.2拟合模型
①对同一月中每个城镇的蔬菜需求量进行累加。
(4-7)
其中i代表第i个城镇;j代表第j月。
②多项式拟合。
设f(x)是一个定义在实轴上的的函数,f(x1),f(x2),…,f(xm)在点x1<x2<…<xm上的函数值。又设pi是某个空间函数。所谓曲线拟合就是要在函数空间pi中找到一个函数φ(x),使得< p="">
(4-8)
其中是权重系数,特殊地,可取就称为在点x1<x2<…<xm上对f(x)的拟合曲线。< p="">
通常,我们要求函数空间π是一个有限维的线性空间,即存在着一组线性无关的函数 使得
(4-9)
于是,求拟合曲线φ(x)就可以化简为求一组实数a1,a2,…,aN,使得下面目标函数达到最小值。
(4-10)
通过偏导数可以得到线性方程组:
(4-11)
其中,
曲线拟合的系数a1,a2, ... ,an 可通过求解正规方程得到。
只考虑多项式拟合,也即
并令,于是,
(4-12)
这时拟合曲线,也即为拟合多项式。
预测出拟合多项式的峰值,就是全省鲜蔬菜需求量的峰值。
(1)分析个城镇需求特征。利用excel求解各城镇需求量的平均值、标准差、变异系数并将变异系数降序排序,选取大于10的结果,见表4。
表 4 部分城镇需求量数值分析
城镇 | 平均值 | 标准差 | 变异系数 |
城镇150 | 1661.204833 | 526.3203535 | 31.683049730 |
城镇94 | 107.8893333 | 33.49600421 | 31.046631926 |
城镇51 | 147.7895 | 44.27108723 | 29.955502404 |
城镇40 | 396.6691667 | 113.6814905 | 28.659018662 |
城镇68 | 1433.024167 | 408.6252685 | 28.514890255 |
城镇86 | 866.582 | 233.8058651 | 26.980235580 |
城镇41 | 456.5551667 | 116.8982483 | 25.604408156 |
城镇124 | 922.238 | 236.0944807 | 25.600168361 |
城镇138 | 577.2691667 | 130.7325862 | 22.646729416 |
城镇67 | 187.0373333 | 41.66267894 | 22.275060384 |
城镇47 | 122.8 | 26.09496113 | 21.249968351 |
城镇110 | 1445.038167 | 285.8264396 | 19.779854000 |
城镇2 | 120.9895 | 22.99937991 | 19.009401564 |
城镇83 | 659.1536667 | 120.4067504 | 18.266871059 |
城镇55 | 176.7381667 | 32.13486883 | 18.182189755 |
城镇107 | 967.2781667 | 171.5059282 | 17.730776329 |
城镇108 | 406.9078333 | 71.33112578 | 17.530044875 |
城镇61 | 405.7856667 | 69.3340611 | 17.086375097 |
城镇76 | 1478.902333 | 243.3522839 | 16.454925954 |
城镇109 | 101.6151667 | 16.68467972 | 16.419477785 |
城镇1 | 1084.089333 | 175.7851084 | 16.215002121 |
城镇39 | 358.741 | 56.99400449 | 15.887229085 |
城镇104 | 2141.91 | 339.4398559 | 15.847531218 |
城镇101 | 2097.4935 | 313.1979511 | 14.932010569 |
城镇129 | 99.2745 | 14.40360765 | 14.508869497 |
城镇62 | 185.4725 | 24.99926742 | 13.478692217 |
城镇10 | 688.0568333 | 91.71353976 | 13.329355267 |
城镇69 | 810.145 | 100.1209229 | 12.358395464 |
城镇130 | 290.7853333 | 33.0892288 | 11.379263328 |
城镇141 | 834.4271667 | 89.32731585 | 10.705226222 |
城镇131 | 147.7898333 | 15.69078608 | 10.616959046 |
分析表4,变异系数反映各城镇每个月份蔬菜需求量的离散程序,取变异系数大于10的这32个城镇,因为这32个城镇蔬菜需求量变化较大,波动较大;而其他城镇变异系数较小,各月份蔬菜需求量变化平稳。
(2)蔬菜需求量预测。将每个月各个城镇需求数加和,对加和后的数据运用MATLAB进行二次拟合,程序见附录3。二次拟合后的公式为:
y = - 2.5726 x2 + 373.2062 x + 106296.9867 (4-13)
得到的拟合趋势图如图2。
图 2 全省蔬菜需求量变化趋势
根据图2得,拟合趋势线为一元二次方程,在第60-90个月之间达到峰值,因而对之后30个月的蔬菜需求量进行预测。
根据拟合后的公式预测出2016年12月之后30个月的月度总需求量,如表5。
表 5 60-90月月度总需求量预测
4.2.4 模型检验
4.2.4.1 检验模型建立
4.2.4.1.1 和方差SSE
该统计参数计算的是拟合数据和原始数据对应点的误差的平方和,计算公式如下:
(4-14)
4.2.4.1.2 SSR
预测数据与原始数据均值之差的平方和,公式如下:
(4-15)
4.2.4.1.3 SST
原始数据和均值之差的平方和,公式如下:
(4-16)
4.2.4.1.4 确定系数R-square
(4-17)
4.2.4.2 模型检验
求解得到R-square=0.9621,由于R-square越接近1,说明模型对数据拟合的越好,因此本次建立的模型较好。
4.3问题三模型的建立与求解
4.3.1模型建立
为公司设计增设销售连锁方案,使全省销售能力达到最大。
根据题中已知条件,设f(i)表示新增的所有连锁店的销售能力;Y(j)表示未来j号城镇的需求;z(j)表示在不足10公里的其他城镇的销售连锁店购买,则这一部分需求量只能实现一半;x(j)表示超过10公里的其他城镇的销售连锁店购买,销售量只能达到需求量的三成;s(j)表示表示原有的23家店在新增设连锁店的销售能力。
设增设的连锁店
(6.1)
目标函数
(6.2)
约束条件
(6.3)
编写Lingo程序对模型进行求解,得到的增设销售店序号及各销售量如表7所示。程序见附录7,具体结果见附录8。
表 7 增设销售店及销售量
增设销售店及销售量 | |
城镇 | 销售量/(kg/日) |
6 | 23453 |
8 | 27491 |
10 | 30000 |
18 | 20673 |
31 | 30000 |
33 | 20304.6 |
50 | 21972.8 |
54 | 22079.7 |
56 | 25328.2 |
62 | 23053.5 |
64 | 25866.1 |
68 | 25122 |
76 | 21909.5 |
100 | 21208.5 |
101 | 27736.1 |
104 | 20478 |
110 | 22389 |
116 | 20382 |
120 | 30000 |
121 | 30000 |
123 | 30000 |
125 | 29010 |
150 | 27089.3 |
154 | 22002.2 |
根据表7可得,一共需要增加24个销售连锁店,增加的销售连锁店的销售额为597548.5kg/日。
得到的原销售店增长之后的销售量如表8所示。
表 8 原有销售店及销售量
原有销售店及销售量 | |
城镇 | 销售量/(kg/日) |
1 | 17692.8 |
10 | 5406 |
11 | 7323.6 |
16 | 1761.1 |
22 | 7650 |
24 | 3901.2 |
27 | 11118 |
31 | 25474.1 |
34 | 541.2 |
36 | 11459 |
42 | 5330.2 |
63 | 46489 |
65 | 7035.5 |
79 | 19577.5 |
94 | 5134.6 |
106 | 42360 |
120 | 64147.4 |
123 | 15370 |
141 | 10270.8 |
145 | 3832 |
根据表8可得,原有销售店增长销售额之后的总量为311874kg/日。
根据上述表7和表8的销售连锁店方案,全省销售量达到最大,最大为909422.5kg/日。
本文建立多个模型,对蔬菜配送问题进行了求解,这些模型在求解过程中既有优点也有缺点。
优点:
(1)Floyd算法容易理解,可以算出任意两个节点之间的最短距离,求解最短距离是灵活可用,适合本文问题。
(2)拟合模型较易求解,较好的拟合蔬菜需求量并预测未来最大值。
(3)非线性规划求解较为简单。
缺点:
拟合模型没有对曲线进行三次及以上方程拟合,可能二次方程不是最优解。 只考虑了非常理想的情况,依然存在很多不稳定因素。
本文由职坐标整理并发布,希望对同学们有所帮助。了解更多详情请关注职坐标大数据云计算大数据安全频道!
喜欢 | 0
不喜欢 | 0
您输入的评论内容中包含违禁敏感词
我知道了

请输入正确的手机号码
请输入正确的验证码
您今天的短信下发次数太多了,明天再试试吧!
我们会在第一时间安排职业规划师联系您!
您也可以联系我们的职业规划师咨询:
版权所有 职坐标-一站式AI+学习就业服务平台 沪ICP备13042190号-4
上海海同信息科技有限公司 Copyright ©2015 www.zhizuobiao.com,All Rights Reserved.
沪公网安备 31011502005948号