气象统计 谱分析课件.pptVIP

谱分析;Spectrum Analysis; 尽管上述变换只是把所研究的资料按频率排列，从表面上看似乎并没有增加什么新东西。但实质上，正因为按频率排列，将物理量随时间的变化作为波动来处理，可以研究其能量的周期分布，并对不同记录进行比较。 实际上大多数大气物理现象在物理上都表现出依赖频率的形式，因而谱表示的意义是极为明显的。;谱表示的优点：;谱分析; 本节侧重在时间序列波动分析上．空间序列的波动分析完全是类似的。;谱的概念;谱分析;ck称为x(t)的复谱．因为ck 是一个复数，可以把ck 表示成振幅 Ak 和位相?k 的形式 Ak 称为振幅谱． ?k 称为位相谱，并可用k 为横坐标， Ak 或位相?k 为纵坐标作图,这种图称为振幅谱或位相谱图．由于k 是整数，谱图上的曲线是不连续的，称这种谱为离散谱或线谱．;如果令周期T为无限大，则有;Power Spectrum; 为功率谱密度;功率谱估计;以波数为横轴，离散功率谱估计值 为纵坐标作谱图．;(2)连续功率谱估计;连续功率谱估计;利用傅里叶反变换得到：;从上式可见，连续功率谱函数在正、负半轴上对称．因而类似离散谱曲作法，只绘出正半轴的连续功率谱即可．;求连续功率谱估计方法步骤： ①计算样本落后自相关系数r(?)(?＝0,1,2,…m);m为最大步长，或最大落后时间长度． ②求粗谱估计，利用近似积分中的梯形公式，变为有限项求和公式;③ 计算平滑功率谱，为消除粗谱估计的抽样误差，还要对粗谱估计作平滑处理，作为功率谱最后估计．常用平滑公式为;最后得到谱估计公式为：; ④作谱图，以波数 l 为横轴，平滑功率谱密度估计值Sl为纵坐标作图．在谱图横轴上，也可同时标上对应波数的周期或频率，以备分析时用．波数 l 与周期 T 关系为; 气象中，对已知序列样本容量为n的情况下，功率谱估计随抽样点m不同而不同，m愈大用来估计谱的采样点就愈多，但并不表明估计就准确些，因为m太大时，相关系数的计算由于样本容量太小，估计自相关函数就差，功率谱估计亦带来较大误差．通常m宜选择在n／3 ～ n／10之间．;利用功率谱作周期分析; 由于功率谱图已经显示不同频率振动的功率大小，同时也是方差贡献的大小，因而可以从谱曲线中的谱值最大来确认主要振动及其对应的周期．但这一确认是否有统计意义还需作显著性检验．;离散功率谱的显著性检验;k;连续功率谱检验;白色噪音过程的功率谱;滤 波;滤波种类; 一种方法可用两个简单的低通滤波器(或高通滤波器)来构成一个带通滤波器．例如先作5年滑动平均，除去5年以下的波，再把原序列做9年滑动平均，除去9年以下的波，然后把两个过滤后的序列相减，则可得到5～9年周期波动的序列．;交 叉 谱;称实部谱为协谱，虚部谱为正交谱．;交叉谱的估计;凝聚谱与位相差谱;位相差谱可用交叉谱的实部与虚部谱求得;对凝聚值的显著性检验，可在两时间函数振动的凝聚期望为0的原假设下，使用统计量;End Of the Class