Numpy库练习

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]=