时间序列的预报幻灯片.ppt

《时间序列的预报幻灯片.ppt》由会员分享，可在线阅读，更多相关《时间序列的预报幻灯片.ppt（82页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、时间序列的预报第1页，共82页，编辑于2022年，星期六 第一节第一节 最小方差估计最小方差估计一一.最小方差估计准则最小方差估计准则n对未知随机变量或未知随机向量进行估计是时间序列分析中预测理论的重要内容。n令 表示依据量测Z对X所得的某种估计，称 为X的估计，它是与X的维数相同的随机向量，而且是Z的函数.第2页，共82页，编辑于2022年，星期六 n记估计的误差为 ，称 为估计的均方误差阵，若存在某种估计，使得估计的均方误差阵比任何其他的均方误差阵都小，则称这个估计是最优的。n注：这里所讲的估计的均方误差阵大小，是指矩阵的大小。设A,B为同阶对称阵，若A-B是非负定阵，则称A不小于B，记为

2、 ；若A-B是正定阵，则A大于B，记为AB。第3页，共82页，编辑于2022年，星期六 二二.条件期望与条件方差条件期望与条件方差n对于两个二阶矩有穷的随机向量，定义在给定Y=y的条件下，X的条件期望和条件方差分别记为E(X|y)和Var(X|y)其中p(x|y)是给定Y=y的条件下X的条件密度函数。第4页，共82页，编辑于2022年，星期六 n条件期望的性质：(1)E(X|Y)是Y的函数向量并与X的维数相同，E(X|Y)使随机向量(2)E(X|Y)具有线性性，即对k个相同维数随机向量 及k个常数 有，(3)EE(X|Y)=EX第5页，共82页，编辑于2022年，星期六(4)设X与Y为相互独立

3、的随机向量，则 Ef(X)|Y=f(X)其中f()为适当的函数向量，仅要求f(x)为随机向量。(5)Ef(Y)|Y=f(Y)第6页，共82页，编辑于2022年，星期六(6)(7)当X与Y的联合分布为正态分布时，则第7页，共82页，编辑于2022年，星期六第8页，共82页，编辑于2022年，星期六第9页，共82页，编辑于2022年，星期六 三三.最小方差估计最小方差估计n寻找最小方差估计就是寻找一个适当的Z的函数 为随机函数向量，较之其它估计 相应估计的均方误差阵达最小，即n可以证明：第10页，共82页，编辑于2022年，星期六 n注：（1）最小方差估计是无偏估计 （2）最小方差估计的均方误差阵

4、为第11页，共82页，编辑于2022年，星期六 四四.线性最小方差估计线性最小方差估计n若估计 为量测随机向量Z的线性函数其中 为n维随机向量，Z为m维随机向量，a为与 同维数的常值向量，A为非随机矩阵，选择a,A是估计的均方误差阵达最小的估计称为线性最小方差估计。第12页，共82页，编辑于2022年，星期六 nX的线性最小方差估计为 n 的性质：(1)线性最小方差估计具有无偏性(2)线性最小方差估计的估计误差为 且 与Z正交。第13页，共82页，编辑于2022年，星期六 n定义1.1 设X与Z分别是二阶矩的n维与m维随机向量，如果存在一个与X同维数随机向量 ，且具有下列性质：(1)可由Z的线

5、性表示，即 ,其中a和A分别为常值向量和常值矩阵(2)(3)与Z正交，即则称 是X在Z上的投影，记为第14页，共82页，编辑于2022年，星期六 n注：(1)线性最小方差估计 是X在Z上的投影(2)X在Z上的投影只能是线性最小方差估计，即投影是唯一的(3)当X与Z为联合正态时,X关于Z的条件期望和X在Z上的投影相等，即第15页，共82页，编辑于2022年，星期六第二节第二节 平稳线性最小均方平稳线性最小均方 误差预报误差预报n预报是根据现在和过去的观察资料，对未来时刻的取值进行估计。n设 为零均值平稳序列，为 的长度为k的样本，根据 对 (为正整数)做出估计 ，取估计的优劣标准为使估计误差 (

6、2.1)的方差 达最小，则 第16页，共82页，编辑于2022年，星期六(1)(2.2)称 为 步最小均方误差预报。(2)如果函数f是 的线性函数,那么(2.2)为正交投影，即 其中 ，称 为 步线性最小均方误差预报。第17页，共82页，编辑于2022年，星期六 n注：(1)当 为正态序列时，最小均方误差预报与线性最小均方误差预报是一致的。(2)当 为非正态序列时，最小均方误差预报要优于线性最小均方误差预报。第18页，共82页，编辑于2022年，星期六 例2.1:设 为相互独立的随机变量,令 分别求Y关于X的线性最小均方误差估计和最小均方误差估计。第19页，共82页，编辑于2022年，星期六第

7、20页，共82页，编辑于2022年，星期六 n线性最小均方误差估计具体的表述：选择使得达最小。由于因此，令第21页，共82页，编辑于2022年，星期六 则有，将理论自协方差函数 换成样本自协方差函数 ，得到，(1.5)第22页，共82页，编辑于2022年，星期六 于是 步线性最小均方误差预报为存在的问题：1.样本自协方差函数 由 计算出来，k很大，故(1.5)的计算量极大2.部分样本自协方差函数 无法计算。第23页，共82页，编辑于2022年，星期六 n为克服上述困难,不妨假设已获得的观测资料为所有的历史资料,即为 ，令 表示由k时刻和它之前的所有历史数据对 所作的 步线性最小均方误差预报，即

8、 其中 是使 步预报误差 的均方达到最小，则有第24页，共82页，编辑于2022年，星期六 1.1 1.1 步预报和预报误差方差步预报和预报误差方差n设 为ARMA(p,q)序列，其传递形式和逆转形式分别为且令 可知，第25页，共82页，编辑于2022年，星期六 n则序列 的 步线性最小均方误差预报为：(1.11)步预报误差为：(1.12)步预报误差方差为：(1.13)第26页，共82页，编辑于2022年，星期六 n注：(1)由(1.11)可得 (1.14)(1.14)表明，在k+1时刻的 步预报等于k时刻的 步预报加上k时刻的一步预报误差得加权修正项。(2)由两部分组成，第一部分 为 步预报

9、误差，第二部分 为 步预报第27页，共82页，编辑于2022年，星期六(3)等式(1.13)表明：在线性最小均方误差预报意义下，步预报误差方差仅与预报步数有关，而与预报起点k无关，并且步数愈大，预报误差也就愈大，即预报精度愈差。第28页，共82页，编辑于2022年，星期六 n对于ARMA(p,q)模型参数 与Green函数 的关系式：其中 当模型建立后，由模型参数 逐步递推得到，再由公式(1.12),(1.13),(1.14)计算出 步预报 和预报误差，误差方差。第29页，共82页，编辑于2022年，星期六 1.2 ARMA(p,q)1.2 ARMA(p,q)序列的平稳线性最小均方误差预报序列

10、的平稳线性最小均方误差预报一一.AR(p).AR(p)序列的预报序列的预报模型：n 步预报为n当 时，有 第30页，共82页，编辑于2022年，星期六 n定理1.1 设 为AR(p)序列，则观测到时刻k为止的各步预报有如下递推公式:(1.16)第31页，共82页，编辑于2022年，星期六 n注：定理1.1表明，对AR(p)序列，只要知道 这p个数据，就可递推求得AR(p)序列的任意步的平稳线性最小均方误差预报，比k-p+1时刻更早的历史数据对预报不起作用。第32页，共82页，编辑于2022年，星期六 n例1.2 试求AR(1)序列的预报和预报误差方差。第33页，共82页，编辑于2022年，星期

11、六 注:若 为正态噪声，则 的95%的置信区间的置信上、下限为第34页，共82页，编辑于2022年，星期六 n例1.3 试求AR(2)序列的预报。步预报为：(1)递推法：由初值 ,可递推出任意 步预报 。(2)差分方程法：满足二阶差分方程 其中推移算子B作用于 ，即第35页，共82页，编辑于2022年，星期六 n设 为特征方程 的特征根，当k固定时，由差分方程理论知：(1)为不同实根，则 步预报为 其中 满足 则，第36页，共82页，编辑于2022年，星期六(2)当 ，则 步预报为 其中 满足 则，第37页，共82页，编辑于2022年，星期六(3)当 ，记 则 步预报为 第38页，共82页，编

12、辑于2022年，星期六 n例1.4 设某地区年平均降水量 为540毫米，其偏差序列 为AR(2)序列其中 ，已知近5年的实测降水量(单位：毫米)为 求今后3年各年降水量的预测值。第39页，共82页，编辑于2022年，星期六第40页，共82页，编辑于2022年，星期六 二二.MA(q).MA(q)序列的预报序列的预报(1)(1)逆函数法预报逆函数法预报nMA(q)模型：(2.1)n其逆转形式为 (2.2)从逆函数出发，给出的预报为逆函数法预报。第41页，共82页，编辑于2022年，星期六 n定理1.2 设 为MA(q)序列，则观测到时刻k为止的各步预报为 (2.3)其中 由下式递推得到，(2.4

13、)第42页，共82页，编辑于2022年，星期六第43页，共82页，编辑于2022年，星期六第44页，共82页，编辑于2022年，星期六 n注：(1)公式（2.4）不仅适用于MA(q)序列，也适用于ARMA(p,q)序列。(2)MA(q)序列和AR(p)序列两者预报有一个根本区别：对MA(q)序列的预报 依赖于k时刻和k时刻以前的全部历史数据。第45页，共82页，编辑于2022年，星期六 n在实际应用中，我们采取有穷和(取求和项适当的大)代替(2.3)式中的无穷和，从而近似的线性最小均方误差预报为 (2.7)n对于M的选择，可根据 收敛于零的速度适当地选择，以保证预报的精度。第46页，共82页，

14、编辑于2022年，星期六(2).(2).向量递推预报向量递推预报n引理1.1 设 为ARMA(p,q)序列，则有如下预报递推公式：(2.8)第47页，共82页，编辑于2022年，星期六第48页，共82页，编辑于2022年，星期六 n定理1.3 设 为MA(q)序列，则线性最小均方误差预报向量 满足如下递推公式 (2.9)其中，第49页，共82页，编辑于2022年，星期六第50页，共82页，编辑于2022年，星期六第51页，共82页，编辑于2022年，星期六 n注：递推公式(2.9)中的初值 可由逆函数公式直接得到，当 很小时，也可取 。n例1.4 求MA(2)序列：的逆函数法预报。解：逆函数为

15、：由系数公式，第52页，共82页，编辑于2022年，星期六 故其逆函数法预报为：当 时，取M=13,其近似预报为：第53页，共82页，编辑于2022年，星期六第54页，共82页，编辑于2022年，星期六 n例1.5求MA(2)序列：的向量递推预报。第55页，共82页，编辑于2022年，星期六 三三.ARMA(p,q).ARMA(p,q)序列的预报序列的预报(1)(1)逆函数法预报逆函数法预报n预报公式为：(2.10)其中 这里，第56页，共82页，编辑于2022年，星期六 n注：在实际应用中，我们采取有穷和(取求和项适当的大)代替(2.10)式中的无穷和，从而近似的线性最小均方误差预报为，(2

16、.11)第57页，共82页，编辑于2022年，星期六(2)(2)向量递推预报向量递推预报nARMA(p,q)模型：(2.12)n其传递形式为：(2.13)n注意到，当 时，(2.14)其中，,故当 时，由(2.14)递推求得。第58页，共82页，编辑于2022年，星期六 n定理1.4 设 为ARMA(p,q)序列，则线性最小均方误差预报向量 满足如下递推公式 (2.15)其中，第59页，共82页，编辑于2022年，星期六第60页，共82页，编辑于2022年，星期六 n注：递推公式中的初值 的选取方法与MA(q)序列选取是一致的。第61页，共82页，编辑于2022年，星期六 n例1.6 求ARM

17、A(1,2)序列的逆函数法预报和向量递推预报。解：(1)模型的逆转形式为：于是，逆函数为：第62页，共82页，编辑于2022年，星期六 逆函数法预报为选取M=30,于是相应的近似预报为第63页，共82页，编辑于2022年，星期六(2)模型的传递形式为：于是，其Green函数为：故，向量递推预报为：第64页，共82页，编辑于2022年，星期六 第三节 一类非平稳序列的线性 最小均方误差预报一一.ARIMA(p,d,q).ARIMA(p,d,q)序列的预报序列的预报n设 为ARIMA(p,d,q)序列,，令，则，可通过 及 表示为 (3.1)第65页，共82页，编辑于2022年，星期六 当tm=d

18、时，由 这里 为ARMA(p,q)序列，我们有如下定理，n定理1.5 设 为ARIMA(p,d,q)序列，则已知条件下关于 线性最小均方误差预报为,其中，第66页，共82页，编辑于2022年，星期六第67页，共82页，编辑于2022年，星期六 n注：定理1.5表明：通过ARMA(p,q)序列 的线性最小均方误差预报，可给出原序列 的预报。n例1.7 求 为ARIMA(1,1,0)序列的预报。解：的模型为：令 ，于是 为AR(1)序列，则第68页，共82页，编辑于2022年，星期六 二二.季节性乘积模型的预报季节性乘积模型的预报n设 为 模型,即：令 ，则 令 这里，为ARMA(p,q)序列，而

19、 为ARIMA(p,d,q)序列。第69页，共82页，编辑于2022年，星期六 序列，则由 的线性最小均方误差预报，得到 的预报为，其中又因为，于是，原序列 的 步预报公式为，或第70页，共82页，编辑于2022年，星期六 n例3.1 设 为 序列：求 的 步预报。第71页，共82页，编辑于2022年，星期六 三三.非平稳序列的叠合模型的预报非平稳序列的叠合模型的预报n具有确定性的趋势分量和周期分量的非平稳序列为 (3.2)其中 为趋势分量，为周期分量，为平稳ARMA(p,q)序列。n将(3.2)式中的三部分各自进行 步预报，之和作为的 步预报，即，第72页，共82页，编辑于2022年，星期六

20、 n具体地，(1)仅含线性趋势：，则(2)仅含多项式趋势：,则(3)仅含指数趋势：实数,则第73页，共82页，编辑于2022年，星期六 n(4)仅含周期分量：,则n(5)既具有趋势又有周期分量的非平稳序列的预报为，第74页，共82页，编辑于2022年，星期六 例3.1 考虑如下叠合模型：故 可表为模型的两部分 步预报之和，即第75页，共82页，编辑于2022年，星期六 n例3.2：化学衰减数据的建模。采用一个温度跃变系统来研究化学反应中浓度变化的衰减过程。具体办法是产生电火花温度甚高，使溶液温度急骤升高而逐渐衰减，趋于平稳状态。它可用一个吸收光谱的测量系统，直接测量到反应溶液的浓度。取采样间隔

21、为0.02秒的100个测点记录数据。第76页，共82页，编辑于2022年，星期六第77页，共82页，编辑于2022年，星期六 n由上图可知，数据变化存在着负指数衰减，且趋于稳态值32左右，故可考虑如下叠合模型去拟合：其中，Wt为平稳序列。第78页，共82页，编辑于2022年，星期六第79页，共82页，编辑于2022年，星期六 n通过运算得，Xt的叠合模型为其中第80页，共82页，编辑于2022年，星期六 n预报结果第81页，共82页，编辑于2022年，星期六 n 随堂练习：设 为平稳正态序列，且为AR(1)序列：已知 的95%的置信区间为(9.5,15.5)试求：（1）模型参数（2）一步预报第82页，共82页，编辑于2022年，星期六