专业和 班级 课程名 称 实 验 目 的 和 要 求 实 验 内 容 和 步 骤 福州大学数学与计算机科学学院 计算机上机实验报告 姓名 学号 成绩 数值计算方法 实验名称 雅可比迭代 高斯迭代法 实验目的：分别用雅可比迭代法、高斯-赛德尔迭代法及逐次超松弛 迭代法解此方程组 雅可比迭代法 m 文件： function x=jacobi(A,b,P,delta,n) N=length(b); for k=1:n for j=1:N x(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j); end err=abs(norm(x'-P)); P=x'; if(err

b=[-12,20,3]'; P=[0,0,0]'; x=jacobi(A,b,P,1e-4,20) 结果： （其中 B 为迭代矩阵） 

高斯-塞德尔迭代法 m 文件： function x=Seidel(A,b,P,delta,n) N=length(b); for k=1:n x(1)=(b(1)-A(1,2:N)*P(2:N))/A(1,1); x(N)=(b(N)-A(N,1:N-1)*(x(1:N-1))')/A(N,N); x(j)=(b(j)-A(j,1:j-1)*x(1:j-1)-A(j,j+1:N)*P(j+1:N))/A(j,j); for j=1:N if j==1 elseif j==N else end end err=abs(norm(x'-P)); P=x'; if(err

（其中 B 为迭代矩阵） 研 究 与 探 讨