3 对周期分析进行统计检验(即F检验),确定你站逐日降水资料是否存在可用周期。
① 首先,“预报方法”选择“方差周期叠加”,“计算项目”选择“逐日降水”,“结果显示”选择“数值列表”,之后点击“计算”按钮,进行计算。瞬间,计算完成后,可显示计算结果,其中上部即为方差分析表,见图-1。方差分析表中,有四个周期,下面我们仅以第1周期为例,说明一下F检验的过程。
首先,“预报方法”选择“方差周期叠加”,“计算项目”选择“逐日降水”,“结果显示”选择“数值列表”,之后点击“计算”按钮,进行计算。瞬间,计算完成后,可显示计算结果,其中上部即为方差分析表,见图-1。方差分析表中,有四个周期,下面我们仅以第1周期为例,说明一下F检验的过程。
图-1 计算结果中的方差分析表
② 第1周期的F检验
第1周期为年周期(这里近似372天),F(372)=1.6341575374552, 自由度n1(372)=371, 自由度n2(372)=1503。
图-2 给出了F检验表在主控页面中的位置。
图-2 F检验表所在主控页面上的位置图
③ 点击图-2 右侧F表中的0.01,调出F=0.01的检验表(见图-3)。表中顶部第一行为自由度 n1,左部第一列为自由度 n2,表中主体数值就为F 0.01 的标准值。当自由度n1(372)=371, 自由度n2(372)=1503时,可查到F 0.01 近似等于1.19(表中,n1没有371,但与500最接近,n2没有1503,但与1000最接近。从n1=500向下引线,从n2=1000向右引线,相交处的数值为1.19)。
由于F(372)=1.6341575374552 大于了F 0.01 的值1.19,达到极显著(**)标准,可以认为第一周期是存在的,可以用于预报。
图-3 F=0.01 的检验表
④ 用同样的方法,我们还可以调出F=0.05和0.10的检验表查出F0.05和F0.10 。
经查,第2周期,F(14)=1.903312881606,大于F0.05的值,达到显著(*)标准,也可以认为第2周期是存在的,可以用于预报。第3周期,F(9)=1.7191466452073,大于 F0.10 的值,达到稍显著(#)标准,可以认为第3周期是可能存在的,可以参考使用。第4周期,F(13)=1.5238143057059,也大于 F0.10 的值,达到稍显著(#)标准,也可以认为第4周期是可能存在的,可以参考使用。
这样,我们已经证明该站逐日降水资料确实存在可用周期,则可输出预报结果,用来预报未来的逐日降水。
*
“农业气象信息及预报工具平台”可在下面两个地址找到:
http://www.nkyzws.com/v6/bencandy.php?fid=67&id=1170
或
http://lidm3.iiskj.com/lxyb2
本“平台”专供气象台站内部使用,不对外开放。研究团队每年对“平台功能”都有更新。有意参加此项研究或使用此平台的气象台站,可与我们联系。
联系办法:
通信地址:黑龙江省佳木斯市安庆街798号,黑龙江省农垦科学院作物所 联系人:李德明
邮编:154007
电子信箱:391039782@qq.com
联系电话:15590919407 QQ:391039782