所以目前我正在使用如下代码：
Q,R = np.linalg.qr(matrix)
Qb = np.dot(Q.T, new_mu[b][n])
x_qr = np.linalg.solve(R, Qb)
mu.append(x_qr)
只要我的矩阵是方阵，代码就可以正常工作，但一旦不是方阵，系统就无法求解，并且会出现错误。如果我理解正确的话，我不能linalg.solve在非满秩矩阵上使用，但是有没有办法让我在不使用解决方案的情况下克服这个障碍lstsquare？