当前位置:首页 >> 数学 >>

预测模型与案例










最近几年,在全国大学生数学建模竞赛常常出现预测模型或是 与预测有关的题目,例如疾病的传播,雨量的预报等。什么是预测模 型?如何预测?有那些方法?对此下面作些介绍。 预测作为一种探索未来的活动早在古代已经出现,但作为一门 科学的预测学,是在科学技术高度发达的当今才产生的。 “预测”是 来自古希腊的术语。我国也有两句古语: “凡事预则立,不预则废” , “人无远虑,必有近忧” 。卜卦、算命都是一种预测。中国古代著名 著作 “易经” 就是一种专门研究预测的书, 现在研究易经的人也不少。 古代的预测主要靠预言家,即先知们的直观判断,或是借助于某些先 兆,缺乏科学根据。预测技术的发展源于社会的需求和实践。20 世 纪初期风行一时的巴布生图表就是早期的市场预测资料, 哈佛大学的 每月指数图表为商品市场、证券市场和货币市场预测提供了依据。然 而这些预测都未能揭示 1929-1930 年经济危期的突然暴发,使工商 界深感失望。尔后,经济学家们从挫折中吸取了教训,采用趋势和循 环技术对商业进行分析和预测,科学预测也因此开始萌生。20 世纪 30 年代凯思斯提出政府干预和市场机制相结合的经济模型,1937 年 诺依曼又提出了扩展经济模型,对近代经济模型产生重要的影响,科 学的经济和商业预测也就步入发展阶段。 技术预测开始于二次世界大战后的 20 世纪 40 年代,直到 20 世 纪 50 年代未才广泛应用于工农业和军事部门。由于社会、科学技术
1

和经济的大量需求,预测技求才成为一门真正的科学,预测未来是当 代科学的重要任务。 20 世纪以来,预测技术所以得以长足进步,一方面,与社会需 求有很大关系,另一方面通过社会实践和长期历史验证,表明事物的 发展是可以预测的。而且借助可靠的数据和科学的方法,以及预测技 术人员的努力,预测结果的可靠性和准确性可以达到很高的程度,这 也是预测技术迅速发展的另一个重要原因。 科学技术、经济和社会预测的应验率也是很高的。维聂尔曾预言 20 世纪是电子时代,法国思想家迈希尔 18 世纪末到 19 世纪初对巴 黎未来几百年的发展进行了预测。从 1950 年的实际情况分析,他的 预测中有 36%得到证实,28%接近实现,只有 36%是错误的。法国 哲学家和数学家冠道塞在法国大革命时期曾采用外推法进行了一系 列社会预测,其中 75%得到证实。沙杰尔莱特 1901 年在《二十世纪 的发明》 一书中的一些预测, 其中 64%得到证实。 凯木弗尔特在 1910 年和 1915 年公布的 25 项预测中,到 1941 年只有 3 项未被证实,3 项是错误的。我国明朝开国功臣刘基就预测将来是天上铁鸟飞,地上 铁马跑,那时还没有火车、飞机。 预测的目的在于认识自然和社会发展规律, 以及在不同历史条件 下各种规律的相互作用,揭示事物发展的方向和趋势,分析事物发展 的途径和条件,使人们尽早地预知未来的状况和将要发生的事情,并 能动地控制其发展,使其为人类和社会进步服务。因而预测是决策的 重要的前期工作。决策是指导未来的,未来既是决策的依据,又是决
2

策的对象,研究未来和预测未来是实现决策科学化的重要前提。预测 和决策是过程的两个方面,预测为决策提供依据,而预测的目的是为 决策服务,所以不能把预测模型和决策模型截然分开,有时也把预测 模型称为决策模型。



预测的前期准备工作

为保证预测结果的精确度,预测之前必须做一系列的准备工作:

(一)数据的准备
数据是预测工作的前提和重要依据,预测不能是臆造和空想,任 何事物的发展都有一定的规律, 认真研究预测对象并充分考察预测对 象所处的环境,以系统分析的方法对过去和现在的数据进行总结,从 中找出规律,便可科学地推断未来。 数据在预测中主要有两个作用: (1) 、用于确定由某些历史观察点 组成的行为模型; (2) 、在因果模型预测中确定自变量的未来值。 预测的初始阶段,首先是从事数据的收集、整理、加工和分析, 为建模创造良好的条件。

(Ⅰ)数据的收集和整理
按时态分,数据可分为历史数据和现实数据;按预测对象分,可 分为内部数据和外部数据;就收集的手段分,可分为第一手数据和第 二手数据。 第一手数据,包括以各种形式初次收集的数据。收集第一手数据 的途径包括:抽样调查,连续调查,或全面调查。在预测的定性方法 中常常需要第一手数据, 例如特尔斐法的第一个阶段就是收集第一手
3

数据。由于获取第一手数据的费用较高,时间较长,所以定量方法常 采用第二手数据。 第二手数据多为已经公布和发表的资料,易于获取,代价低,数 据精度也有一定的保证。其缺点是数据可能不能直接适用于预测情 况。因此,常常需要对已公布的数据进行修正和处理,使其适应于预 测需要。 无论是第一手数据还是第二手数据,都可能是混乱的、无序的、 彼此间孤立的。预测人员都应将原始数据按“单元”或“类别”整理 和集中,以便使其成为内容上完整、有序、系统,形式上简明统一的 数据。

(Ⅱ)数据的分析和处理
建模不仅需要大量的数据,同时数据必须可靠,并适合建模的要 求。这些数据虽然是历史的客观写照,但有可能是失真的数据。对于 失真的数据, 以及不符合建模的数据, 必须通过分析, 加以适当处理。

1.处理的原则
(1)

准确,处理后的数据能正确反映事物发展的未来趋势和状

况; (2)及时,数据的处理要及时; (3)适用,处理的数据能满足建模的需要; (4)经济,要尽量减少数据处理的费用,以降低预测成本; (5) 一致,指处理的数据在整个使用期间内必须是一致的,具 有可比较性。
4

2.处理方法
(1)判别法
通过对历史数据的判断,选择其中可代表整个预测过程中很可能 发生的模式的数据作为建模数据;

(2)剔除法
如果数据量比较大,且非必须具备连续的数据量,这时可剔除数 据中受随机干扰的异常值;

(3)平均值法
在数据比较少或需要连续数据时,则可采取平均值法对数据进行 处理。 对于时间序列数据,可用异常值前后两期数据的算术平均值或几 何平均值对异常值进行修正,即

xt ?1 ? xt ?1 2 xt ? xt ?1 ?xt ?1 xt ?



通常当历史数据的发展趋势呈线性时,取算求平均值,当发展趋 势呈非线性时,取几何平均值。 在利用因果关系建立数学模型时,为去掉偶然因素对建立模型的 影响,可采用下面的计算方法对统计数据中的异常数据加以修正: 当 x 与 y 之间为线性因果关系时,取

yk ?

yl ?xl ? ym ?xm 2 xk
5

当 x 与 y 之间为非线性因果关系时,取

yk ?
式中

yl ?xl ?ym ?xm xk

yk 为有随机因素影响时期因变量的估计值,xk 是与之对应的自

变量; xl , xm 是与 xk 在数值上相差最小的两个自变量,且

xl ? xk ? xm
yl , ym 分别是与 xl , xm 相对应的因变量统计值
(4)拉平法
由于条件发生变化,常常使一些厉史数据不能反映现时的情况, 例如,大型钢铁厂、化肥厂、或油气田的建成投产或开发,可以使产 量猛增,这时历史数据将发生突变,出现一个转折,如用这类数据建 模,则需要处理。这时拉平法是一种较好的方法。它的原理是对转折 点前的数据加一个适当的量值,使其与折点后的数据走向一致。

(5)比例法
销售条件与环境的变化常常会引起一个企业产品市场销售比例的 改变。当比例变化较大时,说明销售条件与环境对销售的影响己超过 其他因素对销售的影响, 也说明以前的销售统计数据所体现出的销售 发展规律不再适用之于目前的情况了。 如果仍然利用这些数据建立预 测模型,将无法体现销售条件和环境变化后的销售量变化的规律,用 这样的模型进行预测,将会造成较大的误差。因此,如果还想利用这 些数据建立模型,进行预测,就应该把它们处理成能体现条件与环境 发生变化之后的情况的数据。对于这类数据,比例法就是一种比较有
6

效的处理方法。 例如,某一生产生产资料的大型企业,80 年代中期前销售额一直 呈递增趋势,而 80 年代中期后,受压缩基建规模的影响,销售量突 然下降。又如轿车在 80 年代中期以前一直是紧俏商品,后因国家实 行控购政策,销售量一度急剧下降。这时,对上述某一生产资料销售 量或对轿车销售量进行预测,都要考虑政策因素的影响,对于前期数 据采用比例法进行适当修正(当时是计划经济,私人买不起轿车。买 轿车的都国家机关、企事业单位。 ) 当然比例法不仅仅限于对数值向下调,也适合向上调。比例法数 据处理公式为

yt ?i ? yt ?i

ut 其中: ut ?i

yt ?i ? t ? i年修正后的数 yt ?i ? t ? i年实际数据 ut ? t年的市场占有率 ut ?i ? t ? i年的市场占有率
(6)移动平均和指数平滑法
如果原始数据总体走向具有一定规律性,但因受随机因素干扰, 数据离散度很大, 采用平均值法也难以处理。 这时可采用一次、 二次、 甚至三次移动平均和指数平滑对数据进行平滑,用平滑的数据建模。 在分解预测时,为处理季节数据,则必须采用高次幂的移动平均 法,对数据平滑。
7

(7)差分法
有些模型,例如鲍克斯-詹金斯模型只能处理平稳数据,如果原始 数据为非平稳数据,则需釆取差分处理。差分有三种主要类型:前向 差分、后向差分、中心差分。 前向差分:在处理时间数列时,一阶前向差分定义为

x ? xt ?1 ? xt
' t
一阶前向差分是当时间由 t 变到 t+1 时, 二阶前向差分定义为

xt 的改变量。

xt'' ? xt' ?1 ? xt? ? xt ?2 ? 2xt ?1 ? xt
同样,可以定义高阶差分。 后向差分:在处理时间数列时,一阶后向差分定义为

xt' ? x t ? xt ?1
一阶后向差分是当时间由 t 递推到 t-1 时, xt 的改变量。 二阶后向差分定义为

xt'' ? xt' ? xt' ?1 ? xt ? 2xt ?1 ? xt ?2
同理可以定义高阶后向差分 中心差分:在处理时间数列时,一阶中心差分定义为

xt' ? xt ? 1 ? xt ? 1
2

2

二阶中心差分定义为

8

xt'' ? xt' ? 1 ? xt' ? 1 ? xt ?1 ? 2 xt ? xt ?1
2 2

同理可以定义高阶中心差分。 在处理时间数列时,主要应用后向差分。一次多项式数据通过一 阶差分就可转换为平稳数据, 二次多项式和三次多项式数据分别通过 二阶和三阶差分可转换为平稳数据, 而三次以上的高次多项式在应用 中很少采用。

(Ⅲ)数据的内涵及数量
在预测过程中,由于预测对象不同,预测内容不同,以及预测期 限不同,所需的数据内涵及数量也不同。

经济预测的数据主要包括:
(1)国民经济总产值及各部类的分配情况; (2)各行业的生产规模和生产能力以及技术水平; (3)政府的经济政策及产业政策; (4)生产力布局; (5)人口发展趋势及就业情况; (6)国民经济投资及分配; (7)国际环境及变化趋势。

市场需求预测需要的数据主要有:
(1)人口及人均收入; (2)国民收入的增长及分配情况; ; (3)与产品消费直接有关的政府政策和法规,如进口限制、进口 税、销售稅和其它税费、信贷管理及外费管理等。
9

(4)一段时期内产量和产值的生产能力; (5)一段时期内的产品的进口量; (6)代用品或近似代用品的产量和进口量; (7)与有关新投入的产品前后关联度高的产品的产量; (8)国家计划规定的产品或代用品的生产指标; (9)产品出口量; (10)个人或集体消费者们的实贯或嗜好; (11)法律方面的资料。

二 专家的选择和专家组的组成
在现实生活中,有时不得不在不确定的条件下作出决策,这是因 为或者决策的制约因素过多,或者其中某些因素无法度量。我们常称 之为定性因素。为这类决策提供预测,因为没有严格的理论依据,定 量方法无法采用。在这种情况下,借助专家的经验判断则有可能作出 定量方法难以得到的科学预测。专家的素质取决于他的知识、经验、 智慧和对未来的预测能力,以及其他一些因素。实践表明,在当今如 此复杂多变的情况下,任何个人或一个专家都难于作出较精确的预 测。必须集中多方专家的意见才能作出科学的预测。因此选择专家组 成员是预测能否成功的重要环节,是预测要做的首要工作。应邀的专 家要具有广泛的知识,对预测所涉及主题的各领域应有较深的造诣。 选择专家不能简单从事, 不能事先未经征得同意就将调查表发给拟邀 请的专家。因为有的专家可能不愿意参加这项预测。那么选择专家应 如何进行呢?
10

(一)什么叫专家
在组织专家预测时,专家是个广义的概念,拟选的专家不能仅仅 局限于一个领域的权威,因为权威人数是有限的。特尔斐法拟选的专 家是指在该领域从事 10 年以上工作的专业干部。

(二)怎样选择专家
怎样选择专家是由预测任务决定的。如果要求比较深入地了解部 门的历史情况和技术政策,或涉及到本部门的机密问题,则最好从本 部门选择专家。从本部门选择专家比较简单,既有档可查,又熟悉干 部的现实情况。如果预测任务仅仅关系到事物的发展,则最好同时从 部门内外挑选。从外部选择专家,大体按以下顺序进行: (1)编制征求专家应答问题一览表; (2)根据预测问题,编制所需专家类型一览表; (3) 将问题一览表发给每个专家, 询问他们能否坚持参加规定问 题的预测。 (4)确定每个专家从事预测所消耗的时间和经费。 从外部选择专家比较困难,一般要经过几轮。首先要收集本部门 职工比较熟悉的专家名单, 而后再从有关期刊和出版物中物色一批知 名专家。以这两部分专家为基础,将调查表发给他们,征求意见,同 时要求他们再推荐 1-2 名有关专家。预测领导小组从推荐的专家名 单中,再选择一批有 2 人以上推荐的专家。

(三)选择什么样的专家
在选专家的过程中,不仅要注意选择精通技术、有一定名望、有
11

学科代表性的专家,同时还需要选择相关学科、边缘学科、社会学和 经济学等方面的专家。选择承担领导职务的专家固然重要,但要考虑 他们是否有足够的时间认真填写调查表。经验表明,一个身居要职的 专家匆忙填写的调查表, 其参考价值还不如一个专事某项工作的一般 专家认真填写的调查表。再者,乐于承担任务,并坚持始终,也是选 择专家时要注意的。

(四)专家组人数
预测小组人数视预测向题规模而定。人数太少,限制学科代表性, 并缺乏权威;人数太多,难于组织,对结果处理也比较复杂。预测的 精度与人数的函数关系是,当人数较少时,随着人数的增加预测精度 很快提高。但人数接近 15 时,进一步增加人数对预测精度影响不大。 小组人数一般以 15-50 人为宜。当然对于一些重大问题,专家人数 也可扩大到 100 名以上。在确定专家人数时,值得注意的是,有的专 家即使同意参加预测,因为种种原因也不见得每轮必答,有时甚至中 途退出,因而预选人数要多于规定人数。

12

预 测 标 准 误 差

专家人数与预测误差的关 系

15

16

17

18

19 人数

定性预测方法尽管有时并不需要外界输入数据, 即使有数据要求, 精度要求也不严格, 但是这并不意味着定性预测方法的精度不如定量 方法。有时定性预测方法所得结果,其精度还高于定量方法。这是因 为每一个专家都是一个数据库,都存贮着大量与预测有关的数据,而 其中相当部分还是社会未发表的数据。 预测小组就可利用专家提供的 数据创造一本脚本,用来描述过去发生了什么事情,未来将发生什么 事情。脚本不仅可以真实地反映一组完整的描述真实事件的数据,同 时这组数据可以同传统形式经常采用的定量变量数据媲美。



预测的数学准备

在预测过程中需要很多数学知识,主要有微分方程、概率与数理 统计、线性规划和非线性规划等等。但使用最多的是统计学的相关知 识:常用的统计量、参数的估算、假设检验、区间估计等。这些我们 就不做介绍了。



实用预测方法
13

(一)定性预测方法
预测方法很多,多达 200 多种,但常用的不过 30 多种,最常用的 只有 10 多种。预测方法的分类没有统一的标准和体系。前苏联的专 家把预测方法分为两类:启发式预测(专家预测)和数学模型预测。 而美国有的专家把预测分为定性方法和定量方法, 有的专家把预测方 法分为定性预测、定量预测、定时预测、概率预测四类。我国多把预 测分为定性和定量两种。下面是我国目前常用几种预测方法: 定性预测方法:主要有特尔斐法、目标预测法; 定量预测方法:主要有时间序列模型,因果关系模型。 而时间序列模型包含移动平均法、指数平滑法、分解预测法、鲍 克斯-詹金斯模型。 因果关系模型包含趋势外推法、回归分析法、数量经济模型、投 入产出模型、灰色模型、系统模型。 每种方法都有它的适用范围和特点,预测程序,预测模型。 下面重点介绍使用最多、应用最广的特尔斐法,至于其它各种方 法请大家自行查看相关的书籍。

特尔斐法
特尔斐法是在专家会议预测法的基础上发展起来的,由美国兰德 公司于 1964 年发明并首先用于技术预测。 专家会议法虽然可以通过会议使专家之间广泛交流意见,互相启 发, 为重大决策提供预测依据, 但专家会议法也有三个重大缺点, 即: (1)易于屈服于权威或多数人的意见;
14

(2)易受劝说性意见的影响; (3)会出现因自尊心影响而不愿公开修正已发表的、然而是不完 全正确、甚至是错误的意见。这就使专家会议作出的预测有时是片面 的,甚至有可能是错误的。特尔斐法克服了以上缺点,它是将所要预 测的问题以信函的方式寄给专家,将回函的意见综合、整理,又匿名 反馈给专家征求意见,如此反复多次,最后得出预测结果。

(一)

特点及适用范围

特尔斐法有三个特点: 1. 匿名性 测 由于特尔斐法采用匿名函询征求意见,应邀参加预

15

? ?特尔斐法 ?定性预测方法 ? ?目杯预测法 ? ? ? ? ?一次移动平均 ? 移动平均法 ? ? ? ? ?二次移动平均 ? ? ? ? ? ? ?一次指数平滑 ? ? ? ? ? ?指数平滑法 ?二次指数平滑 ? ? ? ?三次指数平滑 ? ? ? ? ? ?时间序列模型 ? ? ?分解预测法 ? ? ? ? ?自回归模型 ? ? ? ? ? ? ? ?移动平均模型 ? 鲍克斯 詹金斯模型 ? ? ? ? ? ?鲍克斯-詹金斯模型 ? ? ? ?季节性ARIMA模型 ? ? ? ? ? ? ? ? ?多项式模型 ? ? ? ? ? ? ? ?趋势外推 ?指数模型 ? ?定量预测方法 ? ? ? ?生长曲线 ? ? ? ?包络曲线 ? ? ? ? ? ? ? ?一元线性回归 ? ? ? ? ?回归分析 ? ? ?多元线性回归 ? ? ? ?非线性回归 ? ? ? ? ? ? ? 因果关系模型 ?数量经济模型 ? ? ? ?投入产出模型 ? ? ? ? ? ?灰色系统模型 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?

预测方法分类表 的专家互不相见,可消除心理因素的影响,专家可参照前一轮预 测结果修改自己的意见,而元需作公开说明。 2. 轮间反馈可沟通性 特尔斐法一般要经过四轮,每一轮的汇总意见又匿名反馈给专家,便
16

于互相沟通和启发。 3. 预测结果的统计特性 特尔斐法采用统计方法对结果进行定量处理, 能科学地综合专家们的 预测意见。 特尔斐法是传统定性分析的一个飞跃, 它突破了单纯的定性或定 量分析的界限,为科学、合理地制定决策开阔了思路。由于它能够对 未来发展中可能出现的前景作出概率描述, 因而为决策者提供了多方 案选择的可能性。 采用特尔斐法不仅可以从事技术预测,同时可以从事经济、社会 预测;不仅可以从事短期预测,同时可以从事长期预测;不仅可以预 测事物的量变过程, 同时可以预测事物的质变过程。 因而近几十年来, 特尔斐法已经成为一种广为应用的预测方法。 在长远规划者和决策者 心目中,特尔斐法享有很高威望,并逐渐成为一种重要的规划决策工 具。

(二)预测程序
应用特尔斐法进行预测,主要包括四个阶段: (1)建立预测领导 小组,编制预测日程计划; (2)选择专家; (3)轮间反馈; (4)编写 预测报告。 由于预测结果的准确程度在很大的程度上依赖于专家的知 识广度、深度和经验,因此,如何选择专家是很重要的。具体的预测 程序见下图

17

建立预测领导小组 编制预测日程计划

选择专家









编写预测报告

确定预测 主题

确定预测 事件

设 计 调 查 表,从事第 一轮预测

反馈第一轮 结果,进行 第二轮预测 实现时间的回 答分布:中位 数,上下四分 点

反馈第二轮 结果,进行 第三轮预测 实现时间的回 答分布:中位 数,上下四分 点 回答结界的统 计处理 分析持异端意 见专家的意见 修正脚本

预测结果报告

选定目标

企业环境的调 查(现状和未 来)

需要可行性 实现的时间

修正脚本,编 写企业的未来 趋势脚本

选定期限范围

编制脚本 (1)现状外推 脚本 (2)悲观脚本 (3)乐观脚本 趋势外推法、 目标树法、专

实现的附带 条件 对企业的冲 击程度(重 要性、紧迫 性) 应答者的权 威程度 数据处理

实现的概率及 相对重要性的 统计数据 根据第二轮的 统计数据,制 订第三轮调查 表 研究持异端意 见专家的意见 及其提出的理 由

目标具体化

今后 10 年新 兴领域、新 产品和新技 术预测

利分析法 脚本:10 年 后的企业

特尔斐法预测程序

(三)预测模型
应用特尔斐法需要作归纳、整理等很多工作,无需建立复杂的数 学模型。 在采用特尔斐法进行时间预测时, 一般用中位数代表专家集中意 见,用上下四分点代表专家意见的离散程度。中位数受项目多少的影 响,如果将专家的预测结果在水平轴上按时间的先后顺序排列,则位
18

居中央将全变量分为二等分的年份为中位数。变量的项数为奇数时, 第(n+1)/2 项为中位数。项数为偶数时,位居中央两项的平均数为 中位数。 计算中位数的公式为

? X n?1 ? n为奇数? ? 2 X中 ? ? X n ? X n ?1 / 2 ? n为偶数? ? 2 ? 2

?

?

其中: n 是专家预测的数据个数(该数列是按从小到大顺序排列 的) 用上下四分点表示预测区间时,公式为

X上

X下

? X m?1 ? m为奇数 ? ? 2 ?? X m ? X m ?1 ? / 2 ? m为偶数 ? ? ?? 2 2 ? X 3 m?3 ? m为奇数 ? ? 2 ? X 3 m ?1 ? X 3 m ? 2)/2(m为偶数) ( ? 2 2 ?? ? m为奇数 ? ?X 3m+1 2 ? X 3 m ? X 3 m ?1 ? / 2 ? m为偶数 ? ? ?? 2 2
当 n 为偶数时,m=n/2

上式中:当 n 为奇数时,m=(n-1)/2

人们常常用组合距,即最大预测值与最小预测值之差表示预测值 的变化幅度,而多数用上下四分点的间距表示预测值的变化幅度。 例 1 某单位釆用特尔斐法预测我国“九五”期间的轿车需求量, 经函询,20 位专家的预测值如下表(1) ,试分析该预测结果。
19

1 50 11 60

2 65 12 53

3 53 13 56

4 56 14 54


5 64.5 15 61
(1)

6 50 16 59

7 56 17 59

8 58 18 60

9 59 19 50

10 62 20 51

首先,将专家预测结果按从小到大顺序排列: A(1)=50 A(2)=50 A(3)=50 A(4)=51 A(5)=53 A(6)=53 A(7)=54 A(8)=56 A(9)=56 A(10)=56 A(11)=58 A(12)=59 A(13)=59 A(14)=59 A(15)=60 A(16)=60 A(17)=61 A(18)=62 A(19)=64.5 A(20)=65 n=20 为偶数,则

A中 ? m?

An ? An ?1
2 2

2

?

A10 ? A11 56 ? 58 ? ? 57 2 2

n ? 10为偶数,则 2 Am ? Am ?1 A ? A 53 ? 53 6 2 A上 ? 2 ? 5 ? ? 53 2 2 2 A3 m ? A3 m ?1 A ? A 60 ? 60 16 2 A下 ? 2 ? 15 ? ? 60 2 2 2

组合距 ? Amax ? Amin ? 65 ? 50 ? 15 上下四分点间距 ? A下 ? A上 ? 60 ? 53 ? 7
多方案相对重要性预测近年来应用较广,其专家集中意见用算术 平均值表示,公式为

20

1 Mj ? mj

?C
j ?1

mj

ij

式中:M j-j方案的算术平均值; m j-参加j方案评价的专家数 Cij-i专家对j方案的评分值
专家意见的离散程度用变异系数表示,公式为

Vj ?
式中:

Sj Mj

V j-j方案的变异系数; S j-j方案的标准差。
例 2 用特尔斐法聘请五位专家对 4 种科研方案进行优选排序, 评价结果如下表(2) 科研方案优选预测值 方案 专家 1 2 3 4 5 A B C D

分值 等级 分值 等级 分值 等级 分值 等级 50 60 50 70 90 2 3 2 3 1 100 90 100 80 70
表(2)

1 1 1 2 2

30 70 20 90 50

3 2 4 1 3

20 10 40 20 10

4 4 3 4 4

21

用统计方法分析的意见。 (1) 首先求各方案分数的算术平均值。

1 5 1 M A ? ? Cij ? ? 50 ? 60 ? 50 ? 70 ? 90 ? ? 64 5 j ?1 5 1 ?100 ? 90 ? 100 ? 80 ? 70 ? ? 88 5 1 M C ? ? 30 ? 70 ? 20 ? 90 ? 50 ? ? 52 5 1 M D ? ? 20 ? 10 ? 40 ? 20 ? 10 ? ? 20 5 MB ?
根据各方案的平均值, 可得到按重要程度排序的结论为 B、 A、 C、 D。 (2) 计算各方案专家意见的变异系数。为求得各方案变异系数, 需先计算各方案的方差及标准差,以反映专家意见的离散程度:

1 5 S ? ? ? CiA ? M A ? 5 ? 1 i ?1 1 2 2 2 2 2 ? ?? 50 ? 64 ? ? ? 60 ? 64 ? ? ? 50 ? 64 ? ? ? 70 ? 64 ? ? ? 90 ? 64 ? ? ? 200 ? 4?
2 A 2 SA ? SA ? 14.14

2

同理

1 2 2 2 2 2 2 SB ? ??100 ? 88? ? ? 90 ? 88? ? ?100 ? 88? ? ?80 ? 88? ? 70 ? 88? ? ? 4? ? 170 S B ? 13.04
1? 2 2 2 2 2 30 ? 52 ? ? ? 70 ? 52 ? ? ? 20 ? 52 ? ? ? 90 ? 52 ? ? 50 ? 52 ? ? ? ? 4? ? 820
2 SC ?

SC ? 28.64
22

1 2 2 2 2 2 2 SD ? ?? 20 ? 20 ? ? ? 20 ? 10 ? ? ? 20 ? 40 ? ? ? 20 ? 20 ? ? 20 ? 10 ? ? ? 4? ? 150 S D ? 12.25
各方案的变异系数为

VA ? S A / M A ? 14.14 / 64 ? 0.22 VB ? S B / M B ? 13.04 / 88 ? 0.15 VC ? SC / M C ? 28.64 / 52 ? 0.55 VD ? S D / M D ? 12.25 / 20 ? 0.25
可见专家意见对 B 方案协调程度最高,其余依次为 A、D、C 方 案。 其次常用的定性方法还有目标预测法。目标预测法又称规范性预 测法,是美国霍尼维尔公司首先开发并投入使用的,具有重要实用价 值,但目前在我固应用的示例还很少。 目标预测法的突出特点是:它不是探索在什么时间将达到什么目 标,而是在目标已定的情况下,研究如何实现既定目标。美国霍尼维 尔公司最早成功地利用这种方法建立了一个 PATTERN 模型,用于研 究阿波罗登月课题。阿波罗登月日期,美国政府在考虑多种因素情况 下,经过多次预测和评估己经确定。所以 PATTERN 模型的任务不是 探索阿波罗计划能否实现或何时实现, 而是如何实现, 以便按期登月。 这种方法对我国尤为实用。这是因为,很多领域的宏观或中观目 标,都是由各级主管部门确定的,例如党的十二大确定了 2000 年国 民生产总值比 1980 年翻两番的宏伟目标,同时,粮食、钢材、原煤 产量和发电量也都有了明确目标。有的微观目标也是主管部门决定
23

的,如开发什么产品,获得多少收益等。在这种情况下,探索实现目 标的最佳途径是极为重要的。由于目标预测法是研究实现目标的步 骤、措施,因而目标预测也可称为目标决策,可以说是建立决策实施 保证体系的一种重要方法。由于时间关系,这里就不做介绍了。 (二)定量预测法 定量预测法主要有时间序列模型、因果关系模型。

时间序列模型
时间序列模型主要研究事物的自身发展规律,借以预测事物的未 来趋势。主要方法有移动平均、指数平滑、分解预测、鲍克斯--詹金 斯模型、多变量模型以及历史类推法等。 一. 特点和应用范围

时间序列一般指一组按时间顺序排列的数据, 展示了研究对象 在一定时期的发展变化过程。时间序列模型,就是根据预测对 象时间序列的变化特征,研究事物自身的发展规律,探讨未来 发展趋势,是一种重要的定量预测方法,包括多种模型,主要 适用于经济预测、商业预测、需求预测、库存预测等,预测期 限主要为中、短期,不适用于有拐点的长期预测。 二 . 预测程序

此处介绍的预测程序为定量方法预测程序,既适用于时间序列模 型,又适用于因果关系模型,详见下图。

24

第一阶段:设计

预测开始

明确问题 可能影响预测对象的因素 预测对象的背景、环境

备选方案 无价值继续修改

数据源的分析 内部的 外部的 列出预测方案

第二阶段: 祥细说明 选择新变量 选择一种预测方法 有关变量 函数形式 初步选择预测方案

收集、核实、准备和分析数据

时间序列模型

选择预测模型

因果关系模型

25

估计模型参数

选择新模 型

验证模型合理 性

评价合理性的标 准

第三阶段: 评价

把前段结果文件化

产生模型预测值和置信区间

对预测值作主观判断

向管理部门 未被接受 介绍预测值 被接受 批准预测值

评价预测值的标准

把前面的结果文件化

监视预测的精确度 并确定显著差别的原因

否 差别是否显著 评价差别显著的标准



26

需要时及时修改预测值



运用的参数 是否改变


返回到预测开始

定量预测程序图

三 .

预测方法及数学模型

(一)移动平均值模型 移动平均法是一种最简单的适应模型,是在算术平均的基础上发 展起来的一种预测方法。 算术平均虽能代表一组数据的平均水平,但它不能反映数据的变 化趋势,而原始数据虽然存在某种趋势,但数据可能是零散的或杂乱 无章的,无法直接加以分析。移动平均法克服了上述弱点,其基本方 法是,选一个固定的周期数,对数据进行平均,每递推一个周期就加 上后一个数据,舍去初始数据,依次类推,直至把数据处理完毕。以 n=5 为例:
?? M5 ?
1

?1? M6

Y1 ? Y2 ? Y3 ? Y4 ? Y5 5 Y ? Y3 ? Y4 ? Y5 ? Y6 ? 2 5
27

?? ?? M5 、M6 表示第五、第六个周期的一次移动平均值,依次类
1 1

推。若移动平均的周期为 N,则可得到计算移动平均值的一般公式:

M t?1? ?

Yt ? Yt ?1 ? … ? Yt ? N ?1 N
?1?

其中, M t 表示第 t 期的一次移动平均值 可见,移动平均法实际上是对于某一 t 期数据,取前 N 个数据进 行平均,N 个数权数相同,而其它数据的权数等于零。这样,经过移 动平均, 将消除数据列中异常的因素, 对数据进行修匀。 一般情况下, 如果数据没有明显的周期变化和趋势变化,可用第 t 期的一次移动平 均值作为 t+1 期的预测值,即

? t ?1 ? M ?1? Y t
其一般公式为
1? M t?1? ? M t?? 1?

Yt ? Yt ? N N

表(3)中的第一列和第二列,即是原始数据与一次移动平均值的对 比。如取 N=3 的 3 期移动平均,则第三期数据的移动值为 5766.33, 是由(5600+5796+5930)/3 得到的。如用于预测,它可以作为第 4 期的预测值。 在一次移动平均值的基础上,应用移动平均的原理,还可以进行 二次甚至多次的移动平均,二次移动平均,就是以一次移动平均值为 原始数据,再进行一次移动平均,仍以 N=5 为例:其公式为

M9

? 2?

?1? ?1? ?1? ?1? ?1? M5 ? M6 ? M7 ? M8 ? M9 ? 5

28

? 2? 式中: M 9 表示第 9 期的二次移动平均值,其一般公式为

Mt

? 2?

1? M t?1? ? M t?? N ? M t ?1 ? N

? 2?

二次移动平均使原始数据得到了进一步修匀, 使其显现线性趋势。 表(3)中的第三列数据为 N=3 的二次移动平均值。

序 列 1 2 3 4 5 6 7 8 9 10 11

原 始 数 据 5600 5796 5930 6092 6257 6567 6851 7141 7436 7738 8045

一次移动平均值 M ?1?

二次移动平均值 M ? 2?

5775.3 5939.3 6093.0 6305.3 6558.3 6853.0 7142.7 7438.3 7739.7

5935.9 6112.5 6318.9 6572.2 6851.3 7144.7 7440.2

移动平均值

表(3)

在二次移动平均值的基础上,可建立线性模型:

? t ?? ? a ? b ? Y t t
式中:τ —预测超前期数 通过查表(多项式模型参数估算公式)可知:

at ? 2M t? ? ? M t?
1

2?

bt ?

2 M t?1? ? M t? 2? N ?1

?

?

对于表(3)中的数据,如以 11 期数据预测 12 期值,当取 N=3 时, 则有:
29

?1? ? 2? a11 ? 2M11 ? M11 ? 2 ? 7739.7 ? 7440.2 ? 8039.2

b11 ?

2 ?1? ? 2? ( M11 ? M11 ) ? 7739.7 ? 7440.2 ? 299.5 3 ?1

预测方程为

? 11?? ? 8039.2 ? 299.5? Y ? 12 ? 809.2 ? 299.5 ?1 ? 8338.7 Y
使用移动平均法,最重要的是移动周期 N 的选择。因为

S ?

2 a

S ?? ? N
2 a ——移动平均值方差
2

2

式中: S

S ?? ? ——原始数据点方差
N——数据点数 也就是说,移动平均修匀后的方差,随着 N 的加大而减少。也就是 N 越大,对原始数据修匀能力越强。下表(4)数据可清楚反映这一规律。
(1) 月 份 1 2 3 4 5 6 7 8 9 10 11 12 (2) 期 数 1 2 3 4 5 6 7 8 9 10 11 12 (3) 实际销售额 200.0 135.0 195.0 197.5 310.0 175.0 155.0 130.0 220.0 277.0 235.0 —— (4) 三个月移动平均值 —— —— —— 176.7 175.8 234.2 227.5 213.3 153.3 168.3 209.2 244.2 (5) 五个月移动平均值 —— —— —— —— —— 207.5 202.5 206.5 193.5 198.0 191.4 203.5

30

(某日用品电器销售额的移动平均预测)表(4) 然而修匀能力与对外界变化的反映速度是互相矛盾的,两者不能兼 得。因此,对于 N 值一般应视具体情况,采用折衷办法确定。根据 过程的实际发展趋势,N 值大体有如下四种选择方法: (1)水平式 与 N 值无关。 (2)脉冲式 趋势仅在某一段时间突然增加或减少,随后又保持
2

也就是趋势保持不变,移动平均值是无编差的,M 值

不变,N 取得越大,M 的误差 S a 越小,因此 N 应取得较大些。 (3)阶梯式 趋势仅在开始一段时间保持不变,然后增加或减少
2

到一个新的水平后又保持不变,N 取得越小,M 的误差 S a 越小,因 此 N 应取得较小。 (4)斜坡式 趋势周期的递增或递减,M 总是比实际趋向落后,

因此 N 应取得越小越好。 一般情况下,如欲加大原始数据的修正力度,则 N 宜取大些,如 果希望加大对外界变化的反映力度,则 N 宜取小些。N 的取值范围 一般为 3~20。 例 3 我国 1980~1990 年工业劳动人数见表,用二次移动平均数 法预测 1991~1994 年的劳动人数。
年 份 人 数

1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 5600 5796 5930 6092 6257 6567 6851 7141 7436 7738 8045
1980~1990 年我国工业劳动人数(万人) (表 5)
31

首先,选择移动平均周期 N。本例中数据趋势较明显,呈直线趋势, 为尽量反映近期变化动向,可取 N=3。利用移动平均公式,首先计算 一次移动平均数:
?1? M3 ? ? 5600 ? 5796 ? 5930 ? / 3 ? 5775.33 ?1? M4 ? ? 5796 ? 5930 ? 6092 ? / 3 ? 5939.3
?????????????????????

?1? M11 ? ? 7436 ? 7738 ? 8045? / 3 ? 7739.7

在此基础上再计算二次移动平均数:
? 2? M5 ? ? 5775.3 ? 5939.3 ? 6039 ? / 3 ? 5935.8 ? 2? M6 ? ? 5939.3 ? 6039 ? 6305.3? / 3 ? 6112.5
???????????????????????

? 2? M11 ? ? 7142.7 ? 7438.3 ? 7739.7 ? / 3 ? 7440.2

计算结果见表(4) 。根据表(4)的数据可建立线性趋势模型:

? 11?? ? a ? b ? Y 11 11
前已计算得: a11 ? 8039.2,b11 ? 299.5 则:

? 12 ? 8039.2 ? 299.5 ?1 ? 8338.7 Y ? 13 ? 8039.2 ? 299.5 ? 2 ? 8638.2 Y ? 14 ? 8039.2 ? 299.5 ? 3 ? 8937.7 Y ? 15 ? 8039.2 ? 299.5 ? 4 ? 9236.2 Y
由此得 1991~1994 年劳动人数分别为 8338.7,8638.2,8937.7,和 9236.2 万人。 (二) 指数平滑模型 在时间序列预测过程中,一般来说历史数据对未来发展的影响是
32

不等价的,数据由近及远对未来的影响价值递减。如果这种递减遵循 指数规律,并以此进行预测,则可采用指数平滑法。指数平滑法比移 动平均法需要的数据量少,计算更为方便。 一次指数平滑公式为
1? St?1? ? aYt ? ?1 ? a ? St?? 1

其中 St

?1?

——t

期数据的指数平滑值,α —平滑常数,0<α <1,

Yt -现期数据值。
对上式递推展开 则得
1 ?1? ? St? ? ? aYt ? ?1 ? a ? ? aY ? 1 ? a S ? ? t ? 1 t ?2 ? ? 1? ? aYt ? a ?1 ? a ? Yt ?1 ? ?1 ? a ? St?? 2 2

依此类推可得一次指数平滑的一般公式为

St?1? ? aYt ? a ?1 ? a ? Yt ?1 ? ?1 ? a ? Yt ? 2 ? ……
2

? a ?1 ? a ? Yt ? L ? …… ? a ?1 ? a ? Yt ??t ?1? ? ?1 ? a ? S0?1?
L t ?1 t

? a ? ?1 ? a ? Yt ? k ? ?1 ? a ? S0?1?
k t k ?0

t ?1

上式表明,数据列 Yt , Yt ?1 , Yt ?2 ,……, 的权数分别是α ,α (1-α ),

a ?1 ? a ? ,…… 即离 t 时刻越远的数据, 权数越小, 而且权数的变化
2

呈指数几何级数。 用一次指数平滑法进行预测时, 将 t 期的平滑值作为 t+1 期的预测 值,即

? t ?1 ? S ?1? ? aY ? ?1 ? a ? S ?1? Y t t t ?1
用指数平滑法进行预测时,将会遇到两个影响预测结果的因素,
33

一是初始值 S0 的选取, 这是计算其它平滑值的基础, 如果数据较多, 根据指数平滑的原理,初始值的影响极小,则可用第一个数据代替; 如数据较少,可分析数据的发展趋势给定一个估计值,或采用最初几 个数据的平均值。 二是平滑常数α 的选择,α 对平滑效果影响很大。α 越大,平滑 效果越差,反之,α 越小,平滑效果越好。从推导得知:

?1?

S
其中

2

? St ? ?

a S 2 ?Yt ? 2?a

S

2

? St ? ——指数平滑值方差;

S ? Yt ? ——原始数据方差;
2

α ——平滑常数。 因为 0<α <1,故 S
2

? S1 ? ? S ?Yt ? ,这样 St 的指数平滑值与时间
2

数列 Yt 有相同的均值,即 E

? St ? ? E ?Yt ? ? at ,但方差前者小于后

者。从公式中可以看到,在α 值较大时,指数平滑值的方差与时间数 列 Yt 的方差差别不大。α 越小,指数平滑值方差减少程度越大。因而 指数平滑法是一种滤波器。原始数列各项以一定顺序输入滤波器,而 滤波器的输出数据即为指数平滑值的现值。α 越小,滤波能力越强, 对原始数列修匀程度越好。 与移动平均的 N 值选择相似,α 值的选择亦应采取折衷方式。 如果我们认为初始值选择比较正确, 意欲充分反映初始值对预测值的 影响,α 宜选择小些,亦即参与平滑的数据量多些。如果我们认为初 始值选择不正确,意欲尽快减少初始值影响,α 宣迭择大些,亦即参
34

与平滑数据少些。再者,如果从事长期预测,α 宜选择小些,使更多 数据参与平滑。如果从事短期预测,则α 宜取较大值,使少量数据参 与平滑,以加大对近期数据反映的力度。因为α 值与预测精度和预测 期限间,有一个函数关系,如下图所示。

平滑常数的选择除上述规律外, 平滑常数α 的选择主要还是依靠 经验,视具体问题分别而定。如下几条准则可供参考。 (1)如对初始值的正确性有疑问时,应取较大的α 值,以便扩 大近期数据的作用,而迅速减少初始值的影响。 (2)如果多项式模型中仅有某一段时间的数据为较优估计值, 则需取较大α 值,以便减少较早数据的影响。 (3)如时间数列虽有不同规则变动,但长期趋势接近某一稳定 常数时,则需取较小的α 值(一般为 0.05~0.20) ,使各观察值在现 时指数平滑中具有大小接近的权数。 (4)如果时间数列具有迅速且明显的变动趋势,则α 宜取较大 值(一般取 0.3~0.5) ,使新近数据对于现时的指数平滑值具有较大 价值,从而使新近变动趋势能强烈反映在预测值中。
35

(5)如遇变化甚小的时间数列,则α 宜取稍小些(一般取 0.1~ 0.4) ,使较早的观察值亦能充分反映于指数平滑值中。 据经验统计,α 取值范围较大,α =0.1~0.9 均可得到较好预测 结果。 对于α 取值不同,对消除初始值的影响程度、对数据的平滑能力 以及对外界变化的反映力度,可以从表(6)和表(7)的数据反映出 来。 1
Yt

2 3.14 3.65

3 5.54 3.83

4 2.68 3.81

5 2.93 3.72

6 3.30 3.67

7 2.7 3.57

8 4.30 3.64

9 3.24 3.59

10 2.04 3.42

3.46 3.34

St

表 (6) 1
Yt

2 3.14 3.78

3 5.45 4.79

4 2.68 3.53

5 2.93 3.17

6 3.30 3.25

7 2.71 2.93

8 4.30 3.75

9 3.24 3.45

10 2.04 2.60

3.46 4.74

St

表 (7)

, 假设取 St ?1 ? 3.325 并认为初始值是正确的,则取α =0.1,这
时随着不断引进新的数据,得到表(6)的平滑数值。

St ? 0.1Yt ? ?1? 0.1? St ?1 ? 0.1Yt ? 0.9St ?1
如果另取 St ?1

? 6.65, 并认为初始值不甚正确,取α =0.6,则平滑值

如表(7)所示。

St ? 0.6Yt ? ?1? 0.6? St ?1 ? 0.6Yt ? 0.4St ?1
36

表 (8)是用指数平表滑法处理例 3 数据的结果。表(8)中,对α =0.5,α =0.3 分别进行了指数平滑,一次指数平滑值见表第二、第 五列。以α =0.5 为例,取初始值 S0 ? 5500 ,则:

?? S1? ? ? aY1 ? ?1 ? a ? S0 ? 0.5 ? 5600 ? ?1 ? 0.5 ? ? 5500 ? 5550
1 1 1 1

?? S2 ? aY2 ? ?1 ? a ? S1? ? ? 0.5 ? 5796 ? ?1 ? 0.5 ? ? 5550 ? 5659.5
1

?? ?? S3 ? aY3 ? ?1 ? a ? S2 ? 0.5 ? 5930 ? ?1 ? 0.5 ? 5659.5 ? 5794.8
1

………………
?? ?? S11 ? aY11 ? ?1 ? a ? S10 ? 0.5 ? 8045 ? ?1 ? 0.5 ? ? 7444.7 ? 7744.9
1 1

与移动平均法一样,指数平滑法也可进行二次、三次或更多次 的平滑。 二次指数平滑是以一次指数平滑数据作为原始数据再平滑一 次,其公式为
2? St? 2? ? aSt?1? ? ?1 ? a ? St?? 1

三次指数平滑是二次指数平滑值作原始数据再平滑一次,其公 式为

St3 ? aSt? 2? ? ?1 ? a ? St3?1
二次、三次指数平滑的例子见表(8) ,其中以α =0.5 为例,则 二次指数平滑为
? ? S1? ? ? aS1? ? ? ?1 ? a ? S0 ? 0.5 ? 5550 ? ?1 ? 0.5 ? ? 5500 ? 5525
2 2 1 1 2 2

? ? ?? S2 ? aS2 ? ?1 ? a ? S1? ? ? 0.5 ? 5659.5 ? ?1 ? 0.5 ? ? 5525 ? 5592.3

…………
? 2? ?1? ? 2? S11 ? aS11 ? ?1 ? a ? S10 ? 0.5 ? 7744.9 ? ?1 ? 0.5 ? ? 7164.8 ? 7454.9

三次指数平滑为
37

? ? S1? ? ? aS1? ? ? ?1 ? a ? S0 ? 0.5 ? 5525 ? ?1 ? 0.5 ? ? 5500 ? 5512.5
3 3 2 2 3 3

? ? ? ? S2 ? aS2 ? ?1 ? a ? S1? ? ? 0.5 ? 5592.3 ? ?1 ? 0.5 ? ? 5512.5 ? 5552.4

…………
? 3? ? 2? ? 3? S11 ? aS11 ? ?1 ? a ? S10 ? 0.5 ? 7454.9 ? ?1 ? 0.5 ? ? 6903.6 ? 7179.3

指数平滑法的目的不仅在于对数据进行修匀,同时可用平滑数

Yt
1 2 3 4 5 6 7 8 9 10 11 5600 5796 5930 6092 6257 6567 6851 7141 7436 7738 8045

α =0.5

α =0.3

St? ?
1

St?

2?

St?

3?

St? ?
1

St?

2?

St?

3?

5550 5659.5 5794.8 5943.4 6100.2 6333.6 6592.3 6866.7 7151.4 7444.7 7744.9

5525 5592.3 5693.6 5818.5 5959.4 6146.5 6369.4 6618.1 6884.8 7164.8 7454.9

5512.5 5552.4 5623 5720.8 5840.1 5993.3 6181.3 6399.7 6642.3 6903.6 7179.3

5530 5601.7 5700.2 5817.7 5949.5 6134.8 6349.7 6587.2 6841.8 7110.7 7391

5509 5536.8 5585.8 5655.4 5743.6 5861 6007.6 6181.5 6379.6 6598.9 6836.5

5502.7 5512.9 5534.8 5571.0 5622.8 5694.3 5788.3 5609.3 6048.3 6213.5 6400.4

例 3 指数平滑结果表(8)

据建立多项式模型,二次指数平滑值可建立线性模型,三次指数平滑 值可建立二次抛物线模型。多项式模型的通式为

? t ?? ? a ? b ? ? ct ? 2 ? …+ gt ? n Y t t 2 ! n!
其中
? t ?? ——t+τ 期的预测值 Y

τ —预测超前时间

at、bt、ct ……gt ——多项式参数。
预测中常用的模型为 常数模型(零次多项式) :
38

? t ?? ? a Y t
一次多项式(线性模型) :

? t ?? ? a ? b ? Y t t
二次多项式(二次抛物线模型) :

? t ?? ? a ? b ? ? 1 c ? 2 Y t t t 2
利用指数平滑法和移动平均法计算多项式待定系数的公式,参 看多项式模型参数估算公式(这里省略掉) 。以表 8 中的数据为例, α =0.5 时用三次指数平滑值建立抛物线模型,则:
?? ? ? ? ? a11 ? 3S11 ? 3S11 ? S11 ? 3 ? 7744.9 ? 3 ? 7454.9 ? 7179.3 ? 8063.9
1 2 3

b11 ? ?

a 2 ?1 ? a ? 0.5
2 2

?1? ? 2? ? 3? ? ?? 6 ? 5a ? S11 ? 2 5 ? 4 a S ? 4 ? 3 a S ? ? ? ? 11 11 ? ?

2 ?1 ? 0.5 ? a2

? ?? 6 ? 5 ? 0.5 ? ? 7744.9 ? 2 ? ? 5 ? 4 ? 0.5 ? ? 7179.3? ?

? 350.002 c11 ? ?

?1 ? a ?
0.52
2

2

?1? ? 2? ? 3? ? ? S11 ? 2 S ? S 11 11 ? ?

?1 ? 0.5?

? ? 7744.9 ? 2 ? 7454.9 ? 7179.3?

? 21.74
则可建立预测模型为

? t ?? ? 8063.9 ? 350? ? 1 ? 21.74? 2 Y 2
若预测第 13 期的值,则:

39

? 13 ? 8063.9 ? 350 ? 2 ? 1 ? 21.74 ? 22 ? 8807.4 Y 2
指数平滑法的主要优点是运算简单,费用低,当需要预测的数 据量大时,特别是预测库存时,这种方法比其它方法具有明显优势。 例 4 已知某新产品前 12 个月的销售额如表 9, 试用指数平滑法 建立二次多项式模型,并预测其第 13、14 个月的销售额。 1 562.67 7 781.73 2 586.93 8 791.44 3 599.65 9 765.38 4 668.52 10 832.50 5 597.80 11 880.11 6 670.23 12 934.60

某产品销售额表(9) 首先取α =0.3,初始值 S0 ? 570 ,通过平滑模型计算得:
?? S1? ? ? aY1 ? ?1 ? a ? S0 ? 0.3 ? 562.67 ? ?1 ? 0.3? ? 570 ? 567.8
1 1

?1? S2 ? aY2 ? ?1 ? a ? S1?1? ? 0.3 ? 586.93 ? ?1 ? 0.3? ? 567.8 ? 573.5

…………
?? ?? S12 ? aY12 ? ?1 ? a ? S11 ? 0.3 ? 934.60 ? ?1 ? 0.3? ? 794.1 ? 836.3
1 1

? ? S1? ? ? aS1? ? ? ?1 ? a ? S0 ? 0.3 ? 5678 ? ?1 ? 0.3? ? 570 ? 569.3
2 1 2

…………
? 2? ?1? ? 2? S12 ? aS12 ? ?1 ? a ? S11 ? 0.3 ? 836.3 ? ?1 ? 0.3? ? 725.9 ? 759 ? 3? S1? 3? ? aS1? 2? ? ?1 ? a ? S0 ? 0.3 ? 569.3 ? ?1 ? 0.3? ? 570 ? 569.8

…………
? ? ? ? ? ? S12 ? aS12 ? ?1 ? a ? S11 ? 0.3 ? 759 ? ?1 ? 0.3? ? 672.6 ? 698.5
3 2 3

计算结果列于表(10)
40

由表(10)中数据与原始数据对比可见,数据经平滑修匀后, 与原始数据编离较大。为此,将α 改为 0.5,再进行三次指数平滑, 结果仍列于表(10)中。平滑结果有所改善,可利用三次平滑数据, 建立预测模型:

? 12?? ? a ? b ? ? 1 c ? 2 Y 其中: 12 12 12 2
? ? a1 2 ? 3S ?11 ? S? 2 ? S ? ? 3?1 29 3 3 . 2 2 3 1 2

b1 2 ? c1 2 ?

a 2 ?1 ? a ? a2
2

? ? ?? 6 ? 5 ?? a ? S ?11 ? 2 5 a? 4S ? 2 ? 3? ? 3 ? ? ? ? 4a ? S 2 1 2 ? ?12

53.9

?1 ? a ?

2

? ? 2 ? ? S ?11 ?? ? ? S? ?3 2 2S 1 2 ? ? 1 24.6

从而得:

? 12?? ? 933.2 ? 53.9? ? .3? 2 Y Y13 ? 933.2 ? 53.9 ? 2.3 ?12 ? 989.4 Y14 ? 933.2 ? 53.9 ? 2.3 ? 22 ? 1050.2
亦即第 13 和 14 个月的新产品销售额预测值分别为 989.4 和 1050.2。

Yt
1 2 3 4 5 6 7 8 9 562.67 586.93 599.65 668.52 597.80 670.23 781.73 791.44 765.38

α =0.3

α =0.5

St? ?
1

St? 2?
569.3 570.6 573.8 583.9 590.1 600.4 621.7 647.4 670.7

St?3?
569.8 570 571.1 574.9 579.5 585.8 596.6 611.8 629.5
41

St? ?
1

St? 2?
568.2 572.4 580.3 604.3 608.7 625.2 688.5 710.1 734.3

St?3?
569.1 570.8 575.6 589.9 599.3 612.3 640.4 675.3 704.8

567.8 573.5 581.3 607.5 604.6 624.3 6711.5 707.5 724.9

566.3 576.6 588.1 628.3 613.1 641.7 711.7 751.6 758.5

10 11 12

832.50 880.11 934.60

757.2 794.1 836.3

696.7 725.9 759

649.7 672.6 698.5

795.5 837.8 886.2

764.9 801.4 843.8

734.9 768.2 806

例 4 指数平滑结果表 (10) 时间序列模型还有分解预测模型,鲍克斯-詹金斯模型。定量预测模 型另一类就是因果关系模型。现简要介绍如下:

因果关系模型
因果关系模型是定量预测模型的主要方法之一, 主要用于研究不 同变量之间的相关关系, 用一个或多个自变量的变化来描述因变量的 变化 因果关系模型主要包括:趋势外推、回归分析、数量经济模型、 投入产出模型、灰色系统模型、系统动力学等。 一. 特点和适用范围 事物的发展不仅取决于自身的发展规律, 同时受多种外界因素的 影响,如果把预测值作因变量,那么影响预测对象发展的各变量则称 作自变量。研究因变量与自变量的关系,则是因果关系模型的任务。 因果关系模型在预测中应用最广,它因时间序列模型不同,不仅可以 从事短期预测,而且还可以从事中、长期预测,也可以预测宏观、中 观、微观问题。 二. 预测程序(略) 三. 预测方法及模型 (一)趋势外推 趋势外推法是一种常用的利用事物过去发展的规律, 推导未来趋 势的方法,这种方法简单适用,应用面广。在预测方法分类中,有的
42

将其划归为因果关系模型。有的将其划归为时间序列模型,有的将其 单列为一类。我们将其划归为因果关系模型。因为趋势外推的模型和 预测过程与囬归分析类同,可以作为回归分析的特例,即以时间为自 变量的回归分析。运用趋势外推法,要注意它有两个基本假设: (1) 事物是在同一条件或相近条件下发展的,即决定过去事物发展的原 因,也是决定未来事物发展的原因; (2)事物发展的过程是渐进的, 而不是跳跃的。 趋势外推模型种类很多, 实用预测中最常用的是一些比较简单的 函数模型,如多项式模型、指数曲线、生长曲线和包络曲线等。 1. 多项式模型 很多事物的发展的模型可用多项式表示, 下面举几个常用的多项 式模型。 一次多项式模型(线性模型) : Yt ? a0 ? a1t 二次多项式模型(二次抛物线模型) : Yt 三次多项式模型(三次抛物线模型) :

? a0 ? a1t ? a2t 2

Yt ? a0 ? a1t ? a2t 2 ? a3t 3
n 次多项式模型(n 次抛物线模型) :

Yt ? a0 ? a1t ? a2t 2 ? ……? ant n
多项式的系数一般采用最小二乘法计算。 下面举例说明如何用多 项式进行预测。 例 1 某厂自动车床 1981~1991 年的产量如表(1)所示,试拟 合其发展趋势,并预测 1992 年自动车床的产量。
43

年 份 产 量

1981 1982 1983 1984 1985 1986 1987 1988 405 695 898

1989 1990 1991

1005 1323 1708 2200 24485 2499 2788 3160

表(1) 通过图示分析,可知原始数据大致呈抛物线趋势。 首先将数据拟合成二次抛物线:

? ? a ? a t ? a t2 Y t 0 1 2
将原始数据变为 t 与 Yt 的对应关系如表(2) 。这里将年份改用时 间序号代表。 这里数据个数为 11 个,即 n=11,拟合多项式次数为 2。根据多

a0 =81.8264,a1 =275.002,a2 =0.4066。 元回归参数估算方法可得:
因此

? ? 81.8264 ? 275.002t ? 0.4066t 2 Y t
t

Yt

0 1 2 3 4 5 6 7 8 9 10 405 695 898 1005 1323 1798 2200 2485 2499 2788 3160 表(2)

用此模型预测 1992 年的产量,即 t=11 得

Y1992 =3156 台
即 1992 年的自动车床产量的预测数为 3156 台 2. 指数模型 很多研究表明,大量事物的发展,其定量特征表现为随时间按指 数或接近指数规律增长,其公式为

Yt ? Y0ekt 其中 Y0 是 t=0 时 Yt 的值。
很多技术特性是符合指数增长规律,如飞机速度、电站容量
44

等。经济发展的很多增长过程也表现为指数增长,如产值、投资等。 但是任何事物都不能按指数曲线无限制的发展,因而不能无限制外 推,否则将出现不可能结果。因为任何事物的发展都有一个极限值, 当接近极限值时,事物将不再按指数规律增长。
kt 指数曲线的参数估算可采用最小二乘法,对于曲线 Yt ? Y0e

取自然对数得:

ln Yt ? ln Y0 ? kt
令 ln Yt ? Y ,ln Y0 ? A, k ? B 则变为线性方程 Y=A+Bt 3. 生长(s)曲线 如考虑极值的影响, 发现事物的发展过程如同生物的生长过程 一样,经历发生、发展和成熟三个阶段,而每一阶段的发展速度又不 相同,如绘在坐标图上形同 s,所以生长曲线又称 s 曲线。生长曲线 与指数曲线不同,指数曲线的相对增量是一个常数,例如用微分方程 表示为

dYt ? kYt dt dYt / Yt ? k dt
而生长曲线的相对增量是的一个 Yt 的函数

dYt ? kYt ? l ? Yt ? dt dYt / Yt ? k ? l ? Yt ? dt
45

生长曲线有很多种模型,皮尔模型、龚珀兹模型、替代模型等 都是生长曲线模型。 (1) 皮尔模型为 其中 L-变量的极值 a, b-常数 t-时间 (2) 龚珀兹模型为

Yt ?

L 1 ? ae?bt

Yt ? Le

? be? kt

其拐点为 t=(lnb)/k,

Yt ? L / e

Yt ? e x p a2 ? t? 0t? (3) 替代模型为 1 ? Yt
其中

Yt -待替代的百分比
a -每年替代的百分数的一半 ( 即替代开始以后的替代初始速

度)

t 0 -当 Yt ? 50% 时的时间
t-时间 趋势外推模型还有包络曲线,这里就不作介绍了。

(二)

回归分析

与趋势外推法不同, 这种方法主要用研究不同变量之间的相关 关系。回归分析不仅是一种应用范围极广的预测方法,同时也是建立 数量经济模型的重要基础。回归分析主要包括一元线性回归、多元线 性回归、和非线性回归三种,而非线性回归又可通过一定的变换,转 变为线性回归形或。
46

1.一元线性回归
也称直线回归,这种方法可确定两个变量之间的直线关系,简 单、适用,即可用于短期预测,也可用于长期预测。其公式为

Yt ? a ? bX

这里 Yt 是预测值,a、b 是回归系数, X t 是自变量

例 2 已知 1977~1986 年的钢材消耗量与国民收入的关系大致 呈直线趋势,数据见表(3),用一元回归法拟合其趋势并预测 1987 年 的钢材消耗量。
年份 国民收入 钢材消耗量 1977 289.3 20.13 1978 324.9 24.46 1979 347.6 27.36 1980 1981 368.8 386.9 28.25 29.58 表 (3) 1982 419 31.54 1983 460 35.15 1984 522.1 40.06 1985 588.5 49.83 1986 632.0 51.93

设国民收入为自变量,钢材消耗量为因变量,则回归方程为

? ? a ? bX Y t t
由 a、b 的计算公式可得 a、b 的值

b?

? X Y ? n ? ? X ?? ? Y ?
t t t t 2 1 ? X ? n ?? Xt ? 2 t

1

?

157058.8 ?

1 ? 4339.1? 337.39 10 1 2001291 ? 4339.12 10

=0.08996

a?

1 1 Yt ? ? X t ? b ? n n

=

1 1 ? 337.39 ? ? 4339.1? 0.08996 ? ?5.297 所以一元回归方程为 10 10

? ? ?5.297 ? 0.08996 X 经检验,模型符合实际,给 Y t t
定自变量 X 一个值,如 1987 年国民收入为 675 亿元,则

Y ? ?5.297 ? 0.08996 ? 675 ? 55.429
47

即 1987 年的钢材消费量的预测值为 55.429 万吨。 2. 多元线性回归 一元回归法虽在预测中比较常用,但它毕竟是一种理想化的形式,在 事物发展的局部或某一个侧面、某一个阶段可以用一个自变量反映。 但任何事物都不是孤立存在的,都会受到多种因素的影响,技术的发 展如此,经济系统更是如此。这时仅依据一个自变量进行预测将难以 求得准确的因变量值。 在这种情况下, 就必须采用多元回归进行分析。 多元回归虽然具有预测精度较高的优点,但亦需注意,多个自变量之 间容易出现共相关现象, 而且计算工作量大, 一般需借助计算机进行。 多元回归也有多种方法,这里只讲多元线性回归。 多元回归法比一元线性回归法复杂, 但应用原理和一元回归方法 是一样的。 多元线性回归的基本模型为

Y ? b0 ? b1 X1 ? b2 X 2 ? ……? bm X m ? et
其中 Y 是因变量,X 是自变量,m 是自变量个数;

b1 , b2 ,??, bm 是回归系数; et 是随机误差。
例 3 某科研单位的年收入受每年的科研经费和科研人数的影 响,1980~1989 年的数据见表(4),试建立多元回归模型并预测科研经 费为 380 万元、研究人数为 200 人时的年收入。 年代 收入(万元) 科研经费(万元) 科研人数(人) 1980 1981 235 238 254 257
48

160 163

1982 1983 1984 1985 1986 1987 1988 1989

256 254 271 273 289 298 304 318 表

275 290 295 296 311 318 327 341 (4)

166 169 172 175 178 181 184 187

经考察分析,数据呈线性关系,其回归方程为

? ? b ?b X b X Y 0 1 1 2 2
经计算得: b0

? ?80.2305, b1 ? 0.6852, b2 ? 0.8746

所以回归方程为

? ? ?80.2305 ? 0.6852 X ? 0.8746 X Y t 1 2
当状然还要进行检验,我们不作介绍了。经过相关检验和假设检验, 说明所求模型是成立的,将 X 1 =380, X 2 =200 代入得:

? ? ?80.2305 ? 0.6852 ? 380 ? 0.8746 ? 200 ? 347.84 Y t
即科研经费为 380 万元、研究人数为 200 人时年收入的预测值为 347.84 万元。

49


赞助商链接
相关文章:
时间序列及其它预测案例_图文
时间序列的分解方法时间序列的分解方法比较简单,下面用一个实例加以说明。表 4-...时间序列分解预测法的应用 在求出时间序列各因素之后,即可根据时间序列分解模型...
AR,MA,ARIMA模型介绍及案例分析
AR,MA,ARIMA模型介绍及案例分析 - BOX-JENKINS 预测法 1 适用于平稳时序的三种基本模型 (1) AR( p) 模型(Auto regression Model)——自回...
ARIMA模型预测案例
第5 讲 ARIMA 模型预测案例 【例 1】(1120070693)中国公路客运量 ARIMA 模型(缺中间项的自回归模型) 中国公路客运量数据(1950?2005)序列与差分序列见图。 2,...
数学建模spss 时间预测,心得总结及实例_图文
数学建模spss 时间预测,心得总结及实例_学习总结_总结/汇报_实用文档。数学建模spss 时间预测,心得总结及实例 《一周总结,底稿供参考》我们通过案例来说明: 假设...
案例二 ARMA模型建模与预测指导
案例二 ARMA 模型建模与预测指导 一、实验目的 学会通过各种手段检验序列的平稳性; 学会根据自相关系数和偏自相关系数来初步判断 ARMA 模型的阶数 p 和 q,学会...
R语言文本挖掘预测模型案例分析报告 附代码数据
R语言文本挖掘预测模型案例分析报告 附代码数据 - 【原创】附代码数据 有问题到淘宝找“大数据部落”就可以了 R 语言文本挖掘预测模型案例分析报告 knitr::opts_...
德尔菲法的预测实例
德尔菲预测实例_经管营销_专业资料。德尔菲预测实例某公司开发了一种新产品,现聘请了 9 位专家对新产品投放市场 1 年的销售额进行预测。在专家作出 预测...
指数平滑法应用案例
指数平滑应用案例_数学_自然科学_专业资料。Excel 应用案例指数平滑 移动平均预测值实质上是以前观测值的加权和,且对不同时期的数据给予相同 的加权。这往...
外应占法预测实例
外应占法预测实例_法律资料_人文社科_专业资料。外应预测实例(一车子难以找回 1999 年正月十三下午 2 点多钟,朋友柳万中从大连给我打来电话,说中午几个朋友在...
多元线性回归预测模型论文
多元线性回归预测模型论文 - 伊犁师范学院数学与统计学院 2012 届本科毕业论文 多元线性回归统计预测模型 摘要:本文以多元统计分析为理论基础,在对数据进行统计分析的...
更多相关标签: