求解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也适用。