NetLogo中patch variable空间相关长度计算及FFT扩展问询

我来帮你梳理下两种可行的实现路径，包括你提到的相关函数积分法，以及适合周期性系统的FFT方法：

直接计算空间相关函数并积分 #

这种方法直观易懂，适合中小规模的世界：

• 首先确定你的目标patch变量（比如my-variable），先计算变量的全局均值：

  let mean-val mean [my-variable] of patches
  

• 由于是周期性世界，处理patch平移时要用mod实现边界循环。比如，对于任意滞后量dx（横向偏移）和dy（纵向偏移），某个patch p对应的偏移patch是patch (pxcor + dx) mod world-width (pycor + dy) mod world-height。
• 计算每个滞后（dx, dy）对应的相关系数：

  ; 先计算变量的方差
  let var sum [ (my-variable - mean-val) ^ 2 ] of patches / count patches
  ; 定义计算指定滞后相关系数的过程
  to-report corr-at-lag [dx dy]
    let mean-val mean [my-variable] of patches
    let cov sum [ (my-variable - mean-val) * (([my-variable] of patch (pxcor + dx) mod world-width (pycor + dy) mod world-height) - mean-val) ] of patches
    report cov / (var * count patches)
  end
  

• 得到所有滞后的相关函数值后，就可以积分计算相关长度：通常是按距离r = sqrt(dx^2 + dy^2)分组，对每个r取平均相关值，然后从r=0开始积分，直到相关值降到阈值（比如1/e），此时的r就是相关长度。

利用FFT计算（NetLogo扩展支持） #

对于周期性系统，FFT是更高效的方法，NetLogo确实有第三方扩展支持：

• NetLogo本身没有内置FFT，但你可以使用fft扩展——它提供了一维和二维的FFT、逆FFT功能，完全适配周期性场景。
• 实现步骤大概是：
  1. 先把patch变量转换成二维数组：比如创建一个列表，每个子列表对应一行的patch变量值：

     let patch-data []
     foreach (range world-height) [ y ->
       let row []
       foreach (range world-width) [ x ->
         set row lput ([my-variable] of patch x y) row
       ]
       set patch-data lput row patch-data
     ]
     

  2. 加载扩展后，使用fft:fft-2d对二维数组做傅里叶变换，得到复数频谱；接着计算功率谱（频谱模的平方）。
  3. 对功率谱做逆FFT（fft:ifft-2d），就能得到空间相关函数，记得归一化（除以patch总数）。
  4. 最后同样通过积分相关函数得到相关长度。
• 注意：需要先在NetLogo的扩展管理器中安装fft扩展，然后在代码开头加上extensions [fft]来启用。

实用提示 #

• 直接计算法适合小世界（比如100x100以内），但计算量随世界大小平方增长；FFT法是O(N log N)复杂度，大世界下速度优势明显。
• 处理周期性时，务必用mod确保偏移后的patch不会越界，比如(pxcor + dx) mod world-width能让最右侧的patch偏移后回到左侧。
• 积分相关函数时，要考虑二维空间的权重：相同距离r的滞后方向有很多，计算平均相关值时要按出现次数加权。

内容的提问来源于stack exchange，提问作者Enrico Calzavarini