Scipy练习

from scipy.optimize import least_squares  
import numpy as np  
  
n = 2  
m = 3  
  
A = [[np.random.randint(1, 100) for i in range(m)] for j in range(n)]
print("A:")  
print(A)  
  
b = [np.random.randint(1, 100) for i in range(n)]
  
print("b:")  
print(b)  
  
x = [0, 0, 0]  
  
  
def func(x):  
    return np.linalg.norm(np.matmul(A, x)-np.array(b))  
  
  
solutions = least_squares(func, x)  
print("Solution:")  
print(solutions.x)  
print("Residual:")  
print(func(solutions.x))  

某次运行结果：

import math  
from scipy.optimize import fmin  
  
  
def func(x):  
    return -((math.sin(x-2))**2*math.e**(-x**2))  
  
x = fmin(func, 0)
print("max value:")
print(-func(x))

结果：

import numpy  
import scipy.spatial  
  
m = 3
n = 5  
X = numpy.random.randint(100, size = (m, n))  
print(scipy.spatial.distance.pdist(X))