Python scipy 部分习题

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