蔬菜价格变动分析及采购计划的制定

1、2016年数学建模论文第 3 套 论文题目：蔬菜价格变动分析及采购计划的制定组 别：姓 名：提交日期： 蔬菜价格变动分析及采购计划的制定 摘 要 食品价格是居民消费价格指数的重要组成部分，食品价格波动直接影响居民生活成本和农民收入，是关系国计民生的重要战略问题。蔬菜作为老百姓餐桌上必不可少的元素，其价格的变化关系到千家万户的日常生活。本文应用ARMA和AR模型来研究蔬菜价格的变动以及蔬菜价格指数的编制问题，并运用所构建的模型来进行蔬菜价格的短期预测。针对问题一，首先建立了AR模型。经合理性检验后确定为AR(6)模型，用此模型来研究这6种蔬菜价格随月份的变化规律，并预测2016年后半年以及201

2、7年前半年这6种蔬菜每月的价格，预测结果见表2。为了解决蔬菜因季节性波动和周期性波动引起的AR模型预测与实际趋势有偏差的问题，本文建立了ARMA(1,13)模型来对AR(6)模型进行改进。结果表明，两种模型都能预测出蔬菜价格在年度之内呈上涨趋势。而对于季节性不是很明显的蔬菜，ARMA模型相对于AR模型ARMA模型能更好的预测价格。针对问题二，本文首先利用SPSS软件对17种蔬菜进行了系统聚类，将17种蔬菜分为三类，通过分别计算三类蔬菜价格的平均值来给各类蔬菜对价格指数的影响程度赋予不同的权重值。然后考虑人们的消费习惯对价格指数的影响，本文查找网上资料，按销量将17种蔬菜分为五类，用各类蔬菜的销

3、量在一定程度上反映人们的消费习惯。通过各类蔬菜的销量来给各类蔬菜对价格指数的影响程度赋予不同的权重值。最后对于上述两种因素，本文凭借生活经验，人为的对两种因素赋予不同的权重值，进而计算每月蔬菜价格的加权平均价格，求出每月的定基价格指数。通过检验发现价格指数仍是一平稳的时间序列，因此同第一问一样建立ARMA模型进行研究。针对问题三，本文对问题二所得到的蔬菜价格指数进行回归分析，利用SPSS软件绘制散点图，发现在95%的置信区间内可以进行线性回归分析。然后利用SPSS软件做线性回归，得到显著性水平为0.05时，线性回归模型整体显著。由回归方程可知近几年蔬菜价格总体升高，结合蔬菜价格指数的变动情况可

4、知西安市每年一月至四月蔬菜价格总体处于高位。针对问题四，本文根据题目要求，在满足所有约束条件的情况下，以采购蔬菜的最大重量为目标函数，分别对四个蔬菜批发市场建立整数规划模型。通过LINGO软件进行求解，得出到胡家庙蔬菜批发市场进行一次采购可以使得当天采购蔬菜的总重量最大。关键词：一、 问题重述食品价格是居民消费价格指数的重要组成部分，食品价格波动直接影响居民生活成本和农民收入，是关系国计民生的重要战略问题。在收入增长缓慢的情况下，食品价格上涨将使人民群众明显感到生活成本增加，特别是蔬菜价格的变化关系到千家万户的日常生活，菜价的上涨将严重影响城市低收入群体的生活质量。为监测食品价格的实际变化情况

5、，西安市物价局对食品价格一直进行着严密的监测，每周都会在其官方网站上公布食品价格监测数据（网址：/ptl/def/def/index_1285_3890.html），为了跟踪研究西安市农副产品价格变动的规律，请从该网站下载查阅相关监测数据，建立数学模型解决如下问题： 1）请从监测的17种蔬菜数据中任意选取6种蔬菜，并根据这6种蔬菜近几年的价格数据，建立数学模型研究这6种蔬菜价格随月份的变化规律，并预测2016年后半年以及2017年前半年这6种蔬菜每月的价格。要求至少选择两种以上方法进行对比分析。2）监测17种蔬菜的价格数据给监测人员带来了很大的工作负担

6、，为了综合评价蔬菜价格的总体水平，请建立一个蔬菜价格指数模型，使这个指数的升降能够从总体上较为准确地反映蔬菜价格的水平。注意蔬菜的类型、人们的消费习惯、以及其它因素都可能与这个蔬菜价格指数有关，并说明你所建立指数的合理性。再根据你建立的蔬菜价格指数模型，研究一下近几年蔬菜价格总体的变化趋势，说明一下对于西安市每年什么时候蔬菜价格总体处于高位。3）假设你是一家饭店的采购员，每一天都要根据当天西安市四个主要蔬菜批发市场的蔬菜牌价（如表1所示），选择到其中一家市场进行一次采购。为满足饭店营业需求，饭店制定了采购单，并对采购量做出如下要求：每天必须购买的蔬菜有：豆角至少50公斤，青椒至少30公斤，土豆

7、至少20公斤，西红柿至少100公斤，莲花白至少100公斤，胡萝卜至少20公斤，茄子至少10公斤，其余蔬菜品种采购员可以酌情选择购买至少5种，如果购买这种蔬菜至少购买10公斤。假设采购所用的汽车是一辆载重量不超过1.5吨的小型三轮货车，因为车辆保养的原因，要求每天车辆的公里吨数（即路程载重吨数，空载不计）不得超过8（公里吨），每天采购额不得超过4000元，且采购单中要求必须购买蔬菜的采购额至少要达到实际采购总额的80%以上。请建立数学模型针对表1中蔬菜价格制定某天采购计划的最佳方案，即在满足上述所有条件都满足的情况下，到哪一个蔬菜市场区购买蔬菜，购买哪些蔬菜，各多少公斤，使得当天采购蔬菜的总重量

8、最大。假定蔬菜批发的最小单位为1公斤。表1. 西安市主要蔬菜批发市场某日蔬菜牌价（单位：元/公斤）序号蔬菜品种胡家庙蔬菜批发市场（与饭店距离7km）朱雀路蔬菜批发市场（与饭店距离12km）新北城蔬菜批发市场（与饭店距离11km）西三环蔬菜批发市场（与饭店距离14km）1豆角5.765.85.32茄子4.33青椒44西红柿5青菜3.93.744.26冬瓜7白萝卜28韭菜4.59黄瓜32.43.1310尖椒4.811蒜薹2.612土豆

9、2.913菜花1.414白菜1.915莲花白1.416芹菜17胡萝卜4）如果要在2016年10月1日完成3）中的任务，采购计划的最佳方案又是怎样的？ 二、问题分析 题目要求建立明确的数学模型，分别用来研究蔬菜价格随月份变化的规律，并对蔬菜价格进行预测；研究近几年蔬菜价格总体变化趋势；以及确定蔬菜采购的最优方案。问题一，要根据所选的6种蔬菜近几年的价格数据，建立数学模型研究这6种蔬菜价格随月份的变化规律，并预测2016年后半年和2017年前半年这6种蔬菜每月的价格。需要绘制6种蔬菜价格随月份变化的折线图

10、，仔细观察蔬菜价格随着月份有着什么样的变化。由于6种蔬菜价格分别是6个时间序列，因此可以考虑运用时间序列构建模型来研究蔬菜价格随月份变化的规律，并考虑比较常用的ARMA模型建模。按照建立模型的步骤，首先对序列进行稳定性检验，确定为平稳非白噪声序列后，计算自相关和偏自相关系数，进而进行ARMA模型识别；确定相应的模型后，估计模型中未知参数的值，然后对所得模型进行检验，验证模型是否有效；最后根据所得的模型预测时间序列将来的走势，从而对所选的6种蔬菜价格进行预测。 问题二，考虑到蔬菜的类型、人们的消费习惯、季节性变化等多种因素都会对蔬菜的价格指数造成影响，各种影响是一灰色系统，很难建立确定的数学关系

11、。因此本文选取了季节性变化对蔬菜价格指数的影响进行研究。由于题目中没有各蔬菜的销售量，因此我们采用简单平均法求解价格指数模型，以2011年6月各蔬菜的月均价为基期，研究近五年的蔬菜价格变化趋势，从而得出西安市每年什么时候蔬菜价格总体处于高位。问题三，要求满足题目所有的约束条件的情况下确定采购方案，使得采购蔬菜的总重量最大。显然，这是一个线性规划问题，只需根据题目要求分别对三个蔬菜批发市场建立整数规划模型，确定目标函数和约束条件，利用软件进行求解。最后比较四个结果，选择最优方案即可。问题四，要求完成2016年10月1日的采购蔬菜方案。显然，需要问题一中预测的2016年10月1日的蔬菜价格数据，建

12、立整数规划模型，确定目标函数和约束条件，利用软件进行求解。最后比较四个结果，选择最优方案即可。三、问题假设 1、假设影响蔬菜价格的主要因素为供需关系与各时期产量，与其他变量无关；2、假设在拟合过程中，蔬菜的生产成本与运输成本不发生改变；3、假设在拟合过程中，国家现有关于农产品的计划政策没有发生改变；4、假设所统计的数据全部真实可靠。四、符号说明收购点A,B,C收购点1，2，3x(i,j)从收购点i到菜市场j运送蔬菜的数量b(j)菜市场j每天对蔬菜的需求量c(j)菜市场j的短缺损失d(i)收购点i每天的蔬菜收购量A(i,j)收购点i到菜市场j的最短路程Z目标函数总费用P蔬菜调运费Q各市场供给量小

13、于需求量的短缺损失五、问题一模型建立与求解5.1 蔬菜价格随月份的变化趋势本文采用了西安市物价局公布的蔬菜价格监测数据，具体收集了豆角、茄子、芹菜、西红柿、土豆和胡萝卜6种蔬菜2011年7月至2016年6月每月的价格，并将其整理为表2。表2. 6种蔬菜价格变动（单位：元/千克）年/月价（元）均菜蔬豆角茄子芹菜西红柿土豆胡萝卜2011/73.251.592.202.801.862.082011/83.721.891.853.031.932.272011/94.472.661.513.541.782.202011/104.483.061.533.681.562.312011/114.352.881

14、.013.561.432.012011/125.343.501.074.261.451.842012/16.654.731.434.571.471.682012/27.184.731.744.531.561.882012/37.814.982.245.321.562.012012/47.484.441.975.261.502.352012/55.963.772.534.301.742.782012/64.352.691.782.832.022.932012/73.882.081.512.451.572.572012/33.341.702.692012/94.092.46

15、1.894.091.762.362012/103.702.651.263.741.781.962012/114.272.631.153.442.031.382012/125.903.111.663.832.141.372013/16.463.662.284.242.361.612013/27.883.632.704.032.401.902013/36.613.951.643.952.331.652013/46.924.392.004.832.532.252013/55.813.432.173.642.742.382013/64.402.591.913.062.462.302013/75.582

16、.382.303.252.362.472013/85.832.262.492.762.502.762013/96.392.262.513.522.432.692013/105.862.642.644.342.502.732013/115.823.262.414.502.532.472013/125.923.112.364.312.512.312014/18.953.802.175.752.582.682014/29.644.282.295.882.753.022014/38.074.862.035.292.722.982014/47.433.831.734.942.313.112014/55.

17、282.851.462014/65.062.321.513.381.892.342014/74.282.101.412.751.892.112014/84.922.081.293.451.912.142014/95.452.251.473.321.912.322014/104.842.351.523.591.802.302014/114.682.881.334.131.912.242014/125.383.891.524.682.012.282015/16.204.431.784.832.022.132015/26.884.912.075.992.152.432015/

18、36.293.831.784.862.152.302015/45.353.721.962015/54.783.261.913.831.932.472015/64.222.602.023.351.802.372015/73.931.981.973.401.752.462015/84.132.291.953.931.972.672015/94.432.531.814.232.042.552015/104.092.771.324.441.952.402015/114.483.561.434.692.402.552015/127.244.731.725.102.832.5720

19、16/16.64.832.085.533.002.742016/28.104.982.585.823.223.152016/37.224.642.765.623.333.302016/45.844.432.385.023.783.962016/56.174.392.484.543.914.412016/65.363.892.123.733.424.19绘制6种蔬菜价格随时间变化的折线图，如图1所示。由图可以看出6种蔬菜价格均存在明显的季节变动，变动的总体趋势是冬春价格较高，夏季价格较低。其中豆角、西红柿、茄子每年的价格波动较大，土豆、芹菜和胡萝卜白每年的价格波动较小。从每月价格变动情况看，由于

20、14月天气较寒冷，蔬菜供应大部分为温室蔬菜，其成本较高，所以是一年中蔬菜价格较高时期。随着天气逐渐变暖，59月是蔬菜生产旺季，蔬菜价格逐渐下降。在每年的8月左右，蔬菜价格降到最低值。随着天气逐渐转冷，1012月蔬菜供应逐渐减少，加上北方主要依靠温室种植蔬菜，其成本逐渐提高，蔬菜价格也随之上升。在每年的2月左右，蔬菜价格升到最高值。图1 6种蔬菜月平均价格变动图5.2 AR模型的建立与求解5.2.1 AR模型的建立设是零均值平稳序列，满足下列模型： （1）其中是零值、方差是的平稳白噪声，则称是阶数为p的自回归序列，简记为AR(p)序列，而 称为自回归参数向量，其分量 称为自回归系数。引进算子描述

21、（1）较为方便，算子B的定义如下：记算子多项式 则式（1）可以改写为： 5.2.2 AR模型的求解 由预报差分方程又 给出预报的递推公式： 由此可见， 仅仅依赖于的k时刻以前的p个时刻的值 。这是AR(p)预测的特点。5.2.3 结果分析以2011年6月到2016年6月的数据作为样本，采用AR模型对2016年7月到2017年6月进行预报。首先对AR模型的参数进行确定，经计算采AR(6) 通过合理性检验，所以采用AR模型进行预测。预测结果如图2所示。由图3预测结果可看出，豆角的最高价格出现在2月份，最低价格出现在10月份。芹菜和茄子的最高价格与最低价格变化与豆角类似，只是变化比较平缓。土豆和胡萝

22、卜的价格变化不是很明显。其总体的价格趋势与2011年7月到2016年6月的趋势是相符的。蔬菜整体的价格趋势处于上涨趋势。其中豆角、茄子、芹菜、西红柿这四种蔬菜呈现明显的季节性变动，而土豆和胡萝卜的季节性不是很明显，也可能是由于蔬菜生长受季节性的影响。与此同时，蔬菜价格受季节性影响的同时又各不相同，各有各的特点。 图2 AR模型预测结果5.3 ARMA模型的建立与求解5.3.1 ARMA模型的建立设是零均值平稳序列，满足下面模型： (2)其中是零均值、方差是的平稳白噪声，则称是阶数为p,q的自回归滑动平均序列，简记为ARMA(p,q)序列。当q=0时，它是AR(p)序列；当p=0时，它是MA(q

23、)模型。应用算子多项式，式(2)可以写为 对于一般的平稳序列 ，设其均值，满足下列模型： （3）利用后移算子，式（3）可以表示为 时间序列的m步预报，是根据的取值对未来k+m时刻的随机变量 做出估计。估计量记作，它是的线性组合。引进估计空间 现设是零均值正态ARMA(p,q)序列。所谓平稳线性最小均方预报是指：，且使得 在正态性条件下是对于的条件期望，即 因为是一般相关的，因此，将用的元素的线性表示并不好意思很方便。考虑到对于ARMA(p,q)序列，总存在传递形式，令 则可以证明，这就是说，可以用中的元素来表示。5.3.2 ARMA模型的求解平稳线性最小均方预报具有下面性质。如果是零均值的AR

24、MA(p,q)序列，则下面的预报差分方程成立 （4）事实上，因 当时对于ARMA(p,q)序列，由预报差分方程（4）可知，只需要知道 就可以递推出。定义预报向量 令 可证下面递推预报公式 （5） 式（5）中第三项当时为0。由可逆性条件可知，当较小时，可令初值。 在实际中，参数的模型是未知的。若已建立了时间序列模型，则理论模型中得参数可用估计参数代替，再进行预报。5.3.3 结果分析以2011年6月到2016年6月的数据作为样本，采用ARMA模型对2016年7月到2017年6月进行预报。首先对ARMA模型的参数进行确定，经计算采ARMA(1,13) 通过合理性检验，所以采用ARMA模型进行预测。

25、预测结果如图3所示。 图3 ARMA模型预测结果由图3预测结果可看出，豆角的最高价格出现在2月份，最低价格出现在10月份。茄子的最高价格比豆角迟两个月，出现在12月份，最低价格也是出现在10月份。芹菜的最高价格与最低价格变化与豆角类似，只是变化比较平缓。土豆和胡萝卜的价格变化不是很明显。其总体的价格趋势与2011年7月到2016年6月的趋势是一致的。蔬菜整体的价格趋势处于上涨趋势。其中豆角、茄子、芹菜、西红柿这四种蔬菜呈现明显的季节性变动，而土豆和胡萝卜的季节性不是很明显，也可能是由于蔬菜生长受季节性的影响。与此同时，蔬菜价格受季节性影响的同时又各不相同，各有各的特点。5.4 两种模型预测结果

26、对比分析本文采用了AR模型和ARMA模型来预测蔬菜的价格。由于蔬菜的价格呈现季节性波动、周期性波动和随机性波动。针对季节性波动本文选取了季节性蔬菜豆角和非季节性蔬菜土豆作为典型的研究对象。图4是两种模型下季节性蔬菜豆角和非季节性蔬菜土豆在两种不同的模型下预测的结果。 图4 季节性与非季节性蔬菜在两种模型下的对比由图4可以看出，季节性蔬菜豆角在两种模型下的预测结果差异不是很大，主要原因是ARMA模型和AR模型对于波动性比较大的能够准确预测。对于波动性不是很大、价格趋势不是很明显的土豆来说，相对于AR模型，ARMA模型能更好的预测价格趋势，而AR模型预测的价格偏低。 总体来说，两种模型都能预测出蔬

27、菜价格在年度之内呈上涨趋势，在年内季节性蔬菜呈季节性变动，能反映出不同的蔬菜有不同的特点。除此之外，两种模型也能够较为准确地预测出蔬菜价格的最低点和最高点，反映出蔬菜价格变动的周期性规律。六、问题二模型建立与求解6.1 蔬菜价格指数相关知识在经济统计领域的研究中，统计指数是一种相对数，用来表示某种社会经济现象的变动情况。从广义上来说，一切用来表示社会经济现象数量变动的相对数都是指数，不仅可以反映某一变量随时间变化的相对数，还能够反映某一变量随地域的变化。从狭义上来说，统计指数是一种以相对数（数值比数）的形式综合反映多个变量在数值上总变动情况的方法。价格指数作为一种统计指数，反映着某种商品的价格

28、水平在不同时期变化的方向、程度和趋势。根据统计指数的概念，可以把狭义的价格指数k简单的理解为某一个变量在某一时期（记为时刻t=1）内的价格与该变量在另一个作为比较标准的时期（记为时刻t=0）内的价格的比数。其中作为比较标准的时期即t=0被称为基期；当前时期，即t=1被称为报告期。也就是说把基期的价格指数看成100，以基期水平为标准，计算比值，得到报告期的指数，以此来判断价格在不同时期的变动情况。根据指数计算所选用的基期不同，可以把指数分为定期指数和环比指数两种：定期指数是始终以某一固定的时期作为基期来计算各期的指数；环比指数是始终以前一期作为基期来计算各期的指数。环比价格指数反映相邻两期的价格

29、变化，基期不固定，在长期内没有可比性，因此本文只对蔬菜价格指数的月定基价格指数进行研究，不研究对月环比价格指数。6.2 问题二模型建立及求解6.3.1 蔬菜价格指数模型建立 本文采用简单的平均法来求蔬菜价格指数，以2011年6月作为基期，应用表2中6种蔬菜的月均价数据，建立价格指数模型。 （6）蔬菜种类豆角茄子芹菜西红柿土豆胡萝卜（元/公斤）3.251.592.202.801.862.08表3 2011年7月6种蔬菜的月均价6.3.2 模型求解由表2可查出的值，对于第i种蔬菜月份为j的权重的计算，即同理；接着通过式（1）对2011年7月至2016年6月蔬菜的价格指数I进行计算。 表4 2011

30、年7月-2016年6月蔬菜的价格指数2011/72011/82011/92011/102011/112011/122012/12012/211.091.271.311.271.541.891.962012/32012/42012/52012/62012/72012/82012/92012/102.162.051.661.241.311.221.281.182013/112013/122013/12013/22013/32013/42013/52013/61.221.561.721.971.751.901.561.242013/72013/82013/92013/102013/112013/12

31、2014/12014/21.441.471.611.601.621.612.312.472014/32014/42014/52014/62014/72014/82014/92014/102.171.971.491.341.141.311.401.322014/112014/122015/12015/12015/32015/42015/52015/61.381.611.792.051.761.531.391.222015/72015/82015/92015/102015/112015/122016/12016/21.161.271.351.341.482.011.962.222016/32016

32、/42016/52016/62016/72016/82016/92016/102.061.851.881.65 6.3 结果分析由蔬菜价格指数趋势图可知，此模型很好的反应出了蔬菜价格在近5年的价格随季节变化的趋势，指数的升降代表了一年中不同季节蔬菜价格的上升和下降，说明了季节因素对蔬菜价格指数的影响较大。通过指数模型研究发现，近五年蔬菜价格指数的总体趋势是：2、3月份各蔬菜价格达到最高。4、5月份各蔬菜价格逐渐变低至6、7月份各蔬菜价格达到最低，而从8月份以后蔬菜价格又逐渐增高。本文分析蔬菜价格波动的原因：（1）气候的影响。西安地区四季分明，冬天寒冷，需要在大棚里面种植的蔬菜，这就增加了种植蔬

33、菜的成本。所以冬季的蔬菜价格明显要高于夏季。（2）节假日的影响。前一年的年底月到后一年的2月为过年影响时间段，人们大量采购蔬菜，供小于求，因此也会涨价。（3）国家的调控。国家每年都会制定各种价格调控政策，这也会在某一时间段影响价格的变化。因此西安市每年2月份的时候蔬菜价格总体处于高位。图5 2011年7月-2016年6月蔬菜价格指数图八、问题四模型建立与求解8.1 预测各菜市场蔬菜价格本题需要应用问题二的模型预测2016年10月1日四个菜市场17种蔬菜价格的数据。具体步骤：（1）根据问题二结果可知2016年10月1日6种蔬菜价格，并应用问题二所建模型和表一中的数据预测胡家庙蔬菜批发市场的其他1

34、1种蔬菜价格，从而得到2016年10月1日胡家庙蔬菜批发市场17种蔬菜价格；（2）根据表一中其他三个菜市场各个蔬菜价格与胡家庙蔬菜批发市场中与之对应的蔬菜价格之比，并且考虑到受季节变化影响预测出2016年10月1日另外三个菜市场的蔬菜价格，可得表5为西安市主要蔬菜批发市场2016年10月1日的各蔬菜价格。表5. 西安市主要蔬菜批发市场10.1日蔬菜牌价（单位：元/公斤）序号蔬菜品种胡家庙蔬菜批发市场（与饭店距离7km）朱雀路蔬菜批发市场（与饭店距离12km）新北城蔬菜批发市场（与饭店距离11km）西三环蔬菜批发市场（与饭店距离14km）1豆角3.32茄子3.

35、53青椒34西红柿5青菜6冬瓜2.32.6227白萝卜1.68韭菜3.89黄瓜2.310尖椒3.611蒜薹21.81.8212土豆3.533.22.713菜花1.114白菜1.615莲花白1.116芹菜1.817胡萝卜8.2 模型建立本题所用模型为问题三线性规划模型，其目标函数和约束条件都一致，只有蔬菜价格数据不同。目标函数： 约束条件： 8.3 模型求解根据约束条件和表5的

36、数据，求解得出的具体采购方案如表6所示；表5 某饭店2016年10月1日蔬菜采购计划 蔬菜品种胡家庙蔬菜批发市场（与饭店距离7km）朱雀路蔬菜批发市场（与饭店距离12km）新北城蔬菜批发市场（与饭店距离11km）西三环蔬菜批发市场（与饭店距离14km）豆角50778267茄子10101010青椒30373834西红柿415208227167青菜10101010冬瓜10101010白萝卜10101010韭菜10101010黄瓜15101010尖椒14101010蒜薹10101010土豆39272824菜花10101010白菜10101010莲花白418209227166芹菜10101010胡萝卜

37、40272824总质量11116957405928.3 结果分析由表5可知，在满足所有的约束条件和所预测的数据情况下，到胡家庙需购蔬菜1111公斤，到朱雀路需购蔬菜695公斤，到新城需购蔬菜740公斤，到西三环需购蔬菜为592公斤。可见由于胡家庙距离饭店最近，所以购买的蔬菜的重量最大，而西三环距离饭店距离最远，所以购买的蔬菜重量最少，所以九、模型评价与改进9.1模型优点1、本文在建立模型前，采用均值、标准化、剔除等方法对原始数据进行了处理，使得模型更为准确。2、本文所建立的ARMA模型可以较好的预测蔬菜价格以及价格指数的变化，相对误差较小。3、问题二采用简单平均法说明了蔬菜价格指数的总体变化趋

38、势，方法简单，结果可靠。4、问题三采用线性规划模型制定了采购蔬菜的最佳方案，模型简单，结果可靠。9.2模型缺点1、在得出相关ARMA模型参数后，只对模型进行了检验，以确定所建立的模型正确。但是，在确定模型正确后没有对模型进行进一步的优化，因此本文所建立的模型并不一定是最优的ARMA模型。3、问题二在人们消费习惯和蔬菜种类对蔬菜的影响没有考虑进去，使得蔬菜价格指数缺乏真实性。参考文献1 张志涌，杨祖樱. MATLAB教程M. 北京：北京航空航天大学出版社, 2011.2 姜启源，谢金星，叶俊.数学建模(第三版)习题参考解答M,北京，高等教育出版社，2002.3 齐欢.数学模型方法M.武汉：华中理

39、工大学出版社，2005.4 谢金星，薛毅.优化建模与LINDO/LINGO软件M.北京，清华大学出版社，20055 李工农，阮晓，青徐晨.经济预测与决策及其MATLAB实现M.北京：清华大学出版社，2007.附录1.问题一程序AR模型：clc,clearload data.txt %原始数据以列向量的方式存放在纯文本文件中yt=data; n=length(yt);alpha=0.2 0.5 0.8;m=length(alpha);yhat(1,1:m)=(yt(1)+yt(2)/2;for i=2:nyhat(i,:)=alpha*yt(i-1)+(1-alpha).*yhat(i-1,:)

40、;endyhaterr=sqrt(mean(repmat(yt,1,m)-yhat).2)xlswrite(dianqi.xls,yhat)yhat1=alpha*yt(n)+(1-alpha).*yhat(n,:) ARMA模型：doujiao.mclcclearload doujiao.txt; %把原始数据按照表中的格式存放在纯文本文件 doujiao.txtdoujiao=doujiao;x=doujiao(:);s=12; %周期 s=12n=12; %预报数据的个数m1=length(x); %原始数据的个数for i=s+1:m1 y(i-s)=x(i)-x(i-s);endm2

41、=length(y); %周期差分后数据的个数w=diff(y); %消除趋势性的差分运算m3=length(w); %计算最终差分后数据的个数for i=0:3 for j=0:s+1 spec= garchset(R,i,M,j,Display,off); %指定模型的结构 coeffX,errorsX,LLFX = garchfit(spec,w); %拟合参数 num=garchcount(coeffX); %计算拟合参数的个数 %compute Akaike and Bayesian Information Criteria aic,bic=aicbic(LLFX,num,m3);

42、fprintf(R=%d,M=%d,AIC=%f,BIC=%fn,i,j,aic,bic); %显示计算结果 endendsave bdata x y w n m1 m2 s yuce.mclcclearload bdataspec2= garchset(R,1,M,13,Display,off); %指定模型的结构coeffX,errorsX,LLFX = garchfit(spec2,w); %拟合参数sigmaForecast,w_Forecast = garchpred(coeffX,w,n) %求 w 的预报值yhat=y(m2)+cumsum(w_Forecast) %求 y 的预

43、报值for j=1:n x(m1+j)=yhat(j)+x(m1+j-s);endx_hat=x(m1+1:end) %复原到原始数据的预报值 2.问题二程序clcclearp0=5.06 2.32 1.51 3.381.82.34;p1=8.953.8 2.175.752.582.68; 9.644.282.295.882.753.02; 8.074.862.035.292.722.98; 7.433.831.734.942.313.11; 5.282.851.46; 5.062.321.513.381.892.34; 4.282.1 1.412.751.892.11

44、; 4.922.081.293.451.912.14; 5.452.251.473.321.912.32; 4.842.351.523.591.8 2.3; 4.682.881.334.131.912.24; 5.383.891.524.682.012.28; sum1=0; sum0=0; for j=1:1:12 for i=1:1:6 sum1=sum1+p1(j,i)*p1(j,i); sum0=sum0+p0(i)*p1(j,i); end I=sum1/sum0 sum1=0; sum0=0; end I=1.6659,1.7835,1.6147,1.4539,1.1261,1.0

45、032,0.8706,0.9720,1.0340,0.9940,1.0487,1.2245; t=1:12; h=plot(t,I,-cd); legend(h,蔬菜价格指数);3.问题三、四程序胡家庙蔬菜批发市场采购方案：MODEL: TITLE胡家庙蔬菜批发市场采购方案; SETS: VEGETABLE/1.17/:C,WEIGHT;ENDSETS DATA: C=5.7,4.5,4.6,6,3.9,3.1,2.2,5.5,3,4.7,2.6,3.7,1.6,2.6,1.6,2.2,5;ENDDATA MAX=SUM(VEGETABLE(I):WEIGHT(I); WEIGHT(1)50;

46、WEIGHT(2)10;WEIGHT(3)30; WEIGHT(4)100; WEIGHT(12)20; WEIGHT(15)100; WEIGHT(17)20; W510;W610;W710;W810;W910;W1010;W1110;W1310;W1410;W1610; WEIGHT(5)=IF(X5#EQ#0,0,W5); WEIGHT(6)=IF(X6#EQ#0,0,W6); WEIGHT(7)=IF(X7#EQ#0,0,W7); WEIGHT(8)=IF(X8#EQ#0,0,W8); WEIGHT(9)=IF(X9#EQ#0,0,W9); WEIGHT(10)=IF(X10#EQ#0

47、,0,W10); WEIGHT(11)=IF(X11#EQ#0,0,W11); WEIGHT(13)=IF(X13#EQ#0,0,W13); WEIGHT(14)=IF(X14#EQ#0,0,W14); WEIGHT(16)=IF(X16#EQ#0,0,W16);7.2*(SUM(VEGETABLE(I):WEIGHT(I)8000; SUM(VEGETABLE(I):C(I)*WEIGHT(I)5; (WEIGHT(1)*C(1)+WEIGHT(2)*C(2)+WEIGHT(3)*C(3)+WEIGHT(4)*C(4)+WEIGHT(12)*C(12)+WEIGHT(15)*C(15)+WE

48、IGHT(17)*C(17)-0.8*(SUM(VEGETABLE(I):C(I)*WEIGHT(I)0; BIN(X5);BIN(X6);BIN(X7);BIN(X8);BIN(X9);BIN(X10);BIN(X11);BIN(X13);BIN(X14);BIN(X16); FOR(VEGETABLE(I):GIN(WEIGHT(I);END 朱雀蔬菜批发市场采购方案：MODEL: TITLE朱雀蔬菜批发市场采购方案; SETS: VEGETABLE/1.17/:C,WEIGHT; ENDSETS DATA: C=6,4.3,4.5,5.4,3.7,3,1.9,4.6,2.4,4.9,2.4,3.2,1.8,2.1,1.1,2,4.9; ENDDATA MAX=SUM(VEGETABLE(I):WEIGHT(I); WEIGHT(1)50; WEIGHT(2)10; WEIGHT(3)30; WEIGHT(4)100; WEIGHT(12)20; WEIGHT(15)100; WEIGHT(17)20;W510;W610;W710;W810;W910;W1010;W