#E9.1
import numpy as np
from scipy.linalg import toeplitz
n = 200
m=500
A = np.matrix(np.round(np.random.normal(4,6,n*m),0).reshape(n,m),'int',True)
aList = []
for i in range(1,m+1):
aList.append(i)
B = np.matrix(toeplitz(aList,aList).reshape(m,m),'int',True)
print(A)
print(B)
print(A+A)
print(A.dot(A.T))
print(A.T.dot(A))
print(A.dot(B))
def computes(A,B,lamda):
I = np.eye(m)
temp = B - lamda*I
result = A.dot(temp)
computes(A,B,10)
#E9.2
import numpy as np
from scipy.linalg import toeplitz
m=500
aList = []
for i in range(1,m+1):
aList.append(i)
B = np.matrix(toeplitz(aList,aList).reshape(m,m),'int',True)
b_list = []
for i in range(1,m+1):
b_list.append(i)
b = np.array(b_list)
A_abs = np.linalg.det(B)
B_r = np.linalg.inv(B)
result = B_r.dot(b)
print(B)
print(B_r)
print(result)
#E9.3
import numpy as np
from scipy.linalg import toeplitz
n=200
m=500
aList = []
for i in range(1,m+1):
aList.append(i)
A = np.matrix(np.round(np.random.normal(4,6,n*m),0).reshape(n,m),'int',True)
B = np.matrix(toeplitz(aList,aList).reshape(m,m),'int',True)
B_inf = np.linalg.norm(B, ord=np.inf)
A_fro = np.linalg.norm(A,ord='fro')
print(A_fro)
print(B_inf)# infinity norm
#E9.4
import numpy as np
from scipy.linalg import toeplitz
import time
t0 = time.clock()
n=200
A = np.matrix(np.round(np.random.normal(4,6,n*n),0).reshape(n,n),'int',True)
eigenvalues, eigenvectors = np.linalg.eig(A)
print(str(time.clock()-t0)+'s')
print(eigenvalues)
print(eigenvectors)
pos = []
max = -9999999
t0 = time.clock()
for x in range(0,n):
if eigenvalues[x] > max:
if len(pos) != 0:
pos.clear()
pos.append(x)
elif eigenvalues[x] == max:
pos.append(x)
else:
pass
print(str(time.clock()-t0)+'s')
#n个迭代就可以了
for p in pos:
print("展示最大特征值的位置:")
print(str(p)+' : ')
print("展示最大特征值:")
print(str(eigenvalues[p]))
print("展示最大特征向量:")
print(str(eigenvectors[x]))
#E9.5
import numpy as np
import scipy
n = 200
p = 0.4
C = np.random.binomial(1, p, n*n).reshape(n,n)
print(C.size)
U,S,Vh = np.linalg.svd(C)
print(max(S))
#relationship:n*p 约等于 max(奇异值)
#E9.6
import numpy as np
import math
def func(z,A):
#closest = 0
A_matrix = np.matrix(A,'int',True)
z_list = [z for i in range(0,len(A))]
Z_matrix = np.matrix(z_list,'int',True)
pos = np.argmin(abs(A_matrix-Z_matrix))
return A[pos]
A = [1,2,3,4,5]
z = -9
print(func(z,A))