搜索
气象统计方法实习报告注：本人上传实习报告，并不想为大家抄袭提供方便，只是想给大家一个参考和验证，我认为每一个用心之人都会做得比我更好，本文档只起抛砖之效。专业：大气科学班级：2010级x班学号：xxxxxxxxxx姓名：xx1目录实习一气候场、距平场和均方差场的计算...............................2实习二相关系数及自相关系数的计算.......................................6实习三一元线性回归...................................................................9实习四滑动平均.........................................................................10实习五多元线性回归.................................................................12实习六ＥＯＦ分析.....................................................................15实习心得.......................................................................................262实习一气候场、距平场和均方差场的计算要求：有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N.时段：1982.1~1985.12共48个月。水平分辨率：2.5*2.5，格点数：37*17。编fortran程序，求500hPa高度场的（1）气候场；（2）距平场；（3）均方差场。并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单说明。注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.编写的fortran程序如下：programoneimplicitnoneintegernx,ny,mo,yrparameter(nx=37,ny=17,mo=12,yr=4)realdatas(nx,ny,mo,yr)reala(nx,ny,mo),aa(nx,ny,mo,yr),bb(nx,ny,mo)integeri,j,t,it,ii,mi,irecopen(5,file='e:\h500.dat')doii=1,4domi=1,12read(5,1000)read(5,3000)((datas(i,j,mi,ii),i=1,nx),j=1,ny)enddoenddoclose(5)!气候场dot=1,12doj=1,nydoi=1,nxa(i,j,t)=0doit=1,4a(i,j,t)=a(i,j,t)+datas(i,j,t,it)enddo3a(i,j,t)=a(i,j,t)/4enddoenddoenddo!求距平场dot=1,12doj=1,nydoi=1,nxaa(i,j,t,1)=0doit=1,4aa(i,j,t,it)=datas(i,j,t,it)-a(i,j,t)enddoenddoenddoenddo!求均方差场dot=1,12doj=1,nydoi=1,nxbb(i,j,t)=0doit=1,4bb(i,j,t)=bb(i,j,t)+(datas(i,j,t,it)-a(i,j,t))**2enddobb(i,j,t)=sqrt(bb(i,j,t)/4)enddoenddoenddoopen(11,file='d:\jp.grd',form='unformatted',access='direct',recl=nx*ny)irec=0doit=1,4dot=1,12irec=irec+1write(11,rec=irec)((aa(i,j,t,it),i=1,nx),j=1,ny)enddoenddoclose(11)open(10,file='d:\qhc.grd',form='unformatted',access='direct',recl=nx*ny)irec=0dot=1,12irec=irec+1write(10,rec=irec)((a(i,j,t),i=1,nx),j=1,ny)enddoclose(10)4open(12,file='d:\jfc.grd',form='unformatted',access='direct',recl=nx*ny)irec=0dot=1,12irec=irec+1write(12,rec=irec)((bb(i,j,t),i=1,nx),j=1,ny)enddoclose(12)1000format(2i7)2000FORMAT(37F6.2)3000format(37f8.1)4000format(37f7.2)endprogramone生成图片的gs文件这里只给出生成气候场的gs文件'reinit''open:e\cli.ctl''enableprinte:\cli.gmf''sett1'*'sett7''dh''print''disableprint';图11月气候场5图27月气候场图31982年1月距平场图41982年7月距平场6图51月均方差场图67月均方差场实习二相关系数及自相关系数的计算根据所给年平均气温和冬季平均气温的等级数据（数据略）进行下列计算：年份年平均气温1970~19793.403.303.202.903.402.803.603.002.803.001980~19893.103.002.902.703.503.203.102.802.902.90年份冬季平均气温1970~19793.243.143.262.383.322.712.843.942.751.831980~19892.802.812.633.203.603.403.071.872.632.471）计算两个气温之间的相关系数。72）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度。要求：实习报告中附出简单相关系数或自相关系数程序。求相关系数r的fortran程序：programtwoimplicitnoneintegerirealx(20),y(20)realsum1,sum2,av1,av2,f1,f2,f3,rxydata(x(i),i=1,20)/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(i),i=1,20)/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/!计算均值sum1=0sum2=0doi=1,20sum1=sum1+x(i)sum2=sum2+y(i)enddoav1=sum1/20av2=sum2/20!计算方差的相关系数rf1=0doi=1,20f1=f1+(x(i)-av1)*(y(i)-av2)enddof2=0f3=0doi=1,20f2=f2+(x(i)-av1)**2f3=f3+(y(i)-av2)**2enddorxy=f1/(sqrt(f2*f3))print*,'R=',rxy!!computetheRaboveendprogramtwo运行程序，得到结果如下图所示：图7相关系数R的计算结果8计算自相关系数的程序其中j取1~7：programtwo_2implicitnoneintegeri,jrealx(20),y(20),s1(7),s2(7),r1(7),r2(7)realsum1,sum2,av1,av2,f1,f2,f3,rxydata(x(i),i=1,20)/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(i),i=1,20)/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/!open(5,file='d:\sedata.txt')!read(5,100)x!read(5,100)y!doi=1,20!print*,x(i)!print*,y(i)!enddosum1=0sum2=0doi=1,20sum1=sum1+x(i)sum2=sum2+y(i)enddoav1=sum1/20av2=sum2/20!print*,av1,av2f2=0f3=0doi=1,20f2=f2+(x(i)-av1)**2f3=f3+(y(i)-av2)**2enddof2=f2/20f3=f3/20!print*,f2,f3doj=1,7s1(j)=0s2(j)=0doi=1,20-js1(j)=s1(j)+(x(i)-av1)*(x(i+j)-av1)s2(j)=s2(j)+(y(i)-av2)*(y(i+j)-av2)!print*,s1(j),s2(j)enddo9enddodoj=1,7s1(j)=s1(j)/(20-j)s2(j)=s2(j)/(20-j)!print*,s1(i),s2(i)enddodoj=1,7r1(j)=s1(j)/f2r2(j)=s2(j)/f3enddoopen(10,file='d:\result.txt')!doi=1,7write(10,'(7f10.5)')(r1(i),i=1,7)write(10,'(7f10.5)')(r2(i),i=1,7)!enddo!write(10,*)r2(i)close(10)write(*,'(7f10.5)')(r1(i),i=1,7)write(*,'(7f10.5)')(r2(i),i=1,7)end运行结果如下：图8自相关系数的计算结果由结果可知，年平均气温在滞后长度j=7、冬季序列在j=4最大。实习三一元线性回归根据下列数据，编写一元线性回归程序，当x=14.5时，y的值应该为多少？X012345678910111213141516Y3029.128.428.12827.727.527.22726.826.526.326.125.725.324.824要求：实习报告中附出一元线性回归程序，并给出相应原数据及回归数据图形（一张图）线性回归程序如下：programlinearimplicitnoneintegerirealx(17),y(17)realf1,f2,b,avx,avy,b0,y1data(x(i),i=1,17)/0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16/10data(y(i),i=1,17)/30,29.1,28.4,28.1,28,27.7,27.5,27.2,27,26.8,26.5,26.3,26.1,25.7,25.3,24.8,24/f1=0doi=1,17f1=f1+x(i)*y(i)enddoavx=sum(x)/17!求均值avy=sum(y)/17f2=0doi=1,17f2=f2+x(i)**2enddob=(17*f1-sum(x)*sum(y))/(17*f2-sum(x)**2)!求bb0=avy-b*avxy1=b*14.5+b0print*,'b=',bprint*,'b0=',b0print*,'y=',b0,b,'x'!!computethebabo