成信工统计天气预报实验指导03功率谱分析和滤波.docVIP

实验三 功率谱分析和滤波 一、目的和要求： 周期是大气运动的基本理念之一，周期分析也是统计天气分析中的基础理论。功率谱分析的模型和滤波思想是最基本的周期分析统计模型，是必须理解和掌握的基本理论方法，是后续课程中许多天气分析事实理解的基础，也是毕业设计和论文中的基本分析方法。该方法作为天气分析模型，在实际业务工作和科研工作中有极强的实用意义。通过该实验，能使学生深刻理解气象时间序列周期分析的意义和周期分析的方法，为实际的业务和科研工作打下一定的基础。 二、实验的主要内容： 对原序列滤波前（或后）利用功率谱检测序列的周期，并进行显著性检验。 本实验利用功率谱方法检测原序列的周期，并进行显著性检验；利用带通滤波器提取所需要的周期分量，再利用功率谱检验所提取周期的显著性。 三、步骤： 3.1 熟悉资料方法 3.1.1 资料 本次实验使用国家气候中心提供 58年（1951年-2008年）夏季（6-8月）逐月降水资料。取得是长江中下游及江淮地区降水场第一主要因子的时间序列为例。 3.1.2 功率谱分析方法 谱分析是时间序列在频域上进行分析的方法，亦称频谱分析或波谱分析。主要侧重在时间序列波动分析上。 功率谱的概念是针对功率有限信号的，所表现的是单位频带内信号功率随频率的变换情况。功率谱密度是指对于具有连续频谱和有限平均功率的信号或噪声，表示其频谱分量的单位带宽功率的频率函数。气象上常利用功率谱作周期分析，功率谱图显示不同频率振动的功率大小，同时也是方差贡献的大小，因而可以从谱曲线中的谱值最大来确认主要振动及其对应的周期。 功率谱估计是指对一实测时间序列现实（t=1，2，…，n)的估计，可有两种估计功率谱方法：离散功率谱估计和连续功率谱估计，本实验选择连续功率谱估计。具体求连续功率谱估计方法步骤如下： (1) 计算样本落后自相关系数 ，m为最大步长，或最大落后时间长度，本实验m取18，通常m宜选择在n/3～n/10之间。 （为时间落后步长） (2) 求粗谱估计 波数与圆频率的关系在取样点数为m时有： 波数k最大可取到m/2。习惯上令l=2k，波数l取值点可从0到m。 (3) 计算平滑功率谱 也可 （式中l=0，1，…，m） (4) 作谱图，以波数l为横轴，平滑功率谱密度估计值为纵坐标作图。波数l与周期关系为 连续功率谱检验包括红色噪音过程和白色噪音过程，本实验选择白色噪音过程。 3.1.3 滤波方法 对气象要素作谱分析过程中，一些规则周期占有很大的分量，但这种周期是众所周知的，它的存在，压低了其它周期的表现，一旦把它去掉之后，则可突出地表现其它周期的成分，即把感兴趣的周期成分从原来的序列中识别和提取出来，这一过程就是滤波过程。实际上是原始序列经一定的变换转化为另一序列的过程。 滤波的总体过程： 输入 输出 x(t) 滤波器 y(t) （数字处理） 其中x(t)为原始序列，y(t)为新的时间序列。 要求：过滤系统具有时间不变性和稳定性。 滤波器可分为三类：用于从原序列中滤出低于某个频率成分的滤波器，称为低通滤波器；用于滤出高于某个频率成分的滤波器称为高通滤波器；滤出两个频率之间一个频带成分的滤波器称为带通滤波器。根据本实验的具体情况，选择带通滤波器，滤波范围为2～8年。 3.2 编写程序 要求编写主程序，其中包括资料读入，范围截取，子程序调用。 Cccccccccccccccccccccccccccc（附程序，对关键部分标志出）ccccccccccccccccccccccccc 连续功率谱程序 parameter(n=50,m=n/3-1,ta=1.645,xx=12.59,pi=3.141592653589793) real x(n),r(0:m),ss(0:m),sr(0:m),sss(0:m),ssss(0:m),vv,ssr,sum,rc tta,ave open(1,file=soi-1fil.txt) open(3,file=soi-1filsl1.txt) open(5,file=soi-1filtest1.txt) do i=1,n read(1,*) x(i) end do close(1) **********************求方差***************** ex=0.0 do i=1,n ex=ex+x(i) end do ex=e