Exercise 9.1: Matrix operations
Calculate A + A, AA>,A>A and AB. Write a function that computes A(B−λI) for any λ.
import numpy
def calc(A,B,x):
I=numpy.eye(m)
return numpy.dot(A,B-I*x)
n=200
m=500
A=numpy.random.normal(size=(n,m))
B=numpy.random.normal(size=(m,m))
for i in range(m):
for j in range(m):
if ((i<m-1) and (j<m-1)):
B[i+1][j+1]=B[i][j]
C=A+A
D=numpy.dot(A,A.T)
E=numpy.dot(A.T,A)
F=numpy.dot(A,B)
G=calc(A,B,1)
Exercise 9.2: Solving a linear system
Generate a vector b with m entries and solve Bx = b.
import numpy
n=200
m=500
A=numpy.random.normal(size=(n,m))
B=numpy.random.normal(size=(m,m))
for i in range(m):
for j in range(m):
if ((i<m-1) and (j<m-1)):
B[i+1][j+1]=