计算机数值分析
线性方程组的迭代法:雅可比迭代公式:
书上的流程图如下:
这里直接把输入输出弄到了文件之中
具体代码如下:
里面使用了python的numpy库
#coding=gbk
import numpy as np;
def Jacobi(n,minner,datas,x):
count=1;
numVarity=datas.shape[0];
cols=datas.shape[1];
while count<=n: #第几次迭代
x_new=np.zeros(data.shape[0]);
for j in range(0,numVarity): #变量的个数 第J个变量 也就是第j行
for k in range(0,numVarity):
if k==j:
continue;
x_new[j]+=x[k]*data[j][k];
x_new[j]=(data[j][cols-1]-x_new[j])/data[j][j]
if np.max(x_new-x)<minner:
x=x_new.copy();
return x;
else:
count+=1;
x=x_new.copy();
return x;
data=np.loadtxt("Jacobi.txt",comments="#",delimiter=",");
limit=np.loadtxt("limits.txt",comments="#",delimiter=",");
N=int(limit[0])
Minner=float(limit[1]);
X=np.zeros(data.shape[0])
X=Jacobi(N,Minner,data,X);
print(X)
代码运行的结果如下:
书上的结果如下:
经比较:和书上的结果相同,所以是正确的。
遇事不决,可问春风。