numpy作业

9.1，Matrix operations

都是一些基本操作

n = 200
m = 500
A = numpy.array([numpy.random.normal(0, 1, m) for _ in range(n)])
B = scipy.linalg.toeplitz([1., -.25, *numpy.zeros(m-2)])

# exercise 9.1
print('\033[1;35m A+A: \033[0m')
print(A+A)
print('\033[1;35m AA^T: \033[0m')
print(A.dot(A.T))
print('\033[1;35m A^TA: \033[0m')
print(A.T.dot(A))
print('\033[1;35m AB: \033[0m')
print(A.dot(B))
def compute(Lambda, A, B):
    return A.dot(B-Lambda*numpy.eye(m, m))
print('\033[1;35m A(B-λI): \033[0m')
print(compute(1, A, B))

9.2，Solving a linear system