牛奶价格非线性规划.docx

1、牛奶价格 非线性规划目录目录 1摘要 4关键词：牛奶价格 非线性规划 交叉伸缩性 4一、问题的重述与分析 4二、模型的假设 5三、符号说明 5四、模型的建立与求解 64.1不同奶制品价格的分析： 64.1.1 奶制品价格伸缩性 64.1.2 奶制品价格伸缩性 64.1.3奶制品价格的求解： 64.2模型的建立 64.2.1各种奶制品的脂肪含量 64.2.2各种奶制品的奶粉含量 64.3脂肪含量和奶粉含量满足的约束条件 64.3.1消费量满足的约束条件 74.3.2模型建立 74.4模型求解 84.4.1 奶制品价格伸缩性 84.4.2 奶制品价格伸缩性 84.4.3奶制品价格的求解： 84.4

2、.5建立目标函数 84.4.6脂肪和奶粉的约束条件 8五模型的评价 14六附录 14七参考文献 17摘要 奶牛养殖者与乳制品企业间稳定牛奶价格的存在按照价格规律，商品价格应围绕价值上下波动；按照供求关系影响价格的经济规律，当供给大于需求时，商品价格下降，反之，当供给小于需求时，商品价格上涨。因此，几乎所有商品的市场价格都是变化波动的，例如鸡蛋、猪肉价格的变动等等。而在中国奶牛养殖业中，奶源提供者（奶牛养殖者）与乳制品企业之间的牛奶价格却相对稳定。在变化不定的市场条件下，却有“稳定牛奶价格”，这一现象既违背价格规律（按照价格规律，当饲料价格上涨，牛奶单位价值增加，牛奶的价格也应上涨；反之，当饲料

3、价格下降，牛奶单位价值减少，牛奶的价格则应下降）；又违背供求关系影响价格的规律（按照供求关系影响价格的规律，当牛奶的供给大于需求时，牛奶价格下降，反之，当牛奶的供给小于需求时，牛奶价格上涨那么。某国政府要为其牛奶、奶油和奶酪等奶制品定价。所有这些产品都直接或间接的来自国家的原奶生产，原奶首先要分离成脂肪和奶粉两种组合，去掉生产出口产品和农场消费的产品的部分后，余下的共有60万吨脂肪和70万吨奶粉，可用于生产牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费。需制定合理的价格来促进价格的增长。针对问题1：牛奶、奶油、奶酪1、奶酪2的E值依次为0.4，2.7，1.1和0.4及E12=0.1, E21=

4、0.4试求出4种产品的价格，使所导致的需求使销售总收入为最大。本文建立了一个非线性规划的模型，利用matlab软件求解。针对问题2：然而政策不允许某种产品价格上升过快，假定各种奶制品的价格伸缩性 E满足,交叉伸缩性,试求出4种产品的价格，使所导致的需求使销售总收入为最大。因此，作为对问题的一个特别重要的附加要求，是对这一政策限制的经济代价，给出数量表示。我们在第一问的基础上添加了新的约束条件，让其每种牛奶、奶油和两种奶酪1和奶酪2的总的消费总额不增加。关键词：牛奶价格 非线性规划 交叉伸缩性 一、问题的重述与分析 某国政府要为其牛奶、奶油和奶酪等奶制品定价。所有这些产品都直接或间接的来自国家的

5、原奶生产，原奶首先要分离成脂肪和奶粉两种组合，去掉生产出口产品和农场消费的产品的部分后，余下的共有60万吨脂肪和70万吨奶粉，可用于生产牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费。 表5-17 各种奶制品的百分比组成见下表：产品/成分脂肪奶粉水牛奶4987奶油80218奶酪1353035奶酪2254035 表5-18 往年的国内消费量和价格如下表：产品牛奶奶油1奶酪1奶酪2消费量（千吨）482032021070价格（元/吨）2977201050815价格的变化会影响消费需求。为表现这方面的规律，定义需求的价格伸缩性 E：E=需求降低百分数/价格提高百分数，各种产品的E值，可依据往年的价格和

6、需求变化情况的统计数据，用数理统计方法求出。另外，两种奶酪的需求，随它们价格的相对变化，在某种程度上可以相互替代，这一规律可以用需求关于价格的交叉伸缩性EAB定义作：EAB=A需求提高百分数/B价格提高百分数。奶酪1到奶酪2的交叉伸缩性E12、奶酪2到奶酪1的交叉伸缩性E21、同样可依据往年的数据用统计方法求出。请解决下列问题：1.牛奶、奶油、奶酪1、奶酪2的E值依次为0.4，2.7，1.1和0.4及E12=0.1, E21=0.4试求出4种产品的价格，使所导致的需求使销售总收入为最大。2.然而政策不允许某种产品价格上升过快，假定各种奶制品的价格伸缩性 E满足,交叉伸缩性,试求出4种产品的价格

7、，使所导致的需求使销售总收入为最大。因此，作为对问题的一个特别重要的附加要求，是对这一政策限制的经济代价，给出数量表示。二、模型的假设1.提供的数据真实有效；2.4种产品都能销售出去；3.利用非线性规划有最优解。三、符号说明表示奶制品价格的提高百分比；表示原先奶制品价格；表示现在奶制品价格；表示原先的奶制品消费量；表示现在奶制品消费量；表示脂肪在各种奶制品的比例；表示奶粉在各种奶制品的比例；表示所有奶制品的销售额；其余未说明的符号会在文章中进行说明四、模型的建立与求解4.1不同奶制品价格的分析：4.1.1 奶制品价格伸缩性 价格的变化会影响消费需求。为表现这方面的规律，定义需求的价格伸缩性 E

8、：E=需求降低百分数/价格提高百分数，各种产品的E值，可依据往年的价格和需求变化情况的统计数据，用数理统计方法求出。牛奶、奶油、奶酪1、奶酪2的E值依次为0.4，2.7，1.1和0.4。可以得到，即牛奶、奶油的销量为： 4.1.2 奶制品价格伸缩性 另外，两种奶酪的需求，随它们价格的相对变化，在某种程度上可以相互替代，这一规律可以用需求关于价格的交叉伸缩性EAB定义作：EAB=A需求提高百分数/B价格提高百分数。奶酪1到奶酪2的交叉伸缩性E12、奶酪2到奶酪1的交叉伸缩性E21、同样可依据往年的数据用统计方法求出。奶酪1、奶酪2的E值依次为E12=0.1, E21=0.4。可以得到，即奶酪1、

9、奶酪2的销量为： 4.1.3奶制品价格的求解： 由表示奶制品价格的提高百分比；可以得到现在各种奶制品价格分别为：牛奶、奶油、奶酪1、奶酪2价格分别为：4.2模型的建立4.2.1各种奶制品的脂肪含量牛奶的脂肪含量： 奶油的脂肪含量： 奶酪1的脂肪含量： 奶酪2的脂肪含量： 4.2.2各种奶制品的奶粉含量牛奶的奶粉含量： 奶油的奶粉含量： 奶酪1的奶粉含量： 奶酪2的奶粉含量： 4.3脂肪含量和奶粉含量满足的约束条件 所有这些产品都直接或间接的来自国家的原奶生产，原奶首先要分离成脂肪和奶粉两种组合，去掉生产出口产品和农场消费的产品的部分后，余下的共有60万吨脂肪和70万吨奶粉，可用于生产牛奶、奶油

10、和两种奶酪1和奶酪2，供国内全年消费。得到以下约束条件：4.3.1消费量满足的约束条件牛奶的消费量满足的约束条件：奶油的消费量满足的约束条件：奶酪1的消费量满足的约束条件：奶酪2的消费量满足的约束条件：4.3.2模型建立 所有这些产品都直接或间接的来自国家的原奶生产，原奶首先要分离成脂肪和奶粉两种组合，去掉生产出口产品和农场消费的产品的部分后，余下的共有60万吨脂肪和70万吨奶粉，可用于生产牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费。价格的变化会影响消费需求。为表现这方面的规律，定义需求的价格伸缩性 E：E=需求降低百分数/价格提高百分数，各种产品的E值，可依据往年的价格和需求变化情况的统

11、计数据，用数理统计方法求出。另外，两种奶酪的需求，随它们价格的相对变化，在某种程度上可以相互替代，这一规律可以用需求关于价格的交叉伸缩性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种产品的价格，使所导致的需求使销售总收入为最大。由上面的分析后得到要求解销售总收入满足如下的函数关系： 销售总收入的目标函数是：；约束条件是： 4.4模型求解4.4.1 奶制品价格伸缩性 价

12、格的变化会影响消费需求。为表现这方面的规律，定义需求的价格伸缩性 E：E=需求降低百分数/价格提高百分数，各种产品的E值，可依据往年的价格和需求变化情况的统计数据，用数理统计方法求出。牛奶、奶油、奶酪1、奶酪2的E值依次为0.4，2.7，1.1和0.4。可以得到，即牛奶、奶油的销量为： 4.4.2 奶制品价格伸缩性 另外，两种奶酪的需求，随它们价格的相对变化，在某种程度上可以相互替代，这一规律可以用需求关于价格的交叉伸缩性EAB定义作：EAB=A需求提高百分数/B价格提高百分数。奶酪1到奶酪2的交叉伸缩性E12、奶酪2到奶酪1的交叉伸缩性E21、同样可依据往年的数据用统计方法求出。奶酪1、奶酪

13、2的E值依次为E12=0.1, E21=0.4。可以得到，即奶酪1、奶酪2的销量为： 4.4.3奶制品价格的求解：由表示奶制品价格的提高百分比；可以得到现在各种奶制品价格分别为：牛奶、奶油、奶酪1、奶酪2价格分别为：4.4.5建立目标函数由上面的分析后得到要求解销售总收入满足如下的函数关系：销售总收入的目标函数是： 4.4.6脂肪和奶粉的约束条件由于所有这些产品都直接或间接的来自国家的原奶生产，原奶首先要分离成脂肪和奶粉两种组合，去掉生产出口产品和农场消费的产品的部分后，余下的共有60万吨脂肪和70万吨奶粉，可用于生产牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费。即脂肪和奶粉的满足以下约束条

14、件：可用于生产牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费。其牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费，全年消费不能为负。得到以下不等式：即需要求解的模型是： 所有这些产品都直接或间接的来自国家的原奶生产，原奶首先要分离成脂肪和奶粉两种组合，去掉生产出口产品和农场消费的产品的部分后，余下的共有60万吨脂肪和70万吨奶粉，可用于生产牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费。价格的变化会影响消费需求。为表现这方面的规律，定义需求的价格伸缩性 E：E=需求降低百分数/价格提高百分数，各种产品的E值，可依据往年的价格和需求变化情况的统计数据，用数理统计方法求出。另外，两种奶酪的需求，随它

15、们价格的相对变化，在某种程度上可以相互替代，这一规律可以用需求关于价格的交叉伸缩性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种产品的价格，使所导致的需求使销售总收入为最大。由上面的分析后得到要求解销售总收入满足如下的函数关系：销售总收入的目标函数是：函数满足的约束条件是：在约束条件当中：第一与第二式子，所代表的为原材料的限制，分别为：第一个式子所代表的为“脂肪资源”是

16、有限的，上限为600千吨；第二个式子所代表的为“奶粉资源”是有限的，上限为700千吨；第三、四、五个式子所代表的为牛奶、奶油和两种奶酪1和奶酪2，供国内全年消费不能为负值。模型求解见下文。对于上面的不等式约束条件即目标函数我们利用matlab中的fmincon函数求解程序如下：%求解过程目标函数：function f=fun(x);f=-(297*(1+x(1)*4820*(1-0.4*x(1)+720*(1+x(2)*320*(1-2.7*x(2)+1050*(1+x(3)*210*(1-1.1*x(3)+0.1*x(4)+815*(1+x(4)*70*(1-0.4*x(4)+0.4*x(3

17、)%线性约束：x0=1;0.1;0.01;0.02; A=-1928/25 -3456/5 -1477/20 7/20;-4338/25 -432/25 -581/10 -49/10;1 0 0 0;0 1 0 0;0 0 1.1 -0.1;0 0 -0.4 0.4;b=301/5;844/5;2.5;1/2.7;1;1; Aeq=;beq=; VLB=; VUB=2.5;1/2.7; ; ; x,fval=fmincon(fun,x0,A,b,Aeq,beq,VLB,VUB); xmin=x fmin=fval%得到的最大利润：fmin = -2.3513e+006 %X的取值为：xmin

18、= xmin = 0.7652 -0.1894 0.1649 1.3935可以得到：价格的提高百分率的最优解为，即在保障销售总收入最大的前提下，牛奶的最优价格为524元/吨，奶油的最优价格为483元/吨，奶酪1的最优价格为 201元/吨，奶酪2的最优价格为413元/吨。可获得最高销售总收入为元。此时的奶制品价格（元/吨）：牛奶的价格为： 524奶油的价格为： 583奶酪1的价格为： 1223奶酪2的价格为： 1952奶制品销量（千吨）： 牛奶的销量为： 3344奶油的销量为： 483奶酪1的销量为： 201奶酪2的销量为： 413得到的最大利润为：牛奶奶油奶酪1奶酪2总和价格提高百分率%76.

19、52-18.9416.49139.35销量减少的百分率%30.61-51.1418.1455.82奶制品价格（元/吨）5245831223 1952奶制品销量（千吨）3344483 201413最大利润（千元）1753500282270 2461208078902351300由表格可以发现奶酪2价格最高，牛奶的价格最低。牛奶销售量最大，奶酪1销售量最小。发现了价格与销售量之间存在一定程度上的反相关。从表格中难以直观看出各产品的消费量，为观察各奶制品之间消费量的比例，绘制扇形图。牛奶奶油奶酪1奶酪2总和奶制品原先销量（千吨）4820320210815奶制品销量（千吨）3344483 201413

20、得到奶制品原先销量扇形图（千吨）：从扇形图中可以发现，在奶制品中，牛奶所占的比例最大占72%，奶酪1所占比例最小占3%，且奶酪2比奶酪1的销售量多。得到现在奶制品销量扇形图（千吨）：从扇形图中可以发现，在奶制品中，牛奶所占的比例最大占75%，奶酪1所占比例最小占5%，且奶酪2比奶酪1的销售量多。且可以发现当销售额达到最大时，牛奶、奶酪1、奶酪2的销售较上一年高，而奶油销售量较上一年低。牛奶奶油奶酪1奶酪2总和奶制品原先价格（元/吨）2977201050815奶制品价格（元/吨）5245831223 1952奶制品蓝色为原先价格（元/吨）的柱形图奶制品红色为现在价格（元/吨）的柱形图：A=297

21、 524;720 583;1050 1223;815 1952; h=bar(A)从柱状图中较为直观地可以发现，奶油的价格略微降低，牛奶、奶酪1、奶酪2的价格均比往年有小幅上涨。分析往年的销售量与现在的销售量发现，奶油类地销售量均下降，奶酪1、奶酪2销售量有所升高。牛奶奶油奶酪1奶酪2总和奶制品原先价格（元/吨）2977201050815奶制品价格（元/吨）5245831223 1952奶制品原先销量（千吨）4820320210815奶制品销量（千吨）3344483 201413从表格中可以发现，销售量与价格的增减性往往相反，这是受市场的价格与价值关系规律所决定的。4.4.7问题二的求解针对问

22、题二，由于政策不允许某种产品价格上升过快，我们根据题中现有的数据，即往年各奶制品的消费量和价格得到总的消费费用，以此作为新的限制条件。同时保留第一问中的最终目标及其他约束条件不变，用matlab软件进行求解。再根据新的奶制品的价格及需求得到相对应的销售收入，将销售最大化下销售价格与费用构成的销售总收入作为对比，获得相应的经济收益。牛奶奶油奶酪1奶酪2总和价格提高百分率%010.6618.971.3935奶制品价格（元/吨） 297643 1249 1969奶制品销量（千吨） 4820 412195 414最大利润（千元） 2011600为了保证市场奶制品供应量的相对稳定，新的价格有某种价格不上

23、升。以此作为新的附加约束条件，得到新的奶制品价格，使得销售总收入为最大。价格（元/吨）的最优解分别为，即在保障销售总收入最大的前提下，牛奶的最优价格为297元/吨，奶油的最优价格为643元/吨，奶酪1的最优价格为1249元/吨，奶酪2的最优价格为1969元/吨。销售量（千吨）的最优解为,，即在保障销售总收入最大的前提下，牛奶的最优售量为4820千吨，奶油的最优售量为412千吨，奶酪1的最优售量为195吨，奶酪2的最优售量为414千吨。在价格和售量最优的情况下，可获得最高销售总收入为元。五模型的评价优点：综合考虑了原材料的限制及价格伸缩性与交叉伸缩性 基于单目标最优的方法对该问题建立线性规划模型

24、，为决策有效化提供了良好的参考办法 模型简便，准确，有相当的灵活性与实用性，具有广泛的实际用途缺点：未考虑其他因素对需求量及价格的影响 将需求量近似作为生产量 没有考虑多次价格变动对市场需求的变化 未将市场自身存在的缺点如滞后性等因素考虑在内 六附录%得到的最大利润：fmin = -2.3513e+006 %X的取值为：xmin = xmin = 0.7652 -0.1894 0.1649 1.3935对数据的初步处理：syms x1 x2 x3 x4 t1 t2 t3 t4;w=297*(1+x1)*4820*(1-0.4*x1)+720*(1+x2)*320*(1-2.7*x2)+1050

25、*(1+x3)*210*(1-1.1*x3+0.1*x4)+815*(1+x4)*70*(1-0.4*x4+0.4*x3)w1=297*(1+x1)*4820-297*(1+x1)*4820*0.4*x1+720*(1+x2)*320-720*(1+x2)*320*2.7*x2+1050*(1+x3)*210-1050*(1+x3)*210*1.1*x3+1050*(1+x3)*210*0.1*x4+815*(1+x4)*70*1-815*(1+x4)*70*0.4*x4+815*(1+x4)*70*0.4*x3w2=1431540*x1 + 230400*x2 + 220500*x3 + 5

26、7050*x4 + x4*22050*x3 + x4*22050 + x3*22820*x4 + x3*22820 - x4*22820*x4+x4*22820-x3*242550*x3 + x3*242550 - x1*572616*x1 + x1*572616 - x2*622080*x2 + x2*622080 + 1939490w3=-572616*x12 + 2004156*x1 - 622080*x22 + 852480*x2 - 242550*x32 + 44870*x3*x4 + 485870*x3 - 22820*x42 + 101920*x4 + 1939490p1=482

27、0*(1-0.4*t1)*0.04+320*(1-2.7*t2)*0.8+210*(1-1.1*t3+0.1*t4)*0.35+70*(1-0.4*t4+0.4*t3)*0.25p2=4820*(1-0.4*t1)*0.09+320*(1-2.7*t2)*0.02+210*(1-1.1*t3+0.1*t4)*0.3+70*(1-0.4*t4+0.4*t3)*0.4p1 =2699/5-1928/25*t1-3456/5*t2-1477/20*t3+7/20*t4p2 =2656/5-4338/25*t1-432/25*t2-581/10*t3-49/10*t4p11=p1-600p22=p2-

28、700%求解过程目标函数：function f=fun(x);f=-(297*(1+x(1)*4820*(1-0.4*x(1)+720*(1+x(2)*320*(1-2.7*x(2)+1050*(1+x(3)*210*(1-1.1*x(3)+0.1*x(4)+815*(1+x(4)*70*(1-0.4*x(4)+0.4*x(3)%非线性约束：functiong,ceq=mycon(x)g=(1-0.4*x(1)*(1+x(1)-1;(1+x(2)*(1-2.7*x(2)-1;(1+x(3)*(1-1.1*x(3)+0.1*x(4)-1;(1+x(4)*(1-0.4*x(4)+0.4*x(3)-

29、1;ceq=;%线性约束：x0=1;0.1;0.01;0.02; A=-1928/25 -3456/5 -1477/20 7/20;-4338/25 -432/25 -581/10 -49/10;1 0 0 0;0 1 0 0;0 0 1.1 -0.1;0 0 -0.4 0.4;b=301/5;844/5;2.5;1/2.7;1;1; Aeq=;beq=; VLB=; VUB=2.5;1/2.7; ; ; x,fval=fmincon(fun,x0,A,b,Aeq,beq,VLB,VUB); xmin=x fmin=fval%得到的结果：%得到的最大利润：fmin = -2.3513e+006 %X的取值为：xmin = xmin = 0.7652 -0.1894 0.16491.3935%求解现在的价格，销量：x1=0.7652;x2=-0.1894;x3=0.1649;x4=1.3956;m1=297*(1+x1)m2=720*(1+x2)m3=1050*(1+x