Python-scipy第十三周作业

Exercise 10.1

import numpy as np

m = 10
n = 5
A = np.mat(np.random.randn(m,n))
b = np.mat(np.random.randn(m,1))
x = np.array(np.linalg.lstsq(A,b,rcond=-1)[0])
res = np.array(np.linalg.lstsq(A,b,rcond=-1)[1])
nor = np.linalg.norm(res)

各参数为：

残差值：

norm值：

Exercise 10.2:

from scipy import optimize
f = lambda x:-1 * (np.power(np.sin(x-2),2)) * (np.exp((-1)*x*x))
temp = optimize.minimize_scalar(f)
result = -1 * temp.fun

Exercise 10.3:

from scipy.spatial.distance import pdist 
m = 5
n = 2
city = ["A", "B", "C", "D", "E"]
X = np.mat(np.random.randn(m,n))
Y = pdist(X, 'euclidean')
count = 0
for i in range(m):
	for j in range(i+1, m):
		print("city " + city[i] + " to city " + city[j] +": " + str(Y[count]))
		count += 1