求解超定方程的最小二乘解,求解公式如下:
import scipy.linalg
import numpy as np
m=np.random.randint(20,30)
n=np.random.randint(10,20)
A=np.random.normal(0,1,[m,n])
b=np.random.randn(m)
#直接根据公式求解即可
x=scipy.linalg.inv(A.T.dot(A)) .dot(A.t) .dot(b)
求解方程的最大值,scipy 提供的函数可求解函数值最小时自变量的值,故求函数相反数的最小值即可得到最大值。
import scipy.optimize
import numpy as np
def f_inv(x):
return -np.sin(x-2)**2 * np.exp(-x*x)
#scipy提供求最小值所在位置的函数,通过计算相反数来计算最大值
max_x=scipy.optimize.fmin( f_inv ,0)
print( -f_inv(max_x) )
求解各行代表坐标间的距离,使用scipy中的向量距离计算函数即可。
import scipy.spatial.distance
import numpy as np
m=np.random.randint(20,30)
n=np.random.randint(10,20)
A=np.random.normal(0,1,[m,n])
X=scipy.spatial.distance.pdist(A)