《天气学诊断分析》第1章有限差分方法和客观分析.pptx

天气学诊断分析;学时安排和要求;参考书籍;课程简介;诊断场： 物理量在某时刻的空间分布是该物理量在某时刻的分布实况，在计算该物理量的方程中不含有对时间的微商项。 预报场： 对未来某时刻某物理量的预报结果，在预报方程中含有该物理量对时间的微商项。 诊断方程： 反映各气象要素场之间关系的，不含有对时间微商项的方程。 数据： 观测数据是离散的，不规则的； 网格数据的获得：主观分析和客观分析(今后不做说明，都表示是网格点资料) ;;;天气诊断分析定量客观地展示某些物理量空间分布和时间变化情况，做定量的计算和解释，确定主要和次要的物理因素，给出天气现象发生的空间结构，不涉及这些因素随时间的演变。 将方程组离散化，便于计算机处理： 谱方法，ＦＦＴ，有限差分等等。 差异的原因： 1. 直接参加计算的温、压、湿、风等气象要素值有观测误差，且这种误差在计算过程中会累积放大； 2. 用差分代替微分所产生的截断误差； 因此，在诊断方程求解过程中，采用好的计算方案，使计算误差尽可能小些，数值解和实际场量值的近似程度也就好些。 ;内容;;;;授课内容;第一章 有限差分方法和客观分析;第一章 有限差分方法和客观分析;;;;;;§1 简单有限差分公式;;;1.2 一阶微商的几种差分方案;;1.2.1 两点差分方案及其精度;;1.2.2 三点差分方案及其精度;;1.2.3 两点式和三点式的几何意义;;1.2.4 五点式差分方案及其精度;;;;;1.3 二阶微商的几种差分方案;;1.3.2 五点式差分方案;;;;;;§2 拉普拉斯算子的差分格式;2.1 二维函数的泰勒展开式 ——有限差分方案的出发点;;;;;;;2.2 几种拉普拉斯算子的差分格式;;;;;;(1.1.14)和(1.1.15)式都是用5个点的值来表示▽2，都称为拉普拉斯算子的五点差分格式，都具有二阶精度，其不同之处在于（1.1.14）式是用对角线上的点来表示;2.2.2 拉普拉斯算子的九点差分格式;;;§3 雅可比算子;;3.1 雅可比算子的差分形式;;3.2 Arakawa的雅可比差分方案;;;雅可比算子的差分形式之一;;;;§4 差分方程的精确性问题;以一阶中心差为例;;讨论;§1.5 客观分析;;;垂直方向上插值;(一)线性插值 线性插值也叫一次插值，这种方法比较简单。 设x0，x1两点上有函数y0，y1，要求多项式y＝P1(x)，P1(x0)＝y0，P1(x1)＝y1其几何意义是通过两点式直线方程，即已知两点求直线上任一点插值函数。 已知两点A(x0,y0)，B(x1,y1)求任一点M(x，y) 两点式直线方程：;;（二)二次插值 二次插值也叫抛物插值。 已知三点(x0,y0)，(x1,y1)，(x2,y2)求点(x,y) 插值点 ;(三)拉格朗日插值 拉格朗日插值公式是二次插值公式向n次插值推广的结果，因此也叫n次插值。 任一点的函数可用它本身x0,x1,x2，…,xn上的y0,y1, y2,…,yn，乘以适当的系数来表示，可把插值公式写成 ;;;; 由于气压随高度呈对数变化，应用插值公式时常取气压的对数形式lnP。 例如，已知探空资料500hPa、700hPa、850hPa的高度值分别为5500、3000、1500（单位：位势米），用拉格朗日二次插值公式求600hPa的高度。 ;小结