气象统计实习报告专业：大气科学班级：xxxx级x班学号：2012130xxxx姓名：*实习一求500hPa高度场气候场、距平场和均方差场一、实习结果1981.1距平场1981.1 500hpa高度场在欧亚大陆为正距平，在印度洋和太平洋为负距平。

1981.1气候场1982.1的气候场呈现明显的沿纬度的上升而下降的趋势。

1981.1均方差场1981.1的均方差场在欧亚大陆的南侧有最小值，在10N~20N的南侧岁纬度而降低，在其北侧随纬度升高。

二、相关的fortran程序，gs文件（1）Fortran程序:program mainparameter(nx=37,ny=17,mo=12,yr=4)real var(nx,ny,mo,4)!数据real vars(nx,ny,mo)!4年气候态real jp(nx,ny, mo,4)!距平real fc(nx,ny,mo)!方差real jfc(nx,ny,mo)!均方差integer i,j,m,y,irecreal::summ=0.0!*********************求均方差*********************do m=1,12do i=1,37do j=1,17do y=1,4summ=summ+jp(i,j,m,y)**2enddofc(i,j,m)=summ/4!方差jfc(i,j,m)=sqrt(fc(i,j,m))!均方差summ =0.0enddoenddoenddo!******************************************open(11,file='d:\shixi\vars.grd',form='binary' )open(12,file='d:\shixi\jp.grd',form='binary' )open(13,file='d:\shixi\jfc.grd',form='binary' )do m=1,12write(11)((vars(i,j,m),i=1,nx),j=1,ny)enddodo m=1,12write(13)((jfc(i,j,m),i=1,nx),j=1,ny)enddodo y=1,4do m=1,12write(12)((jp(i,j,m,y),i=1,nx),j=1,ny)enddoenddoclose(11)close(12)close(13)End（2）Gs程序1.距平场'open e:\anomaly.ctl''set gxout shaded''d xiaobao''set gxout contour''d xiaobao''draw title 1982年1月''enable print e:\anomaly.gmf''print''disable print'2.气候场'open e:\climate.ctl''set gxout shaded''d xiaobao''set gxout contour''d xiaobao''draw title 1982.1''enable print e:\climate.gmf''print''disable print'3.均方差场'open e:\deviation.ctl''set gxout shaded''d xiaolu''set gxout contour''d xiaolu''draw title 1982.1 deviation''enable print e:\deviation.gmf''print''disable print'*实习二计算给定数据资料的简单相关系数和自相关系数根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。

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

（滞后长度τ最大取10）一、实习结果（1）计算简单相关系数计算出相关系数为r=0.4685170（2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。

（滞后长度τ最大取10）可以知道，年平均气温在滞后长度j=7，冬季j=4最大二、相关的fortran程序（部分）implicit nonereal x(20)!年平均气温real y(20)!冬季平均气温real :: zx(10)=(/0,0,0,0,0,0,0,0,0,0/)real :: rzx(10)=(/0,0,0,0,0,0,0,0,0,0/)real ::s=0.0 !协方差real ::jx=0.0real ::jy=0.0real ::jfx=0.0real ::jfy=0.0,rinteger i,t,jreal ::m=0.0data x/3.40,3.30,3.20,2.90,3.40,2.80,3.60,3.00,2.80,3.00,&&3.10,3.00,2.90,2.70,3.50,3.20,3.10,2.80,2.90,2.90/data y/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/do i=1,20jx=jx+x(i)jy=jy+y(i)end dodo i=1,20s=s+(x(i)-jx/20)*(y(i)-jy/20)jfx=jfx+(x(i)-jx/20)**2jfy=jfy+(y(i)-jy/20)**2end dor=s/20/sqrt(jfx/20*jfy/20)print *,"r=",rdo t=1,10do i=1,20-tzx(t)=zx(t)+(x(i)-jx/20)*(x(i+t)-jx/20)end dorzx(t)=zx(t)/(20-t)/(jfx/20)print *,t,rzx(t)if(abs(rzx(t))>m) thenm=abs(rzx(t))j=tend ifend doprint *,'年平均温度的自相关系数绝对值最大的滞后时间长度'print *,j实习三计算给定数据的落后交叉相关系数和偏相关系数根据下表北京冬季（12月~2月）气温资料计算：12月气温与1月和2月气温的落后交叉相关系数（滞后长度τ最大取10）和偏相关系数。

在实习报告中给出程序。

一、实验结果偏相关系数计算出的 12 月气温与 1 月气温的偏相关系数为 0.327，12 月气温与 2 月气温的偏相关系数为 0.290。

二、相关的fortran程序Funtion area(x,y)implicit nonereal x(30)real y(30)real ::jx=0.0real ::jy=0.0real ::jfx=0.0real ::jfy=0.0,rinteger i,real ::s=0.0do i=1,30jx=jx+x(i)jy=jy+y(i)end dodo i=1,30s=s+(x(i)-jx/30)*(y(i)-jy/30)jfx=jfx+(x(i)-jx/30)**2jfy=jfy+(y(i)-jy/30)**2end dor=s/30/sqrt(jfx/30*jfy/30)print *,"r=",rend funtion area*实习四求给定数据的一元线性回归方程利用下表数据，以环流指标为预报因子，气温为预报量，计算气温和环流指标之间的一元线性回归方程，并对回归方程进行检验。

一、实习结果（1）用excel制作的气温-环流的医院线性回归方程回归方程为：ˆ7.5-0.23 y x（2）回归方程的检验检验结果：F=20.18>Fα=4.41，回归方程显著二、Fortran程序（部分）（1）回归方程的检验!回归方程显著性检验!计算两数组的距平及均方差do i=1,mvar(i)=0do j=1,nxdiff(j,i)=dat(j,i)-ave(i)var(i)=var(i)+diff(j,i)**2end dovar(i)=sqrt(var(i)/nx)end do!计算协方差i=1;t=0do j=1,20t=t+diff(j,i+1)*diff(j,i)end doE=t/20!计算相关系数r=E/(var(1)*var(2))print*,'the relative value r is:'， rF=r**2*(nx-2)/(1-r**2)实习五 求给定数据的多元线性回归方程说明：x 1-x 4为四个预报因子，y 为预报量；样本个数n =13 要求：选取预报因子1、2、4，求预报量的标准化回归方程。

一、用excel 制作的医院线性回归方程标准化变量回归方程：124ˆ =0.5679 +0.43230.2613yx x x*实习七计算给定数据的11年滑动平均和累积距平利用数据ma.dat，编写11点滑动平均的程序，ma.for给出了阅读资料的fortran程序。

数据在文件夹中单独给出。

要求：实习报告中附出程序，并给出原数据和滑动后数据的图形（1张图）和累积距平数据图形（1张图）一、实习结果累积距平数据图形滑动后数据二、部分fortran程序：program ma! dimension x(1000),x1(1000),nny1(1000)! write(*,10)! 10format(5x,'n=?,ih=?,nyear=?')! Read(*,*)n,ih,nyear! **********************************************! * n: sample size of the time series *! * ih:moving length *! * nyear: first year of the series *! * x(n): oroginal time series *! * x1(n-ih+1): moved series *! **********************************************integer i,n,ih,nyearparameter (n=85,ih=11,nyear=1922)30real x(n),x1(n-ih+1)open(2,file='d:\ma.dat ')read(2,*)(x(i),i=1,n)close(2)do i=1,n-ih+1x1(i)=sum(x(i:i+ih-1))/ihend doprint*,x1open(10,file='d:\moveaverage.txt')write(10,'(f10.6)')(x1(i),i=1,n-ih+1)close(10)end*实习八对给定的海温数据进行EOF分析给出海表温度距平数据资料sstpx.grd，以及相应的数据描述文件sstpx.ctl，对其进行EOF分析，资料的时空范围可以根据sstpx.ctl获知。