时间戳

首页 > 解决方案 > 对称矩阵乘积不是对称的

matlab - 对称矩阵乘积不是对称的

问题描述

当我们有两个对称矩阵x和y时，矩阵z = xyx理论上也是对称的。然而，这在 Matlab 中并不完全正确：

x = randn(3);
y = randn(3);
x = x*x';
y = x*x';
z = x*y*x;
issymetric(z)

为什么会发生这种情况，我该怎么办？既然不想做

z = .5*(z+z')

这个stackoverflow问题的所有答案都不令人满意。

标签： matlabmatrixsymmetric

解决方案

问题是 Matlab 不能精确地执行计算，因此浮点计算的有限精度会引入一些截断误差。

如果您使用符号数学（没有截断，精确）运行示例，您会看到它z实际上是对称的。

x = sym(randn(3));
y = sym(randn(3));
x = x*x';
y = y*y';
z = x*y*x;

% issymmetric does not take symbolic expressions as argument, convert to double
issymmetric(double(z))

由于您现在在最后而不是在中间步骤中进行双重转换，因此矩阵将保持对称。

推荐阅读