R数组中的元素矩阵代数

Question

我有两个尺寸为arry1[2,2,n]和array2[2,2,n]的数组。也就是说，n两两个矩阵。

我想生成第三个数组，它是array1和array2的n元素矩阵乘法（即％*％）。生成另一个维度为array3[2,2,n]的数组。

令人沮丧的是，我无法弄清楚如何使用%*%来解决此问题，以下似乎无法正常工作

array3 <- array1[1:2,1:2,]%*%array2[1:2,1:2,]

此外，apply()系列函数似乎无法启用我的操作。 任何帮助将不胜感激。

Answer 1

如果我理解正确，这将有效

array(sapply(1:5,function(x) a1[,,x]%*%a2[,,x]),dim = c(2,2,5))

数据

a1=array(outer(outer(1:2,1:2),1:5),dim = c(2,2,5))
a2=array(outer(outer(1:2,1:2),1:5),dim = c(2,2,5))