Scipy部分习题

题目10.1：

代码：

import numpy as np  
from scipy.linalg import * 
m=10
n=5
A = np.random.randint(0,30,size=(m,n))  
b = np.random.randint(0,20,size=(m,1))  
x, res, rnk, s = lstsq(A, b)  
b1 = norm(b-np.dot(A,x),ord=2)    
print("A:\n",A)  
print("b:\n",b)  
print("b1=",b1)   

输出：

A:
 [[ 7 19 16  1  5]
 [20 20 29 22 13]
 [19 19 19 14 28]
 [ 8  1 11 21 11]
 [20 10  4 13 11]
 [10 12 15 29  9]
 [ 4 28 16 28 14]
 [ 4 14 21  9  6]
 [17  7 13 14 21]
 [15 11  0  9 17]]
b:
 [[ 2]
 [ 5]
 [ 7]
 [ 6]
 [19]
 [ 6]
 [ 7]
 [12]
 [ 5]
 [16]]
b1= 15.309997441988713

题目10.2：

代码：

import numpy as np  
from scipy.optimize import *
import math  
def fx(x):  
    return(-(math.sin(x-2)**2)*math.exp(-(x**2)))  
n=minimize_scalar(fx)
print(-n.fun)  

输出：

0.9116854118471548

题目10.3：

代码：

import numpy as np  
from scipy.spatial.distance import * 
import math  
m=2
n=2
X = np.random.rand(m, n)  
print(X)  
print(cdist(X,X,metric="euclidean"))  

输出：

[[0.52627919 0.99969535]
 [0.82253715 0.25840462]]
[[0.         0.79829864]
 [0.79829864 0.        ]]