R求解多个线性系统,其中每个RHS是矩阵的一列
R solving multiple linear systems where each RHS is a column of a matrix
我知道我可以使用 R 中的 backsolve
求解系统 Lx = b
其中 L
是上三角。但是,现在我有一个矩阵 B
其中每个B
的列是 Lx = b
的右侧。基本上我想解决 Lx = b
,其中 b
是矩阵 B
的每一列。有没有比循环更有效的方法?
我没有发布 MWE,因为我很确定应该有一个函数可以做到这一点。
或许你可以试试下面的方法
x <- ginv(L)%*%b
我知道我可以使用 R 中的 backsolve
求解系统 Lx = b
其中 L
是上三角。但是,现在我有一个矩阵 B
其中每个B
的列是 Lx = b
的右侧。基本上我想解决 Lx = b
,其中 b
是矩阵 B
的每一列。有没有比循环更有效的方法?
我没有发布 MWE,因为我很确定应该有一个函数可以做到这一点。
或许你可以试试下面的方法
x <- ginv(L)%*%b