【云计算】蔬菜配送方案制定与实施
小标 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


 

 

σ


 

 

各城镇蔬菜需求量标准差


 


四、模型的建立与求解


4.1问题一模型的建立与求解


本问题要求设计运输成本最低的生产与配送方案,建立Floyd模型画出全省交通网络拓扑图。首先求解出所有城镇之间的最短路径,再建立0-1矩阵,选取距离销售连锁店最近的生产基地,即得到各个生产基地的配送方案,将距离加和,之后求解出运输成本即为最低运输成本。


4.1.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


 


4.2 问题二模型的建立与求解


本问题要求分析各城镇需求特征并预测全省鲜蔬菜何时达到峰值,以及达到峰值时排序前5位和后5位的城镇,在分析特征时,拟求解各城镇蔬菜需求量的平均值、标准差,由于标准差不能完全反映需求特征,因而再求解变异系数,分析这三个数据体现各城镇需求量特征。拟建立拟合模型,对全省各月总需求量进行拟合,求解出何时达到峰值,再对每一个城镇需求量进行拟合,求解出在全省需求量达到峰值的月份各城镇蔬菜需求量,之后对各城镇需求量排序,得到排名前五和后五的城镇。


4.2.1需求特征


通过均值和方差以及变异系数的计算,可以看出每个镇对蔬菜的需求量及其波动情况。


 (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)


这时拟合曲线,也即为拟合多项式。


预测出拟合多项式的峰值,就是全省鲜蔬菜需求量的峰值。


 


4.2.3模型求解


(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)


4.3.2模型求解


编写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
看完这篇文章有何感觉?已经有0人表态,0%的人喜欢 快给朋友分享吧~
评论(0)
后参与评论

您输入的评论内容中包含违禁敏感词

我知道了

助您圆梦职场 匹配合适岗位
验证码手机号,获得海同独家IT培训资料
选择就业方向:
人工智能物联网
大数据开发/分析
人工智能Python
Java全栈开发
WEB前端+H5

请输入正确的手机号码

请输入正确的验证码

获取验证码

您今天的短信下发次数太多了,明天再试试吧!

提交

我们会在第一时间安排职业规划师联系您!

您也可以联系我们的职业规划师咨询:

小职老师的微信号:z_zhizuobiao
小职老师的微信号:z_zhizuobiao

版权所有 职坐标-一站式AI+学习就业服务平台 沪ICP备13042190号-4
上海海同信息科技有限公司 Copyright ©2015 www.zhizuobiao.com,All Rights Reserved.
 沪公网安备 31011502005948号    

©2015 www.zhizuobiao.com All Rights Reserved