本文共 654 字,大约阅读时间需要 2 分钟。
求解超定方程的最小二乘解,求解公式如下:
import scipy.linalgimport numpy as npm=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.optimizeimport numpy as npdef 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.distanceimport numpy as npm=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)
转载地址:http://sywzi.baihongyu.com/