大气海洋数据同化方法Kalman滤波课件.ppt
大气/海洋数据同化方法（2）Kalman滤波,Kalman滤波基本思想：,预报场,初始场,分析场=预报场+K（观测-模式解）,最小方差估计 P，K,1 表示求逆T 表示转置,一个简单例子：,信号源运动方程为,系统误差:方差为0.1的高斯白噪声,问题：请估计状态,状态矩阵,状态转换矩阵A,观测矩阵H,观测误差:方差为0.25的高斯白噪声,结果分析,简单推广应用：,动力模型要简单,气象/海洋MODEL,统计模型,统计预报,回归方程,应用猜想1）气温旬预报2）Elnino预报,一个实际例子一维大气边界层动力初值化模式如下:,h,x,Ug,Vg,PBL,问题：如何消化观测数据形成风廓线？,三种思路：1）：直接利用观测数据进行插值2）：运用松弛逼近法3）：运用Kalman filters哪些可取？哪些不可取？,运用松弛逼近法,运用Kalman filters首先将方程写成 x(t+1)=A*x(t)的形式,湍流垂直交换系数方案采用OBrien(1970)的廓线系数公式,状态向量：,离散化之后的1-D边界层方程可以写为：,(2N,2N),一维标准Kalman滤波实验结果,KF,Without KF,观测所在垂直层次2468,Kalman滤波,Kalman滤波的几个发展阶段,标准（standard）Kalman滤波 扩展（extended）Kalman滤波简化（simplified）Kalman滤波集合（ensemble）Kalman滤波,Kalman filters,标准Kalman滤波给出的状态估计是最优的 扩展（extended）Kalman滤波是针对非线性系统提出的。在非线性情形下，模式状态转换矩阵可以是模式状态的函数，观测转换矩阵可以是模式状态与观测的函数保留1阶导数项（即获取非线性方程的切向线性方程），类似标准Kalman滤波的形式同样可以得到，但此时的Kalman滤波给出的估计不再是最优的，而是次优的。为了提高精度，每进行一次Kalman滤波之后保持不变并使用相同的观测，重复滤波过程对模式状态进行更替，这就是叠代（iterated）扩展Kalman滤波。,简化Kalman滤波与集合Kalman滤波的诞生Kalman滤波原理并不复杂，但具体实施起来有相当的难度。其一：模式只是对真实物理过程的一种近似，近似的好坏程度即误差很难确定，因而Kalman滤波中的模式误差 较难给定。其二：计算模式预报误差协方差矩阵 是Kalman滤波的核心关键算法，需要 倍的模式积分时间，而且还需要存储 维矩阵的空间，给出转换矩阵 的形式；对于当前原始方程数值预报模式，通常达到，而且很难给出 的形式。这些意味着要付出昂贵的计算时间并占用大量的机器内存。,基于以上原因，针对 对Kalman滤波进行简化成为必然。目前基本上有两种简化途径：一种是简化（simplified）Kalman滤波，减少 的维数；一种是集合（ensemble）Kalman滤波，改变 的算法。,简化Kalman滤波有多种思路。一种简化Kalman滤波认为(Phillips,1986;Dee,1991;Cohn,1991)，模式预报误差协方差的传播方程向前传播的信息过多，其实并没有必要传播有关每一个状态变量的信息；它主张通过了解模式里物理过程的内在规律和抓住模式的动力学特征，用一小部分参量去代表整个模式误差；相应地，只计算与这一小部分参量所对应的模式误差相关的模式预报误差协方差，对预报误差协方差的其余部分作静态处理即认为它们不随时间传播，从而达到减小 维数的目的。,集合Kalman滤波,MonteCarlo法（集合积分法）来计算模式预报误差协方差。将模式状态预报看成近似随机动态预报，用一个模式状态集合（设数目为M）去代表随机动态预报中的概率密度函数，通过向前积分模式状态集合很容易计算不同时间的概率密度函数所对应的统计特性（如均值与协方差）。,假定 代表模式预报集合，则 当集合数目增加时，解概率密度的误差以 速率趋近0；对于一个实用的集合数目，如100左右，此误差将会被统计噪音而非动态随机预报所控制；此时，所消耗的计算时间约为100倍模式积分时间。,标准Kalman滤波与集合Kalman滤波实验结果对比,En-KF,S-KF,Kalman滤波的优点与缺点优点：模式预报状态协方差 不断更新，不需要伴随模式，可以考虑非线性问题（ensemble Kalman filters)。缺点:计算量太大！,Kalman滤波的应用进展 Kalman滤波在大气数据同化中正蓬勃发展。在简单、低维的数据同化问题中，标准Kalman滤波和扩展Kalman滤波发挥着重要作用，算法已经成熟。但对复杂、高维和强非线性数据同化问题，扮演主要角色的是简化Kalman滤波和集合Kalman滤波。海洋数据同化方面，Evensen（1996）借助一个二层准地转海洋模式用集合Kalman滤波同化卫星高度计数据来模拟海流）。,欧洲中期天气预报中心（ECMWF）1993年就在Courtier的建议下开始发展自己的简化Kalman滤波，并将它的4DVAR系统向简化Kalman滤波扩展；1997年初次得到了简化Kalman滤波结果，发现4DVAR中运用简化Kalman滤波可以考虑更多的动力影响。1998年又用一个低分辨率的扩展Kalman滤波作为工具诊断和测试了简化Kalman滤波中的几个理论和实践方面的问题，将简化Kalman滤波朝业务使用方面又推进了一步。最近，Mitchell等（2000）发展了一种基于集合Kalman滤波的估计模式误差的方法。可以说，Kalman滤波技术正在与伴随变分技术展开“竞争”，瓜分或者共同联手占有大气数据同化“市场”。,从Kalman Filters得出Optimal interpolationdynamic system:x(t+1)=A*x(t)+u(k)E(u*u)=Qmeasurement system:y(t)=M*x(t)+w(k)E(w*w)=ROptimal interpolation无dynamic systemA=1,Kalman Filters,Optimal interpolation,谢谢！,
