Scipy 习题

    求解超定方程的最小二乘解，求解公式如下：

                                             

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)