【python】第二十一次课作业

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))