Jacobi迭代法 Gauss-Seidel迭代法【DOC精选】.docVIP

Matlab线性方程组的迭代解法（Jacobi迭代法 Gauss-Seidel迭代法）实验报告 2008年11月09日 星期日 12:49 1.熟悉Jacobi迭代法，并编写Matlab程序matlab程序 按照算法(Jacobi迭代法)编写Matlab程序(Jacobi.m) function [x, k, index]=Jacobi(A, b, ep, it_max) %求解线性方程组的Jacobi迭代法，其中 % A ---方程组的系数矩阵 % b ---方程组的右端项 % ep ---精度要求。省缺为1e-5 % it_max ---最大迭代次数，省缺为100 % x ---方程组的解 % k ---迭代次数 % index --- index=1表示迭代收敛到指定要求； % index=0表示迭代失败 if nargin 4 it_max=100; end if nargin 3 ep=1e-5; end n=length(A); k=0; x=zeros(n,1); y=zeros(n,1); index=1; while 1 for i=1:n y(i)=b(i); for j=1:n if j~=i y(i)=y(i)-A(