1、氣象統(tǒng)計實習(xí)報告 專業(yè)：大氣科學(xué) 班級：xxxx級x班 學(xué)號：xxxx 姓名：*實習(xí)一 求500hPa高度場氣候場、距平場和均方差場一、實習(xí)結(jié)果1981.1距平場1981.1 500hpa高度場在歐亞大陸為正距平，在印度洋和太平洋為負(fù)距平。1981.1氣候場1982.1的氣候場呈現(xiàn)明顯的沿緯度的上升而下降的趨勢。1981.1均方差場1981.1的均方差場在歐亞大陸的南側(cè)有最小值，在10N20N的南側(cè)歲緯度而降低，在其北側(cè)隨緯度升高。二、相關(guān)的fortran程序，gs文件（1）Fortran程序: program main parameter(nx=37,ny=17,mo=12,yr=4) re

2、al var(nx,ny,mo,4)!數(shù)據(jù) real vars(nx,ny,mo)!4年氣候態(tài) real jp(nx,ny, mo,4)!距平 real fc(nx,ny,mo)!方差 real jfc(nx,ny,mo)!均方差 integer i,j,m,y,irec real:summ=0.0!*求均方差*do m=1,12 do i=1,37 do j=1,17 do y=1,4 summ=summ+jp(i,j,m,y)*2 enddo fc(i,j,m)=summ/4!方差 jfc(i,j,m)=sqrt(fc(i,j,m)!均方差 summ =0.0 enddo enddoend

3、do!*open(11,file=d:shixivars.grd,form=binary )open(12,file=d:shixijp.grd,form=binary )open(13,file=d:shixijfc.grd,form=binary )do m=1,12 write(11)(vars(i,j,m),i=1,nx),j=1,ny)enddodo m=1,12 write(13)(jfc(i,j,m),i=1,nx),j=1,ny)enddodo y=1,4 do m=1,12 write(12)(jp(i,j,m,y),i=1,nx),j=1,ny) enddoenddo cl

4、ose(11) close(12) close(13)End（2）Gs程序1. 距平場open e:anomaly.ctlset gxout shadedd xiaobaoset gxout contourd xiaobaodraw title 1982年1月enable print e:anomaly.gmfprintdisable print2.氣候場open e:climate.ctlset gxout shadedd xiaobaoset gxout contourd xiaobaodraw title 1982.1enable print e:climate.gmfprintdisa

5、ble print3. 均方差場open e:deviation.ctlset gxout shadedd xiaoluset gxout contourd xiaoludraw title 1982.1 deviationenable print e:deviation.gmfprintdisable print*實習(xí)二 計算給定數(shù)據(jù)資料的簡單相關(guān)系數(shù)和自相關(guān)系數(shù) 根據(jù)下表中年平均氣溫和冬季平均氣溫的等級數(shù)據(jù)進(jìn)行下列計算：1）計算兩個氣溫之間的簡單相關(guān)系數(shù)。2）分別找出兩個氣溫數(shù)據(jù)自相關(guān)系數(shù)絕對值最大的滯后時間長度。（滯后長度最大取10）一、實習(xí)結(jié)果（1） 計算簡單相關(guān)系數(shù)計算出相關(guān)系數(shù)為r

6、=0.（2）分別找出兩個氣溫數(shù)據(jù)自相關(guān)系數(shù)絕對值最大的滯后時間長度。（滯后長度最大取10）可以知道，年平均氣溫在滯后長度j=7，冬季j=4最大二、相關(guān)的fortran程序（部分） implicit none real 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 !協(xié)方差 real :jx=0.0 real :jy=0.0 real :jfx=0.0 real :jfy=0.0,r integer i,t

7、,j real :m=0.0 data 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,20 jx=jx+x(i) jy=jy+y(i) end do do i=1,20 s=s+(x(i)-jx/2

8、0)*(y(i)-jy/20) jfx=jfx+(x(i)-jx/20)*2 jfy=jfy+(y(i)-jy/20)*2 end do r=s/20/sqrt(jfx/20*jfy/20) print *,r=,r do t=1,10 do i=1,20-t zx(t)=zx(t)+(x(i)-jx/20)*(x(i+t)-jx/20) end do rzx(t)=zx(t)/(20-t)/(jfx/20) print *,t,rzx(t) if(abs(rzx(t)m) then m=abs(rzx(t) j=t end if end do print *,年平均溫度的自相關(guān)系數(shù)絕對值最大

9、的滯后時間長度 print *,j實習(xí)三 計算給定數(shù)據(jù)的落后交叉相關(guān)系數(shù)和偏相關(guān)系數(shù) 根據(jù)下表北京冬季（12月2月）氣溫資料計算：12月氣溫與1月和2月氣溫的落后交叉相關(guān)系數(shù)（滯后長度最大取10）和偏相關(guān)系數(shù)。在實習(xí)報告中給出程序。年份12月1月2月19511.0-2.7-4.31952-5.3-5.9-3.51953-2.0-3.4-0.81954-5.7-4.7-1.11955-0.9-3.8-3.11956-5.7-5.3-5.91957-2.1-5.0-1.619580.6-4.30.21959-1.7-5.72.01960-3.6-3.61.31961-3.0-3.1-0.81962

10、0.1-3.9-1.11963-2.6-3.0-5.21964-1.4-4.9-1.71965-3.9-5.7-2.51966-4.7-4.8-3.31967-6.0-5.6-4.91968-1.7-6.4-5.11969-3.4-5.6-2.01970-3.1-4.2-2.91971-3.8-4.9-3.91972-2.0-4.1-2.41973-1.7-4.2-2.01974-3.6-3.3-2.01975-2.7-3.70.11976-2.4-7.6-2.21977-0.9-3.5-2.31978-2.7-4.2-0.51979-1.6-4.5-2.91980-3.9-4.8-1.41、

11、 實驗結(jié)果偏相關(guān)系數(shù)計算出的 12 月氣溫與 1 月氣溫的偏相關(guān)系數(shù)為 0.327，12 月氣溫與 2 月氣溫的偏相關(guān)系數(shù)為 0.290。二、相關(guān)的fortran程序Funtion area(x,y) implicit none real x(30) real y(30) real :jx=0.0 real :jy=0.0 real :jfx=0.0 real :jfy=0.0,r integer i, real :s=0.0 do i=1,30 jx=jx+x(i) jy=jy+y(i) end do do i=1,30 s=s+(x(i)-jx/30)*(y(i)-jy/30) jfx=j

12、fx+(x(i)-jx/30)*2 jfy=jfy+(y(i)-jy/30)*2 end do r=s/30/sqrt(jfx/30*jfy/30) print *,r=,rend funtion area*實習(xí)四 求給定數(shù)據(jù)的一元線性回歸方程利用下表數(shù)據(jù)，以環(huán)流指標(biāo)為預(yù)報因子，氣溫為預(yù)報量，計算氣溫和環(huán)流指標(biāo)之間的一元線性回歸方程，并對回歸方程進(jìn)行檢驗。年份氣溫T環(huán)流指標(biāo)19510.93219521.22519532.22019542.4261955-0.52719562.5241957-1.128195802419596.21519602.71619613.2241962-1.130196

13、32.52219641.23019651.82419660.63319672.42619682.52019691.2321970-0.835一、實習(xí)結(jié)果（1）用excel制作的氣溫-環(huán)流的醫(yī)院線性回歸方程回歸方程為：（2）回歸方程的檢驗檢驗結(jié)果：F=20.18F=4.41，回歸方程顯著2、 Fortran程序（部分）（1） 回歸方程的檢驗!回歸方程顯著性檢驗!計算兩數(shù)組的距平及均方差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

14、do!計算協(xié)方差i=1;t=0do j=1,20t=t+diff(j,i+1)*diff(j,i)end doE=t/20!計算相關(guān)系數(shù)r=E/(var(1)*var(2)print*,the relative value r is:， rF=r*2*(nx-2)/(1-r*2)實習(xí)五 求給定數(shù)據(jù)的多元線性回歸方程說明：x1-x4為四個預(yù)報因子，y為預(yù)報量；樣本個數(shù)n=13要求：選取預(yù)報因子1、2、4，求預(yù)報量的標(biāo)準(zhǔn)化回歸方程。i12345678910111213x17111117113122111110x226295631525571315447406668x3615886917221842

15、398x46052204733226442226341212y78.574.3104.387.695.9109.2102.772.593.1115.983.8113.3109.4一、用excel制作的醫(yī)院線性回歸方程標(biāo)準(zhǔn)化變量回歸方程：*實習(xí)七 計算給定數(shù)據(jù)的11年滑動平均和累積距平 利用數(shù)據(jù)ma.dat，編寫11點滑動平均的程序，ma.for給出了閱讀資料的fortran程序。數(shù)據(jù)在文件夾中單獨給出。要求：實習(xí)報告中附出程序，并給出原數(shù)據(jù)和滑動后數(shù)據(jù)的圖形（1張圖）和累積距平數(shù)據(jù)圖形（1張圖）一、實習(xí)結(jié)果累積距平數(shù)據(jù)圖形滑動后數(shù)據(jù)二、部分fortran程序：program ma! dimen

16、sion 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,i

17、h=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+1 x1(i)=sum(x(i:i+ih-1)/ih end doprint*,x1 open(10,file=d:moveaverage.txt) write(10,(f10.6)(x1(i),i=1,n-ih+1) close(10)end*實習(xí)八 對給定的海溫數(shù)據(jù)進(jìn)行EOF分析 給出海表溫度距平數(shù)據(jù)資料sstpx.grd，以及相應(yīng)的數(shù)據(jù)描述文件sstpx.ctl，對其進(jìn)行EOF分析，資料的時空范圍可以根據(jù)sstpx.ctl獲知。 數(shù)據(jù)在文件夾中單獨給出，距平或者標(biāo)準(zhǔn)化距平處理后再進(jìn)行EOF。Zhunsst.for給出了如何讀取資料，Ssteof.for為對距平或者標(biāo)準(zhǔn)化距平處理后的資料進(jìn)行EOF分析。要求：實習(xí)