求解 R 中的冗余线性方程组

Question

我有一个不可逆矩阵A 和一个向量b，我相信有一个解决方案x 到Ax = b。我想找到一个这样的例子x。当我在 R 中尝试 solve(A,b) 时，它会产生错误，因为 A 是单数。有什么办法可以让 R 给我一个随机的解决方案？

Answer 1

我最终尝试了lm(b ~ 0 + A)，它有效。它会将某些列的估计值保留为 NA，您可以将其替换为 0 以获得示例解决方案。例如

A = matrix(c(1,1,0,0),nr=2,byrow=F)

b = c(2,2)

lm(b ~ A)

将为A 的两列生成系数 2 和 NA。 solve(A,b)、lsfit(A,b) 和 qr.solve(A,b) 不起作用。

编辑：MASS::ginv(A) %*% b 也可以。