光明市的菜篮子工程Word下载.docx
- 文档编号:7760771
- 上传时间:2023-05-09
- 格式:DOCX
- 页数:16
- 大小:53.36KB
光明市的菜篮子工程Word下载.docx
《光明市的菜篮子工程Word下载.docx》由会员分享,可在线阅读,更多相关《光明市的菜篮子工程Word下载.docx(16页珍藏版)》请在冰点文库上搜索。
(3)假设各蔬菜市场的蔬菜只来源于A、B、C3个收购站,不包含外来送货状况。
(4)假设该市经济保持相对稳定,3个收购站每年的收购量相对稳定。
(5)假设蔬菜价格一致,不存在恶意竞争。
(6)假设各收购站可作为中转站。
(7)不计算该市的新建市场,只考虑该题固定的8个菜市场。
(8)不考虑该市的新增路段,只在题中路段进行选择。
3.2符号说明
xij:
第i个收购点向j市场供给的数量
cij:
第i个收购点向j市场供给的单位运费
ai:
第i个收购点供应量
bj:
第j个市场需求量
dj:
第j个市场因供给量小于需求量的单位短缺损失
4.问题分析、模型建立及求解
4.1模型(a)的分析及建模
4.1.1模型分析
目标函数总费用Z,包括两项:
蔬菜调运费Q,各市场供给量小于需求量的短缺损失P
Z=P+Q
其中:
P=
Q=
约束条件为
3个收购点的蔬菜全部供给给8个市场
(i=1,2,3)
3个收购点分别向每个市场供应的总量不超过每个市场的需求量
(j=1,…,8)
③变量非负性限制
(i=1,2,3,j=1,…,8)
从而得出问题(a)的数学模型如下:
minZ=
+
s.t.
4.1.2模型求解
(一)为了求解模型,必须求出系数(
)
,其中每一
表示第i个收购点向j市场供给单位量蔬菜的运费,但因为从收购点至各菜市场单位量蔬菜单位路程的调运费用为1元/(100kg*100m),而蔬菜的单位量为100kg,单位距离为100m,则可求出第i个收购点到第j市场每单位蔬菜的单位距离运费为1元/(100m*100kg)*100m*100kg=1元。
因而
在数值上等于第i个收购点到第j市场的距离值,从而等价于一个求最短路的问题,
(1)标志距离:
将图中15个点标号,分别为A,B,C,o,p,q,r,1,2,3,4,5,6,7,8.并由此构成15*15的权矩阵W15*15,其中Wij表示第I个点到第j个点的距离,若第I个点和第j个点不相邻,则wij=
。
(2)列出矩阵:
对得到的W,使用弗洛依德算法,得到最短距离,也就是单位最小费用矩阵。
从中抽取出第i(i=1,2,3)行和第j(j=8,...,15)列的子矩阵W
,其中的值w
ij即对应为第i个收购站到第j个市场的单位最小费用。
表1单位最小运费
1
2
3
4
5
6
7
8
A
19
11
22
20
B
14
16
12
23
17
C
15
10
(3)结合上述分析,根据建立的模型,利用LINGO软件,输入目标函数和约束条件,求解模型的最优解,解如表2。
表2各收购点向市场供应量分配表
75
70
55
60
80
30
90
40
总计费用:
4610(元)
已知各市场每单位短缺损失(表3)
表3各市场每单位短缺损失(元/100kg)
市场
5
6
短缺
损失
(二)比较表1中每个收购点到市场的单位蔬菜的运价cij和表3每个市场的单位蔬菜短缺的损失价格dj,若cij<
dj,即运费小于短缺损失价格,则表明可运,可以减少宏观经济的损失。
若cij=dj,即运费等于短缺损失,对整个宏观经济来说,即可以运也可以不运。
若cij>
dj,即运费大于短缺损失,则不运,否则增加宏观经济的损失。
由此,我们得出
表5
可运
运或不运
而表2中B---3,B----4,A----5,C---8的路线上发生了运输往来,不利于整个宏观经济值增加。
(1)考虑到如C收购点到8市场的单位量蔬菜的运输费用大于8市场单位量蔬菜的短缺损失等情况,模型2修改模型1的
假设,为允许3个收购点分别向每个市场供应的总量可超过每个市场的需求量。
即改变约束条件2,此时模型为
minZ=
s.t.
根据建立的上述模型,利用LINGO软件,输入目标函数和约束条件,求解模型(b)的最优解。
表6各收购点向市场供应量分配表
145
100
4460(元)
比较表2和表6,从第二个模型所求得分配方式中可以看到A—5,C---8两条不合理的运输路线已被取消,同时最终的运费也有所下降,下降了150元。
(2)仍然考虑到如C收购点到8市场的单位量蔬菜的运输费用大于8市场单位量蔬菜的短缺损失等情况,在模型1的基础上,对模型1的
假设做出了另一种修改,为允许每个收购点的蔬菜可以只运部分。
即改变约束条件1,可得模型
(j=1,…,7)
(i=1,2,3,j=1,…,7)
根据建立的模型,利用LINGO软件,输入目标函数和约束条件,求解模型的最优解。
表7各收购点向市场供应量分配表
B
C
3840(元)
(3)比较表2与表7,从第三个模型所求得分配方式中可以看到B—3,B—4,A—5,C---8四条不合理路线都被取消,同时总运费减少了545元。
(三)在市场经济下,模型c,随着市场的调节,最终A只愿供应12000千克,B只愿供应6000千克,C只愿供应10000千克,大大小于各收购点常年的每天收购量20000千克,17000千克,16000千克。
4.2模型(b)的分析及建模
4.2.1模型分析
按题中问题(b)规定各菜市场短缺量一律不超过需求量的20%的条件,我们对需求量的约束条件进行了修改。
s.t.
(j=1,…,8)
(j=1,…,8)
4.2.2模型求解
(1)根据建立的模型,利用LINGO软件,输入目标函数和约束条件,求解模型的最优解。
表8各收购点向市场供应量分配表
50
64
56
24
72
4806(元)
(2)比较表2和表8,主要是对3,4,7,8市场的供应量作出了调整。
其中的主要原因是对于3,4,7,8市场,从收购点到其的单位量蔬菜的运输费用大于该市场单位量蔬菜的短缺损失,所以,当加入各菜市场短缺量一律不超过需求量的20%的约束条件后,为了保证4,8市场的需求,在考虑到3,7市场相对其他市场运输代价较高的情况下,在这四个市场之间做出平衡供给量的调整。
4.3模型(c)的分析及建模
4.3.1模型分析
为满足城市居民的蔬菜供应,光明市的领导规划增加蔬菜种植面积,即模型为
s.t.
t
0(i=1,2,3)
4.3.2模型求解
(1)根据建立的模型,利用LINGO软件,输入目标函数和约束条件,求解模型的最优解。
表9各收购点向市场供应量分配表
4770(元)
各收购点增加的蔬菜收购量如下表:
t
(2)比较表2和表9,对于4,8市场做了比较大的调整,主要是因为从收购点到其的单位量蔬菜的运输费用大大高于这两个市场单位量蔬菜的短缺损失,供应这两个市场并不能获得收益,反而会受到较大损失。
在市场机制的主导下,无法满足这两个市场的需求量。
但是,光明市的领导为保证城市居民的蔬菜供应,规划增加蔬菜种植面积,提高收购点的蔬菜收购量,在承担一定损失的情况下,满足了这两个市场的蔬菜需求。
5模型评价
(1)以上各模型的优点:
所建立的模型简洁明了,便于使用数学工具。
如Lingo,降低了编程求解的难度,缩短了运行时间,提高了工作效率。
②对同一个问题从不同的角度进行了考虑,建立了多个模型,并进行了结果的比较分析,既结合题目要求,又考虑了实际意义。
从社会效益和经济效益对问题进行了分析,也表现出现实生活中政府在寻求两者之间的平衡中做出的努力。
(2)不足之处:
以上模型均只考虑在降低运输费用和短缺费用的目标下的优化方案,并未涉及到市场上蔬菜供过于求和收购点蔬菜积压而导致的存储费用等,而使所建立的模型不能很好地符合实际情况,还有待改进。
参考文献
[1]姜启源,谢金星,叶俊.数学建模[M],北京:
高等教育出版社,2004.
[2]谢金星,薛毅.优化建模与LINDO/LINGO软件[M],北京:
清华大学出版社,2005.7.
[3]胡运权.运筹学基础及应用[M],北京:
高等教育出版社,2008.6.
附录部分
1.数据预处理部分:
求最小费用的LINGO的文件如下:
model:
SETS:
NODES/A,B,C,o,p,q,r,1,2,3,4,5,6,7,8/;
ROADS(NODES,NODES)/A,oA,pA,1A,2A,3A,6B,oB,rB,2B,3C,qC,5C,7C,8o,2o,3p,3p,5p,6q,5q,6q,7r,3r,4r,5r,81,21,63,57,8/:
W0;
LINK(NODES,NODES):
W,D;
NNN(Nodes,nodes,nodes):
U;
ENDSETS
DATA:
BIG=1000;
W0=7448866117786510354756710653675511;
@TEXT(FinalCost.txt)=@writefor(nodes(i)|i#le#3:
@writefor(nodes(j)|j#ge#8#and#j#le#15:
@format(D(i,j),'
5.0f'
)));
ENDDATA
CALC:
@FOR(LINK(i,j)|@IN(ROADS,i,j):
W(i,j)=W0(i,j);
W(j,i)=W0(i,j);
);
@FOR(LINK(i,j)|i#eq#j:
W(i,j)=0);
@FOR(LINK(i,j)|i#ne#j#and##not#@IN(ROADS,i,j)#and#
#not#@IN(ROADS,j,i):
W(i,j)=BIG;
W(j,i)=BIG;
@FOR(NNN(i,j,k)|k#eq#1:
U(i,j,k)=W(i,j));
@For(nodes(k)|k#lt#@size(nodes):
@FOR(LINK(i,j):
U(i,j,k+1)=
@if(U(i,j,k)#le#U(i,k,k)+U(k,j,k),
U(i,j,k),U(i,k,k)+U(k,j,k))));
@FOR(NNN(i,j,k)|k#eq#@size(nodes):
D(i,j)=
U(i,j,k),U(i,k,k)+U(k,j,k)));
ENDCALC
End
求得最短路部分结果如下:
VariableValue
D(A,B)13.00000
D(A,C)17.00000
D(A,O)7.000000
D(A,P)4.000000
D(A,Q)13.00000
D(A,R)14.00000
D(A,1)4.000000
D(A,2)8.000000
D(A,3)8.000000
D(A,4)19.00000
D(A,5)11.00000
D(A,6)6.000000
D(A,7)22.00000
D(A,8)20.00000
2.模型主体LINGO程序如下(只取其中一个,其他类同):
MODEL:
SUPPLY/A,B,C/:
S;
NEED/1..8/:
B,P;
LINK(Supply,need):
C,X;
ENDSETS
S=200170160;
B=75608070100559080;
P=10851010858;
C=@TEXT(FINALCOST.txt);
@TEXT(FinalResult.txt)=@writefor(supply(i):
@writefor(need(j):
@format(x(i,j),'
)),@newline
(1));
[obj]MIN=@sum(link(i,j):
c(i,j)*x(i,j))
+@sum(need(j):
P(j)*(b(j)-@sum(supply(i):
x(i,j))));
@for(supply(i):
[con1]@sum(need(j):
x(i,j))=S(i));
@for(need(j):
[con2]@sum(supply(i):
x(i,j))<
=b(j));
END
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 光明 菜篮子工程