我想找到矩阵的所有等价矩阵,假设它是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分解,但这不是我想要的。
这个问题有什么关系???
谢谢大家!