Scipy作业

10.1：

 

import numpy
import scipy.linalg
n = 3
m = 3
 
A = numpy.random.rand(n, m)
b = numpy.random.rand(m, 1)
x, res , rk, s = scipy.linalg.lstsq(A, b)
 
print(b)
print(A.dot(x))

 

10.2：

import scipy.optimize
import numpy
 
def f(x):
    return -(numpy.sin(x - 2) ** 2) * numpy.e ** (-x ** 2)
 
mint = scipy.optimize.fmin(f, 0)
print(mint, -f(mint))

 

10.3：

import numpy
import scipy.spatial.distance
 
n = 3
m = 3
X = numpy.random.rand(n, m)
Y = scipy.spatial.distance_matrix(X, X)
print(X, '\n')
print(Y)