R中多個矩陣的加法

#第一種算法，採用循環的方法

 matrices <- vector(mode="list", length=10000)#生成一個空的列表，用於存儲多個矩陣

for (i in 1:10000) {matrices[[i]] <- matrix(rnorm(10000), 100, 100)}#將生成的多個矩陣存儲到列表中

  system.time({

   S <- matrix(0, 100, 100)

               for (i in 1:10000) S <- S + matrices[[i]]

 })  #對多個矩陣求和，並計算所用時間

system.time(S1 <- Reduce("+", matrices))#第二種算法，Reduce函數將函數'+'依次鏈接列表中的矩陣

range(S1 - S)#檢驗兩種算法的結果是否一致，即求差值的最大值和最小值是否均爲0