方程组如下:
我们使用赛得尔迭代法求解改方程组,运行的误差在0.00001。代码如下:
import numpy as np
ef g(x,x0,n,d):#赛德尔迭代法
x1=np.ones(n);y=np.ones(n);y[:-1]=-1;qd=0;x2=np.ones(n);x3=np.ones(n);x3[:]=x0[:]
while np.any(abs(x1-x3)>d):
qd+=1
for u in range(n-1):
x1[u]=x3[u]
for j in range(n-1):
tq=0
for m in range(n):
if m!=j:
tq+=(x[j,:]/x[j,j]*y)[m]*x3[m]
x3[j]=tq;x2[j]=tq
for re in range(n-1):
x3[re]=x2[re]
return x3,qd
运行结果:
print('赛德尔迭代法,最终的解X_i(角标从小到大):{},及迭代次数K:{}'.format(g(x,x0,n,d)[0][:-1],g(x,x0,n,d)[1]))
赛德尔迭代法,最终的解X_i(角标从小到大):[1.0000023 1.99999984 2.99998957 0.99998176],及迭代次数K:28