matlab - 使用ode45后如何在matlab中查找矩阵错误？

Question

目前，我已将 ode45 与一个输出 3945x9 的矩阵 (C) 的函数一起使用。这些应该是 3x3 矩阵，所以我做了 C = reshape(C.',3,3,[])。现在，我有一个 3x3x3945 矩阵。我想要做的是找到每个 3x3 矩阵的误差。这是使用 C*C 完成的。- 眼睛（3）。但是，我现在不知道如何处理我的所有数据。它适用于 C(:,:,1)*C(:,:,1)。- 眼睛（3）。但是，不是 C(:,:,:)*C(:,:,:) - eye(3)。

标签： matlab

Answer 1

使用 nD 矩阵乘法例程 pagemtimes。例如，

pagemtimes(C,'none',C,'transpose') - eye(3)