python测试代码:
import numpy as np
'方阵A'
A = np.array([[1,2],[3,4]])
A
array([[1, 2],
[3, 4]])
'逆矩阵...非奇异矩阵python测试 :
import numpy as np
'方阵A'
A = np.array([[1,2],[3,4]])
A
array([[1, 2],
[3, 4]]...'病态矩阵测试 Ax = b'
'x = la.inv(A).dot(b) 如下python代码所示:'
A = np.array([[400,-201],[-800,401]])
b = np.array...([200,-200])
x = la.inv(A).dot(b)
x
array([-100., -200.])...条件数的定义:K(A)= ‖inv(A)‖ * ‖A‖ 的大小。
其中,‖‖ 表示对矩阵取某一种范数。
接下来测试上面提到的病态矩阵的条件数,和一个良好的矩阵的条件数,看看它们的大小。