import numpy as np
import scipy.linalg as lig
m = 10
n = 5
A = np.random.random((m, n))
b = np.random.random(m)
x = lig.lstsq(A, b)[0]
norm = np.linalg.norm(np.dot(A, x)-b, ord=2)
print(norm)
Result
0.811680468476
Exercise 10.2: Optimization
Solution
from scipy.optimize import fmin
import numpy as np
import matplotlib.pyplot as plt
deff(x):return -1 * np.sin(x-2)**2 * np.e**(-1*x**2)
opt = fmin(f, 0)
X = np.linspace(-2, 2, 1000)
plt.plot(X, -f(X))
print('maxmium:', -f(opt)[0])
plt.show()
import numpy as np
import scipy.linalg as lig
import scipy.spatial.distance as dis
m = 4
n = 5
X = np.mat(np.random.rand(m, n))
print(X)
cities = dis.pdist(X)
ans = dis.squareform(cities)
print(ans)