利用ncl计算两个自相关时间序列的有效自由度

AXD · 发表于 2020-10-6 09:16:20

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 AXD 于 2020-10-6 09:21 编辑

; Calculate the effective degrees of freedom
; Reference: Jianping Li et al. (2013, GRL: 40)
N = dimsizes(y_reg)
temp = new(N,float)
do i = 0,N-1
  mxlag = i
  acr1 := esacr(y_reg,mxlag)
  acr2 := esacr(SATI,mxlag)
  temp(i)  = ((N-i)/N)*acr1(i)*acr2(i)
end do
Neff = N/(1+2*dim_sum(temp))
print(Neff)

其中，y_reg代表序列1，SATI代表序列2。附件为计算方法参考文献。

夏无且 · 发表于 2023-10-18 15:04:28

参考原公式后，感觉 i 应该是从1开始循环样

niuyingli · 发表于 2020-12-13 21:12:32

我也用了这个，为什么会出现这样的错误呢
fatal:Loop end must be scalar, can't execute loop

niuyingli · 发表于 2020-12-13 21:24:12

其中一个自相关算出来是三维，另一个自相关是3维，就导致错误，该怎么办呢

JGD飘 · 发表于 2021-11-9 15:18:04

如果是要计算时间序列和空间场之间的有效自由度，该如何修改呢

xxr · 发表于 2023-9-20 20:40:26

请问i=N-1的时候，算出来的相关居然大于1，这是什么问题啊

		自动登录	找回密码
密码			立即注册

[其他] 利用ncl计算两个自相关时间序列的有效自由度

登录后查看更多精彩内容~