自然正交分解在1996年降水场预测中的应用研究.pdfVIP

自然正交分解在1996 年降水场预测中的应用 张书俊 武汉中心气象台 武汉430074 摘要 利用自然正交分解将降水场分解为特征向量和时间系数 利用周期分析方法延长时间系数 再利用延长后的时间系数和 特征向量求出预测的降水场 关键词 自然正交分解 降水场预测 特征向量 时间系数 1 引言 自然正交分解多用于对要素场分布特征的分析 或对要素场随时间变化的特征分析 它还能用来对要 素场进行预测 它可以根据不同的精度要求 利用较少的数据去逼近原始资料分布 预测时只要对几个特 [1] 征向量的时间系数作出预报 然后利用特征向量和预报出的时间系数 便可做出要素场的分布预报 该 预报方法的关键是如何预报时间系数 有的选取因子建多元线性回归方程预报时间系数 该方法若选取的 时间系数项数较多 则需要建立多个多元线性回归方程 这样所要选用的因子太多 计算的工作量太大 也不便于计算操作 本文采用周期分析方法延长时间系数 它操作简单 自动化程度高 运算速度快 特 别对于时间系数选项较多时更为适用 该方法不仅适用于对空间场的预报 也可以将时间序列当作空间场 处理进行计算预测 例如将一年12个月的降水量序列 当作12个点的空间场 利用上述方法可以预报下 一年12个月的降水量分布 该方法经十几年的业务应用效果不错 有一定的实用价值 2 样本资料及处理 用1953-1995年湖北省40站5-9月降水量资料作为样本资料 根据国家局制定的6级制化分标准 将月降水量化为6级 平均值为1961-1990年的30年平均 这样 处理后的资料和预报值都是6级制级别值 应用方便 另外6级资料将降水资料缩小在1-6之间 在某种 程度上起到标准化的作用 该资料计算速度快 收敛快 3 基本思路及原理 这里简要介绍利用自然正交分解和周期分析方法进行降水场预报的思路和基本原理 首先建立湖北省 40个代表站降水资料 用矩阵形式F表示 其中i为样本数 i=12 n J为站点数 J=12 m ij F F F 11 12 1m = F (1) ij F F F n1 n2 nm 利用雅可比旋转法进行自然正交分解 将矩阵F分解为时间系数和特征向量两个矩阵相乘 ij F T X ij ij ij 根据预报精度的要求 可选P个特征向量和时间系数进行计算 由自然正交分解的性质上式可近似的 写为 F T X ij ip pj 式中T为时间系数 它只与年序有关 X为给定资料的特征向量 它不随年代变化 上式用标量表示 ip pj 可写成 p F T X T X T X T X 4 ij ip pj i 1 1j i 2 2 j ip pj p 1 根据 1 式 预报场可写成 p F T X T X T X T X (5) i lj i 1,p pj i 1,1 1j i 1,2 2 j i 1,p pj p 1 式中F 为降水场的预报场 i 1j 周期分析是人们很熟悉的统计方法 其计算方法和过程不再详述 经多年应