没有合适的资源?快使用搜索试试~ 我知道了~
线性方程数值解中jacobi迭代matlab实现
5星 · 超过95%的资源 需积分: 10 11 下载量 81 浏览量
2009-04-28
10:33:31
上传
评论
收藏 506B TXT 举报
温馨提示
试读
2页
线性方程数值解中jacobi迭代matlab实现
资源推荐
资源详情
资源评论
%jacobi迭代
%author:野渡无人
%create date:2007-11-23
function res=jacobi(a,b,x0,eps)
%a为系数矩阵,b为方程组常数向量,x0为初值
if nargin==3
eps=1.0e-6;
elseif nargin<3
error
return;
end
D=diag(diag(a));%求矩阵对角阵
D=inv(D);
L=tril(a,-1);%求严格下三角矩阵
U=triu(a,1);%求严格上三角矩阵
B=-D*(L+U);
if norm(B,inf)>=1
fprintf('B范数>=1.');
return;
end
f=D*b;
res=B*x0+f;
i=0;
while norm(res-x0)>=eps
x0=res;
res=B*x0+f;
i=i+1;
%author:野渡无人
%create date:2007-11-23
function res=jacobi(a,b,x0,eps)
%a为系数矩阵,b为方程组常数向量,x0为初值
if nargin==3
eps=1.0e-6;
elseif nargin<3
error
return;
end
D=diag(diag(a));%求矩阵对角阵
D=inv(D);
L=tril(a,-1);%求严格下三角矩阵
U=triu(a,1);%求严格上三角矩阵
B=-D*(L+U);
if norm(B,inf)>=1
fprintf('B范数>=1.');
return;
end
f=D*b;
res=B*x0+f;
i=0;
while norm(res-x0)>=eps
x0=res;
res=B*x0+f;
i=i+1;
资源评论
- SIMIRMORNING2012-11-15不错,真实好东西,看了过后顿时明白了··············好资源,你值得拥有。
wanjin0
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功