Question

我一直试图解决两个 m -by- n 矩阵 A 和 B 是相当于通过

[a,ja] = rref(A,tol)
[b,jb] = rref(B,tol)

然后比较

isequal(a,b) & isequal(ja,jb)

首先，我真的不明白ja和jb是什么。我的问题是，对于 A 和 B ，减少行的梯形形式非常简单，并且在所有情况下都相同。我不知道这是否有意。例如，我得到等同于

A = rand(40,3)
B = rand(40,3)

我不确定是否是胡说八道。

Answer 1

看起来你正试图检查两个矩阵的reduced row echelon forms是否与元素相等。从您定义A和B的方式来看，它们是（它实际上是一个超定系统，我相信）。但是，我认为您可能已经翻转了行和列。相反，如果您创建A和B，则列数多于行数（即未确定系统的augmented matrix）：

A = rand(3,40)
B = rand(3,40)

然后当你运行rref时，你会看到一个非常不同的输出，你的比较会像你预期的那样返回false。

另外，我认为使用以下内容就足够了，因为两个元素相等的矩阵肯定会共享相同的等级（或其近似值）：

a = rref(A,tol);
b = rref(B,tol);
isequal(a,b)

减少行梯形矩阵等价

1 个答案: