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

南京信息工程大学实验（实习）报告 实验（实习）名称 短期气候预测实习 日期11.29得分 指导教师 系 ~ 专业 ~ 年级 ~ 班次~姓名Trichtu 学号~ 实习目的： 掌握短期气候预测中物理统计预测的基本步骤。 实习要求： 能运用提供的资料和方法子程序，编写或补充完成程序当中的部分片断，了解区域降水的预测方法及其建立过程，输出实验要求的相应结果，并就方法对区域降水的拟合及试验预测效果进行分析。 实习内容： 利用前期1月的海温关键区（Nino3.4指数）和环流特征量（西太平洋副高脊线、西太平洋副高西伸脊点、亚洲极涡面积、南方涛动指数）等前期冬季预测因子；运用多元回归方法，对1952-2001年华北夏季降水建立预测方程；对2002-2008年进行多元回归预测试验； 实习资料： 前期1月的Nino3.4指数（来自CPC） 西太平洋副高脊线、西太平洋副高西伸脊点、亚洲极涡面积、南方涛动指数（来自中国气象局整编的74个环流指数）、夏季华北区域10站的降水量距平百分率。 实习方法： 回归分析(mregrssion.for)是用来寻找若干变量之间统计关系的一种方法，利用所找到的统计关系对某一变量作出未来时刻的估计，称为回归预报值。 效果分析—— 回归拟合效果的参数分析 残差平方和 (SSR)误差方差标准差 (2)标准差 (3)复相关系数 当 R 近似等于1，则相对误差将近似0，说明回归效果很好。 回归方差 反映回归拟合的程度,其值越大反映效果越好 (5)总离差平方和 它遵从自由度为m和n-m-1的F分布。 实习步骤： 编写程序：（原文件的子例行程序不计入在内） PROGRAM MAIN INTEGER,PARAMETER::N=50 INTEGER,PARAMETER::K=5 REAL,DIMENSION(K,N)::X REAL,DIMENSION(N)::Y REAL,DIMENSION(K+1)::A REAL,DIMENSION(K+1,K+1)::B REAL,DIMENSION(K)::V REAL Q,S,R,U,ind(6,60),year(N),expect(50) C OPEN THE INPUT DATA FILE open(8,FILE=e:\copy6\weight.txt) open(9,FILE=e:\copy6\compare.grd,form=binary) OPEN(10,FILE=e:\copy6\shixi.txt) C READ THE DATA and give data to X and Y read(10,*) do i=1,N read(10,*) year(i),y(i),x(1,i),x(2,i),x(3,i),x(4,i),x(5,i) end do MM=K+1 call DYHG(X,Y,K,MM,N,A,Q,S,R,V,U,B,DYY) !CCCCCCCCCCCCCCCCCCC预测1952-2001 do i=1,N expect(i)=a(1)+x(1,i)*a(2)+x(2,i)*a(3)+x(3,i)*a(4) expect(i)=expect(i)+x(4,i)*a(5)+x(5,i)*a(6) enddo write(9) ((expect(i),y(i)),i=1,50) write(*,88) A(1) 88 format(/1x,b 0=,f19.5) do 89 j=2,MM 89 write(*,100) j-1,A(j) 100 format(1x,b,i2,=,f9.5) cccccccccccccccccccccccccccccccccccccccccccccccccccccc write(*,20)Q,S,R 20 format(1x,Q=,f13.6,3x,S=,f13.6,3x,R=,f13.6) write(*,22)U,DYY 22 format(1x,U=,f13.6,3x,DYY=,f13.6) write(*,30)(i,V(i),i=1,K) 30 format(1x,V(,i2,)=,f13.6) write(*,40)U 40 format(1x,U=,f13.6) open(6,file=e:\copy6\table.txt) ! output dat