高斯赛德尔迭代法matlab编程(共2页).doc
《高斯赛德尔迭代法matlab编程(共2页).doc》由会员分享,可在线阅读,更多相关《高斯赛德尔迭代法matlab编程(共2页).doc(2页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
精选优质文档-倾情为你奉上functionx,k=GaussSeidel(A,b,x0,eps,M)%高斯赛德尔迭代法求方程组的解(矩阵公式求解)%A为方程组的系数矩阵;b为方程组的右端项%x为线性方程组的解了;x0为迭代初值%eps为误差限;M为迭代的最大次数if nargin=3 eps= 1.0e-6;%默认精度 M = 10000;%参数不足时默认后两个条件elseif nargin =4 M = 10000;%参数的默认值elseif nargin=1 error(迭代矩阵谱半径大于1迭代法不收敛); return;endk=0;tol=1;while tol=eps x = B*x0+g; k = k+1; %迭代步数 tol = norm(x-x0);%前后两步迭代结果的误差 x0 = x; if(k=M) disp(Warning: 迭代次数太多,可能不收敛!); return; endend专心-专注-专业
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高斯赛德尔 迭代法 matlab 编程
限制150内