matlab幻方变换_矩阵和幻方矩阵 - MATLAB & Simulink - MathWorks 中国

矩阵求和、转置和对角矩阵

您可能已经注意到，幻方矩阵的特殊属性与元素的不同求和方法相关。如果沿任何行或列求和，或者沿两条主对角线中的任意一条求和，您将始终得到相同数字。让我们使用 MATLAB 来验证这一点。尝试的第一个语句是

sum(A)

MATLAB 返回的结果为

ans =

34 34 34 34

如果未指定输出变量，MATLAB 将使用变量 ans(answer 的缩略形式)来存储计算结果。您已经计算包含 A 的列总和的行向量。每个列的总和都相同，即幻数和 34。

行总和如何处理？MATLAB 会优先处理矩阵的列，因此获取行总和的一种方法是转置矩阵，计算转置的列总和，然后转置结果。

MATLAB 具有两个转置运算符。撇号运算符(例如，A')执行复共轭转置。它会围绕主对角线翻转矩阵，并且还会更改矩阵的任何复数元素的虚部符号。点撇号运算符 (A.') 转置矩阵，但不会影响复数元素的符号。对于包含所有实数元素的矩阵，这两个运算符返回相同结果。

因此

A'

生成

ans =

16 5 9 4

3 10 6 15

2 11 7 14

13 8 12 1

而

sum(A')'

生成包含行总和的列向量

ans =

34

34

34

34

有关避免双重转置的其他方法，请在 sum 函数中使用维度参数：

sum(A,2)

生成

ans =

34

34

34

34

使用 sum 和 diag 函数可以获取主对角线上的元素的总和：

diag(A)

生成

ans =

16

10

7

1

而

sum(diag(A))

生成

ans =

34

从数学上讲，另一条对角线(即所谓的反对角线)并不是十分重要，因此 MATLAB 没有对此提供现成的函数。但原本用于图形的函数 fliplr 可以从左往右地翻转矩阵：

sum(diag(fliplr(A)))

ans =

34

您已经验证丢勒雕刻中的矩阵确实是一个幻方矩阵，同时在验证过程中，您已经尝试了几个 MATLAB 矩阵运算。下面各部分继续使用此矩阵来演示 MATLAB 的其他功能。