收敛性:
对应的矩阵J=A-D谱小于1
import numpy as np print('start') N=2 M=20 #A=np.array([[2,1],[1,3]]) D=np.array([[2,0],[0,3]]) F=np.array([[0,-1],[-1,0]]) X=np.array([1,1]) B=np.array([10,10]) B=np.matmul(np.linalg.inv(D),B) for i in range(M): X=np.matmul(np.linalg.inv(D),np.matmul(F,X)) X=X+B print(X)