《气象统计方法课程实习》学生姓名 x x学号 xxxxxxxxxxxx院系大气科学专业大气科学任课教师 x x二Ｏ一四年十二月二十日实习一求500hPa高度场气候场、距平场和均方差场（1）气候场二月份高度场的气候场呈现南高北低的状态，陆地上的高度场比较稀疏，而在西太平洋上高度场比较密集。
7月份高度场的气候场总体呈现东高西低的状态，在印度半岛出现低压中心，而在赤道西太平洋地区出现高压中心，位置在130°E，25°N附近。
35°N以北高度分布很密集，而35°N以南比较稀疏。
（2）距平场1982年5月距平场在我国华东地区出现负距平,在亚洲西南部也出现低压中心，在青藏高原处为正距平。
1984年4月距平场在日本东部海洋地区形成低压中心，印度半岛的西部有一低压中心，在35°N-40°N基本都为正距平。
（3）均方差场三月份高度的均方差场整体呈现南小北大的状态。
说明低纬地区高度的波动幅度比较小，而中高纬地区高度的波动比较大。
在太平洋北部波动最大。
十月份高度的均方差场在西太平洋有极大值，其余地区波动都较小。
实习二计算给定数据资料的简单相关系数和自相关系数单相关系数和自相关系数程序：program mainparameter(n=20,m=10)integer i,j,t,max1,max2real r,s1,s2real a(n),b(n),ano1(n),ano2(n),bzh1(n),bzh2(n),r1(m),r2(m)real ave1,ave2,sum12,sum11,sum22data a/,,,,,,,,,,,,,,,,,,,data b/,,,,,,,,,,,,,,,,,,,!求平均ave1=ave2=do i=1,nave1=ave1+a(i)ave2=ave2+b(i)enddoave1=ave1/nave2=ave2/n!求距平ano1(n)=ano2(n)=do i=1,nano1(i)=a(i)-ave1ano2(i)=b(i)-ave2Enddo!求标准差s1=s2=do i=1,ns1=s1+ano1(i)*ano1(i)s2=s2+ano2(i)*ano2(i)enddos1=sqrt(s1/n)s2=sqrt(s2/n)!标准化bzh1(n)=bzh2(n)=do i=1,nbzh1(i)=ano1(i)/s1bzh2(i)=ano2(i)/s2enddo!求相关系数sum12=sum11=sum22=do i=1,nsum12=sum12+ano1(i)*ano2(i)sum11=sum11+ano1(i)*ano1(i)sum22=sum22+ano2(i)*ano2(i)enddor=sum12/sqrt(sum11*sum22)print*print*,'中国1970-1989年年平均和冬季平均气温的相关系数为r=',r print*!求自相关系数r1(m)=r2(m)=do t=1,mdo j=1,n-tr1(t)=bzh1(j)*bzh1(j+t)+r1(t)r2(t)=bzh2(j)*bzh2(j+t)+r2(t)enddor1(t)=r1(t)/(n-t)r2(t)=r2(t)/(n-t)enddo!比较自相关系数绝对值大小max1=1max2=1do t=2,mif(abs(r1(t))>abs(r1(max1)))max1=tif(abs(r2(t))>abs(r2(max2)))max2=tenddoprint*,'年平均气温自相关系数绝对值最大的滞后时间长度t为:',max1,r1(max1) print*print*,'冬季平均气温自相关系数绝对值最大的滞后时间长度t 为:',max2,r2(max2)print*end分析：中国1970-1989年年平均和冬季平均气温相关系数为，为正相关;年平均气温自相关系数绝对值最大的滞后时间长度为7，自相关系数为负，呈负相关;冬季平均气温自相关系数绝对值最大的滞后时间长度为4，自相关系数为负，呈负相关实习三（附加）计算给定数据的落后交叉相关系数和偏相关系数程序：program mainparameter(n=30,m=10)integer i,j,treal ave1,ave2,ave3,r12,r13,r23,ry1,ry2,ry3real a(n),b(n),c(n),ano1(n),ano2(n),ano3(n),bzh1(n),bzh2(n),bzh3(n) real rt12(m),rt13(m)!a-12月；b-1月；c-2月--(30个数据)data a/,,,,,,,,,,,,,,,,,,,,,,,,,,,,,data b/,,,,,,,,,,,,,,,,,,,,,,,,,,,,,data c/,,,,,,,,,,,,,,,,,,,,,,,,,,,,,!求平均ave1=ave2=ave3=do i=1,nave1=ave1+a(i)ave2=ave2+b(i)ave3=ave3+c(i)enddoave1=ave1/nave2=ave2/nave3=ave3/n!求距平ano1(n)=ano2(n)=ano3(n)=do i=1,nano1(i)=a(i)-ave1ano2(i)=b(i)-ave2ano3(i)=c(i)-ave3enddo!求标准差s1=s2=s3=do i=1,ns1=s1+ano1(i)*ano1(i)s2=s2+ano2(i)*ano2(i)s3=s3+ano3(i)*ano3(i)enddos1=sqrt(s1/n)s2=sqrt(s2/n)s3=sqrt(s3/n)!标准化bzh1(n)=bzh2(n)=bzh3(n)=do i=1,nbzh1(i)=ano1(i)/s1bzh2(i)=ano2(i)/s2bzh3(i)=ano3(i)/s3enddo!求落后交叉相关系数（滞后长度τ最大取10）12月与1月rt12；12月与2月rt13 rt12(m)=rt13(m)=do t=1,mdo i=1,n-trt13(t)=bzh1(i)*bzh3(i+t)+rt13(t)rt12(t)=bzh1(i)*bzh2(i+t)+rt12(t)enddort12(t)=rt12(t)/(n-t)rt13(t)=rt13(t)/(n-t)enddoprint*,'12月气温与1月气温的落后交叉相关系数依次为(1-10年):'print '',rt12print*print*,'12月气温与2月气温的落后交叉相关系数依次为(1-10年):'print '',rt13print*!求相关系数，12月和1月r12；12月和2月r13，1月和2月r23r12=r13=r23=do i=1,nr12=r12+bzh1(i)*bzh2(i)r13=r13+bzh1(i)*bzh3(i)r23=r23+bzh2(i)*bzh3(i)enddor12=r12/nr13=r13/nr23=r23/n!求偏相关系数，12月和2月(消除1月)ry1；1月和2月(消除12月)ry2；12月和1月(消除2月)ry3ry1=(r13-r12*r23)/sqrt((1-r23*r23)*(1-r12*r12))ry2=(r23-r12*r13)/sqrt((1-r13*r13)*(1-r12*r12))ry3=(r12-r13*r23)/sqrt((1-r23*r23)*(1-r13*r13))print*,'消除1月影响，12月与2月气温的偏相关系数:',ry1print*print '(a,',' 消除12月影响，1月与2月气温的偏相关系数:',ry2print*print*,'消除2月影响，12月与1月气温的偏相关系数:',ry3print*end分析：消除1月影响，12月与2月气温的偏相关系数为正，呈正相关;消除12月影响，1月与2月气温的偏相关系数为正，呈正相关;消除2月影响，12月与1月气温的偏相关系数为正，呈正相关实习四求给定数据的一元线性回归方程程序program mainparameter(n=20)integer i! x为环流指标（预报因子），y为气温（预报量）real x(n),y(n)real ave1,ave2,s12,s1,s2,b,b0,r,Fdata x/32,25,20,26,27,24,28,24,15,16,24,30,22,30,24,33,26,20,32,35/ data y/,,,,,,,0,,,,,,,,,,,,!求平均ave1=ave2=do i=1,nave1=ave1+x(i)ave2=ave2+y(i)enddoave1=ave1/nave2=ave2/n!求协方差、预报因子\预报量的方差s12=s1=s2=do i=1,ns12=s12+(x(i)-ave1)*(y(i)-ave2)s1=s1+(x(i)-ave1)*(x(i)-ave1)s2=s2+(y(i)-ave2)*(y(i)-ave2)enddos12=s12/ns1=s1/ns2=s2/n!求b，b0b=s12/s1b0=ave2-b*ave1!求回归方程print*,'气温和环流指标之间的一元线性回归方程为：'print'(a,,,a)',' y=',b0,b,'x'print*!检验Fr=sqrt(s1/s2)*bF=r*r/((1-r*r)/(n-2))print'(a,',' F =',Fend分析：F=>Fα=，回归方程显着实习五（附加）求给定数据的多元线性回归方程实习六（附加）分析中国夏季降水线性趋势的分布特征程序：program mainparameter(m=160,n=25)integer i,t(n),avetinteger sta(m) !站号real lon(m),lat(m),f(m,n) !经，维，记录real ave(m),ano(m,n),anot(n)real b(m),sxy(m),streal timeinteger level1000 format(3a,25i)2000 format!读数据open(5,file='d:\qxtj\6\')read(5,1000)do i=1,mread(5,*),sta(i),lon(i),lat(i),(f(i,j),j=1,n) enddo!计算数据平均，距平，得到距平数组ano(m,n)ave(m)=ano(m,n)=do i=1,mdo j=1,nave(i)=ave(i)+f(i,j)end doave(i)=ave(i)/ndo j=1,nano(i,j)=f(i,j)-ave(i)end doenddo!计算时间距平anot(n)t(n)=0avet=0anot(n)=0do i=1,nt(i)=1981+iavet=t(i)+avetenddoavet=avet/ndo i=1,nanot(i)=t(i)-avetenddo!计算b(m)(160个)b(m)=0sxy(m)=st=do j=1,nst=anot(j)*anot(j)+stEnddodo i=1,mdo j=1,nsxy(i)=ano(i,j)*anot(j)+sxy(i)enddob(i)=sxy(i)/stenddoprint*print*,'160站夏季降水线倾向率:'print '',(b(i),i=1,m)print*End分析：b(m)为正时，降水有随时间增多的趋势；b(m)为负时，降水有随时间减小的趋势实习七计算给定数据的11年滑动平均和累积距平程序：program mainparameter(n=85,k=11,nyear=1922)real dat(n),ano(n),h(n-k+1),l(n)real ave1000 format2000 format3000 format4000 format!读文件open(5,file='d:\qxtj\7\')do i=1,nread(5,*) dat(i)enddo!求距平ave=ano(n)=do i=1,nave=ave+dat(i)enddoave=ave/ndo i=1,nano(i)=dat(i)-aveenddo!滑动平均h(n-k+1)h(n-k+1)=0do i=1,n-k+1do j=i,i-1+kh(i)=h(i)+dat(j)enddoh(i)=h(i)/kenddo!累计距平l(n)l(n)=do i=1,ndo j=1,il(i)=l(i)+ano(j)enddoEnddo!输出open(6,file='d:\qxtj\7\')open(7,file='d:\qxtj\7\')write(6,1000) (h(i),i=1,n-k+1)write(7,2000) (l(i),i=1,n)close(6)close(7)write(*,'("11年滑动距平为")')write(*,3000) (h(i),i=1,n-k+1)print*write(*,'("累计距平为")')write(*,4000) (l(i),i=1,n)print*End分析：数据从1922年到2006年共85年。
