气象统计方法实习报告

目录TOC\o"1-5"\h\z\o"CurrentDocument"实习一求500hPa高度场气候场、距平场和均方差场 3\o"CurrentDocument"1、资料介绍 3\o"CurrentDocument"2.要求 3\o"CurrentDocument"3、实习结果 4\o"CurrentDocument"1)、FORTRAN源程序 4、grads文件 6、实习结果 6\o"CurrentDocument"实习二计算给定数据资料的简单相关系数和自相关系数 4-1\o"CurrentDocument"1、 资料介绍 U-2、 要求 1--1--\o"CurrentDocument"3、 实习结果 投-\o"CurrentDocument"、Fortran源程序 -R\o"CurrentDocument"、程序运行结果: -4\o"CurrentDocument"实习三分析中国夏季降水线性趋势的分布特征- -5\o"CurrentDocument"资料介绍及要求： 15-\o"CurrentDocument"2。实习结果 1-5-\o"CurrentDocument".matlab程序 -1-5\o"CurrentDocument"(2).程序运行结果 16\o"CurrentDocument"实习四求给定数据的一元线性回归方程 17\o"CurrentDocument"1、 资料介绍及要求 47-\o"CurrentDocument"2、 实习结果 1-8-\o"CurrentDocument"（1）、MATLAB程序 18-（2） 、程序运行结果- 18\o"CurrentDocument"（3） 、结果分析 -19-\o"CurrentDocument"实习五对给定的海温数据进行EOF分析 -21\o"CurrentDocument"1、 资料介绍 21-\o"CurrentDocument"2、 要求 -24--\o"CurrentDocument"3、 实习结果： 24-（1） 、FORTRAN源程序 24\o"CurrentDocument"（2） 空间场和时间序列的ctl文件 -23（3） 运行结果 24-\o"CurrentDocument"（4） 分析 26-\o"CurrentDocument"实习三（附加）计算给定数据的11年滑动平均和累积距平 -28\o"CurrentDocument"1、 资料介绍 28-2、 要求 2--8--\o"CurrentDocument"3、 实习结果 28-实习四（附加）求给定数据的多元线性回归方程 3-1\o"CurrentDocument"1、 说明 -31■-2、 要求 3--1--\o"CurrentDocument"3、 实习结果： 3-1--（1） Matlab源程序 34（2） 运行结果 36-实习一求500hPa高度场气候场、距平场和均方差场1、资料介绍有一500hPa高度场资料，文件名h500。dat,范围：60〜150E,0〜40N。时段：1982.1〜1985。12共48个月。水平分辨率：2.5大2.5,格点数：37大17。2,要求编fortran程序，求500hPa高度场的气候场；距平场；均方差场.并能用Grads做出图形，实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。注：h500。For给出了如何用fortran读取ASCII码资料h500。dat.3、实习结果1）、FORTRAN源程序programex_gradsimplicitnoneinteger,parameter::nx=37,ny=17,nz=4,nt=12integeri，j，iz，itrealvar（nx,ny,nz,nt）,cl（nx，ny，nt）,sum，jp（nx,ny，nz，nt），jfc（nx,ny，nt）！Openingfileopen（10，file=’g：\gradsdata\h500。dat’）doiz=1,nzdoit=1，ntread（10，1000）read（10,3000）（（var（i，j，iz，it），i=1，nx），j=1,ny）enddoenddo1000format（2i7）2000format（37f6.2）3000format（37f8.1）4000format（37f7。2）close（10）！Outputopen（16，file='g:\gradsdata\h500。grd’，form='binary'）doiz=1,nzdoit=1，ntwrite（16） （（var（i，j，iz,it）,i=1，nx），j=1，ny）enddoenddo!CalculatingtheClimatologicalFielddoit=1，ntdoi=1，nxdoj=1,nysum=0doiz=1,nzsum=sum+var（i,j，iz，it）enddocl（i，j,it）=sum/4enddoenddoenddo!Outputclimate-fileopen(12,file='g:\gradsdata\climate.grd',form=’binary’)doit=1，ntwrite(12)((cl(i,j,it),i=1,nx)，j=1，ny)enddo!CalculatingtheAnomalydoiz=1，nzdoit=1,ntdoi=1，nxdoj=1，nyjp(i，j,iz，it)=var(i,j，iz，it)-cl(i，j,it)enddoenddoenddoenddoopen(13，file='g:\gradsdata\anomaly。grd’，form='binary')！Outputanomaly—filedoiz=1,nzdoit=1，ntwrite(13)((jp(i，j,iz，it)，i=1,nx)，j=1,ny)enddoenddo!CalculatingtheMean—squareDeviationdoit=1，ntdoi=1,nxdoj=1,nysum=0doiz=1，nzsum=sum+(jp(i，j，iz,it))*大2enddojfc(i，j,it)=sqrt(sum/4)enddoenddoenddo!Outputmean—squaredeviation-fileopen(14，file=’g:\gradsdata\deviationogrd’，form=’binary’)doit=1，ntwrite(14) ((jfc(i，j,it)，i=1,nx),j=1，ny)

enddoend(2)、grads文件,openg:\gradsdata\***。ctl’ (***为所求场对应的ctl文件名)'setlat040'’setlon60150’’setlev500''enableprintg:\gradsdata\***.gmf'(***为所求场名称)i=1while(i<=48(或12))'sett’i’dh’’print’'c’i=i+1endwhile’disableprint’；、实习结果、原始场1982年1月,气， /Of sGr 页 iiiiii iThi TihE F E I：lthebCChPciHsiz-,气， /Of sGr 页 iiiiii iThi TihE F E I：l1982年7月

'■-56605B5Q110€'■-56605B5Q110€the500hPaHeightFelc/1382July结果分析：冬季（此处以1月为代表）等高线分布整体平缓，表明高度场分布相对均匀，且北部接近极地位势高度低，赤道地区位势高度高，这与太阳直射点在1月在南半球，使北半球整体的辐射吸收随纬度增加而减小有关。北半球的气块受热随纬度递增而递减，因而膨胀率递减，故南方接近赤道地区的气体膨胀大，位势高，而北方近极地地区气体膨胀小，位势低夏季（此处以7月为代表），大洋上空出现副热带高压（588线位置），在东亚地区存在较为明显的位势高度槽，即东亚大槽.夏季在30°N以北的地区位于西风带中，从图中可看出明显的西风带长波特征。除东亚大槽外，在中亚地区也存在一长波槽，这些槽线发生长波调整时，会在部分地区发生较剧烈的天气变化。此外，从图中可以看出，1982年7月副热带高压脊线的平均位置位于25°N,125°E附近，我国华东地区位于副高北侧西南气流控制下，西南气流为水汽输送的主要通道，该地区发生降水较为频繁。从图中还可以看出冬夏季的位势高度分布存在明显差异，这与太阳直射点的年纪变化密切相关.、气候场2月

8月结果分析：气候场是多年数据中同时段的平均值序列，表征了区域内多年平均的位势高度变化.从气候场图形可以看出多年平均的500hPa高度场中，冬季（此处以多年平均的2月气候场为代表）等高线较平直，大陆上等高线较稀疏，而海洋上等高线较密集，这表明大陆上空的位势高度变化率较海洋小。此外，冬季等位势高度线分布平直，还说明冬季的天气变化显著。夏季（此处以多年平均的8月气候场为代表）在太平洋上有副热带高压，副高是深厚的系统，所以可以一直延伸到500hPa高度。东亚地区存在明显的西风带长波槽，即东亚大槽，东亚大槽的移动和变化配合副热带高压和夏季风的影响，会使我国大部分地区天气造成强烈变

化，如形成大范围降水或强对流天气。同时，孟加拉湾处存在一低涡.由于高原的阻挡作用，这一系统对我国的影响并不显著。图中副高脊线8月的多年平均位置位于30°N以北，且东部长波槽位于110°E附近，故由气候场分析，华北地区位于长波槽前，又槽线受到副高的阻挡作用，因而华北地区容易形成降水。、距平场1982年1月th已A-omcIi.sHcchtI-icrl h"□/'1闩萨亍：i1982年7月.”itheAnomalousHeightFieldof500hPa/1982：i1982年7月.”itheAnomalousHeightFieldof500hPa/1982Julymi：i1H：ll?l：l结果分析：距平场指示了位势高度的震荡趋势，因距平的平均值为0，则大于0的值表明位势高度偏高，小于0的值表示位势高度偏低.从图中可看出冬季（此处以1982年1月距平场做代表）在大陆位势高度为正距平,而在大洋则存在明显的负距平.则由距平场的性质得，冬季在大洋上位势高度偏高，在大陆上位势高度偏低。其原因是，海水的比热大于陆地，则冬季海洋温度比陆地高，所以海洋上气块膨胀更多，位势高度更高。夏季（此处以1982年7月距平场做代表）相反,大洋上位势高度偏低，而陆地偏高。在70°E〜90°E附近位势偏高的原因在于夏季青藏高原接受太阳辐射使之在对流层中层形成热源，位势高度因比大陆其他位置高.由此可见，位势高度的变化不仅与海陆差异有关，同时也与地形有关。在海洋上副热带高压所在的位置存在证据平值。

、均方差场6月theWeanSquareDeviationalHeightFieldof500hPd/June1D0C11QE120CMOE12月1D0C11QE120CMOE12月theMeanSquareDeviationalHeightFieldofSOOhPa/Dec结果分析：均方差场反应同一时段内的位势高度变化幅度的大小。由图可以看出，整体位势高度在大陆上的变化幅度比海洋小，且海洋上冬季的变化幅度比夏季大，而陆地上相反，冬季的变化幅度比夏季小。因为陆地的比热小，所以陆地在夏季白天与夜间的温差大于冬季，对应的高度场震荡就比冬季剧烈。海洋上的位势高度变化幅度的影响方面温差为次要因素，其主要受到副热带高压,西风带长波槽脊影响，位势高度根据天气系统的移动而变化，所以震荡幅度较大。除受到天气系统影响外，海洋上的位势高度场还受到洋流等因素的影响。

实习二计算给定数据资料的简单相关系数和自相关系数1、资料介绍根据下表中年平均气温和冬季平均气温的等级数据进行下列计算:1） 计算两个气温之间的简单相关系数。2） 分别找出两个气温数据自相关系数绝对值最大的滞后时间长度.（滞后长度7最大取10）中国1970-1989年年平•均和冬季平切［温年份年乎均气温1970-19793.403.303.202.903.402.S03.603.002803.00198019893.103.002.902.703.503.203.102.80Z902.90年份冬季"均气温1970-19793.243143.262383.322713.942.751831980-19892.802812.633203.603403.071.872.632、要求实习报告中附出简单相关系数或自相关系数程序。答案：r=0。47年平均气温在滞后长度j=3、冬季序列在j=4最大.3、实习结果、Fortran源程序PROGRAMEXAMIMPLICITNONEINTEGER,PARAMETER：:N=20INTEGERi,j,k,ty,tw,tywREAL::avr_y=0,avr_w=0,sy=0,sw=0,rxy=0,max_y=0,max_w=0,max_yw=0REALy(N),w(N)DATAy/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/DATAw/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/REALsyy(N),sww(N),r(N),rty(N),rtw(N),rtyw(N),rxy_ty(N),rxy_tw(N),rxy_tyw(N)!求两数组平均值DOi=1,Navr_y=avr_y+y(i)avr_w=avr_w+w(i)ENDDOavr_y=avr_y/Navr_w=avr_w/N!简单相关系数DOj=1,Nsyy(j)=(y(j)—avr_y)**2sy=sy+syy(j)sww(j)=(w(j)—avr_w)**2sw=sw+sww(j)ENDDOsy=sqrt(sy/N)sw=sqrt(sw/N)DOj=1,Nr(j)=((y(j)—avr_y)/sy)*((w(j)—avr_w)/sw)rxy=rxy+r(j)ENDDOrxy=rxy/NPRINT”(/’1970-1989年全年平均气温与冬季平均气温的简单相关系数rxy=',f5.2)”,rxyk=0!自相关系数DOty=1,N/2DOi=1,N-tyrty(i)=((y(i)—avr_y)/sy)大((y(i+ty)-avr_y)/sy)rxy_ty(ty)=rxy_ty(ty)+rty(i)ENDDOrxy_ty(ty)=rxy_ty(ty)/(N—ty)rxy_ty(ty)=ABS(rxy_ty(ty))IF(rxy_ty(ty)〉max_y)THENmax_y=rxy_ty(ty)k=tyENDIFENDDOPRINT”('全年平均气温绝对值最大自相关系数rxy_ty=',f7。4,/,’滞后时间长度k=',I2)”,rxy_ty(k)，kk=0DOtw=1,N/2DOi=1,N-twrtw(i)=((w(i)-avr_w)/sw)*((w(i+tw)-avr_w)/sw)rxy_tw(tw)=rxy_tw(tw)+rtw(i)ENDDOrxy_tw(tw)=rxy_tw(tw)/(N-tw)rxy_tw(tw)=ABS(rxy_tw(tw))IF(rxy_tw(tw)>max_w)THENmax_w=rxy_tw(tw)k=twENDIFENDDOPRINT"('冬季平均气温绝对值最大自相关系数rxy_tw=',f7。4,/,’滞后时间长度k=’，I2)”,rxy_tw(k),kk=0!落后交叉相关系数DOtyw=1，N/2DOi=1,N—tywrtyw(i)=((y(i)—avr_y)/sy)*((w(i+tyw)—avr_w)/sw)rxy_tyw(tyw)=rxy_tyw(tyw)+rtyw(i)ENDDOrxy_tyw(tyw)=rxy_tyw(tyw)/(N—tyw)rxy_tyw(tyw)=ABS(rxy_tyw(tyw))IF(rxy_tyw(tyw)>max_yw)THENmax_yw=rxy_tyw(tyw)k=tywENDIFENDDOPRINT”(’全年平均温度与冬季平均气温之间的落后交叉相关系数rxy_tyw=’，f7.4,/,'滞后时间长度k=',I2)”，rxy_tyw(k)，kEND(2)、程序运行结果:实习三分析中国夏季降水线性趋势的分布臃1。 资料介绍及要求：利用数据160zhan-rainfall—summer。txt,编写求1982-2006年中国160站各站夏季降水线性倾向率，给出分布图，并进行简单分析。read—rain.for给出了阅读资料的fortran程序。数据在文件夹中单独给出.2。 实习结果(1)。matlab程序%编写求1982-2006年中国160站各站夏季降水线性倾向率clearallclcfid=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);fori=1：160;j(i,1：25)=data2(i,4:28);n1=1982：1：2006；pp(i,：)=polyfit(n1,j(i,1：25),1)；endb=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)

xlabel（'精度’）；ylabel（’纬度’）；title（’1982—2006年中国160站各站夏季降水线性倾向率分布图'）（2）.程序运行结果1982-2006年中国160站各站夏季降水线性倾向率分布图151050-5-10151050-5-10实习四求给定瞄的一元魄回归方程1、资料介绍及要求利用下表数据，以环流指标为预报因子，气温为预报量，计算气温和环流指标之间的一元线性回归方程，并对回归方程进行检验。年份气温T环流指标19510。93219521。22519532。22019542。4261955—0。52719562。5241957—1。128195802419596。21519602。71619613。2241962—1。13019632。52219641.23019651。82419660。63319672.42619682.52019691。2321970-0.835答案：宁 7.5-0.23x仃=20.18〉Fa=4.41，回归方程显著

2、实习结果(1)、MATLAB程序%实习四求给定数据的一元线性回归方程ClimateData=xlsread(’F:气象统计方法\实验数据\气象统计实验四数据。xls');%从Excel文件读取数据x=ClimateData(:,3);%提取ClimateData的第三列，即环流指标y=ClimateData(:，2);%提取ClimateData的第三列，即气温Txdata=[ones(size(x,1),1),x]； %在原始数据x的左边加一列1,即模型包含常数项[b,bint,r,rint,s]=regress(y,xdata);%调用regress函数作一元线性回归yhat=xdata大b； %计算y的估计值%定义元胞数组，以元胞数组形式显示系数的估计值和估计值得95%置信区间head1={'系数的估计值’，'估计值的95%置信下限'，'估计值的95%置信上限'}；[head1;num2cell([b,bint])]%定义元胞数组，以元胞数组形式显示y的真实值、y的估计值、残差和残差的95%置信区间head2={'y的真实值’,'y的估计值’，'残差'，'残差的95%置信下限’，'残差的95%置信上限'};%同时显示y的真实值，y的估计值、残差和残差的95%置信区间[head2;num2cell([y，yhat,r,rint])]%定义元胞数组，以元胞数组形式显示判定系数、F统计量的观测值、检验的P值和误差方差的估计值head3={'判定系数'，'F统计量的观测值'，’检验的P值’，'误差方差的估计值'}；[head3；num2cell(s)](2)、程序运行结果ans=4。6554] [—0.3433] [10o3637]-0o12534。6554] [—0.3433] [10o3637]-0o1253][ 7.5095] [[ -0o2343] [ans='y的真实值’限’'y的估计值’’残差’'残差的95%置信下限’’残差的95%置信上[ 0。9000][0.0123][0o8877][ —1o5633][ 3.3388][ 1.2000][1.6523][—0.4523][ —3o0331][ 2.1285][2o2000][2o8237][-0.6237][ —3o1171][ 1.8696]

[[2。4000]—0o5000][1o4180][1o1837[0.9820] [] [—1.6837] [—1o5611] [—4o1250][3o5251]0o7576][2.5000][1o8866][0o6134] [—1o9531][3o1800][-1o1000][0o9494][—2.0494] [-4.4072] [0.3084][0][1.8866][—1.8866] [—4o2867] [0o5136][6o2000][3o9952][2o2048] [0o1971][4o2125][2o7000][3o7609][—1o0609] [-3.3535][1o2317][3o2000][1o8866][1.3134] [-1o1840] [3o8108][-1o1000][0o4808][-1o5808] [-3.9959][0o8342][2o5000][2o3552][0.1448] [-2o4137] [2.7034][1.2000][0o4808][0o7192] [-1o8001] [3.2385][1.8000][1o8866][-0o0866] [-2.6717] [2.4986][0o6000][—0.2220][0o8220] [—1o5996] [3.2437][2o4000][1.4180][0o9820] [—1.5611] [3o5251][2.5000][2o8237][-0o3237] [—2o8318] [2o1844][1o2000][0.0123][1o1877] [-1o2301] [3.6056][-0.8000][-0o6906][-0.1094] [—2.4794] [2.2606]ans='判定系数’’F统计量的观测值’’检验的P值’ ’误差方差的估计值’[0.5313] : 20.4045] [2.6673e-004] [ 1。5134]（3）、结果分析从输出的结果看，常数项和回归系数的估计值分别为7.5095和一Go2343,从而可以写出线性回归方程为y=7.5095-0.2343x回归系数估计值的置信区间为[-Go3433,-0.1253].对回归直线进行显著性检验，原假设和对立假设分别为H0:气=0,H1:R1古0检验P的值为2o6673X10-4〈0.01，可知显著性水品0=0.01下应拒绝原假设H0，可认为y（环流指数）与x（气温T）的线性关系是显著的。原始数据散点图与回归直线图＞＞plot（X,y,'k.’，’Markersize’，15）%原始数据散点图＞〉holdon〉＞plot（x,yhat,’linewidth',3）%回归直线图〉〉xlabel（'环流指标（x）’）％标注x轴〉〉ylabel（'气温（y）'）%标注y轴＞〉legend（'原始散点'，'回归直线’）％加标注框实习五对给定的海温数据进行EOF分析1、 资料介绍给出海表温度距平数据资料sstpxogrd,以及相应的数据描述文件sstpx.ctl,对其进行EOF分析，资料的时空范围可以根据sstpxoctl获知.数据在文件夹中单独给出，距平或者标准化距平处理后再进行EOFoZhunsstofor给出了如何读取资料，Ssteofofor为对距平或者标准化距平处理后的资料进行EOF分析。2、 要求实习报告中给出第一特征向量及其时间系数，并分析其时空特征。3、 实习结果：(1)、FORTRAN源程序!preparedataforeofanalysis!theprogramistonormalizeseasurfacetemperature(SST)!mt:thelengthoftimeseries；!mo:themonthnumbers; my：theyearnumbers;!sst:seasurfacetemperaturedata;!sst3:theworkarray； avf:theaverageofSST；!df:thevarianceofSST；programmainparameter(mo=12,my=43,nx=18,ny=12,mt=516)dimensionavf(mo,nx,ny),df(mo,nx,ny)dimensionsst(nx,ny,mt),sst3(nx,ny,mt)open(1,file=’g：\sstpx\sstpx°grd’,form=’unformatted',access='direct；recl=nx*ny)irec=1doit=1,mtread(1,rec=irec)((sst(i,j,it),i=1,nx),j=1,ny)irec=irec+1enddo!averagedoj=1,nydoi=1,nxdok=1,modoit=k,mt,12avf(k,i,j)=avf(k,i,j)+sst(i,j,it)enddoavf(k,i,j)=avf(k,i,j)/myenddoenddoenddo!variancedoj=1,nydoi=1,nxdok=1,modoit=k,mt,12df(k,i,j)=df(k,i,j)+(sst(i,j,it)—avf(k,i,j))*大2enddodf(k,i,j)=sqrt(df(k,i,j)/my)enddoenddoenddo!standardizingdoj=1,nydoi=1,nxdok=1,modoit=k,mt,12if(sst(i,j,it)==—999.0)thensst3(i,j,it)=-999.0elsesst3(i,j,it)=(sst(i,j,it)—avf(k,i,j))/df(k,i,j)endifenddoenddoenddoenddo!outputfileopen(2,file=’g：\sstpx\standard.grd',form='unformatted,access='direct',recl=nx*ny)irec=1doit=1,mtwrite(2,rec=irec)((sst3(i,j,it),i=1,nx),j=1,ny)irec=irec+1enddoclose(2)close(1)end分解后的时间系数写入tcf.grd文件中，空间场写入evf。grd文件中，特征值和分析误差写入sstpx文件夹下的°[8。洒31如且3七文件，特征向量写入。[8。洒360甜且3七文件。由eigenva1ues.dat中的标准特征向量可得出一般特征值的前两个模态有效。用grads打开evf.ct^^tcf。ctl，分别画出海平面气温EOF分解后的空间场和时间序列。(2)空间场和时间序列的ctl文件evf。ctldsetg：\sstpx\evf.grdtitleCoadsSSTAEundef—999o0xdef18linear12010ydef12linear—27o55zdef1linear10001tdef2linear1jan19481monthvars1S099CoadsSSTanomalyinterperatedusingendvarstcf.ctldsetg：\sstpx\tcfogrdtitleCoadsSSTATundef—999.0xdef1linear12010ydef1linear-27.55zdef1linear10001tdef516linear1jan19481monthvars2a099timecoefficient1b099timecoefficient2endvars(3)运行结果第一模态空间场-口.WI1CF:h=SS~SpacialFieldforEOFAnalysis/NO.1I20E时间系数16DWIDOW15N--u.os——L一如 ffqy•-Cij2-01-- -u工•"八-一"、一i"N：SSTTi~eCoef'ic'e^tforEOFA-avsis/Ko.1-012-—'..-HSf O.OS--.20S-12CKI40E..史二-I6OE,_-DI6CWHOWIDOWGOW时间系数I.HuSS-Ti-'-ic? 七’-012-—'..-HSf O.OS--.20S-12CKI40E..史二-I6OE,_-DI6CWHOWIDOWGOW时间系数I.HuSS-Ti-'-ic? 七’-C：1FAn::vxb：j.I9：i0巧E1^15511-9150I9B31W019^0第二模态空间场:heSS-Spcc'I~icld"c-EC"A-icys's/f\0.2TOC\o"1-5"\h\z-n1■?.、 .广drJ -■/ r._--，•-…"二们L■- --0.09第一特征向量Eigenvaluesodat文件给出了EOF分析的第一特征向量值的216个值000 -0.02200.01800.04300.03400.06400.05400o06000o06900o05500o03700.0100-0.0190 —0o0210—0o04600000-0o02300o02100.05000o06100.05400o05000.02600o01300.0070-0.0600—0.0580—0o0690-0.0530-0o05500—0.01300-0.00100o02300.02800o03500.05100.0320 0o0080—0.0550-0o0730-0.1070-0o1180-0.0990—0.0760-0.0680 —0o07800-0o0140-0.01400o00100o03600o02700.0180—0.0010-0o0560-0o0610-0o1050

-0.1230—0.1230-0.1220 -0。1140 -0.0970 —0.0870 -0。11100-0.0130-0.00500.0170 0.0490 -0。0140 -0。0510 —0。0540—0.1050-0。1120—0。0970—0。1010-0.1400—0。1410—0。1270-0.1270-0.1070—0。11800-0.0080 0.0170 0。0490 0。0120-0.0640-0。0950-0。1100—0.1330-0.1250—0.1220 —0。1130 -0。1220 -0.1270-0。1360-0。1190-0。11800 0—0。0200 0。0120 0。0150-0.0010-0.0480-0。1040-0.1030 —0.1270 —0.1300 —0.1160 -0。1180-0。1130—0。1010—0.1000-0.1060—0。1240 —0.1080 0 -0.02800。00700。01400.01200.0020-0。0340—0.0710 —0.0810 -0。08000。1020-0.1150—0。0980-0.0950 —0.0890 -0.1080 -0。1380 -0.11200-0.03200。01200。0250 0。0120 0.0010 —0.0190-0。00800。0440—0。0620—0.0770 -0.0810 —0.0630 —0。0530 -0。0810 -0。0770-0。1310—0。0780-0.0640 —0。0210—0。0090 0.0210 0.03100.02400.0040—0。0090—0。0370 —0。0610 —0。0560 -0。06400。0650-0.0720—0。0840-0.0850 -0。0520 -0.0660 -0。0780 —0.0110-0。0260-0.0100—0.01100。0280 0。0180 0.0240 0 —0。0230-0。0420-0。0660-0.0630—0。0650 -0.0960 0.0160 -0.0340 —0.0480—0.0630-0。0080—0.0110-0。0070-0。0050 0。0130 0。0350 0.04500。06000。05500.0280—0。0230 —0。0590 -0.0950 -0。0630 0.00700.00400。0100-0.0140（4）分析第一模态空间场theSSTSpacialFieldforEOFAnalysis/N0.1miHIDE昉匠1B0 IO1201 IDON 8DW时间系数the:?STTimi?CO'?"'fo-EOF^ralv?^/No.1ifl&o iaaa ig口 a-am iwid i-sna此次试验EOF分析中的前两个特征向量最大限度地表征了海平面温度场的主要结构。第一特征向量所描绘的第一经验正交函数的特征场（即第一模态）具有海表面气温516个样本的最相似的特征。若其可以解释为516个月的标准化特征，它指示出海表温度变化的扰动。其对应的时间系数可以表示为第一模态空间场的时间权重。从第一模态的空间特征场可以看出，受到大尺度环流影响，整场的空间变化基本全为负值。而其值乘以时间权重后均变为负值。也就是大的时间系数乘以空间特征值对应海表温度的低值，而小的时间系数乘以空间特征值则对应高值。海表温度的低值对应了气象上的拉尼娜年，而海表温度的高值对应了厄尔尼诺年。厄尔尼诺现象泛指赤道附近的东部太平洋表层海水温度上升引起的气候异常现象，它是热带海洋洋流与大气互作用的产物。其基本特征是太平洋沿岸的海面水温异常升高，海水水位上涨，并形成一股暖流向南流动。它使原属冷水域的太平洋东部水域变成暖水域，结果引起海啸和暴风骤雨，造成一些地区干旱，另一些地区又降雨过多的异常气候现象。所以，在空间特征场乘以时间系数后的高值表示厄尔尼诺年。拉尼娜现象是指海洋中的赤道的中部和东部太平洋，东西上万公里，南北跨度上千公里的范围内，海洋温度比正常温度东部和中部海面温度偏低0。2摄氏度，并持续半年（与厄尔尼诺现象正好相反），东南信风将表面被太阳晒热的海水吹向太平洋西部，致使西部比东部海平面增高将近60厘米，西部海水温度增高，气压下降，潮湿空气积累形成台风和热带风暴，东部底层海水上翻，致使东太平洋海水变冷的现象。所以，在空间特征场乘以时间系数后的低值表示拉尼娜年。实习三（附加）计算给定瞄的11年滑动平均和累积距平1、 资料介绍利用数据ma.dat,编写11点滑动平均的程序,ma.for给出了阅读资料的fortran程序。数据在文件夹中单独给出。2、 要求实习报告中附出程序，并给出原数据和滑动后数据的图形(1张图)Matlab程序load’g：\MA。DAT'x=MA；year=[1922:1:2006]';year2=year(1+(ih-1)/2：length(x)—(ih—1)/2)；ih=11；fori=1:length(x)-10avr(i)=sum(x(i：i+10))/ih;endplot(year,x,'b:'holdonplot(year2,avr,'r')save('Exam_4_output_data.txt','avr',’-ascii’)3、 实习结果①、FORTRAN程序滑动平均计算值(已导入文件Exam_4outputfile_DATAdat)3.0727272。9909102.9181822。9000002.7818193。0454552。9545462。9363642.8545462。7363642。 7454552.7545462。7272732.7272732.7545462。7363642.7727282。7727282.7636372。7727282。7636372.7272732.6818182。6636372。6636372.6181822.5727272。5363642。5272732.5363642。5727272。6090912.6545462。 6818182.7545452.7909092.8818183.0090913。 0363643.0727273.1181823.0909093.1090913。1090913。1454553。1272733。1090913。0818193。0545463.0818193.1181823.1363643.1818193。1545463.2000013.1727283.2181833.2272733.1818193。1636373.1363643.1000013。0636373。0181823.0727273。0545463.0818183。0090913.0000003。0090912.9727272.9909093。0545463.0727273.045455②、原始数据及实验滑动平均拟合曲线图表（图表由matlab画出）Figure1;Exarri_4_MQ¥rig_^ei-igfC-_[igMomngAvBrageht>aiDalm—11pwMnvAh"gFigure1;Exarri_4_MQ¥rig_^ei-igfC-_[igMomngAvBrageht>aiDalm—11pwMnvAh"g③、分析图中滑动平均值彳艮好的描述了变量x随时间的变化趋势.滑动平均值滤掉了较大的震荡，使趋势更加明显.从图中看出，x的整体趋势被体现而震荡的极大值却在趋势线中没有显著表示。所以，滑动平均在分析数据时可以更好的体现变化趋势但无法体现较大的异常值。在分析异常时，不宜使用滑动平均。

(枷)求给定数据的多元线ffl^^g(枷)求给定数据的多元线ffl^^g1、说明%一%为四个预报因子，J为预报量；样本个数n=132、要求选取预报因子1、2、4,求预报量的标准化回归方程。•i12345678910111213%7111117113122111110%26295631525571315447406668%615886917221842398%46052204733226442226341212y78.574。3104.387.695。9109.2102。772。593。1115.983.8113。3109.4答案：标准化变量回归方程：宁=0.5679气+0.4323气°.2613七附加：利用上表资料，尝试编写逐步回归程序。3、实习结果:(1)Matlab源程序ClimateData=xlsread('F：\—O®\—。㈤伍\气象统计方法'实验数据'气象统计实验四数据(附加题)°xls')；%从Excel文件读取数据X=ClimateData(2:5,:);%提取ClimateData的第二到5行数据，即自变量观测矩阵X

y=ClimateData(6,:)； %提取ClimateData的第六行，即预报量reglm(y,X)%调用reglm函数做4重线性回归，显示回归分析的方差分析表和参数估计表functionstats=reglm(y,X,model,varnames)%多重线性回归分析或广义线性回归分析%%reglm(y,X),产生线性回归分析的方差分析表和参数估计结果，并以表格形式显示在屏幕上.参% 数X是自变量观测值矩阵，它是n行p列的矩阵.y是因变量观测值向量，它是n行1列的列向量.%%stats=reglm(y,X),还返回一个包括了回归分析的所有诊断统计量的结构体变量stats.%stats=reglm(y,X,model)，用可选的%stats=reglm(y,X,model)，用可选的model参数来控制回归模型的类型。model是一个字符串，%%%%%%%标签.%其可用的字符串如下’linear’’interaction'’quadratic’'purequadratic'带有常数项的线性模型(默认情况)带有常数项、线性项和交叉项的模型带有常数项、线性项、交叉项和平方项的模型带有常数项、线性项和平方项的模型stats=reglm(y,Xvarnames可以是字符矩阵或字符串元胞数组，它的每行的字符或每个元胞的字符串是一个变量model,varnames)，用可选的varnames参数指定变量的标签，它的行% 数或元胞数应与X的列数相同.默认情况下，用X1,X2,…作为变量标签。%%例：%x=[215250180250180215180215250215215136。5136.5136。5138.5139。5138.5140.5140.5140。5138。5138.5]'；%y=[6。27.54.85.14。64.62.83。14。34。94.1]'；% reglm(y,x,'quadratic')% —-一方差分析表-— % 方差来源 自由度F值 p值%回归 % 方差来源 自由度F值 p值%回归 5。00007.6122 0.0219%残差 5.00001.9742平方和15.02770。3948均方3.0055%总计10.0000 17.0018%% 均方根误差(RootMSE)0.6284 判定系数(R-Square)0.8839% 因变量均值(DependentMean)4。7273 调整的判定系数(AdjR-Sq)0.7678%% —— — — — ——-参数估计-—%变量估计值标准误t值p值%常数项30.94282011。11170.01540.9883%X10。70400。64051.09920.3218%X2—0 。848729.1537—0.02910.9779% X1*X2-0。00580.0044—1.31320。2461% X1*X10.00030。00030.83840。4400%X2大X20.00520.10550.04920.9626%%Copyright2009-2010xiezhh.% $Revision:1。0.0。0$$Date： 2009/12/2221:41：00$ifnargin〈2error('至少需要两个输入参数');endp=size(