11 气象统计方法课程实践内容2013

实习二 计算给定数据资料的简单相关系数和自相关系数

1、资料介绍

根据下表中年平均气温和冬季平均气温的等级数据进行下列计算： 1）计算两个气温之间的简单相关系数。

2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。（滞后长度τ最大取10）

2、要求

实习报告中附出简单相关系数或自相关系数程序。

答案：r=0.47

年平均气温在滞后长度j=3、冬季序列在j=4最大。

第 11 页 共 36 页

12 气象统计方法课程实践内容2013

3、实习结果

（1）、Fortran源程序

PROGRAM EXAM

IMPLICIT NONE

INTEGER,PARAMETER::N=20 INTEGER i,j,k,ty,tw,tyw

REAL::avr_y=0,avr_w=0,sy=0,sw=0,rxy=0,max_y=0,max_w=0,max_yw=0 REAL y(N),w(N)

DATA y/3.4,3.3,3.2,2.9,3.4,2.8,3.6,3.0,2.8,3.0,3.1,3.0,2.9,2.7,3.5,3.2,3.1,2.8,2.9,2.9/ DATA

w/3.24,3.14,3.26,2.38,3.32,2.71,2.84,3.94,2.75,1.83,2.80,2.81,2.63,3.20,3.60,3.40,3.07,1.87,2.63,2.47/

REAL syy(N),sww(N),r(N),rty(N),rtw(N),rtyw(N),rxy_ty(N),rxy_tw(N),rxy_tyw(N)

!求两数组平均值 DO i=1,N

avr_y=avr_y+y(i) avr_w=avr_w+w(i) END DO

avr_y=avr_y/N avr_w=avr_w/N !简单相关系数 DO j=1,N

syy(j)=(y(j)-avr_y)**2 sy=sy+syy(j)

sww(j)=(w(j)-avr_w)**2 sw=sw+sww(j) END DO

sy=sqrt(sy/N) sw=sqrt(sw/N) DO j=1,N

r(j)=((y(j)-avr_y)/sy)*((w(j)-avr_w)/sw) rxy=rxy+r(j) END DO rxy=rxy/N

PRINT \年全年平均气温与冬季平均气温的简单相关系数rxy=',f5.2)\k=0

!自相关系数 DO ty=1,N/2 DO i=1,N-ty

第 12 页 共 36 页

13 气象统计方法课程实践内容2013

rty(i)=((y(i)-avr_y)/sy)*((y(i+ty)-avr_y)/sy)

rxy_ty(ty)=rxy_ty(ty)+rty(i) END DO

rxy_ty(ty)=rxy_ty(ty)/(N-ty) rxy_ty(ty)=ABS(rxy_ty(ty)) IF(rxy_ty(ty)>max_y) THEN max_y=rxy_ty(ty) k=ty END IF END DO

PRINT \全年平均气温绝对值最大自相关系数rxy_ty=',f7.4,/,'滞后时间长度 k=',I2)\k=0

DO tw=1,N/2 DO i=1,N-tw

rtw(i)=((w(i)-avr_w)/sw)*((w(i+tw)-avr_w)/sw) rxy_tw(tw)=rxy_tw(tw)+rtw(i) END DO

rxy_tw(tw)=rxy_tw(tw)/(N-tw) rxy_tw(tw)=ABS(rxy_tw(tw)) IF(rxy_tw(tw)>max_w) THEN max_w=rxy_tw(tw) k=tw END IF END DO

PRINT \冬季平均气温绝对值最大自相关系数rxy_tw=',f7.4,/,'滞后时间长度 k=',I2)\k=0

!落后交叉相关系数 DO tyw=1,N/2 DO i=1,N-tyw

rtyw(i)=((y(i)-avr_y)/sy)*((w(i+tyw)-avr_w)/sw) rxy_tyw(tyw)=rxy_tyw(tyw)+rtyw(i) END DO

rxy_tyw(tyw)=rxy_tyw(tyw)/(N-tyw) rxy_tyw(tyw)=ABS(rxy_tyw(tyw)) IF(rxy_tyw(tyw)>max_yw) THEN max_yw=rxy_tyw(tyw) k=tyw END IF END DO

PRINT \全年平均温度与冬季平均气温之间的落后交叉相关系数rxy_tyw=',f7.4,/,'滞后时间长度 k=',I2)\

第 13 页 共 36 页

14 气象统计方法课程实践内容2013

END

（2）、程序运行结果：

第 14 页 共 36 页

15 气象统计方法课程实践内容2013

实习三 分析中国夏季降水线性趋势的分布特征

1.资料介绍及要求：

利用数据160zhan-rainfall-summer.txt，编写求1982-2006年中国160站各站夏季降水线性倾向率，给出分布图，并进行简单分析。read-rain.for给出了阅读资料的fortran程序。数据在文件夹中单独给出。 2.实习结果

（1）.matlab程序

%编写求1982-2006年中国160站各站夏季降水线性倾向率

clear all clc

fid=fopen('E:/160zhan-rainfall-summer.txt','rt'); tline=fgets(fid);

data1=fscanf(fid,'%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f',[28,160]); data2=data1'; fclose(fid); for i=1:160;

j(i,1:25)=data2(i,4:28); n1=1982:1:2006;

pp(i,:)=polyfit(n1,j(i,1:25),1); end

b=pp(:,1);

jd=data2(:,3); wd=data2(:,2); jdc=75:0.5:135; wdc=18:.5:55;

bz=griddata(jd,wd,b,jdc,wdc','cubic'); c=contour(jdc,wdc,bz)

第 15 页 共 36 页