从另一个矩阵的项总和创建矩阵

Question

我正在尝试使用矩阵近似另一个矩阵。

我有一个名为 Bpp 的1000x96矩阵，我需要以此形式创建一个新的矩阵 Omega ：

其中b''（z_i）_l是Bpp的第i，l个条目（Bpp是1000x96矩阵）

我有这种尝试：

Omega = matrix(0, 96, 96)

for(k in 1:96){
   for(l in 1:96){
     Omega[k,l] = sum(Bpp[,k]*Bpp[,l]*delta)
   }
}

但是我几乎可以肯定这是不对的，因为它不能产生我以后需要解决的问题。

预先感谢您的帮助/指导。

Answer 1

我认为您想做内部产品。它将Bpp的每一列及其全部96列相乘并求和，这与循环的作用相同。

Omega = t(Bpp)*Bpp*delta