第十一章 并行准高斯高阶递归滤波算法研究.pdf

2014年 5月 计 算 数 学 第36卷第 2期 May．2014 MATHEMATICA NUMERICA SINICA Vo1．36，No．2 并行准高斯高阶递归滤波算法研究冰1) 王玉柱1,2姜金荣 ，)迟学斌 岳天祥。 (1．中国科学院计算机网络信息中心超级计算中心，北京 100190) (2．中国科学院大学，北京 100049) (3．中国科学院地理科学与资源研究所，北京100101) 摘 要 三维变分同化系统中一个重要的问题是背景误差协方差矩阵 B及其逆的求解．背景误差协 方差矩阵的水平变换部分采用递归滤波运算，可以简化矩阵的求解，解决了背景误差协方差矩阵 B及其逆难以求解的问题．本文对准高斯高阶递归滤波的算法原理和过程进行了深入研究．因为 递归滤波并行的低可扩展性制约了高阶递归滤波算法在三维变分同化系统中的应用，所以本文提 出了阶段二维区域剖分并行化方法，实现了并行准高斯高阶递归滤波算法库．数值试验表明，四阶 递归滤波 1次的效果明显优于一阶 4次的滤波效果；并且高阶递归滤波并行算法 64核时能达到 大约 50倍的加速，并行效率高达 78％，具有 良好的加速效果和较强的可扩展性． 关键词：三维变分资料同化；高阶递归滤波；背景误差协方差；并行算法；区域剖分 MR (2000)主题分类：65Y05，68W10 1．引 言 在 20世纪80年代把变分方法引入资料同化后，三维和四维变分资料同化在数值天气预 报中起着越来越重要的作用．在三维变分同化系统中尽可能正确地求解背景误差协方差矩阵 B及其逆是非常重要的．但是背景误差协方差矩阵逆的计算非常困难，此矩阵不仅规模巨大 (阶数超过 10 ×10)，而且还常常是病态的．因此，出现了采用滤波方法来近似求解 B 的相 关研究． 1982年，英国气象中心的Purser和 McQuigg第一次把递归滤波作为气象客观分析的工 具，用于区域尺度的二维分析．递归滤波是一种连续的近似方法，在一系列的连续迭代中，通 过严格的质量控制和逐步降低 地“域性的影响”来修改背景场的信息 [1J．后来，Purser等 【j和 Hayden等 lJ通过运用平滑算子来模拟一种高效的数值和空间适应分析方案，每个空间的平 滑算子都是由简单的递归滤波数值算子合成，并把递归滤波应用在美国国家环境卫星、数据 及信息服务中心 (NESDIS)的三个业务中．另外，他们的研究表明，当迭代次数接近于无穷大 时，递归滤波器会逐渐趋近于高斯滤波器．2003年，Purser等 5【_进一步对高阶递归滤波进行 研究，研究表明：n阶递归滤波运算 1次与一阶递归滤波运算 札次花费的代价相当，但是前者 的滤波效果明显好于后者．2011年，Zhuang等 _6J针对北海／波罗的海设计和实现了预运行的 三维变分资料同化系统，并在背景误差协方差的水平部分采用各向同性递归滤波来处理． 在国内，国家卫星气象中心的孙东联和王宗皓 L7J详细介绍了递归滤波的原理和方法，并 结合个例讨论和描述了它的特点和优点．张华等 Is】在处理背景误差协方差矩阵的逆时，对 目 标函数进行预调节，水平变换部分采用二维水平一阶递归滤波，可以大大减少内存消耗和计算 2013年 7月 19日收到． )基金项 目：国家高技术发展计划 (2010AA012301，2012AA01A309)资助项 目 )通讯作者：姜金荣，E．Mail：jjr@sccas．crt． 180 计 算 数 学 量．曾忠一 9【J在高斯滤波器的基础上，简单讨论了一阶递归滤波器的数学原理和滤波系数的 求解过程． 由于三维变分资料同化系统的计算量和 内存需求量比较大，为了使同化技术能满足业务 系统的高时效要求，必须实现资料变分同化方法的高效并行计算．因此，开展相应高阶递归滤 波并行算法的研究与实现是十分必要的．本文首先研究高阶递归滤波算法的数