MATLAB 块雅克比迭代法求线性方程组Ax=b的解 块高斯-赛德尔迭代法求线性方程组Ax=b的解.docVIP

function [x,N]= BJ(A,b,x0,d,eps,M) %块雅克比迭代法求线性方程组Ax=b的解 if nargin==4??? eps= 1.0e-6;??? M = 10000;elseif nargin4??? error??? returnelseif nargin ==5??? M = 10000;???????????? %参数的默认值end NS = size(A);n = NS(1,1);??????????????? if(sum(d) ~= n)??? disp(分块错误！);??? return;endbnum = length(d);bs = ones(bnum,1);for i=1:(bnum-1)??? bs(i+1,1)=sum(d(1:i))+1; ??? %获得对角线上每个分块矩阵元素索引的起始值end DB = zeros(n,n);??????????????? for i=1:bnum??? endb = bs(i,1)+d(i,1)-1;??? DB(bs(i,1):endb,bs(i,1):endb)=A(bs(i,1):endb,bs(i,1):endb);??? %求A的对角分块矩阵endfor i=1:bnum??? endb = bs(i,1)+d(i,1)-1;??? invDB(bs(i,1):endb,bs(i,1):endb)=inv(DB(bs(i,1):endb,bs(i,1):endb));??? %求A的对角分块矩阵的逆矩阵end N = 0;tol = 1;while tol=eps ??? x = invDB*(DB-A)*x0+invDB*b;??? %由于LB+DB=DB-A??? N = N+1;?????????????????????????????? %迭代步数??? tol = norm(x-x0);?????????????????? %前后两步迭代结果的误差??? x0 = x;??? if(N=M)??????? disp(Warning: 迭代次数太多，可能不收敛！);??????? return;??? endend function [x,N]= BGS(A,b,x0,d,eps,M) %块高斯-赛德尔迭代法求线性方程组Ax=b的解 if nargin==4??? eps= 1.0e-6;??? M = 10000;elseif nargin4??? error??? returnelseif nargin ==5??? M = 10000;end NS = size(A);n = NS(1,1);bnum = length(d);bs = ones(bnum,1);for i=1:(bnum-1)??? bs(i+1,1)=sum(d(1:i))+1;??? %获得对角线上每个分块矩阵元素索引的起始值end DB = zeros(n,n);?????? for i=1:bnum??? endb = bs(i,1)+d(i,1)-1;??? DB(bs(i,1):endb,bs(i,1):endb)=A(bs(i,1):endb,bs(i,1):endb);??? %求A的对角分块矩阵endLB = -tril(A-DB);????? %求A的下三角分块阵UB = -triu(A-DB);????? %求A的上三角分块阵 N = 0;tol = 1;while tol=eps ??? invDL = inv(DB-LB);??? x = invDL*UB*x0+invDL*b;??? %块迭代公式??? N = N+1;??? tol = norm(x-x0);??? x0 = x;??? if(N=M)??????? disp(Warning: 迭代次数太多，可能不收敛！);??????? return;??? endend 类别：matlab 编程 | | 添加到搜藏 | 分享到i贴吧 | 浏览(168) | 评 论?(0) ? 上一篇：MATLAB 共轭梯度法求线性方程组A...????