气象统计方法第一章气象资料及其表示方法课件.ppt

1、气象统计方法气象统计方法主讲：温 娜南京信息工程大学大气科学学院2015年3月本课件主要参考南信大李丽平老师课件 要求：要求：1、掌握气象上常用的一些统计方法，运用这些方法进行资料分析，在此基础上作一些简单的气候预测等2、了解其它方法在业务和科研中的应用，对后期学习起引导作用 气象统计方法：气象统计方法：利用统计方法对气象资料样本进行分析来估计和推测总体的规律性，为气象预报提供依据。 课时安排 总学时：48学时（1-16周） 讲授为主, 课堂练习 考核方式 平时成绩（出勤、课堂作业） 期末考试参考书目：参考书目：1、李湘阁等气象统计方法，2、黄嘉佑著气象统计分析与预报方法，气象出版社，2004

2、.33、施能著气象统计预报，气象出版社，2009.114、吴洪宝等著气候变率诊断和预测方法，气象出版社，2010.65、魏凤英著现代气候统计诊断与预测技术，气象出版社，2009.9教学内容教学内容第一章第一章 气象资料及其表示方法气象资料及其表示方法第二章第二章 气候稳定性检验气候稳定性检验第三章第三章 选择最大信息的预报因子选择最大信息的预报因子第四章第四章 一元线性回归分析一元线性回归分析第五章第五章 多元线性回归分析多元线性回归分析第六章第六章 气候趋势分析气候趋势分析第七章第七章 主分量分析主分量分析第八章 聚类分析1. 气象要素气象要素 大气温度、压力、空气湿度、风向和风速、降水、云

3、、雾、雷暴、辐射、能见度等 还有土壤、陆面植被、海洋等监测要素一、气象资料（研究对象）一、气象资料（研究对象）2. 气象监测全球监测系统全球监测系统ARGO计划气象监测意义：1. 记录天气、气候的实际情况2. 了解气候的基本状况3. 分析研究气候变化规律4. 气候预测（第一张天气图的诞生）（第一张天气图的诞生）拉萨站（29N, 91E）气温观测数据江苏省（116-122E，30-36N）月平均气温数据（1958-2007）江苏省（116-122E，30-36N）月平均气温数据（1958-2007） 气象中单个或多个要素可看成为统计学气象中单个或多个要素可看成为统计学中单个或多个变量。中单个或多

4、个变量。二、气象资料的表示二、气象资料的表示第一节第一节 单个要素的气象资料单个要素的气象资料1 1、数据的表示：数据的表示： 某要素x有n 次观测值，其向量表达形式为或者 n 样本容量 时间序列概念时间序列概念 数据是随时间变化的序列，习惯称为时间序列。数据是随时间变化的序列，习惯称为时间序列。例如：例如：取某要素月平均值的取某要素月平均值的n n年资料年资料 几何意义：几何意义： （1 1）n n 维空间中的一个点维空间中的一个点 （2 2）一维空间（单坐标）中的）一维空间（单坐标）中的n n个点个点 举例：举例： 1.某站点1958至2007年的气温 2.某站点1960至2010年一月份

5、的气温 3.某站点某时段冬季/夏季降水 要素样本中资料分布的特点要素样本中资料分布的特点-用一些统计量用一些统计量表征。表征。 1 1）平均值平均值 含义：平均值是要素总体数学期望的一个估计。含义：平均值是要素总体数学期望的一个估计。反映了该要素的平均（气候）状况。反映了该要素的平均（气候）状况。平均值概念在气象上的应用：平均值概念在气象上的应用： 气象上的月平均气温、年平均气温及某要素多气象上的月平均气温、年平均气温及某要素多年平均值就是这种统计量。年平均值就是这种统计量。例如：某地气候状态 （对于逐月资料，一般分别求各月多年的平均值，所以会有12个月平均值场；逐日资料也是类似）。 气象研究

6、中，不同时间分辨率气象资料的使用（逐日资料、月资料、年资料等） 举例： 江苏气候？江苏省1958-2007年月平均气温江苏省各月气温平均值（气候态）江苏省1958-2007年冬季月平均气温 2 2）距平）距平 含义：含义：反映数据偏离平均值的状况反映数据偏离平均值的状况 ，也是通常所，也是通常所说的说的异常异常。 距平序列距平序列：单要素样本中每个样本资料点的距平：单要素样本中每个样本资料点的距平值组成的序列称为距平序列，也可以记为距平向值组成的序列称为距平序列，也可以记为距平向量。量。气象上的应用：气象上的应用： 中心化的概念：中心化的概念： 把资料处理为距平的方法叫把资料处理为距平的方法叫

7、中心化。中心化。 气象上常用距平值代替原样本中的资料值作气象上常用距平值代替原样本中的资料值作为研究对象。为研究对象。中心化的必要性：中心化的必要性： 因为气象要素的年变化周期影响很大，各月因为气象要素的年变化周期影响很大，各月的平均值不一样，为了使之能在同一水平下比较，的平均值不一样，为了使之能在同一水平下比较，常使用距平值（比如之前的举例）。常使用距平值（比如之前的举例）。特性：特性：距平值的平均值为距平值的平均值为0 0，使用方便；，使用方便； 直接作为预报值，比较直观（偏高直接作为预报值，比较直观（偏高/ /偏偏 低）。低）。 举例：举例： 某地气温的变化情况？ （气温偏高、偏低是相对

8、于它的平均值而言）江苏省冬季气温的异常（1958-2007）如何正确计算异常场？不可取江苏省气温异常（1958-2007）3 3）方差和均方差（标准差）方差和均方差（标准差） 对气象要素对气象要素x,x,资料长度资料长度n,n,其表达式：其表达式： 含义：含义： 是均方差（标准差），描述样本中资料与平均是均方差（标准差），描述样本中资料与平均值差异的值差异的平均平均状况，反映变量围绕平均值的状况，反映变量围绕平均值的平均平均变化程度变化程度（离散程度），（离散程度）， 是方差。是方差。方方 差差向量表示形式：气象上的应用：气象上的应用： 1)1)如果如果1212月份气温标准差比月份气温标准差比

9、1 1月份大，反映了月份大，反映了1212月份气温随时间变化幅度比月份气温随时间变化幅度比1 1月大。月大。 2)2)对于同一个月（例如对于同一个月（例如1212月），月），如果南京气如果南京气温的标准差比温的标准差比拉萨小拉萨小，表表明明拉萨冬季气温的变化拉萨冬季气温的变化幅度比较大。幅度比较大。 ( (内陆日变化较沿海大内陆日变化较沿海大, ,这个日变化大小的比较使这个日变化大小的比较使用的是标准差的比较用的是标准差的比较) )江苏省各月气温标准差3 3）均方差小的要素预报比大的困难还是容易？）均方差小的要素预报比大的困难还是容易？原因？原因？4 4）变量减去某常数后均方差相同。）变量减去

10、某常数后均方差相同。5）标准差与变量值同量纲，）标准差与变量值同量纲，一般用标准差表一般用标准差表示变量取值变化的大小。示变量取值变化的大小。 对气象要素对气象要素x,x,资料长度资料长度n,n,其表达式：其表达式： 特征：特征：1 1）标准化变量的平均值为）标准化变量的平均值为0 0。 2 2）标准化变量的方差为）标准化变量的方差为1 1。1 2,tztxxxxtnS， 为何要进行标准化为何要进行标准化？ 各要素各要素单位不同单位不同、平均值平均值和和标准差标准差也也不同。为使它们在同一水平上比较，采用不同。为使它们在同一水平上比较，采用标准化方法，使它们变成同一水平的无单标准化方法，使它们

11、变成同一水平的无单位的变量位的变量-标准化变量标准化变量。江苏省气温异常及其标准化降水距平百分率降水距平百分率距平距平/ /平均值平均值* *100%100%1 1）计算降水距平，即观测值减去平均值）计算降水距平，即观测值减去平均值2 2）1 1步骤所得结果除以该平均值，乘以步骤所得结果除以该平均值，乘以100100 ，即为降水距平百分比，即为降水距平百分比注意：注意：当观测值序列时间比较长，超过当观测值序列时间比较长，超过3030年，可以年，可以选择选择1980198020092009的平均值，作为步骤的平均值，作为步骤1 1中的平均值中的平均值4 4） 1 1）意义：）意义： 说明变量值变

12、化的大小。说明变量值变化的大小。2 2）变率）变率 绝对变率绝对变率：距平绝对值的平均。：距平绝对值的平均。 相对变率：相对变率： 绝对变率与平均值之比绝对变率与平均值之比11naiiVxxnarVVx3 3）变差系数）变差系数 标准差与平均值之比（标准差与平均值之比（% %） 表示变量的相对变化，表示变量的相对变化， 2111nxpiiSVxxxxn（ - ） 注意：注意： 绝对变率和标准差的数量级与变量平均绝对变率和标准差的数量级与变量平均值的量级有关。值的量级有关。 有些同类型变量，彼此之间有些同类型变量，彼此之间平均值差别平均值差别大大，若要比较它们的变化性用绝对变率和，若要比较它们的

13、变化性用绝对变率和标准差不恰当，应当利用标准差不恰当，应当利用相对变率相对变率或或变差变差系数系数。5 5）频率分布频率分布 累积频率概念的引入：累积频率概念的引入： 平均值和均方差相同，但取值很大区别，区别其平均值和均方差相同，但取值很大区别，区别其特征，就需要引入新的统计量特征，就需要引入新的统计量-累积频率累积频率。 累积频率累积频率：变量小于某上限的次数与：变量小于某上限的次数与总次数之比。（总次数之比。（样本特征样本特征直方图直方图） 总体总体( (母体母体) )：统计分析对象的全体。统计分析对象的全体。 样本：样本：总体中的一部分。总体中的一部分。三、总体和样本三、总体和样本理解与

14、应用：理解与应用：总体的特征是总体的特征是客观客观存在的存在的样本的特征样本的特征随样本而变随样本而变，与其有关，与其有关 的变量均称为随机变量，如平的变量均称为随机变量，如平 均值、均方差等均值、均方差等选取选取有代表性有代表性的样本很重要的样本很重要样本量样本量n=30,n=30,根据数理统计中的大数定理推断得根据数理统计中的大数定理推断得到。到。气象上的总体指无限总体，一组气象资料就是无气象上的总体指无限总体，一组气象资料就是无限总体的样本。总体与样本关系的相对性。限总体的样本。总体与样本关系的相对性。大数定律大数定律 大数定律又称大数法则、大数率。大数定律又称大数法则、大数率。 在一个

15、随机事在一个随机事件中，随着件中，随着试验次数的增加试验次数的增加，事件，事件发生的频率趋发生的频率趋于一个稳定值于一个稳定值；同时，在对物理量的测量实践中，；同时，在对物理量的测量实践中，大量测定值的大量测定值的算术平均也具有稳定性算术平均也具有稳定性。在数理统。在数理统计中，根据计中，根据贝努利定理贝努利定理 辛钦定理辛钦定理：当当n n很大时，很大时，算术平均值接近数学期望；算术平均值接近数学期望；频率以概率收敛于事频率以概率收敛于事件的概率。件的概率。1 1）分布函数）分布函数 无限无限总体总体的累积频率的累积频率 称为概率密度函数，其最常见的形式是称为概率密度函数，其最常见的形式是正

16、态分布正态分布 和和 分别是总体平均值（期望）和标准分别是总体平均值（期望）和标准差，可以用样本平均值和均方差去估计。差，可以用样本平均值和均方差去估计。， xdxxfxPxF)()()()(xf222)(21)(xexf正态分布曲线 标准化变量的平均值为标准化变量的平均值为0 0、方差为方差为1 1。特点特点：1）标准化正态分布随机变量的绝对值大于2.58（1.96）的概率仅为0.01（0.05） 2）不同要素变量作图 3）标准化变量值的取值范围(+/-3之间),大于3的概率仅为0.0027 世界气象组织（WMO）旱涝年确定标准：距平达到或者大于2倍均方差，出现概率不到5%年份。 江苏省月平

17、均气温标准化数据P(|T|=2)=0.0483; P(|T|=3)=0.005 江苏省冬季气温数据分布图2 2）正态分布的统计检验）正态分布的统计检验 大多数气候诊断方法和预测模型是在气候变量呈正态分布假定前提下进行的，所以对气候变量是否呈正态分布形态的检验是十分必要的。正态分布检验不仅可以判断原始变量是否遵从正态分布，还可以检验那些原本不遵从正态分布，但经过数学变换后的变量是否已成为正态分布形式。（1 1）概念概念 峰度系数峰度系数与与偏度系数偏度系数是用来衡量随机变量概是用来衡量随机变量概率密度分布曲线形状的数字特征，描述了气候变率密度分布曲线形状的数字特征，描述了气候变量的分布特征。量的

18、分布特征。 偏度系数偏度系数：表征曲线峰点对期望值（平均值）：表征曲线峰点对期望值（平均值）偏离的程度。偏离的程度。 峰度系数峰度系数：表征曲线分布形态顶峰的凸平度：表征曲线分布形态顶峰的凸平度（即渐进于横轴的陡度）。（即渐进于横轴的陡度）。（2）标准偏度系数和峰度系数的计算公式为：）标准偏度系数和峰度系数的计算公式为：其中，其中， 和和 s 分别为样本均值和标准差。分别为样本均值和标准差。xg2= g1= 偏度系数：偏度系数： 峰度系数：峰度系数： 标准偏度系数的意义标准偏度系数的意义测量数据分布的不对称性情况，测量数据分布的不对称性情况，刻画以平均值为中心的偏向情况，刻画以平均值为中心的偏

19、向情况，g g101010，表示正偏，即均值在峰值的右边；，表示正偏，即均值在峰值的右边；g g1=01=0，表示对称分布，表示对称分布 标准峰度系数的意义标准峰度系数的意义测量分布图像坡度的缓急程度g20, 表明图形坡度偏陡若若g1=0, g2=0g1=0, g2=0时，表明研究变量为理想的正态分布变量时，表明研究变量为理想的正态分布变量（3）显著性检验）显著性检验（4 4）资料的正态化资料的正态化 1 1、正态化的必要性、正态化的必要性 各类统计预报模型和统计检验方法各类统计预报模型和统计检验方法(Ftu (Ftu 检验检验) )要求资料是符合正态分布。要求资料是符合正态分布。 年年 月平

20、均气温月平均气温 气压气压 多雨地区的月降水多雨地区的月降水量符合量符合! ! 日降水日降水和和少雨地区月降水少雨地区月降水通常偏态。旬通常偏态。旬 候降水不一定候降水不一定! !2 2 2、资料正态化处理方法、资料正态化处理方法 1.立方根或四次方根； 2.双曲正切转换（纠正课本公式）-旬降水； 3.化为有序数后的正态化转换（标准化和正态化）。 江苏省全年月降水数据分布图 *也可以理解为同一要素多个格点也可以理解为同一要素多个格点（站点）的资料，下面慢慢体会。（站点）的资料，下面慢慢体会。 多个气象要素的样本如何表示？多个气象要素的样本如何表示？-矩阵矩阵。 设有设有m m个气象要素，每个要

21、素有个气象要素，每个要素有n n次观测值，次观测值，则数据矩阵为：则数据矩阵为：（2.1）11121212221212()mmnmmnnnmxxxxxxxxxXx xx第第t t个样本的资料向量为个样本的资料向量为nt, 2 , 1（2.2）123()tttmtx x xxtx两个方面来研究问题两个方面来研究问题( () )“R R型分析型分析”：研究不同要素或变量研究不同要素或变量( (如同一如同一 要素不同格点之间）的关系。（列）要素不同格点之间）的关系。（列）“Q Q型分析型分析”：研究样本之间的关系（行）。研究样本之间的关系（行）。1 1、n n维空间中的维空间中的 m m个点（列）个

22、点（列） m m个变量（格点、站点）确定了个变量（格点、站点）确定了n n维空间中的维空间中的 m m个点；个点； 用来研究变量（或者不同格点、站点）之间的关系。用来研究变量（或者不同格点、站点）之间的关系。-相关关系，如两个变量之间的相关系数相关关系，如两个变量之间的相关系数2 2、m m维空间中的维空间中的n n个点（行）个点（行）-空间点聚图空间点聚图 一个样本对应一个样本对应m m维空间中的一个点维空间中的一个点； 分析样本之间的关系时用到，如寻找相似个例。分析样本之间的关系时用到，如寻找相似个例。三、均值向量三、均值向量 m m个变量的样本平均值组成的向量。个变量的样本平均值组成的向

23、量。 m m维空间中的维空间中的n n个点的个点的重心重心（各部分受到的（各部分受到的重力作用集中于一点，这一点就是重心重力作用集中于一点，这一点就是重心 ）。）。 12()mx xxxntitixn11xmi, 2 , 1（2.3） 多年平均多年平均1月气温月气温（19712010年）年） 多年平均多年平均7月气温月气温（19712010年）年） 多年平均多年平均1月降水量月降水量（19712010年）年） 多年平均多年平均7月降水量月降水量（19712010年）年）四四、距平向量距平向量1112121222d1d2dm12()ddd mddd mddndndnmxxxxxxxxxXxxx(

24、1,2,1,2,)dijijjxxxin jm，其表达式为其表达式为：1.1.协方差协方差 衡量任意两个气象要素（变量）之间衡量任意两个气象要素（变量）之间关系的统计量关系的统计量( (正、负相关关系正、负相关关系) )，另外一，另外一个统计量叫相关系数（以后讲解）。个统计量叫相关系数（以后讲解）。表达式：表达式：)( )(11jjtntiitijxxxxnsjintjtitxxxxn11mji, 2 , 1,（2.4）距平的内积距平的内积协方差气象意义的进一步理解：协方差气象意义的进一步理解： 1）反映了两个气象要素异常关系）反映了两个气象要素异常关系的平均状况，或的平均状况，或者两个变量的

25、正、负相关关系。者两个变量的正、负相关关系。两变量关系越密切，两变量关系越密切，其协方差的绝对值越大其协方差的绝对值越大， 如理解（气温为例）：如理解（气温为例）： 前冬气温负距平（冷）、后冬正距平（暖）前冬气温负距平（冷）、后冬正距平（暖）-协协方差负值方差负值-反相关反相关 前冬气温正距平（暖）、后冬正距平（暖）前冬气温正距平（暖）、后冬正距平（暖）-协协方差正值方差正值-正相关正相关 2）变量自身的协方差就是方差。）变量自身的协方差就是方差。 距平的乘积距平的乘积X1X1与与x2x2的距平符号相同率高，有相同的变化趋势，的距平符号相同率高，有相同的变化趋势，x2x2与与x3x3的距平符号

26、相反率高，有相反的变化趋势；的距平符号相反率高，有相反的变化趋势；两组变量均有良好的相关关系。两组变量均有良好的相关关系。问题：问题：协方差带单位，不同要素之间不好比较，以后学习相关系协方差带单位，不同要素之间不好比较，以后学习相关系 数可解决这个问题。数可解决这个问题。 2.2.协方差矩阵协方差矩阵 m m阶对称矩阵阶对称矩阵, ,对角线元素是第对角线元素是第i i个变量的个变量的方差方差, ,撇号代表距平。撇号代表距平。)(ijsSmji, 2 , 1,（2.6）1TddnX X ppnnnppppppnppppnnddxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

27、xnXXnS22112222112122111121222222211111211111 ppppppsssssssssS212222111211 协方差的另一种表示：协方差的另一种表示： 用用离差积离差积表示，构成离差矩阵。表示，构成离差矩阵。 总体协方差矩阵的无偏估计总体协方差矩阵的无偏估计)( )(1jjtntiitijxxxxss)(ijssSS（2.7）mji, 2 , 1,SS11Sn（2.8）1月气温标准差月气温标准差7月气温标准差月气温标准差1月降水标准差月降水标准差7月降水标准差月降水标准差状态资料和统计特征状态资料和统计特征1 1、状态资料、状态资料 表征气象要素的各种状态

28、，观测结果无法用数据表示。 如雾、冰雹、霜-用“有”、“无”、“强”、“弱”等表示。 雨的强度-等级表示，如大暴雨、大雨等。 风的强度等级-1级、2、。12 2 2、频率表、分布列、频率表、分布列 如何描述状态资料的统计特征？如何描述状态资料的统计特征？ 列出各个状态出现的频率。列出各个状态出现的频率。 对样本而言是对样本而言是频率表频率表，总体而言就是，总体而言就是分分布列布列。针对气象要素的状态资料；针对气象要素的状态资料；概念：概念：统计多个气象要素（现象）的各种情况统计多个气象要素（现象）的各种情况 下的频率，组成一张多维频率表。下的频率，组成一张多维频率表。优点：优点：对预报有参考价

29、值。对预报有参考价值。 1 1、台风、台风20122012年年8 8月我省先后受到月我省先后受到5 5个近海个近海和和登陆台风登陆台风的直接或外围环流系统影响，台风影响个数偏多。的直接或外围环流系统影响，台风影响个数偏多。5 5个热带气旋分别是个热带气旋分别是9 9号台风号台风“苏拉苏拉”、1010号台风号台风“达维达维”、1111号强台风号强台风“海葵海葵”、1414号台风号台风“天天秤秤”和和1515号超强台风号超强台风“布拉万布拉万”。其中。其中，“海葵海葵”对我省的影响特别严重，其风雨程度超过对我省的影响特别严重，其风雨程度超过20052005年年的台风的台风“麦莎麦莎”，是，是199

30、11991年以来影响江苏最严重年以来影响江苏最严重的台风的台风。 据不完全统计，据不完全统计，20122012年年8 8月江苏月江苏省因强省因强对流天气过程共造成多人伤死，房屋损坏，对流天气过程共造成多人伤死，房屋损坏，农田受灾。农田受灾。8 8月苏州出现月苏州出现8 8个强对流天气日，其中出个强对流天气日，其中出现雷雨大风短时强降水，冰雹，龙卷风袭击，现雷雨大风短时强降水，冰雹，龙卷风袭击， 8 8月月2626日日1717：3030左右，江苏洪泽湖上出左右，江苏洪泽湖上出现了难得一见的现了难得一见的“龙吸水龙吸水”现象现象，巨大的水，巨大的水柱连接水天。灰色水柱从远处的乌云中垂落柱连接水天。

31、灰色水柱从远处的乌云中垂落湖面，像吊在空中晃晃悠悠的一条巨蟒，这湖面，像吊在空中晃晃悠悠的一条巨蟒，这时候的天空依然可见蓝天白云，拖着水柱的时候的天空依然可见蓝天白云，拖着水柱的那片乌云分外显眼。那片乌云分外显眼。 3 3、暴雨洪涝、暴雨洪涝4 4、高温、高温8 8月苏南大部地区出现月苏南大部地区出现2 26 6天高温天气常州出现天高温天气常州出现高温高温6 6天。因高温引发头痛的病人频频增加，其中感天。因高温引发头痛的病人频频增加，其中感冒、发烧引发头痛者最多。冒、发烧引发头痛者最多。5 5、雾霾、雾霾8 8月连云港市共出现区域性大雾天气月连云港市共出现区域性大雾天气2 2次，全市大次，全市

32、大部分地区出现能见度低于部分地区出现能见度低于500500米的浓雾，最强时段能米的浓雾，最强时段能见度低于见度低于200200米，对交通运输和旅客出行造成了一定米，对交通运输和旅客出行造成了一定影响。影响。课堂练习利用该资料，计算相关统计量和分析数据特征江苏省夏季降水资料（1978-2007）江苏省夏季降水的气候态及变化特征？1958-2007年海温(x)和降水(y)资料序列时间时间原始资料原始资料中心化中心化标准化标准化txyxyx*y*195827.3122.70.5050.540.482.07200627.657.70.73-14.580.70-0.60200725.158.7-1.68

33、-13.58-1.61-0.56平均值平均值 均方差均方差 协方差协方差 12.6 1 1、什么叫资料均匀或均一、什么叫资料均匀或均一2 2、什么是可靠的观测序列、什么是可靠的观测序列3 3、资料的订正、资料的订正 插补插补、纠正、延长、纠正、延长资料的订正资料的订正(插补、延长插补、延长)空间均一性空间均一性* *在统计各种气候指标、进行气候分析在统计各种气候指标、进行气候分析之前，对气候资料的质量应该进行之前，对气候资料的质量应该进行审查审查。 从统计意义上讲，各种气象要素观从统计意义上讲，各种气象要素观测值都可以看作随机变量。一地的气候测值都可以看作随机变量。一地的气候状态可以用各种要素

34、的统计特征值来表状态可以用各种要素的统计特征值来表示，如算术平均值、均方差等实质上是示，如算术平均值、均方差等实质上是要素总体数字特征的样本值。要素总体数字特征的样本值。气候指标气候指标是否反映真实的气候状态是否反映真实的气候状态，取决于，取决于: :1)数字特征样本值对总体值的数字特征样本值对总体值的抽样误差抽样误差如样本平均值的抽样误差如样本平均值的抽样误差 )()()()(xDnmxEnmnxnEmxE222221111nxDxDnxnDxD)()()()( 211一般的，样本容量越大，特征值的抽样误差越一般的，样本容量越大，特征值的抽样误差越小小。 m为真实平均数为真实平均数.2 2）

35、要素观测值对真值要素观测值对真值的观测误差的观测误差* *系统误差系统误差: :由于仪器不良、观测方法不完善由于仪器不良、观测方法不完善, ,在各次观测中的大小和符号保持不变。在各次观测中的大小和符号保持不变。 通过校正仪器和改善观测方法可基本消除。通过校正仪器和改善观测方法可基本消除。 温度不超过温度不超过0.05C0.05C，气压不超过，气压不超过0.050.05百帕，降百帕，降水量不超过水量不超过0.05mm0.05mm（根据规定进行器差订正后的根据规定进行器差订正后的观测值对真值的系统误差观测值对真值的系统误差） 如果在整个记录时期，误差无改变，如果在整个记录时期，误差无改变，则由则由

36、N N次次观测求得的算术平均值对真值的系统误差等于每观测求得的算术平均值对真值的系统误差等于每次观测中的系统误差。次观测中的系统误差。而而N N次观测值总和的系统误次观测值总和的系统误差将是每次观测的系统误差的差将是每次观测的系统误差的N N倍。倍。， *偶然误差（过失误差）偶然误差（过失误差）: 操作不慎等，如读错数、操作不慎等，如读错数、计算错误。计算错误。 一般通过对资料的仔细审核、校对一般通过对资料的仔细审核、校对.*随机误差随机误差: 随机因素造成，如四舍五入的小随机因素造成，如四舍五入的小数位数的取舍、仪器刻度的限制。数位数的取舍、仪器刻度的限制。 随机误差各次观测值是相互独立的；

37、所随机误差各次观测值是相互独立的；所以，以， 旬、月、年平均气温旬、月、年平均气温的随机误差一般的随机误差一般可以忽略不计；可以忽略不计；1、质量要求、质量要求 精确性：精确性：指测量得到数据的精确指测量得到数据的精确程度程度 准确性：准确性：指观测值与真值的符合程度指观测值与真值的符合程度 对对准确性的审查准确性的审查：主要检查资料中有无明显过：主要检查资料中有无明显过失误差，如发现有，应将其改正或订正，无法改正失误差，如发现有，应将其改正或订正，无法改正的，则不参加统计分析；的，则不参加统计分析； 精确性的审查：精确性的审查：主要看其记录是否达到规定的主要看其记录是否达到规定的精度，精度，

38、如降水观测应达到如降水观测应达到0.1mm0.1mm。精度不够会影响。精度不够会影响分析结果，但对观测随机误差较大的一些项目，过分析结果，但对观测随机误差较大的一些项目，过高的记录精度也没有实际意义。高的记录精度也没有实际意义。 如果测站气象记录序列仅仅是气候变如果测站气象记录序列仅仅是气候变化的反映，那么资料是均一的。化的反映，那么资料是均一的。 测站位置的迁移，周围环境的改变，测站位置的迁移，周围环境的改变，观测仪器和安装方法的更新，观测时制的观测仪器和安装方法的更新，观测时制的改变改变都可能破坏资料的均一性。分析时应都可能破坏资料的均一性。分析时应注意。注意。 测站所在位置是否对周围一定

39、区域范围内的气候具有代测站所在位置是否对周围一定区域范围内的气候具有代表性，如山谷的测站代表性差。一般说，由于地理环境的差表性，如山谷的测站代表性差。一般说，由于地理环境的差异异各个测站所能代表的范围是很不一样各个测站所能代表的范围是很不一样的。的。 气候记录的代表性气候记录的代表性是指某个测站记录能否反映我们研是指某个测站记录能否反映我们研究地区范围内的气候状况。究地区范围内的气候状况。在不同目的、不同尺度的气候在不同目的、不同尺度的气候分析中，测站代表性当然可以是不同的分析中，测站代表性当然可以是不同的。如河谷中的测。如河谷中的测站由于地形对降水的影响站由于地形对降水的影响, ,常偏大常偏

40、大, ,大范围气候分析中无大范围气候分析中无代表性代表性, ,但研究地形影响时就有代表性但研究地形影响时就有代表性. . 只要观测场地的设置符合规范要求，气象站的记录都只要观测场地的设置符合规范要求，气象站的记录都可代表类似地理环境下相当大范围的气候状况。可代表类似地理环境下相当大范围的气候状况。 不同测站进行气候时空变化分析时，不同测站进行气候时空变化分析时，资料要有可比性。如不同测站气候变化的资料要有可比性。如不同测站气候变化的空间比较时，要求空间比较时，要求各测站的资料都在同一各测站的资料都在同一时期时期；另外是否相同的观测制式，仪器性另外是否相同的观测制式，仪器性能等能等，如，如城郊对

41、比观测时城郊对比观测时，最好观测同时，最好观测同时段，使用相同的设备。段，使用相同的设备。 气候资料的审查是气候资料的审查是项复杂而细致项复杂而细致的工作。的工作。需要具备一定的气象学、天气需要具备一定的气象学、天气学、气候学和气象观测知识学、气候学和气象观测知识，更需要有，更需要有相当相当丰富的气候资料丰富的气候资料工作经验工作经验 通常将审查工作分为两类，通常将审查工作分为两类， 一一 类：类：技术性检查技术性检查 另一类：另一类：合理性检查合理性检查 技术性检查技术性检查主要从下列几方面进行：主要从下列几方面进行： (1)(1)查阅测站历史沿革和资料说明，分查阅测站历史沿革和资料说明，分

42、析析是否存在因测站迁移、仪器和观测方是否存在因测站迁移、仪器和观测方法更新、观测时制改革法更新、观测时制改革等引起的非均一等引起的非均一性。性。 (2)(2)根据观测规范、统计规定、检查观根据观测规范、统计规定、检查观测记录和统计结果测记录和统计结果是否符合规定是否符合规定、校对、校对统计计算有否错误。统计计算有否错误。 (3)(3)检查同一要素的检查同一要素的各个统计项目之间各个统计项目之间是否协调。是否协调。 合理性检查合理性检查主要是以主要是以天气学、天气学、气候学气候学知识为依知识为依据，从气象要素的据，从气象要素的时、空变化规律时、空变化规律和和各要素问的相互各要素问的相互联系规律联

43、系规律出发，分析气候资料是否合理。主要采用以出发，分析气候资料是否合理。主要采用以下几种方法：下几种方法： (1)(1)本站本站前后期资料比审法前后期资料比审法。气候变化是缓慢的、。气候变化是缓慢的、连续的。因此虽然逐年的观测值不一样，但它们应在连续的。因此虽然逐年的观测值不一样，但它们应在一个大致的水平上随机波动。如果通过前后资料的对一个大致的水平上随机波动。如果通过前后资料的对比，比，发现资料序列中存在明显的不连续变化发现资料序列中存在明显的不连续变化，则可能，则可能存在非均一性，应配合测站历史沿革情况作进一步的存在非均一性，应配合测站历史沿革情况作进一步的分析判断。分析判断。 (2)(2

44、)区域资料比审法区域资料比审法。相邻测站由于受同一天。相邻测站由于受同一天气系统的影响，常常有相当好的一致性和相关性。气系统的影响，常常有相当好的一致性和相关性。相邻测站气象要素之间的相互联系规律相邻测站气象要素之间的相互联系规律，可以成为，可以成为我们发现和订正错误记录的重要依据。我们发现和订正错误记录的重要依据。 (3) (3)气象要素相关法。气象要素相关法。各种不同的气象要素各种不同的气象要素从不同的测面描写了一地的天气气候特征以及从不同的测面描写了一地的天气气候特征以及它们之间存在着各种不同程度的相关。在实际它们之间存在着各种不同程度的相关。在实际工作中，工作中，常常用同一测站或者若干

45、个地理环境常常用同一测站或者若干个地理环境相似测站的相关密切的两个要索的观测值作成相似测站的相关密切的两个要索的观测值作成相关图进行审查相关图进行审查。当资料符合要求时，所有点。当资料符合要求时，所有点子应密集地落在一根曲线或直线附近。如果个子应密集地落在一根曲线或直线附近。如果个别点子明显偏离相关线，则这个点子的观测值别点子明显偏离相关线，则这个点子的观测值可能有明显的误差。可能有明显的误差。1 1）回归订正方法（关系密切的站）回归订正方法（关系密切的站） 两站同一要素两站同一要素: :先计算相关系数先计算相关系数, ,再建立再建立一元线性回归一元线性回归 方程方程 对周围的基本站进行逐步筛

46、选来建立对周围的基本站进行逐步筛选来建立多元回归订正公多元回归订正公 式式, ,是一个较好的途径。是一个较好的途径。 屠先生提出屠先生提出: :用基本站不同要素值作为逐步回归因子用基本站不同要素值作为逐步回归因子 建立与订正站某一要素的联系方程建立与订正站某一要素的联系方程。他用上海站。他用上海站 为基本站，用逐步回归方法对为基本站，用逐步回归方法对7 7个因子筛选，建立个因子筛选，建立 南京南京1 1月平均气温多元回归插补公式，月平均气温多元回归插补公式， 因此公式插补抗日战争期间南京缺测资料。这因此公式插补抗日战争期间南京缺测资料。这7 7个因子是上海个因子是上海1 1月平均气压、平均气温

47、、平均最高气温、月平均气压、平均气温、平均最高气温、平均最低气温、月降水量、降水数和日照百分率平均最低气温、月降水量、降水数和日照百分率。2 2）差值订正法）差值订正法（地理环境近似一致的站点，（地理环境近似一致的站点，站与站之间差值基本为常数站与站之间差值基本为常数-气温、气气温、气压）压）3 3）比值订正法）比值订正法 两台站降水比值为准常数两台站降水比值为准常数* *订正适当性的论证订正适当性的论证 * 只有在基本站与订正站距离不很远，只有在基本站与订正站距离不很远，两站气象要素相关相当密切，且平行观测两站气象要素相关相当密切，且平行观测时期比较长时才用回归订正法进行序列延时期比较长时才

48、用回归订正法进行序列延长。如果平行观测时期不长，通常还是以长。如果平行观测时期不长，通常还是以采用计算简便的差值订正法或比值订正法采用计算简便的差值订正法或比值订正法为好。为好。1.理解并掌握统计量平均值、距平、标准差、协理解并掌握统计量平均值、距平、标准差、协方差、变率、变差系数的概念，特别是其在气方差、变率、变差系数的概念，特别是其在气象上的意义。象上的意义。2.思考如何求出北半球思考如何求出北半球1980-2010年逐月年逐月850hPa纬向风场的气候场、异常场和均方差场纬向风场的气候场、异常场和均方差场？3.理解总体、样本、频率表、分布列的概念。理解总体、样本、频率表、分布列的概念。4.何谓中心化、标准化、正态化，这样做资料处何谓中心化、标准化、正态化，这样做资料处理的必要性为何理的必要性为何?如何处理？如何处理？5.如何整理区域资料。如何整理区域资料。