第四届文鼎创杯华中地区大学生数学建模邀请赛
题目：城市居民食品分类及零售价格预测
摘 要
本文第一部分是关于城市居民消费食品的分类问题，根据武汉市 2010 年 3 月～
2011 年 3 月的食品零售价格数据，建立谱系聚类分析模型，即模型一，采用欧式距离平
方统计分析的方法，应用“”统计分析软件，得出了一种食品分类方式，并分析
了每类食品的特点。
本文第二部分是关于居民食品零售价格的预测问题，从时间序列和灰色系统两个视
角解决该问题，建立了两个预测模型。
模型二：这里时间序列理解为食品价格按照时间顺序排列而成的一组观测值，它由
食品价格在不同时间上的观测值和食品价格所属的时间构成，并认为不同年份在同一月
份时间段内的变化率相同，建立逐期增长量模型，用 2010 年 4、5 月份的价格环比增长
率预测出了 2011 年 4、5 月的食品零售价格，并利用 Matlab 程序对模型进行了计算和
分析。
模型三：基于 模型，即一个变量一阶灰色预测模型，引入其基本形式
，是一个近似的差分微分方程，因其具有微分、差分、指数兼容的性
质，将系统看成一个随时间而变化的函数。在此模型的基础上讨论食品零售价格走势，
进行最小二乘估计，根据已知数据算得模型系数，由此预测出了 2011 年 4、5 月份武汉
市居民食品零售价格，并对预测结果与已有数据进行了误差分析。
根据所建立的两个预测模型，对最后得到的两组结果进行了比对，最后对研究合理
的规范性做了充分和必要的论证。
本文第三部分是对于武汉市居民食品零售价格情况的分析，并通过市场价格基本态
势的预测，显示出 4、5 月蔬菜价格受春节节日效应消褪和三月份回落影响，整体呈回
落态势，市场供需情况较为平衡等。根据所得结果，可以帮助有关部门在保证民生和稳
定物价等方面制定出合理的应对方案。
1,1GM
0 1x k az k b
1 问题的提出
消费者物价指数（Consumer Price Index），英文缩写为 CPI，也称消费价格指数，
是反映与居民生活有关的产品及劳务价格统计出来的物价变动指标，通常作为观察通货
膨胀水平的重要指标，是与人民生活密切相关的参考指标。
而城市居民食品零售价格是消费者物价指数的重要组成部分，权威机构研究认
为粮食生产、流通成本上涨一定会带动农产品价格总体上涨，特别是 2011 年异常的气
候情况，导致生产成本大量增加，国际粮价对国内供需的影响，食品价格未来可能发生
上涨。刚公布 3 月份的 CPI 增幅达 %，创 32 个月来的新高，这使得年内的通货膨胀
压力正在增强，对于国民经济的发展非常不利的。
（1）若能根据一定时间段内部分食品的零售价格信息，建立恰当的数学模型对食品
进行分类，并分析每类食品的特点，这对于更好地保障民生、发展经济显得格外重要；
（2）若能根据在给定时间段内部分食品的零售价格信息，准确预测未来几个月城市
居民食品零售价格走势，十分利于政府采取一定的措施，防止因价格的过快增长而产生
的危机；
（3）若能对于城市居民食品零售价格的情况进行分析，同时对有关部门提出相关的
建议，这些对于稳定人民的生活十分迫切。
2 问题的分析
对于问题 1，根据附录（数据为城市居民经常购买的 42 种不同类别食品，在 2010
年 3 月 5 日—2011 年 3 月 25 日期间，每隔 10 天给出的食品的零售价格）所给出的城市
居民食品零售价格数据，进行更近一步分析，发现要得出食品的适当分类，必须先对数
据进行统一的标准化处理，再采用“”统计分析工具[1]，进行谱系聚类分析得到
聚类树形图，最后对于分析的结果进行检算。
对于问题 2，根据附录数据进行的可视化图表分析，发现在给定的时间段内，根据
统计学的时间序列分析[2]，建立相应的预测模型，对于该预测模型采用 Matlab 软件进行
计算，但是这样做有失完整性，故在时间序列预测模型的基础上可以考虑更好的预测模
型，首先想到的预测方法便是灰色预测法，可以继续采用 Matlab 软件[4]进行计算和误差
检验。
对于问题 3，在对附录数据进行更进一步的分析后，在参考前两问的基础上，可以
利用主成分分析对数据进行研究，以及找出影响物价大幅度变化的因素和得出相应的结
论。从而提出对于相关部门提出一定的建议，可以使得有关部门在保证民生和稳定物价
方面做出合理的措施，同时这些建议和意见也必须具有可操作性。
3 模型的假设
(1)各食品价格会受到多种因素的影响；
(2)各食品变量间彼此独立；
(4)两年相应月份影响价格变化率的因素相同；
(5)所用数据稳定，无奇异数据。
4 符号说明
符号 说明
食品
食品
食品 、 之间的距离
ix i
jx j
( , )i jd x y i j
符号 说明
食品 类
食品 类数量
食品 类数量
类 和 之间的距离
每个食品 指标量
食品个数
欧氏距离
欧氏距离平方
两两之间距离构成的对称矩阵
是一个 阶对称阵
表示现象所属的时间
表示现象在不同时间上的观测值
称为最初发展水平
称为最末发展水平
整个观测期内的发展水平与参照基期
观测期内增长的绝对数量
环比增长速度
原始价格序列
为一阶累加序列 的紧邻均值序列
值为发展系数， 反映模拟值的未来发展情况
为灰色作用量
准光滑条件系数
准指数规律系数
相对误差值
5 模型的建立和求解
问题一模型的建立与求解
模型一的建立
本文就食品聚类予以讨论，并为简单起见，以 、 等分别表示食品 ， ,以 表
示食品 、 之间的距离 ，用 和 表示两个类，它们所包含的食品个数分别
记为 和 ，类 和 之间的距离用 表示。
先就食品聚类予以讨论，为简单起见，以食品间的“相近性”度量—距离设每个食品
有 个指标，它们的观测值可表示为：
这时，每个食品 可看成 维空间中的一个点， 个食品就组成 组空间中的 个
点，很自然地用各点之间的距离来衡量各食品之间的靠近程度。
在聚类分析中，有时所用的距离并不满足 (三角不等式)，
在广义的角度上仍称它为距离。
为此，引入 Euclidean distance（欧氏距离）：
qG q
pn p
qn q
( , )p qD G G pG qG
p ix
n
( , )EUCLID x y
( , )SEUCLID x y
D
(1)D 1n
it
iY
1Y
nY
0t
i
iG
(0)X
(1)z (1)X
a
b
( )k
(1) ( )k
i j ix jx ijd
i j ( , )i jd x y pG qG
pn qn pG qG ( , ) min | ,p q ij p qD G G d i G j G
ix p 1 2( , , , ) , 1, 2, , .
T
i i i pix x x x i n
ix p n p n
( , ) ( , ) ( , )i j j i k jd x x d x x d x x
（1）
Squared Euclidean distance（欧氏距离平方）：
（2）
在考虑实际情况后，本文模型一采用 Squared Euclidean distance（欧氏距离平方）
模型，区组间联结为主要影响因素。
模型一的求解步骤
欧氏距离在 型聚类中是比较常见的，但当指标的测量值相差悬殊时，应先对数据
进行标准化，然后再用标准化的数据计算距离，具体求解步骤如下：
令： （3）
（4）
（5）
对数据的标准化采用以下方法：
（6）
（7）
这里不考虑 个指标的相关性问题。
有了食品之间的距离以及类与类之间距离的定义后，便可进行谱系聚类，其基本步
骤归纳如下：
1. 个食品一开始就作为 个类，计算两两之间的距离构成一个对称矩阵 ,
其对角线上的元素全为零(对相似系数矩阵，其对角线上元素全为 1)。显然，此时有
,记 。
2.选择 中对角线元素的下三角部分（或上三角部分）中的最小元素，设其为
，则将 与 合并成一个新类 。在 中划去 与 所对应
的两行与两列，并加入由新类 与剩下的未聚合的各类之间的距离所组成的
一行和一列，得到一个新的距离矩阵 ， 是一个 阶对称阵（若在 中最小元
素不唯一，对其它的最小元素也作如上相同处理，每合并两类，矩阵 则降低一
阶）。
3.由 出发，重复步骤 2 得到对称矩阵 ，从 出发得到 ，依次类推，直
到 个食品聚为一个大类为止，经“” [3]进行分析后，得出表 5-1 食品聚类表。
表 5-1： 食品聚类表
群集成员
案例
18 群
集
17 群
集
16 群
集
15 群
集
14 群
集
13 群
集
12 群
集
11 群
集
10 群
集
9 群
集
8 群
集
7 群
集
6 群
集
1: 菜籽油 1 1 1 1 1 1 1 1 1 1 1 1 1
2: 大豆油 1 1 1 1 1 1 1 1 1 1 1 1 1
3: 花生油 2 2 2 2 2 2 2 2 2 2 2 2 2
2 ( , ) ( )i i
i
EUCLID x y x y
2( , ) ( )i i
i
SEUCLID x y x y
Q
1
1
,
n
i
i
x x
n
1,2, , p
2 2
1
1
( ) ,
1
n
i
i
s x x
n
1,2, , p
1 1
,i ii n i n
r max x min x
1,2, , p
* ,ii
x x
x
s
1,2, , ,i n 1,2, , p
' 1 ,
i ii n
i
x min x
x
r
1,2, , ,i n 1,2, , p
p
n n ( )ij n nD d
( , )p q pqD G G d (0)D D
(0)D
( , )p qD G G pG qG ,r P qG G G (0)D pG qG
,r P qG G G
(1)D (1)D 1n (0)D
(0)D
(1)D (2)D (2)D (3)D
n
4:大豆调和油 3 3 3 3 3 3 3 3 3 3 3 3 3
5:鲜猪肉(精瘦) 4 4 4 4 4 4 4 4 4 4 4 4 4
6:鲜猪肉(肋条) 5 5 5 5 5 5 5 5 5 5 5 5 5
7: 鲜牛肉 6 6 6 6 6 6 6 6 6 6 6 6 6
8:鲜羊肉(去骨) 7 7 7 7 7 7 7 7 7 7 7 6 6
9:鲜羊肉(带骨) 8 8 8 8 8 8 8 8 6 6 6 6 6
10: 活鸡 9 9 9 9 9 9 4 4 4 4 4 4 4
11: 鸡肉 10 10 10 5 5 5 5 5 5 5 5 5 5
12: 鸡蛋 11 11 11 10 10 10 9 9 8 8 8 7 1
13: 带鱼 12 12 12 11 11 11 10 10 9 9 5 5 5
14: 草鱼 1 1 1 1 1 1 1 1 1 1 1 1 1
15: 鲤鱼 1 1 1 1 1 1 1 1 1 1 1 1 1
16: 芹菜 13 13 13 12 12 10 9 9 8 8 8 7 1
17: 大白菜 14 14 14 13 13 12 11 11 10 8 8 7 1
18: 油菜 15 15 14 13 13 12 11 11 10 8 8 7 1
19: 黄瓜 13 13 13 12 12 10 9 9 8 8 8 7 1
20: 萝卜 14 14 14 13 13 12 11 11 10 8 8 7 1
21: 茄子 16 11 11 10 10 10 9 9 8 8 8 7 1
22: 西红柿 13 13 13 12 12 10 9 9 8 8 8 7 1
23: 土豆 13 13 13 12 12 10 9 9 8 8 8 7 1
24: 胡萝卜 13 13 13 12 12 10 9 9 8 8 8 7 1
25: 青椒 13 13 13 12 12 10 9 9 8 8 8 7 1
26: 尖椒 17 16 15 14 14 13 12 1 1 1 1 1 1
27: 圆白菜 14 14 14 13 13 12 11 11 10 8 8 7 1
28: 豆角 17 16 15 14 14 13 12 1 1 1 1 1 1
29: 蒜苔 18 17 16 15 1 1 1 1 1 1 1 1 1
30: 韭菜 13 13 13 12 12 10 9 9 8 8 8 7 1
31: 芦柑 13 13 13 12 12 10 9 9 8 8 8 7 1
32: 苹果 1 1 1 1 1 1 1 1 1 1 1 1 1
33: 香蕉 13 13 13 12 12 10 9 9 8 8 8 7 1
34: 西瓜 13 13 13 12 12 10 9 9 8 8 8 7 1
35: 豆腐 13 13 13 12 12 10 9 9 8 8 8 7 1
36: 食用盐 14 14 14 13 13 12 11 11 10 8 8 7 1
37: 绵白糖 5 5 5 5 5 5 5 5 5 5 5 5 5
38: 白砂糖 1 1 1 1 1 1 1 1 1 1 1 1 1
39: 红糖 18 17 16 15 1 1 1 1 1 1 1 1 1
40: 酱油 1 1 1 1 1 1 1 1 1 1 1 1 1
41: 醋 16 11 11 10 10 10 9 9 8 8 8 7 1
42: 鲜牛奶 16 11 11 10 10 10 9 9 8 8 8 7 1
4.在合并过程中记下两类合并时食品的编号以及合并两类时的距离或相似系数的大
小(我们称之为水平)，并绘成聚类的谱系图，然后可根据市场、气候等背景和要求选定
相应的临界水平确定出类的个数，经过分析得到群集数图 5-1：
图 5-1 谱系聚类聚集数
进一步分析，得到食品聚类谱系图 5-2：
图 5-2 食品聚类谱系图
结论：在分析所有的分类后，发现在分为 7 类后较合适，分类结果为：
第一类：菜籽油，大豆油，草鱼，鲤鱼，尖椒，豆角，蒜苔，苹果，白砂糖，红糖，酱
油；
第二类：花生油；
第三类：大豆调和油
第四类：鲜猪肉（精瘦），活鸡；
第五类：鲜猪肉（肋条），鸡肉，带鱼，绵白糖；
第六类：鲜牛肉，鲜羊肉(去骨)，鲜羊肉(带骨)；
第七类：鸡蛋，芹菜，大白菜，油菜，黄瓜，萝卜，茄子，西红柿，土豆，胡萝卜，青
椒，圆白菜，韭菜，芦柑，苹果，香蕉，西瓜，豆腐，食用盐，醋，鲜牛奶。
问题二模型的建立与求解
模型二的建立
时间序列模型时间序列式同一现象按照时间顺序排列而成的一组观测值，它由两部
分构成：现象在不同时间上的观测值和现象所属的时间。现象的观测值根据表现形式的
不同有绝对数、相对数和平均数等，因此时间序列可以分为绝对数时间序列、相对数时
间序列和平均数时间序列。对于本文，每月的零售商品价格是绝对数时间序列，价格消
费指数就是相对数时间序列，由绝对数时间序列可以派生出相对数时间序列和平均数时
间序列。
时间析序列分的第一步通过作图来观察数据碎时间变化的规律，然后在此基础上展
开分析和建模。
时间序列中，令 表示现象所属的时间， 表示现象在不同时间上的观测值，也称
为发展水平。设 的取值为 ，相应的发展水平为 ，则 称为最初发展
水 平 ， 称 为 最 末 发 展 水 平 。 若 将 整 个 观 测 期 内 的 发 展 水 平 与 参 照 基 期
的发展水平做对比，则 对应的发展水平 称为基期发展水平；而分
析研究的其他时期称为报告期，所对应的指标值称为报告期发展水平。
图 5-3 各食品价格随时间变化
增长量用来描述现象在观测期内增长的绝对数量，由报告期发展水平减去基期发展
水平得到。根据基期的不同选择，增长量分为逐期增长量和累计增长量。设时间序列观
测值为 ，增长量为 ，其计算公式为：
逐期增长量 ， （8）
it iY
it 1 2, , , nt t t 1 2, , , nY Y Y 1Y
nY
0 0 1 2, , , nt t t t t 0t 0Y
0,1, ,iY i n
1i i iY Y 1,2, ,i n
累计增长量 ， （9）
根据基期的不同选择，增长速度可以分为环比增长速度和定基增长速度，这选用环
比增长速度进行计算设时间序列观测值为 ，增长速度为 ，其计算公式为：
环比增长速度 ， （10）
时间序列模型求解
根据在 2010 年 4 月至 5 月的数据利用环比增长速度 可以得到在
2011 年 4 月至 5 月的预测价格。
日期
食品种类
菜籽油 散装
元/500
克
大豆油 散装
元/500
克
花生油
一级桶
装
元/5 升
大豆调
和油
一级桶
装
元/5 升
鲜猪肉 精瘦肉
元/500
克
鲜猪肉 肋条肉
元/500
克
鲜牛肉
新鲜去
骨
元/500
克
鲜羊肉
新鲜去
骨
元/500
克
鲜羊肉
新鲜带
骨
元/500
克
活鸡
活肉鸡
公
斤
元/500
克
鸡肉
白条鸡、
开膛 上
等
元/500
克
鸡蛋
新鲜完
整 鸡场
蛋
元/500
克
带鱼
冻 250
克左右
元/500
克
草鱼
活 1000
克左右
一条
元/500
克
模型三的建立
随着社会的不断进步和现代化建设的日益发展，人们在土地规划、经济决策、大型
工程施工及安全监测、灾情预报等方面都越来越离不开科学的预测方法， 好的预测方
法和科学的预测成果无疑能对上述工作提供帮助， 减少不必要的损失， 而用于预测的
0i iY Y 1,2, ,i n
1 2, , , nY Y Y G
1
1 1
1i i ii
i i
Y Y Y
G
Y Y
1,2, ,i n
1
1 1
1i i ii
i i
Y Y Y
G
Y Y
方法很多， 在众多的预测方法中， 较好的预测方法是灰色预测法[5]， 这已被许多专家
和学者认同。本文在结合附录数据无规律的变化情况下，由于不涉及到其它外生变量，
因此考虑灰色系统理论下的 模型，只从系统的背景序列对未来的食品零售价格
趋势进行预测。
假设原始价格序列 为，将附录中每 10 天给出的均价作为一个值，因此每个商
品有 39 个价格值： ， （11）
则 模型的基本形式为:
， （12）
其中: 为一阶累加序列 的紧邻均值序列; 值为发展系数， 反映模拟值的未
来发展情况; 为灰色作用量，反映数据变化的关系，其确切内涵是灰的。
对于每一种食品，均建立 模型，采用建立的模型进行未来 4 和 5 月的价格
预测，此时 。
模型的求解
对于每一种食品建立的 模型，分别去求 ，在求出时间响应式：
（13）
继续对于求出的 的模拟值，
由： （14）
对于 进行还原，再对于还原的结果采用相对误差进行误差检验，最后根据所求的时
间响应式 ，将 ，代入继续求得 的模拟
值 ，还原求出 的模拟值，
， 即可。
下面取第一种食品菜籽油进行举例分析，其它食品的计算均按照第一种计算即可：
第一步：先取菜籽油生成 序列，
（15）
再对 作 ，得:
（16）
图 5-4 菜籽油价格走势图
1,1GM
(0)X
(0) (0) (0) (0)( (1), (2), ( ))X x x x k L , 1,2,3 39k L
1,1GM
( 0 ) (1 )( ) ( )x k a z k b 1,2,3 39k L
(1)z (1)X a
b
1,1GM
40,41,42,43,44,45k
1,1GM
1,1GM ˆ , Ta a b
1 0 ( 1)ˆ 1 a kb bx k x e
a a
1X
(0) (1) (1)ˆ ˆ( ) ( )x k x k
0X
1 0 ( 1)ˆ 1 a kb bx k x e
a a
40,41,42,43,44,45k (1)X
(1) (1) (1) (1) (1) (1) (1)ˆ ˆ ˆ ˆ ˆ ˆ ˆ( (40), (41), (42) (43) (44) (45))X x x x x x x , , ,
0X
(0) (1) (1) (1) (1)ˆ ˆ ˆ ˆ( ) ( ) ( ) ( 1)x k x k x k x k 40,41,42,43,44,45k
0X
(0) (0) (0) (0) (0)( (1), (2), (3), (39))X x x x x L , (,,, ,) L
0X 1 AGO
(1) (1) (1) (1) (1)( (1), (2), (3) (39))X x x x x L, , (,,, , ) L
第二步：对 作准光滑性检验。由：
（17）
得： ，
当 时，准光滑条件满足。
第三步：检验 是否具有准指数规律。由
（18）
得： ，
当 时， ， ，准指数规律满足，故可对 建立
模型。
第四步：对 作紧邻均值生成。
令： （19）
得：
， （20）
第五步：对参数列 进行最小二乘估计。得：
（21）
第六步：确定模型：
（22）
及时间响应式
，将 代入得
（23）
第七步：求 的模拟值
（24）
第八步：还原求出 的模拟值。由式（14）可得：
第九步：检验误差
（25）
第十步：对 4、5 月份的价格进行预测
中，将 代入得到：
0X
(0)
(1)
( )
( )
( 1)
x k
k
x k
(3) < (4) <
2k
(1)X
(1)
(1)
(1)
( )
( )
( 1)
x k
k
x k
(1) (3) (1) (4)
2k (1) ( ) 1, (1)X
1,1GM
(1)X
(1) (1) (1)( ) ( ) ( 1)z k x k x k
(1) (1) (1)( (2), (39))Z z z L , (,, ,) L
(0)
(0)
(0)
(2)
(3)
(39)
x
x
Y
x
MM
(1)
(1)
(1)
(2) 1 1
1(3) 1
1(39) 1
z
z
B
z
M MM M
ˆ , Ta a b
1 ˆ ( )
T Ta B B B Y
(1)
(1)
dx
x
dt
1 0 ( 1)ˆ 1 a kb bx k x e
a a
ˆ
a
1 ( 1)ˆ k e
(1)X
(1) (1) (1) (1) (1)ˆ ˆ ˆ ˆ ˆ( (1), (2), (3) (39))X x x x x L, , (,,, , ) L
0X
(0) (0) (0) (0) (0)ˆ ˆ ˆ ˆ ˆ( (1), (2), (3) (39))X x x x x L, , (, , , ,) L
39
2
1
38 kk
1 ( 1)ˆ k e 40,41,42,43,44,45k
(1)ˆ (,,,,,)X
根据式（14）： ， 还原，
同时再减去误差积累值，求出在 4、5 月份的价格预测值如表 5-1：
表 5-1 菜籽油 4、5 月份的预测价格
预测价格
食品种类 计量单价
2011/4/5 2011/4/15 2011/4/25 2011/5/5 2011/5/15 2011/5/25
菜籽油 散装
元/500
克
从而利用 Matlab（源代码见附录二）求出的所有食品价格预测：
表 5-2 所有食品 4、5 月份的预测价格
预测价格
食品种类 计量单价
2011/4/5 2011/4/15 2011/4/25 2011/5/5 2011/5/15 2011/5/25
菜籽油 散装 元/500 克
大豆油 散装 元/500 克
花生油 一级桶装 元/5 升
大豆调
和油
一级桶装 元/5 升
鲜猪肉 精瘦肉 元/500 克
鲜猪肉 肋条肉 元/500 克
鲜牛肉 新鲜去骨 元/500 克
鲜羊肉 新鲜去骨 元/500 克
鲜羊肉 新鲜带骨 元/500 克
活鸡
活肉鸡
公斤
元/500 克
鸡肉
白条鸡、开
膛 上等
元/500 克
鸡蛋
新鲜完整
鸡场蛋
元/500 克
带鱼
冻 250 克
左右
元/500 克
草鱼
活 1000 克
左右一条
元/500 克
鲤鱼
活 500 克
以上一条
元/500 克
芹菜 新鲜一级 元/500 克
大白菜 新鲜一级 元/500 克
油菜 新鲜一级 元/500 克
黄瓜 新鲜一级 元/500 克
萝卜 新鲜一级 元/500 克
茄子 新鲜一级 元/500 克
西红柿 新鲜一级 元/500 克
土豆 新鲜一级 元/500 克
胡萝卜 新鲜一级 元/500 克
(0) (1) (1) (1) (1)ˆ ˆ ˆ ˆ( ) ( ) ( ) ( 1)x k x k x k x k 40,41,42,43,44,45k
青椒 新鲜一级 元/500 克
尖椒 新鲜一级 元/500 克
圆白菜 新鲜一级 元/500 克
豆角 新鲜一级 元/500 克
表 5-3 各食品相对误差分析表:
误差及日期
食品种类
2010/3/5 2010/3/15 2010/3/25 2010/4/5 2010/4/15 2010/4/25 2010/5/5
菜籽油 散装
大豆油 散装
花生油 一级桶装
苹果
红富士
一级
豆腐 无包装
酱油
当地主销
瓶装
500ml
醋
当地主销
瓶装
500ml
鲜牛奶
当地主销
(袋装)
问题三
基于城市居民食品零售价格情况给有关部门的建议：
尊敬的食品相关管理机构及部门领导：
经分析，可以看出 4、5 月蔬菜价格受春节节日效应消褪影响，整体呈回落态势，
市场供需情况较为平衡，蔬菜日交易量也较春节 2 月份期间稳定增加，一些大陆蔬菜品
种如小白菜、茄子、黄瓜、莲藕、西红柿等品种降幅都超过一成，4、5 月青椒价格涨幅
较大，主要是薄皮青椒产地收购价格高，青椒市场需求量又较大，因此批零终端都居高
不下。我们预计随着四月气温的稳定上升，蔬菜种植将逐渐走出“春淡传统”效应影响，
本地蔬菜上市量将增加，会拉动蔬菜市场价格呈小幅下降态势。
本月粮油市场价格基本保持稳中温和上行态势，其中粮食价格的上涨主要还是由于
国家 2、3 月份公布提高稻谷收购价的政策效应影响，进入月底，粮食批零价格上涨态
势放缓，市场逐步恢复稳定，但受国际市场粮食价格依然高位运行以及国内通胀预期较
高的影响，今后一段时间内，粮食价格仍然会保持缓和上行态势。食用油在进入传统消
费淡季后，仍然呈高位运行态势，主要由于当前国际油脂期货价格仍然居高不下，而我
国油脂压榨原料对外依存度较高，因此后期食用油价格下降可能性不大，基本在当前价
位小幅波动运行。
猪肉价格在经历了春节期间高位运行期后，4、5 月价格稳中回落，猪肉市场零售价
格跌幅高于生猪收购价跌幅，主要是当前猪肉消费量较春节期间有所下降，而生猪收购
价受饲料以及人工成本等因素影响，仍然保持高位坚挺运行态势，预计后期猪肉价格涨
跌幅不会太大，基本保持当前价位。
为保障 4、5 月份期间市场价格基本稳定，进一步规范市场价格秩序，营造良好的
节日消费环境，市物价局需出台一些政策，应当加强市场价格监管，维护好广大经营者
和消费者的合法权益。
根据我们研究，具体地来说，要做好 4、5 月份市场价格监管，首先要加强价格监
测，掌握市场动态。密切关注粮油、蔬菜、肉禽蛋奶、食用盐等居民基本生活必需品等
价格的变动情况，及时分析变动原因，把握变动趋势。针对节日消费的特点，结合当前
监管工作的重点，加强市场流通领域的价格行为监管。
今年以来，市场物价上涨的压力明显，消费者心理价格承受能力比较脆弱。价格是
国民经济的综合反映，需要多部门综合施策，地方政府在稳定物价上要有所作为。一是
要保障低收入群众的生活，在低收入群众的生活救助上反应要灵敏，措施不要滞后。二
是要发挥价格调节的作用，要充分运用价格调节基金平抑市场，稳控物价。
此致
敬礼！
某数学建模深资研究小组
2011 年 5 月 3 日
6 研究规范合理性的讨论
（1）时间序列分析预测模型是根据市场过去的变化趋势预测未来的发展，它的前
提是假定事物的过去会同样延续到未来。事物的现实是历史发展的结果，而事物的未来
又是现实的延伸，事物的过去和未来是有联系的。市场预测的时间序列分析法，正是根
据客观事物发展的这种连续规律性，运用过去的历史数据，通过统计分析，进一步推测
市场未来的发展趋势。市场预测中，事物的过去会同样延续到未来，其意思是说，市场
未来不会发生突然跳跃式变化，而是渐进变化的。
（2）通过少量的、不完全的信息，建立灰色微分预测模型，对事物发展规律作出
模糊性的长期描述，结合该灰色预测模型的原理特点，是可行的。
7 模型的优缺点分析
（1）针对预测建立的灰色模型，具有原始数据需求量小、对分布规律性要求不严、
预测精度较高等优点。
（2）时间序列预测法因突出时间序列而很少考虑外界因素影响，因而存在着预测
误差的缺陷。
（3）对于模型在预测价格时，可以采用残差 模型，来进一步降低相对误差，
对原来的模型进行修正，以提高精度。
(1,1)GM
参考文献
[1].金勇进，统计学，北京：中国人民出版社，2010。
[2].梅长林，周家良，实用统计方法，北京：科学出版社，2002。.
[3].何源等，SPSS 统计分析，北京：人民邮电出版社，2010。.
[4].周燕等，Matlab 在统计与工程数据分析中的应用，北京：电子工业出版社，2010。
[5].刘思峰，灰色系统与建模，开封：科学出版社，2004。周燕等，Matlab 在统计与工
附 录
一、时间序列 Matlab 源程序
clc;
avg=[];
sum=[];
for j=1:42
for k=1:6
for i=k:k+1
sum(j，k)=0;
sum(j，k)=sum(j，k)+(p10(j，i+1)-p10(j，i))/p10(j，i);
end
end
end
avg=sum./2;
fcast1=base;
for n=1:42
for g=1:6
fcast1(n，g+1)=fcast1(n，g)*(avg(n，g)+1);
end
end
二、单价预测 Matlab 源程序
clc;
x=[];
for i=1:42
x(i，1)=data(i，1);
end
for j=1:42
for i=2:39
x(j，i)=x(j，i-1)+data(j，i);
end
end
z=[];
for k=1:42
for i=1:38
z(k，i)=*(x(k，i)+x(k，i+1));
end
end
one=[];
for i=1:38
one(end+1)=1;
end
B=[];
Y=[];
a=[];
b=[];
for g=1:42
for j=1:38
Y(j，1)=data(g，1+j);
B(j，1)=(-1)*z(g，j);
B(j，2)=1;
end
A=inv((B'*B))*B'*Y;
a(end+1)=A(1，1);
b(end+1)=A(2，1);
end
X=[];
for j=1:42
for i=1:45
X(j，i)=(data(j，1)-b(j)/a(j))*exp(a(j)*(1-i))+b(j)/a(j);
end
end
XX=[];
for i=1:42
for n=2:45
XX(i，1)=X(i，1);
XX(i，n)=X(i，n)-X(i，n-1);
end
end
for i=1:45
XX(36，i)=;
end
avg=[];
for j=1:42
for m=1:39
avg(j，m)=data(j，m)-XX(j，m);
end
end
rate=[];
for i=1:42
for f=1:39
rate(i，f)=avg(i，f)/data(i，f);
end
end
avgrate=[];
arate=[];
for i=1:42
sam=0;
sum=0;
for u=1:39
sum=sum+abs(rate(i，u));
sam=sam+rate(i，u);
end
avgrate(i)=sum/39;
arate(i)=sam/39;
end
price=[];
for s=1:42
for q=40:45
price(s，q)=XX(s，q)*(arate(s)+1);
end
end
三、所有食品的相对误差值
食品种类 平均相对误差 食品种类 平均相对误差 食品种类 平均相对误差
菜籽
油
散装 鲤鱼
活 500 克
以上一条
蒜
苔
新鲜一级
大豆
油
散装 芹菜 新鲜一级
韭
菜
新鲜一级
花生
油
一级桶装
大白
菜
新鲜一级
芦
柑
一级
大豆
调和
油
一级桶装 油菜 新鲜一级
苹
果
红富士 一级
鲜猪
肉
精瘦肉 黄瓜 新鲜一级
香
蕉
国产 一级
鲜猪
肉
肋条肉 萝卜 新鲜一级
西
瓜
当地主销 一
级
鲜牛
肉
新鲜去骨 茄子 新鲜一级
豆
腐
无包装
鲜羊
肉
新鲜去骨
西红
柿
新鲜一级
食
用
盐
精致含碘盐
鲜羊
肉
新鲜带骨 土豆 新鲜一级
绵
白
糖
当地主销(袋
装)
活鸡
活肉鸡
公斤
胡萝
卜
新鲜一级
白
砂
糖
当地主销(袋
装)
鸡肉
白条鸡、开
膛 上等
青椒 新鲜一级
红
糖
当地主销(袋
装)
鸡蛋
新鲜完整
鸡场蛋
尖椒 新鲜一级
酱
油
当地主销 瓶
装 500ml
带鱼
冻 250 克左
右
圆白
菜
新鲜一级 醋
当地主销 瓶
装 500ml
草鱼
活 1000 克
左右一条
豆角 新鲜一级
鲜
牛
奶
当地主销(袋
装)
