数值天气预报(4-6-3)--第6节空间平滑和滤波.pptVIP

第4章 正压原始方程模式 * * * * * * * * * * * * * * * * * § 4.6 空间平滑和滤波 在数值模式积分过程中，非线性相互作用会产生短波扰动以及混淆误差，而短波系统能量的迅速增长，会引起不稳定现象，有时甚至会使数值积分中断。为了滤除数值模式积分过程中的短波扰动，有效地抑制非线性计算不稳定，通常采用空间平滑滤波的方法。 本节将介绍常用的平滑滤波方法。大致包括三个方面 1。所用的计算公式； 2。计算公式的性能； 3。计算公式的应用。 一、一维平滑 1、三点平滑 对于任意一个离散变量 (4.118) 定义一维平滑算子： 式中S为空间平滑系数。由于上述空间平滑公式仅涉及到三个网格点上的F 值，故称为三点平滑公式。 也可写成： (4.119) 空间平滑滤波的性能分析： 首先，将要素表示为谐波形式： 其中C 为一常数， 、k 和 ? 分别为谐波的振 幅、波数和位相。将上式代入空间平滑公式： (4.120) 请同学们思考C代表的物理意义是什么？ (4.121) 可以看出：三点平滑只改变了原来谐波的振幅，而波数和位相均未发生变化。 原谐波形式 定义三点平滑的响应函数R为：平滑后的振幅/原振幅 表示经过平滑后各波长的波的振幅消长程度。 响应函数R与平滑系数S和波长L有关，选择适当的平滑系数可以使平滑运算具有良好的滤波性能。 则有： 如果取 S=1/2 (4.122) (4.123) 经过平滑后： 波长为L=2?x的波动 完全 滤除，L?2?x的波受不同程 度的衰减，例如，L=10 ?x 的波R=0.905,而L很大的波 并未受到影响. 图4.14 三点平滑的响应函数随波长的变化 (S=1/2) 响应函数 R 变化曲线如右图所示。 在数值天气预报的过程中，一般对某一要素场需要进行多次平滑。 2、多重平滑 如果对要素场平滑M 次，则可以求得其响应函数: 例：对波长为L=10?x的波，如果采用5次平滑运算，即M=5，则有： 可见，多次平滑对波长较长的波具有相当严重的衰减作用。 当S取-1/2时，三点平滑响应函数为： 经过平滑运算，各波长的波均有不同程度的增长。 (4.124) 利用空间平滑所具有的以上特性，可以设计出一种所谓的“正逆平滑”运算，既过滤短波，又使长波的衰减很小。 正、逆平滑系数分别取为 对要素场进行连续两次的平滑，其响应函数为： 上式及图4.14均表明，经过这种平滑运算后,波长为L=2?x 的波仍然可以被过滤掉，而对波长较长的波的衰减却可以减弱。例如,同样对 L=10?x的波,如果采用5次正逆平滑运算?,其(R2)5=0.9777. (4.125) 二、二维平滑 对二维变量Fi,j ，则进行 1、五点平滑 (4.126) x方向三点平滑 y方向三点平滑 然后将上述两式相加，取其平均则得： 五点平滑公式 下面来求其对应的响应函数，设如下形式的波动解： (4.127) 式中k和l分别为x方向和y方向的波数 显然，五点平滑公式对应的响应函数为： 将其代入平滑公式，得到： (4.128) (4.129) 应用一维三点平滑公式 (4.119)式，对先进行x 方向的平滑，然后再对其进行y 方向平滑，就可得到九点平滑公式： 2、九点平滑 其对应的响应函数为： 由于九点平滑对于长波的衰减比五点平滑要强的多，所以，通常采用五点平滑公式。 (4.130) 第4章 正压原始方程模式 * * * * * * * * * * * * * * * * *