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