第十三周作业

1.

代码如下：

import numpy as np  
import scipy.linalg  as sla  
  
m = 10 
n = 5  
A = np.random.random((m, n))  
b = np.random.random((m, 1))  
  
x = sla.lstsq(A, b)[0]  
print(np.linalg.norm(np.dot(A, x)-b))  

运行结果：

2.

代码如下：

import numpy as np  
import scipy.optimize as opt  
  
def fun(x):  
    return -((np.sin(x-2)**2) * np.exp(-(x**2)))  
  
max_num = opt.minimize(fun,0)  
print(-max_num.fun)  

运行结果：

3.

代码如下：

import numpy as np  
import scipy.spatial.distance as dist  
  
m = 4 
n = 4 
x = np.random.normal(0, 16, (m, n))  
distance = dist.pdist(x)  
print(dist.squareform(distance))  

运行结果：