如何在R中找到矩阵的所有等效矩阵?

时间:2015-04-27 15:35:57

标签: r matrix permutation equivalent

我想找到矩阵的所有等价矩阵,假设它是M, 但我希望这些等效矩阵与M具有完全相同的对角元素,并且与M具有相同的非对角元素。

例如,如果M是

> M<-matrix(x,4,4);M
     [,1] [,2] [,3] [,4]
[1,]   22   -2   -2    2
[2,]   -2   22    2    2
[3,]   -2    2   22   -2
[4,]    2    2   -2   22 

我想要一个等效矩阵,例如

     [,1] [,2] [,3] [,4]
[1,]   22    2   -2    2
[2,]    2   22    2   -2
[3,]   -2    2   22   -2
[4,]    2   -2   -2   22

正如你不仅可以看到M,而且所有等效矩阵必须是对称的(M =转置(M))。

我尝试了很多方法,比如“permut”和“Matrix”套餐,但我找不到能帮到我的东西。 我使用过LU和QR分解,但这不是我想要的。

这个问题有什么关系???

谢谢大家!

0 个答案:

没有答案