Scipy学习

Exercise 10.1: Least squares

代码：

import scipy.linalg as lina
import numpy as np

m = 10
n = 6
A = np.random.random_sample((m, n))
x = np.random.random_sample((n,))
b = A.dot(x) + np.random.normal(0, 1, (m,))
x_ = lina.inv(A.T.dot(A)).dot(A.T).dot(b)
print(lina.norm(A.dot(x_) - b))

结果：

2.7711555872938862
[Finished in 3.8s]

Exercise 10.2: Optimization

代码：

import scipy.optimize as so
from math import sin, exp

def f(x):