import numpy as np
A = np.random.randint(-20, 20, size = 30 * 20).reshape(30, 20)
b = np.random.randint(-20, 20, size = 30 * 1).reshape(30, 1)
x = np.dot(np.dot(np.linalg.inv(np.dot(A.T, A)), A.T), b)
print(x)
print ("The norm of the residual:", np.linalg.norm(x))
结果:
import numpy as np
from scipy import optimize
def func(x):
return -np.power(np.sin(x - 2), 2) * np.exp(-np.power(x, 2))
max = -optimize.minimize_scalar(func).fun
print(max)
结果:
import numpy as np
from scipy import spatial
X = np.random.randint(-20, 20, size = 10 * 10).reshape(10, 10)
print(X)
print(spatial.distance.cdist(X, X))
结果: