Scipy练习

import numpy as np
from scipy.optimize import leastsq

m = 20
n = 10
A = np.random.rand(m,n)
b = np.random.rand(m,1)


#Now find x = arg minx ||Ax - b||2

#S(b) = ||Xb-y||2  -->_b = arg min(S(b))
#XT * X * _b = XT * y
#_b = (XT * X)-1 * XT * y

a = np.dot(A.T,A)
a = np.linalg.inv(a)
a = np.dot(a,A.T)
x = np.dot(a,b)

print(x)

import numpy as np
import math
from scipy.optimize import minimize as f

x = np.linspace(0,2,500)

def func(x):
	y1 = np.power(np.sin(x-2),2)
	y2 = y1*np.power(math.e,-(x ** 2))
	return y2

answer = f(func, -1)

print(-func(answer.x))

import numpy as np 
from scipy import optimize
import scipy.spatial.distance as dis

X = np.random.randint(10,size = (8,10))
distance = dis.pdist(X,'euclidean')
print(dis.squareform(distance))