当你有一个Ax = b
形式的方程时,np.solve()
工作得很好 我的问题是我实际上有一个xC = D
形式的方程,其中 x 是我想要找出的 2x2 矩阵,C 和 D 是我给出的 2x2 矩阵。
而且因为矩阵乘法通常不是可交换的,所以我不能只是交换两者。
有没有一种有效的方法可以在numpy(或python中的其他库(中解决这个问题?
x @ C = D
与D^-1 @ x @ C @ C^-1 = D^-1 @ D @ C^-1
相同,后者D^-1 @ x = C^-1
的形式为Ax = b,其中 A 为np.linalg.pinv(D)
,b 为np.linalg.pinv(C)
归结为
x = D @ np.linalg.pinv(C)
你可以通过将等式的两边乘以 C 的倒数来得到