第十三周 Scipy练习

#10.1_Least squares 
import numpy as np
import scipy.linalg
m = 30
n = 20

A = np.random.rand(m,n)
b = np.random.rand(m)
x, residues , rank, s = scipy.linalg.lstsq(A,b)

print(x)
print(residues)

#10.2_ Optimization 
import math
import scipy.optimize
def func(x):
    return -(math.sin(x-2) * math.sin(x-2) * math.exp(-x**2))

maxi = scipy.optimize.fmin(func, 1)
print(maxi)

#10.3_Pairwise distances 
import scipy.spatial.distance as dis
m = 12
n = 8
X = np.random.rand(m,n)
Y = dis.pdist(X, 'euclidean')
print(Y)