【科学与工程计算】 实验6 SciPy

【科学与工程计算】

实验6 SciPy

题目1 随机产生一个5行3列的矩阵(元素值的范围1-100)，使用scipy计算行与行的距离

import scipy.spatial.distance 
import numpy as np 
A = np.random.randint(0,101,[5,3]) #随机生成矩阵各个元素
#print(A)
X=scipy.spatial.distance.pdist(A) #计算行与行之间的距离
print(X)

题目2 使用scipy计算如下函数的最大值

import scipy.optimize 
import numpy as np 
def f(x): 
    return -np.sin(x-2)**2 * np.exp(-x*x)
#scipy提供求最小值所在位置的函数，通过计算相反数来计算最大值 
max_x=scipy.optimize.fmin( f ,0) 
print( -f(max_x) )