R语言中如何高效将固定对角矩阵添加到矩阵列表？

你已经创建了包含5个3×3随机矩阵的列表：

l1 = list()
for(i in 1:5){
 l1[[i]] = matrix(rnorm(9),3)
}

同时定义了一个3阶对角矩阵：

My_D = diag(c(1:3))

尝试了l1 + My_D和lapply(l1, function(i) l1[[i]] + My_D)两种方法都无效，下面是正确的高效实现方式：

错误原因分析 #

• l1 + My_D：列表无法直接与矩阵执行元素级加法，R不支持这种操作，会直接报错。
• 第二种lapply写法逻辑错误：匿名函数的参数i是列表中的每个矩阵元素，而非索引，你却用l1[[i]]去取元素，自然会出错。

正确实现方法 #

方法1：修正lapply的基础写法 #

直接用匿名函数参数指代列表中的每个矩阵，执行加法即可：

l1_updated = lapply(l1, function(x) x + My_D)

lapply会遍历列表里的每一个矩阵，完成加法后返回新的更新列表，这是R原生的高效实现方式。

方法2：用purrr包的map简化代码 #

如果你熟悉tidyverse生态，purrr包的map函数写法更简洁：

library(purrr)
l1_updated = map(l1, ~ .x + My_D)

这里~代表匿名函数，.x指代列表中的每个矩阵元素，代码更易读。

方法3：极简的lapply写法 #

利用R的函数式特性，直接把加法运算符作为函数传递：

l1_updated = lapply(l1, `+`, My_D)

这种写法省略了匿名函数，把My_D作为加法的右操作数，执行效率和第一种方法一致，代码更紧凑。

验证效果可以查看更新后列表的任意元素，比如l1_updated[[1]]，能看到原矩阵的对角元素分别增加了1、2、3，其他元素保持不变。

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