成信工统计天气预报课件第9章_谱分析.pptVIP

3·连续的功率谱 称为功率谱密度 或证为 三·功率谱的估计 1·离散功率谱的估计 对于实现时间序列现实 （1）计算不同波数K的傅氏系数估计 其中K的取值为 [ ] 表示取整数部分 1·离散功率谱的估计 （2）计算不同波数K的功率谱值 （3） 波数为横轴，离散功率谱估计值 为纵坐标做谱图。 在横坐标上通常也标上对应的振动周期值或频率值。 2·连续功率谱的估计 通过时间函数的自相关函数作间接估计 （1）自相关函数表示的连续功率谱 为标准化时间函数 即数学期望为0 方差为1 则自相关函数 2·连续功率谱的估计 利用傅里叶反变换 展开有 为奇函数，具有反对称性 为偶函数 ,具有对称性，加上平稳随机过程的自相关函数，也是偶函数，也有对称性。上式容易变成 可见，连续功率谱函数在正负半轴上对称，因而类似离散谱做法，只给出正半轴的连续功率谱即可。 2·连续功率谱的估计 （2）连续功率谱的估计步骤 1）计算样本落后自相关系数 为最大步长，或最大落后时间长度。 2）求粗谱估计 利用近似积分中的梯形公式 K最大可取 令 3）计算平滑功率谱 2·连续功率谱的估计 可用平滑公式 设计一个窗函数 则谱估计公式可写为 另外，为保持权重一致，端点也还要乘上0.5，在此，计算平滑功率谱密度估计公式 式中 2·连续功率谱的估计 4）做谱图 以波数 为横轴，平滑功率谱密度估计值 为纵坐标作图 其中波数 与周期关系 四·利用功率谱做周期分析 1·利用功率谱图显示的不同频率振动的功率 大小就可以确认主要振动及其对应的周期 但是，由此确认的周期是否有统计意义，还需要做显著性检验。 四·利用功率谱做周期分析 2·离散功率谱检验 ∵ 的总方差可以分解为 个波动方差之和 ∴可利用方差检验的方法 假设 统计量 其中 为原序列方差 显著 反之不显著 四·利用功率谱做周期分析 3·连续功率谱检验 （1）连续功率谱检验基本思路 与非周期性随机过程比较来进行 非周期性随机过程：一般指两种过程，白色噪音过程和红色噪音过程。 这里的比较过程可通过谱曲线进行。 四·利用功率谱做周期分析 （2）红色噪音过程的功率谱 红色噪音过程即是一阶马尔柯夫过程或一个自回归模型。 ∴红色噪音功率谱为 四·利用功率谱做周期分析 上面求和式展开可由三部分组成。 当 有 当 有 当 有 四·利用功率谱做周期分析 将上三部分带回原式 四·利用功率谱做周期分析 可见，红噪音谱在 有其最大值，即只有线性直线趋势而无其他周期振荡。 （3）白色噪音过程功率谱 白噪音：自协方差 可见，白色噪音过程是另一种周期过程 （4）显著性检验 五 滤波 1·问题的提出 2·滤波方法的基本思想 3·低通过滤 4·高通过滤 5·带通过滤 1·问题的提出 在前述的谱和功率谱的讨论中，我们可以概括功率谱图上现实的不同频率振动的功率大小,（也即方差贡献的大小）即谱值最大来确认主要振动及其相应的主要周期；并通过显著性检验来确认显著性主要周期的存在，从而获得时间序列中的主要 周期信息。 然而，在气象要素的谱分析过程中，一些规则周期占有很大的分量，例如月均序列分析中，年变化的周期一定是主要周期，当这种周期确实众所周知的，它的存在，压低了其他周期的表示，一旦把它去掉，则可突出的表现其他周期的成分，这在实际的气象分析和研究工作中都是很有意义的需求。 1·问题的提出 那么如何突出提取和某些需要的周期成分？ 解决办法:滤波 （1）含义：去掉某些周期成分，突出所需的周期成分 （2）方法：实际上是将原来序列经过一定的变换转化为另一序列（的过程） 这个过程通常称为滤波过程。 2·滤波方法的基本思想 要求:过滤器系统具有时间不变性和稳定性，如果用一脉冲函数 作为输入，其输出记为 称为脉冲响应。 对任意输入函数 ，可表示成无穷脉冲的 积分形式 根据线性和时间不变性的要