exercises
Scripy
10.1
题目
源代码
import numpy as np import scipy.linalg as sl m = int(input("Please input m: ")) n = int(input("Please input n: ")) if n > m: t = m m = n n = t A = np.random.rand(m, n) b = np.random.rand(m, 1) A = np.mat(A) b = np.mat(b) x = sl.inv(A.T * A) * A.T * b print(x)
10.2
题目
源代码
import numpy as np import scipy.optimize as opt def func(nums): return (np.sin(nums - 2) ** 2) * (np.e ** (-nums ** 2)) func_inv = lambda x: -func(x) res1 = opt.brute(func_inv, ((-5, 5),), full_output=1) res2 = opt.fmin(func_inv, 0, full_output=1) max1, max2 = -res1[1], -res2[1] print(max1, max2)
10.3
题目
源代码
import numpy as np import scipy.spatial.distance as ssd import math m = int(input("Please input m: ")) n = int(input("Please input n: ")) X = np.random.rand(m, n) Y = ssd.pdist(X) print(ssd.squareform(Y))