Jianping Li's Climate Indices

常见问题 Frequently Asked Questions

一维离散功率谱分析 Discrete spectrum analysis of one-dimensional series

dspectrum(n,lw,x,tl,sl,st95)
Q:	在计算滞后相关的模块中，求滞后相关为什么用n而不用n-i？
A:	这是由实际资料的特性和统计学的基本原理决定的（很多统计学书中也谈到这一点）。一是序列的平稳性假设就要求序列的均值和方差保持不变；二是通常资料太短，导致用较短的资料得到的相关结果具有不稳定性，掩盖了事物的真相。因此，用n要好于n-i。一些同学也做了试验证实了这一点（有时n-i还得到错误的结论，需要各自使用者注意的地方）。

旋转经验正交函数分解 Rotated Empirical Orthogonal Functions

reof(m,n,mnl,np,f,ks,er,egvt,ecof,rer,regvt,recof)
Q:	缺省值和奇异点的去除？
A:	程序中给定了缺省值，变量名为undef，可以根据具体数据作相应改变，如本程序中为32767；程序中的变量mg1为原始数据去掉了缺省值之后的有效格点数；而mg2为在去掉了缺省值的基础上，进一步去掉了奇异点之后的有效格点数（即方差等于0的点）。
Q:	关于nd的取值？
A:	若序列为逐月时间序列，则取nd=12；若为逐日时间序列，则取nd=365；若为季节平均的时间序列，则取nd=1；以此类推。
Note:	经验正交函数分解的子程序eof(m,n,mnl,f,ks,er,egvt,ecof) 中，没有包含去缺省值的子程序部分，若原始数据中包含缺省值，则需要预处理。

如有任何问题请联系我：Questions, comments, suggestions?

Last updated: 2010-07-01 (Jianping Li: ljp@lasg.iap.ac.cn)

dspectrum(n,lw,x,tl,sl,st95)