from scipy.optimize import least_squares
import numpy as np
n = 2
m = 3
A = [[np.random.randint(1, 100) for i in range(m)] for j in range(n)]
print("A:")
print(A)
b = [np.random.randint(1, 100) for i in range(n)]
print("b:")
print(b)
x = [0, 0, 0]
def func(x):
return np.linalg.norm(np.matmul(A, x)-np.array(b))
solutions = least_squares(func, x)
print("Solution:")
print(solutions.x)
print("Residual:")
print(func(solutions.x))
某次运行结果:
import math
from scipy.optimize import fmin
def func(x):
return -((math.sin(x-2))**2*math.e**(-x**2))
x = fmin(func, 0)
print("max value:")
print(-func(x))
结果:
import numpy
import scipy.spatial
m = 3
n = 5
X = numpy.random.randint(100, size = (m, n))
print(scipy.spatial.distance.pdist(X))