import numpy as np
from scipy.linalg import lstsq
m = 80
n = 40
A = np.random.normal(1, 2, (m, n))
b = np.random.normal(1, 2, (m, 1))
r_norm = lstsq(A, b)[1]
print(r_norm)
输出:
[203.82746375]
import numpy as np
from scipy.optimize import fmin
def func(x):
return (np.sin((x-2)*(np.exp(-x**2)))**2)*(-1)
xopt = fmin(func, 0)
print(xopt[0], -func(xopt[0]))
输出:
Optimization terminated successfully.
Current function value: -1.000000
Iterations: 21
Function evaluations: 42
0.2906875000000003 0.99999999978