python - numpy.dot 的倒数

Question

我可以很容易地计算出如下内容:

R = numpy.column_stack([A,np.ones(len(A))]) 
M = numpy.dot(R,[k,m0])

其中 A 是一个简单数组，k,m0 是已知值。

我想要一些不同的东西。确定 R、M 和 k 后，我需要获得 m0。有没有办法通过函数 numpy.dot() 的反函数来计算这个？或者只有重新排列矩阵才有可能？

Answer 1

M = numpy.dot(R,[k,m0])

正在执行矩阵乘法。 M = R * x。

所以要计算倒数，您可以使用 np.linalg.lstsq(R, M) :

import numpy as np
A = np.random.random(5)
R = np.column_stack([A,np.ones(len(A))]) 
k = np.random.random()
m0 = np.random.random()
M = R.dot([k,m0])

(k_inferred, m0_inferred), residuals, rank, s = np.linalg.lstsq(R, M)

assert np.allclose(m0, m0_inferred)
assert np.allclose(k, k_inferred)

请注意， k 和 m0 是确定的，给定 M 和 R(假设 len(M) >= 2)。