2014-12-02 42 views
-1

假设矩阵A有5000行6000列,其中每个元素为a_ij,还有另一个矩阵S,它有6000行6000列,它是一个对称矩阵,其中每个元素是s_lk,我想计算m_ij = sum(k=1..6000, l=1..6000)(a_il * a_jk * s_lk)(即this formula),它应该是5000 * 5000的矩阵。做这样的问题有什么好处(多重数字产品是双倍数额)?我们应该先简化它吗?R:处理“三数乘积和”的策略

非常感谢。

+0

等一下。前两天我们没有看到相同的问题吗? – 2014-12-02 21:41:22

+1

我想你正在寻找'a%*%S%*%t(a)' – nicola 2014-12-02 21:44:34

+0

这是两个数字的产品。有一点不同,但对我而言,我仍然无法解决当前的问题。 – BioChemoinformatics 2014-12-02 21:44:41

回答

0

@nicola,尼科拉给了正确的anwser。谢谢。

A %*% S %*% t(A)

在这里,我只是张贴Nicola的答案。如果不适合,请删除它。