如何生成向量的循环邻域向量数组？求高效实现方案

需要将向量v转换为向量数组，数组中每个向量的第i个元素对应原向量第i个元素的k个循环邻域值（包含自身及左右各k个循环移位后的对应元素）。例如：

• 输入v = [1, 2, 3, 4, 5]、k = 1时，期望输出数组：

  M = [
      5, 1, 2, 3, 4;  % 原向量循环左移1位
      1, 2, 3, 4, 5;  % 原向量本身
      2, 3, 4, 5, 1   % 原向量循环右移1位
  ]
  

  此时M的第一列[5;1;2]就是原向量第一个元素1的循环邻域值（左邻、自身、右邻）。

你提供的MATLAB代码存在几处逻辑错误，无法正确生成目标数组：

• 未定义变量r，aux1=[a(:,(end-r+1):end),a(:,1:(end-r))];这一行会直接报错；
• 循环for j=2:k中错误地始终给a(k,:)赋值，实际应为a(j,:)来存储对应移位后的向量；
• 最终生成的I是cell数组，且逻辑上仅生成了k个移位，未覆盖左右两侧的完整邻域（如示例中k=1需要3个移位结果）。

在MATLAB中，向量化操作是这类数组变换的最优选择，避免显式循环能大幅提升大向量场景下的性能。以下是两种高效实现方案：

方案1：基于circshift的循环实现 #

通过circshift生成每个需要的循环移位向量，拼接为结果矩阵：

function M = get_cyclic_neighbors(v, k)
    n = length(v);
    % 定义移位范围：从左移k位到右移k位
    shifts = -k:k;
    M = zeros(length(shifts), n);
    for idx = 1:length(shifts)
        M(idx, :) = circshift(v, shifts(idx));
    end
end

调用示例：

v = [1,2,3,4,5];
k = 1;
M = get_cyclic_neighbors(v, k)

输出符合预期：

M =
     5     1     2     3     4
     1     2     3     4     5
     2     3     4     5     1

方案2：矩阵索引实现（性能最优，无显式循环） #

通过构造循环索引矩阵，直接从原向量中提取元素，完全利用MATLAB底层矩阵运算优化：

function M = get_cyclic_neighbors_fast(v, k)
    n = length(v);
    % 构造循环扩展的索引矩阵，确保索引落在1~n范围内
    extended_idx = mod((1:n) + (-k:k)' - 1, n) + 1;
    % 直接索引原向量得到结果
    M = v(extended_idx);
end

调用方式与输出结果和方案1一致，且对于长度大于1e4的向量，该方案的执行速度远快于方案1。

方案2逻辑解释 #

• (-k:k)'生成列向量[-k; -k+1; ...; 0; ...; k-1; k]，对应每个邻域的移位偏移；
• (1:n) + (-k:k)'生成(2k+1)×n的矩阵，每个元素为原索引加上偏移量；
• mod(..., n) + 1处理循环边界，将超出范围的索引映射回1~n；
• 最终通过索引矩阵直接提取原向量元素，得到邻域值矩阵。

内容的提问来源于stack exchange，提问作者Rodrigo G. Garrastázul