jacob迭代法 gauss-seidel迭代法 matlab命令

1. jacobi迭代法解方程组的MATLAB程序

function x=jacobi(a,b,x0,m)%a:系数矩阵 b:等号右边矩阵 x0:迭代初值 m:迭代次数 x_temp=x0; for k=1:m

for i=1:length(b) sum=0; for j=1:length(b) if(j~=i)

sum=sum+a(i,j)*x_temp(j); end end

x(i)=-(sum-b(i))/a(i,i); end x_temp=x; end 2.

gauss-seidel迭代法解方程组的MATLAB程序

function x=gaussseidel(a,b,x0,m) x=x0; for k=1:m

for i=1:length(b) sum=0; for j=1:length(b) if(j~=i)

sum=sum+a(i,j)*x(j); end end

x(i)=-(sum-b(i))/a(i,i); end

x' end m=inv(a)*d';