Generate matrices A, with random Gaussian entries, B, a Toeplitz matrix, where A ∈ R(n×m) and B ∈ R(m×m), for n = 200, m = 500.
import numpy
import scipy.linalg
n = 200
m = 500
Mu = 0
Sigma = 1
A = numpy.random.normal(Mu, Sigma, (n, m))
print(A)
base = list(range(m))
B = scipy.linalg.toeplitz(base)
print(B)
Exercise 9.1: Matrix operations
Calculate A + A, A(A^T), (A^T)A and AB. Write a function that computes A(B − λI) for any λ.
#计算A + A
result = A + A
print('A + A =\n', result)
#计算A * (A ^ T)
result = numpy.dot(A, A.T)
print('A * (A ^ T) =\n', result)
#计算(A ^ T) * A
result = numpy.dot(A.T, A)
print('(A ^ T) * A =\n', result)
#计算A * B
result = numpy.dot(A, B)
print('A * B =\n', result)
#计算A * (B - λI)
Lambda = int(input(