对矩阵的主对角线进行排序

时间:2015-04-10 13:50:28

标签: sorting matrix wolfram-mathematica

我有一个矩阵,它在主对角线上的元素没有排序,所以我需要一个能在主对角线上返回带有有序元素的新矩阵的函数。我无法理解为什么这不起作用。

   Function[A_] := Module[{res, diagonal = {}, m, n},
   {m, n} = Dimensions[A];
   Table[AppendTo[diagonal, A[[i, i]]], {i, 1, m}];
   dijagonal = SelectionSort[diagonal];
   Table[A[[i, i]] = dijagonal[[i]], {i, 1, m}];
   Return[A // MatrixForm];
   ];

选择排序有效。

这可以是矩阵的一个例子:

   A={{60, 10, 68, 72, 64},{26, 70, 32, 19, 29},{94, 78, 86, 59, 17},
 {77, 13, 34, 39, 0}, {31, 71, 11, 48, 83}}

当我运行它时,会显示:

Set :: setps:{{60,10,68,72,64},{26,70,32,19,29},{94,78,86,59,17},{77,13,零件赋值中的34,39,0},{31,71,11,48,83}}不是符号。 >>

1 个答案:

答案 0 :(得分:1)

主要问题

  • 您无法定义自己的名为Function的函数。 (通常避免使用任何以大写字母开头的符号以避免冲突)
  • 您无法修改输入参数,因此请复制
  • 只使用Sort而不是SelectionSort

我做了一些其他更具风格的变化:

  function[A0_] :=
     Module[{res, diagonal = {}, m, n, A = A0},
     {m, n} = Dimensions[A];
     diagonal = Table[A[[i, i]], {i, 1, m}];
     dijagonal = Sort[diagonal];
     Do[A[[i, i]] = dijagonal[[i]], {i, 1, m}]; A]

 function[A] // MatrixForm

enter image description here

请注意,您可以内联完成所有这些操作:

 ReplacePart[ A, 
      Table[ {i, i} -> (Sort@Diagonal[A])[[i]], {i, Length@A} ]]

 (A + DiagonalMatrix[Sort@# - # &@Diagonal[A]])