如何对矩阵列表按每3个一组批量执行指定计算并生成结果列表？

没问题，我来帮你搞定这个批量分组计算的需求！要实现每3个矩阵一组自动计算并把结果存入新列表，我们可以分两步来做：

第一步：自动生成分组索引 #

首先我们需要把原列表的元素按每3个一组拆分，用split()结合ceiling()就能轻松实现，不管你的列表长度是多少（只要是3的倍数，这里刚好12个）：

# 生成每3个一组的索引分组
groups <- split(seq_along(L), ceiling(seq_along(L)/3))

运行后groups会是一个包含4个元素的列表，每个元素对应一组矩阵的位置索引（比如第一个元素是c(1,2,3)，第二个是c(4,5,6)，以此类推）。

第二步：批量计算每个分组的结果 #

接下来用lapply()遍历每个分组，对每组的矩阵子集执行你需要的计算逻辑。这里要修正你之前的小笔误：应该先对组内矩阵求和，再除以组内所有矩阵行和的总和，完整代码如下：

# 假设你的矩阵列表是L（替换成你实际的列表）
L <- list(M1, M2, M3, M4, M5, M6, M7, M8, M9, M10, M11, M12)

# 生成分组索引
groups <- split(seq_along(L), ceiling(seq_along(L)/3))

# 批量计算并存储结果到新列表
result_list <- lapply(groups, function(idx) {
  # 取出当前分组的矩阵子集
  current_group <- L[idx]
  # 计算组内矩阵的和
  sum_mat <- Reduce(`+`, current_group)
  # 计算组内每个矩阵行和的总和
  total_row_sums <- rowSums(sapply(current_group, rowSums))
  # 最终计算结果
  sum_mat / total_row_sums
})

这样result_list里就会包含4个计算后的矩阵，分别对应你要求的4个分组。而且这个方法是通用的——如果之后你想按其他数量分组，只需要把ceiling(seq_along(L)/3)里的3换成你想要的分组大小就行。

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