农产品问题数学建模.docx
《农产品问题数学建模.docx》由会员分享,可在线阅读,更多相关《农产品问题数学建模.docx(16页珍藏版)》请在冰豆网上搜索。
农产品问题数学建模
摘要
本文针对某国政府奶产品定价和销售总收入的问题,运用最优化理论建立了数学模型,通过LING求解,能迅速定出各产品的价格,具有较强的实用价值。
对问题一,由于受原料数量约束,首先建立目标函数.其次,通过综合分析建立产品价格的优化模型;最后,我们用Lingo软件进行求解。
得出在销售总收入最大的情况下4种产品的价格及所导致的需求,如下表:
牛奶
奶油
奶酪1
奶酪2
销量(万吨)
334.4686
48.36238
20.11823
3.559875
价格(元/吨)
524.2657
583.6468
1223.093
1950.675
销售收入(万元)
175350.4
28226.55
24606.47
6944.160
销售总收入(万元)
235127.6
问题二,由于给出的政策约束,在问题一建立产品价格的优化模型的基础上,分别结合了问题二中的每一种价格指标限制的限定条件,运用Lingo软件可以迅速得出这一政策限制情况下的最大销售总收入,并得出经济代价给出的数量,如下表:
牛奶价格限制
奶油价格限制
奶酪1价格限制
奶酪2价格限制
限制下的销售额(万元)
201159.3
229985.5
231050.6
231050.6
经济代价(万元)
33968.3
5142.1
4077
4077
主要特点在于,所用求解模型的效率十分显著。
在对原始数据做简单预处理的条件下,通过LINGO的求解可以直接求得相应的最优解。
另外,本文所建立的模型简单,LINGO求解,方便模型求解实现,对模型推广有很大帮助。
关键字:
数学建模最优解LINGO原料限制政策限制
一、问题的提出
某国政府要为其牛奶、奶油和奶酪等奶制品定价。
所有这些产品都直接或间接的来自国家的原奶生产。
原奶首先要分离成脂肪和奶粉两种组合,去掉生产出口产品和农场消费的产品后,余下的共有60万吨脂肪和70万吨奶粉,可用于生产牛奶、奶油和两种奶酪,供国内全年消费。
各种产品的百分比组成见下表:
往年的国内消费和价格如下表:
产品\成分
脂肪
奶粉
水
牛奶
4
9
87
奶油
80
2
18
奶酪1
35
30
35
奶酪2
25
40
35
产品
牛奶
奶油
奶酪1
奶酪2
消费量(千吨)
4820
320
210
70
价格(元/吨)
297
720
1050
815
价格的变化会影响消费需求。
为表现这方面的规律,定义需求的价格伸缩性E:
E=需求降低百分数/价格提高百分数;
各种产品的E值,可以据往年的价格和需求变化情况的统计数据,用数理统计方法求出。
另外,两种奶酪的需求,随它们价格的相对变化,在某种程度上可以相互替代。
表现这一规律要用需求关于价格的交叉伸缩性EAB其定义为:
EAB=A需求提高百分数/B价格提高百分数。
奶酪1到奶酪2的E12值和奶酪2到奶酪1的交叉伸缩性E21值,同样可以凭数据用统计方法求出。
现已经求出牛奶、奶油、奶酪1、奶酪2的E值依次为0.4,2.7,1.1和0.4以及E12=0.1,E21=0.4。
问题一:
试求出4种产品的价格及所导致的需求,使销售总收入为最大。
问题二:
如果政策不允许某种价格指标上升,这使得新的价格必须使消费的总费用较上一年度不增加,这是对这个问题的一个特别重要的附加要求。
请对这一政策限制的经济代价给出数量表示。
二、问题的分析
a、题中要求制定一个总费用(决策目标)最大的生产定价方案,属于优化问题,并且对于问题二新的价格必须使消费的总合和单价用较上一年度不增加,产量还受原料的限制(此两点为主要约束条件)
b、故此模型即为基于以上约束条件的整数规划(最优决策目标)问题。
三、模型假设
a、工厂在生产期间不会出现因设备、人员发生意外而导致产量下降。
b、水这种原料可以无限提供使用。
c、产品卫生质量达标。
d、数理统计数据接近实际消费变化。
e、产品正常销售。
四、符号说明
i:
第i种产品
Zi:
第i种产品的销售收入
Pi:
第i种产品的价格
P0i:
第i种产品去年的价格
Si:
第i种产品的销量
S0i:
第i种产品的去年的销量
vii:
第i种产品自身对自身的影响下的销量的变化率(正负号表示提高和降低)
vij:
第i种产品在j的影响下的销量的变化率(正负号表示提高和降低)
Ei:
第i种产品的价格伸缩率
Eij:
第i产品对第j种产品的交叉伸缩率
Yj:
第j种原料的总量(万吨)
yij:
第i种产品j种原料的占有率
五、模型建立
由题意知道,四种产品由三种原料加工,原料一(脂肪)、原料二(奶粉)分别为60、70万吨,原料三(水)是无限提供,所以用于生产各种产品的原料之和不能超过这种原料。
即:
Siyij
Yj
价格的变化会影响消费需求,根据价格伸缩率的定义,因此对于第i种产品的销售量有以下关系:
产品三和产品四的销售量,除了受到自身的价格影响,还要受到除自身外另外一种产品的影响。
对产品三的销量变化率:
市场对产品三和产品四的需求,随着他们价格相对变化,在某种程度上可以相互替代(属于替代品),根据交叉伸缩率的定义,产品四的价格和产品三销售量变化率:
由以上两式得,得到产品三的销售总量和产品三的价格及产品四的价格的关系为:
产品三和产品四的销售量,除了受到自身的价格影响,还要受到除自身外另外一种产品的影响。
对产品四
市场对产品三和产品四的需求,随着他们价格相对变化,在某种程度上可以相互替代(属于替代品),根据交叉伸缩率的定义,产品四销售量与产品三的价格关系如下:
由以上两式得,得到产品三的销售总量和产品三的价格及产品四的价格的关系为:
问题一:
目标函数:
目标函数:
Max=z1+z2+z3+z4zi=pisi
约束条件:
(i=1,2)
Siyij
Yj
Si
0
Pi
0
六、模型求解
将上述模型输入LINGO可得到
Globaloptimalsolutionfound.
Objectivevalue:
235127.6
Objectivebound:
235127.7
Infeasibilities:
0.2910383E-10
Extendedsolversteps:
93
Totalsolveriterations:
19646
VariableValueReducedCost
E10.40000000.000000
E22.7000000.000000
E31.1000000.000000
E40.40000000.000000
E340.10000000.000000
E430.40000000.000000
S01482.00000.000000
S0232.000000.000000
S0321.000000.000000
S047.0000000.000000
P01297.00000.000000
P02720.00000.000000
P031050.0000.000000
P04815.00000.000000
Z1175350.40.000000
Z228226.550.000000
Z324606.470.000000
Z46944.1600.000000
S1334.46860.000000
P1524.26570.000000
S248.362380.000000
P2583.64680.000000
S320.118230.000000
P41950.6750.000000
P31223.0930.000000
S43.5598750.000000
RowSlackorSurplusDualPrice
1235127.61.000000
20.000000515.2343
30.0000001.000000
40.000000403.0198
50.0000001.000000
60.0000001144.069
70.0000001.000000
80.0000001894.229
90.0000001.000000
100.000000225.7837
1131.471160.000000
12334.46860.000000
1348.362380.000000
1420.118230.000000
153.5598750.000000
16524.26570.000000
17583.64680.000000
181223.0930.000000
191950.6750.000000
牛奶
奶油
奶酪1
奶酪2
销量(万吨)
334.4686
48.36238
20.11823
3.559875
价格(元/吨)
524.2657
583.6468
1223.093
1950.675
销售收入(万元)
175350.4
28226.55
24606.47
6944.160
销售总收入(万元)
235127.6
即最大销售利润为Z=235127.6(万元);仅当定牛奶、奶油、奶酪1、奶酪2的价格分别为3524.2657(元)、583.6468(元)、1223.093(元)、1950.675(元)时,它们的销量分别为334.4686(万吨)、48.36238(万吨)、20.11823(万吨)、3.559875(万吨)时,总销售利润达最大
问题二:
由题意可得,不允许某种价格指标上升,这使得新的价格必须使消费的总费用较上一年度不增加,计算该政策的经济代价。
分析可得这个政策有4种情况
即:
(1)以牛奶作为政策对象
(2)以奶油作为政策对象
(3)以奶酪1作为政策对象
(4)以奶酪2作为政策对象
以下进行逐一分析:
(1)以牛奶作为政策对象
目标函数:
Max=z1+z2+z3+z4zi=pisi
约束条件:
(i=1,2)
p1
p01
p1s1
s01p01
si
0
pi
0
siyij
Yj
利用LONGO进行分析计算可得下表
牛奶
奶油
奶酪1
奶酪2
销量(万吨)
482.0000
41.21395
19.59366
3.564245
价格(元/吨)
297.0000
643.2171
1249.158
1969.634
销售收入(万元)
143154.0
26509.52
24475.57
7020.260
销售总收入(万元)
201159.3
经济代价(万元)
201159.3-235127.6=-33968.3
即是,当牛奶价格指标受到限制时,销售总收入降为201159.3万元,经济代价数量为:
-33968.3万元,即销售总收入减少了33968.3万元。
(2)以奶油作为政策对象
目标函数:
Max=z1+z2+z3+z4zi=pisi
约束条件:
(i=1,2)
p2
p02
p2s2
s02p02
si
0
pi
0
siyij
Yj
利用LONGO进行分析计算可得下表
牛奶
奶油
奶酪1
奶酪2
销量(万吨)
337.4000
32.00000
20.91353
3.553249
价格(元/吨)
519.7500
720.0000
1183.577
1921.932
销售收入(万元)
175363.6
23040.00
24752.77
6829.101
销售总收入(万元)
229985.5
经济代价(万元)
229985.5-235127.6=-5142.1
即是,当牛奶价格指标受到限制时,销售总收入降为229985.5万元,经济代价数量为:
-5142.1万元,即销售总收入减少了5142.1万元。
(3)以奶酪1作为政策对象
目标函数:
Max=z1+z2+z3+z4zi=pisi
约束条件:
(i=1,2)
p3
p03
p3s3
s03p03
si
0
pi
0
siyij
Yj
利用LONGO进行分析计算可得下表
牛奶
奶油
奶酪1
奶酪2
销量(万吨)
334.0788
46.92106
21.00000
7.000000
价格(元/吨)
524.8662
595.6578
815.0000
1050.000
销售收入(万元)
175346.7
27948.90
22050.00
5705.000
销售总收入(万元)
231050.6
经济代价(万元)
231050.6-235127.6=-4077
即是,当牛奶价格指标受到限制时,销售总收入降为231050.6万元,经济代价数量为:
-4077万元,即销售总收入减少了4077万元。
(4)以奶酪2作为政策对象
目标函数:
Max=z1+z2+z3+z4zi=pisi
约束条件:
(i=1,2)
p4
p04
p4s4
s04p04
si
0
pi
0
siyij
Yj
利用LONGO进行分析计算可得下表
牛奶
奶油
奶酪1
奶酪2
销量(万吨)
334.0788
46.92106
21.00000
7.000000
价格(元/吨)
524.8662
595.6578
815.0000
1050.000
销售收入(万元)
175346.7
27948.90
22050.00
5705.000
销售总收入(万元)
231050.6
经济代价(万元)
231050.6-235127.6=-4077
即是,当牛奶价格指标受到限制时,销售总收入降为231050.6万元,经济代价数量为:
-4077万元,即销售总收入减少了4077万元。
七、结果分析
将以上讨论的4个问题制成图表:
牛奶价格限制
奶油价格限制
奶酪1价格限制
奶酪2价格限制
限制下的销售额(万元)
201159.3
229985.5
231050.6
231050.6
经济代价(万元)
33968.3
5142.1
4077
4077
画出经济限制关系条形图:
图一:
销售总收入在各个产品价格限制下与无限制情况下的对比
图二:
经济代价在各种产品政策限制下的数量表示
由上表可知,当牛奶价格指标受到限制时,产生的经济代价最大,而其他产品价格指标的限制影响较小,其中限制影响程度(牛奶>奶油>奶酪1=奶酪2)。
八、方案评价
(1)本文把所解决的问题归结为优化问题,建立的数学模型清晰合理。
(2)运用MATLAB和LINGO软件处理数据和进行运算,降低运算量,简单易行,有很大的可操作性,且所得数据较为合理可靠。
(3)、“需求的价格伸缩性”和“价格的交叉伸缩性”是由数理统计方法算出,与真实数据存在一定误差。
(4)、模型与实际紧密联系,符合价格的变化规律,本文最后还对模型进行了推广,使得模型更贴近实际,通用性强。
(5)、对于模型的推广,文中的变量以及约束条件可以扩大,例如:
当有N种产品以及有M种原料约束条件等,通过所建立的模型就能够快速、有效的解决。
九、参考资料
[l]陈东彦主编.数学建模.科学出版社,2008
[2]化存才主编.数学建模应用与实践.云南科技出版社,2008.10.
[3]姜启源,数学模型(第3版),北京:
高等教育出版社,1999.
[4]李伯德,数学建模方法,北京:
高等教育出版社,2006年05月第1版
[5]沈继红施久玉高振滨张晓威编著.数学建模.哈尔滨工程大学出版社,2007.
[6]美)MARKM.MEERSCHAERT著.数学建模方法与分析原书第3版.机械工业出版社,2009.05.
[7]袁新生,邵大宏,郁时炼主编.LINGO和EXCEL在数学建模中的应用.科学出版社,2006.
[8]肖华勇编著.基于MATLAB和LINGO的数学实验.西北工业大学出版社,2009.03.
十、附录
一、部分LONGO代码:
!
数据;
data:
E1=0.4000;E2=2.7000;E3=1.1000;E4=0.4;E34=0.1000;E43=0.4000;
s01=482;s02=32;s03=21;s04=7;
p01=297;p02=720;p03=1050;p04=815;
enddata
!
目标函数;
max=z1+z2+z3+z4;
!
约束;
s1=s01-(p1-p01)*E1*s01/p01;
z1=s1*p1;
s2=s02-(p2-p02)*E2*s02/p02;
z2=s2*p2;
s3=s03+(p4-p04)*E34*s03/p04-(p3-p03)*E3*s03/p03;
z3=s3*p3;
s4=s04+(p3-p03)*E43*s04/p03-(p4-p04)*E4*s04/p04;
z4=s4*p4;
!
约束条件;
0.04*s1+0.80*s2+0.35*s3+0.25*s4<=60;
0.09*s1+0.02*s2+0.30*s3+0.40*s4<=70;
s1>=0;s2>=0;s3>=0;s4>=0;
p1>=0;p2>=0;p3>=0;p4>=0;
end