time_analysis.py 文件源码

def corr(data):  
  ns = data.shape[0];
  nt = data.shape[1];

  pairs = make_pairs(ns);
  npp = len(pairs);

  mean = np.nanmean(data, axis = 0);
  var = np.nanvar(data - mean, axis = 0);

  c = np.zeros(nt);
  for p in pairs:
    c += np.nanmean( (data[p[0]] - mean) * (data[p[1]] - mean), axis = 0) / var;
  c /= npp;

  return c;