《雅可比迭代法的MATLAB程序》由会员分享，可在线阅读，更多相关《雅可比迭代法的MATLAB程序（2页珍藏版）》请在金锄头文库上搜索。
1、雅可比迭代法的MATLAB程序：Functionx，k,index=Jacobi（A，b,ep,itmax）% 求线性方程组的雅可比法; A 为方程组的系数矩阵； b 为方程组的右端项； x 为方程组的解；% ep 为精度要求，缺省值为 le-5；% it_max为最大迭代次数，缺省值为100；% k 为迭代次数;% index为指标变量,index=0表示计算失败，index=1表示计算成功； if nargin4it_max=100;endif nargin3ep=le5;endn=length（A）;k=0； x=zeros（n， 1）;y=zeros（n,1）;index=1; whi
2、le k=it_maxfor i=1： nif abs （A（i， i）le10index=0；return;endy（i）=（b（i）-A（i,1:n）x（1： n）+A（i， i）x（i）/A（i， i）；endif norm（y-x， inf）ep break;endk=k+1;x=y;end高斯一赛德尔迭代的MATLAB程序Functionx， k,index=Gau-seidel（A,b， ep， itmax）% 求线性方程组的高斯-赛德尔迭代法；% A 为方程组的系数矩阵;% b 为方程组的右端项；% x 为方程组的解;% ep 为精度要求,缺省值为 le5;% it_max 为最
3、大迭代次数,缺省值为 100;% k 为迭代次数;% index 为指标变量， index=0 表示计算失败， index=1 表示计算成功; if nargin4it_max=100；endif nargin3ep=le-5;end n=length（A）;k=0；x=zeros（n,1）；y=zeros（n，1）;index=1; while k=it_maxfor i=1:nif abs （A（i，i）le-10index=0； return；Endif i=1 y（i）=（b（i）A（i,i+1：n）x（i+1:n）/A（i,i）; elseif i=n y（i）=（b（i）A（i,1：i-1）y（1:i1）/A（i，i）; elsey（i）=（b（i）A（i,1：i1）y（1：i-1）A（i，i+1:n）*x（i+1:n）/A（i,i） end endif norm（y-x，inf）epbreak;endk=k+1;x=y；end
