《水文随机分析第四章090918》由会员分享,可在线阅读,更多相关《水文随机分析第四章090918(24页珍藏版)》请在金锄头文库上搜索。
1、第四章第四章 水文随机序列的预报水文随机序列的预报概述概述 平稳线性最少方差预报平稳线性最少方差预报 AR(p)序列预报序列预报 门限自回归模型门限自回归模型 2024/7/211水文水资源学院水文水资源学院第一节第一节 概述概述 模机水文模型应用:模机水文模型应用:最简单一个例子:最简单一个例子:对中心化(离均差后)对中心化(离均差后) 系列建立了一阶自回归模型:系列建立了一阶自回归模型:(1) 水文序列水文序列 已知已知1984年年 2024/7/212水文水资源学院水文水资源学院对对(1)式取数学期望得:式取数学期望得:(预报公式)(预报公式) 其实,这是一个期望预报,把随机变量取为其实
2、,这是一个期望预报,把随机变量取为0值(平均值)值(平均值)因此一步预报误差:因此一步预报误差: 已知水文序列作预报。已知水文序列作预报。建立水文随机模型(假如无趋势及周期)建立水文随机模型(假如无趋势及周期)给出预报公式和预报方法。给出预报公式和预报方法。利用现在及过去观测值作预测,必要时作实时利用现在及过去观测值作预测,必要时作实时修正。修正。区间预报,误差分析,对模型作评定及检验。区间预报,误差分析,对模型作评定及检验。2024/7/213水文水资源学院水文水资源学院第二节第二节 平稳线性最小方差预报平稳线性最小方差预报 对正态平稳系列作对正态平稳系列作l步预报步预报 问题:当前时刻问题
3、:当前时刻k和过去时刻和过去时刻 为已知,需对未来时刻随机变量为已知,需对未来时刻随机变量 ,可记为,可记为 。所谓平稳线性最小方差预报定义为:所谓平稳线性最小方差预报定义为: 水文变量水文变量 2024/7/214水文水资源学院水文水资源学院一、一、ARMA(p、q)序列差分形式公式:序列差分形式公式:离均化离均化 ARMA(p、q)模型两端取数学期望(条件)两端取数学期望(条件)当当MA(0,q)模型时,如模型时,如 当为当为AR(p)模型时,将在后面作介绍。模型时,将在后面作介绍。2024/7/215水文水资源学院水文水资源学院二、二、ARMA(p、q)序列传递形式预报公式:序列传递形式
4、预报公式:把上式把上式(2)式中式中t用用 代替代替2024/7/216水文水资源学院水文水资源学院上式两边取条件数学期望上式两边取条件数学期望该公式该公式 2024/7/217水文水资源学院水文水资源学院 三、三、ARMR(p、q)递推预报公式(可用于实时修正)递推预报公式(可用于实时修正)对对(2)式作一些变换,把式作一些变换,把k改为改为k+12024/7/218水文水资源学院水文水资源学院事实上以上公式含义:事实上以上公式含义: 要做要做k+1时刻时刻l步预报步预报= =在在k时刻作时刻作l+1步步预报值预报值 + + (k+1时刻实测值与在时刻实测值与在k时刻对时刻对k+1时刻预测值
5、之差)时刻预测值之差)* *权重系数权重系数 。 相当于对相当于对 时刻作预测时,用误差来校正时刻作预测时,用误差来校正 预测预测值,把这个称为实时校正(现时校正)。值,把这个称为实时校正(现时校正)。预测预测如如 显然在显然在 时刻预测时刻预测 比用比用 预测预测 效果效果要好,用了误差校正。要好,用了误差校正。 2024/7/219水文水资源学院水文水资源学院预报误差(评定与检验)预报误差(评定与检验)合格率(允许相对误差范围内如合格率(允许相对误差范围内如20%)合格率)合格率85%,甲等;,甲等; 70%-84%,乙等;,乙等;60%-69%,丙等,丙等(枯季径流允许误差(枯季径流允许
6、误差30%,每日径流,每日径流20%)2024/7/2110水文水资源学院水文水资源学院如如 区间估计。区间估计。2024/7/2111水文水资源学院水文水资源学院第三节 AR(p)序列预报)序列预报(一般(一般ARMA(p,q)模型有此结果)模型有此结果)一步预报公式一步预报公式 先估计模型参数先估计模型参数 2024/7/2112水文水资源学院水文水资源学院 再作一步预报再作一步预报 ,又由,又由 预报预报 如此递推可以作如此递推可以作l步预测,当然还可以作实时校正预报。步预测,当然还可以作实时校正预报。实例:某站中心化实例:某站中心化 模型已建立(正态分布假模型已建立(正态分布假定)定)
7、求预测求预测1983-1985平均流量。平均流量。2024/7/2113水文水资源学院水文水资源学院以以1982年为年为k时刻进行预报时刻进行预报区间预报区间预报 2024/7/2114水文水资源学院水文水资源学院实时校正:以实时校正:以1983年为年为k+1进行实时校正预报进行实时校正预报2024/7/2115水文水资源学院水文水资源学院门限自回归模型门限自回归模型1978年提出年提出H.Tong)()(Threshold Autoregressive Model) 观察:观察: 这样生成这样生成 不是线性模型生成随机过程,若生成数据不是线性模型生成随机过程,若生成数据100年,用线性模型拟
8、合效果较差即残差相关比较大,不一定年,用线性模型拟合效果较差即残差相关比较大,不一定独立,正态。独立,正态。2024/7/2116水文水资源学院水文水资源学院模型形成(一般)模型形成(一般)一般一般l取取23,不宜取太多了,否则太复杂。,不宜取太多了,否则太复杂。 这个模型实际说明了随机序列是分段线性的,即每个这个模型实际说明了随机序列是分段线性的,即每个 区间内可用自回归模型来描述,但这些模型序数在不同区间是区间内可用自回归模型来描述,但这些模型序数在不同区间是不一样的。不一样的。2024/7/2117水文水资源学院水文水资源学院建模:建模:已知已知 当用传统线性模型建模时发现,独立性差,或
9、预测精度较当用传统线性模型建模时发现,独立性差,或预测精度较差,则可以考虑用门限自回归模型来建模,(当然还应考虑是差,则可以考虑用门限自回归模型来建模,(当然还应考虑是否在成因上做分段线性门限值?)否在成因上做分段线性门限值?)据以上模型据以上模型 采用以上已知数据采用以上已知数据 ,来辩认模型参数步骤如下:,来辩认模型参数步骤如下:1:令:令 。分成段,如确定各。分成段,如确定各段自回归阶数段自回归阶数 ,L为逐段自回归允许最大阶数,这为逐段自回归允许最大阶数,这里都认为各段有相同里都认为各段有相同L值。值。 2024/7/2118水文水资源学院水文水资源学院把实数轴分成把实数轴分成l个区间
10、个区间一般是预报值取为一般是预报值取为 这样可以把这样可以把 取值区间分取值区间分l成个。成个。对于已知对于已知 观察分析动态数据中在不同区间个数观察分析动态数据中在不同区间个数(以排在(以排在 开始统计),如开始统计),如L=3,d=1,从,从 开始统开始统计在不同区间个数。计在不同区间个数。2024/7/2119水文水资源学院水文水资源学院记第一个区间记第一个区间N1个个第二个区间第二个区间N2个个 第第l个区间个区间Nl个个设动态数据中有:设动态数据中有:对于第一段线性自回归模型具体形式对于第一段线性自回归模型具体形式已知:已知:2024/7/2120水文水资源学院水文水资源学院用最小二
11、乘法可求出用最小二乘法可求出 可以想象,给定以上一组系数则可求出一组残差可以想象,给定以上一组系数则可求出一组残差 及其残差平方和及其残差平方和 ,即可得第一段自回归模型,其他,即可得第一段自回归模型,其他 段求法一致。显段求法一致。显然不同然不同 2024/7/2121水文水资源学院水文水资源学院2、对于固定、对于固定 RSS( k1 )为第为第1个模型残差平方和(最小二乘法估计参数)个模型残差平方和(最小二乘法估计参数)显然显然k1 愈大,则残差愈小,但愈大,则残差愈小,但AIC考虑惩罚因子考虑惩罚因子 则第二项大,故找到则第二项大,故找到 。其他段。其他段的自回归模型阶数一样估计。的自回归模型阶数一样估计。 2024/7/2122水文水资源学院水文水资源学院3、固定、固定 达到最少的即为门限值达到最少的即为门限值 4、估计值、估计值AIC( )一般一般 5、 2024/7/2123水文水资源学院水文水资源学院不同不同 那么那么d初值取多少?初值取多少? 若效果不好,可先进行常用若效果不好,可先进行常用 2024/7/2124水文水资源学院水文水资源学院