Python学习——numpy作业

Numpy Exercises

Generate matrices A, with random Gaussian entries, B, a Toeplitz matrix, where A 2 Rn×m and B 2 Rm×m,for n = 200, m = 500.

import numpy as np
from scipy.linalg import toeplitz
n, m = 200, 500
A = np.random.randn(n, m)
b = [i+0.0 for i in range(1,m+1)]
B = toeplitz(b, b)

Exercise 9.1: Matrix operations

Calculate A+A, A*A.t ; A.T*A and A*B. Write a function that computes A(B - λI) for any λ.

print(A+A)
print(np.dot(A,A.T))
print(np.dot(A.T,A))
print(np.dot(A,B))

def func(para):
    I = np.identity(m)
    return np.dot(A, B - para * I)

Exercise 9.2: Solving a linear system

Generate a vector b with m entries and solve Bx = b.

分析：可以直接用numpy.linalg.solve，也可以用B的逆矩阵乘以b，两者等价

b = np.random.random(m)
print(np.linalg.solv