使用matlab实现Jacobi迭代法与Gauss-Seidel迭代法求解线性方程组
jacobi.mfunction [y,n]=jacobi(A,b,x0,ep) %A:系数矩阵;b:常数矩阵;x0:迭代初值;ep:迭代精度;D=diag(diag(A)); %取A的主对角线上的元素建立对角矩阵;L=-tril(A,-1); %取A的主对角线以下(不包括主对角线)的元素建立下三角矩阵;U=-triu(A,1); %取A的主对角线以上(不包括主对角线)的元素建立上三角矩阵;B=D\(L+U);f=D\b;y=B*x0+f;n=1;while norm(y-x0)>
原创
2021-09-21 11:25:59 ·
5275 阅读 ·
1 评论