短期气候预测实习报告六.docx

《短期气候预测实习报告六.docx》由会员分享，可在线阅读，更多相关《短期气候预测实习报告六.docx（5页珍藏版）》请在第一文库网上搜索。

1、南京信息工程大学实验实习报告实验实习名称短期气候预测实习日期11.29得分指导教师系二专业二年级二班次;姓名TriChtU学号实习目的：掌握短期气候预测中物理统计预测的基本步骤.实习要求：能运用提供的资料和方法子程序;编写或补充完成程序当中的部分片断;了解区域降水的预测方法及其建立过程;输出实验要求的相应结果;并就方法对区域降水的拟合及试验预测效果进行分析.实习内容：利用前期1月的海温关键区Nino3.4指数和环流特征量西太平洋副高脊线、西太平洋副高西伸脊点、亚洲极涡面积、南方涛动指数等前期冬季预测因子;运用多元回归方法;对1952-2001年华北夏季降水建立预测方程；对2002-2008年进

2、行多元回归预测试验；实习资料：前期1月的Nino3.4指数来自CPC西太平洋副高脊线、西太平洋副高西伸脊点、亚洲极涡面积、南方涛动指数来自中国气象局整编的74个环流指数、夏季华北区域10站的降水量距平百分率.实习方法：回归分析mregrssion.for是用来寻找若干变量之间统计关系的一种方法;利用所找到的统计关系对某一变量作出未来时刻的估计;称为回归预报值.效果分析一一回归拟合效果的参数分析(I)残差平方和SSR误差方差标准差Q=SSR=f(y_必)2=力比一(4+4%+b2x2+zm)2ii=2标准差S=ySSR/n3复相关系数R=I-ssrj(yi-y)2Vi=1当R近似等于1;则相对误

3、差将近似0;说明回归效果很好.(4)回归方差U=Zu(%+1x1+b2x2+bmxmi)2/=1反映回归拟合的程度;其值越大反映效果越好5总离差平方和DYY=U+QF_U1mQ/(-m-1)它遵从自由度为m和n-m-1的F分布.实习步骤：(1)编写程序：原文件的子例行程序不计入在内PROGRAMMAININTEGERPARAMETER:N=50INTEGERPARAMETER:K=5REA1;DIMENSTONK;N:Xrea1；Dimensionn：yREA1;D1MENS1oNK+1:Area1；dimension+i；k+i：brea1；Dimenstonk：vREA1Q;S;R;U;i

4、nc16;60;yearN;expect50COPENTHEINPUTDATAFI1EOPCn8;FI1E=e:copy6weight.txt,OPCn9;F11E=e:copy6compare.grd,;form=，binary,OPEN1O;FT1E=,e:copy6shixi.txt,CREADTHEDATAandgivedatatoXandYread1;*doi=1;Nread1;*yeari;yi;x1;i;x2;i;x3;i;x4;i;x5;ienddoMM=K+1ca11DYHGX;Y;K;MM;N;A;Q;S;R;V;U;B;DYYCCCCCCCCCCCCCCCCCCC预测19

5、522001doi=1;Nexpecti=a1+x1;i*a2+x2;i*a3+x3;i*a4expccti=expecti+x4;i*a5+x5;i*a6enddowrite9expecti;yi;i=1;50write*;88A188format/1x;,b0=，;f19.5do89j=2;MM89write*;100j-1;Aj100formatIx;,b,;i2=;f9.5ccccccccccccccccccccccccccccccccccccccccccccccccccccccwrite*;20Q;S;R20formatIxQ=;f13.6;3x/S=;f13.6;3x;R=;f13

6、.6write*;22U;DYY22formatIx/U=;f13.6;3x/DYY=;f13.6write*;30i;Vi;i=1;K30format1x/V*;f13.6write*;40U40format1x;,U=，;f13.6open6;fi1e=，e:copy6tab1e.txt,outputdatawrite6J80180formatx/regressioncoefficients/write6;88A1do189j=2;MM189write6;100j-1;Ajwrite6;200200format/1x;,GenericAna1ysisofVarianceTab1eforth

7、eMu1tip1e* 1inearRegression,write6;202202format/1x* ,write6;204204format3x;,SourcedfSSMS,write6;202write6;206NT;DYY206format/1x;,Tota1nT=;i2;SST,;f13.4u2=Urea1Kwrite6;208K;U;U2208format/1x/RegressionK=;i2;SSR=;f13.4;MSR=SSRK,*1f13.4q2=qrea1n-k-1write6;209n-k-1;q;q2209format/1x/Residua1n-k-1=，;i2;,SS

8、E=;f13.4;MSE=SSEn-k-1* =;f13.4f=Urea1KQrea1N-K-Iwritc6;220f220format/1x;F=MSRMSE=;f13.4write6;202c1osc6stopend2结果输出得到回归效果的参数及回归系数;拟合预测与观测资料的对比数据;独立预测试验与观测资料的对比数据.输出反映回归效果的参数及回归系数;并就相关参数分析回归效果；回归系数：h=-196.06274b1=4.40570b2=1.19246b3=-0.00290b4=0.32351b5=0.50310即：Nino3.4指数=-196.06274+4.40570*西太平洋副高脊线+

9、1.19246*西太平洋副高西伸脊点-0.00290*亚洲极涡面积+0.32351*南方涛动指数+0.50310*夏季华北区域IO站的降水量距平百分率GenericAna1ysisOfVarianceTab1efortheMu1tip1e1inearRegressionTota1n-1=49DYY=20596.3984RegressionK=5U=1492.2344MSR=SSRK=298.4469Residua1n-k-1=44Q=19104.5078MSE=SSEn-k-1=434.1934R=0.269137F=MSRMSE=0.6874相关为0.29137;F小于标准值2.40.05且

10、大于0.416;效果不显著;说明拟合程度不好.预测量与回归方程计算的估计值和观测值的历年曲线变化图I9522OO1年;并附简单的说明；Ct1文件：,reinitopene:copy6compare.ct1enab1eprinte:copy6comparation.gmf,setx1sety1,set1150,drea1dexpprintdisab1eprint,Gs文件：,reinitopene:copy6compare.ct1enab1eprinte:copy6comparation.gmf,setx1sety1,set1150),drea1dexpprintdisab1eprint实际值黑线、预测值绿线从图中可以看出;某些年份;预测值和实际值的趋势是反向的;实际的极大值在预测中是极小值.所以预测效果不好.输出独立预测试验的观测与预测值下表.2002-2008夏季降水预测试验结果年份2002200320042005200620072008预测值-9.70-7.80-8.074.18-1.62-10.835.885实际值29.75-17.450.6422.96-22.96-19.403.30