我有一个矩阵,它在主对角线上的元素没有排序,所以我需要一个能在主对角线上返回带有有序元素的新矩阵的函数。我无法理解为什么这不起作用。
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}}不是符号。 >>
答案 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
请注意,您可以内联完成所有这些操作:
ReplacePart[ A,
Table[ {i, i} -> (Sort@Diagonal[A])[[i]], {i, Length@A} ]]
或
(A + DiagonalMatrix[Sort@# - # &@Diagonal[A]])