function [num,x] = Jacobi(A,b,X,Nmax,tol)
%用Jacobi迭代法求解方程组Ax = b
%输入:A为方程组的系数矩阵,b为方程组右端的列向量,X为迭代初值构成的列向量,nm为最大迭代次数,tol为误差精度
%输出:x为求得的方程组的解构成的列向量,Nmax为迭代次数
num = 1;
D = diag(diag(A)); %令A = D-L-U,计算矩阵D
L = tril(-A)+D; %令A = D-L-U,计算矩阵L
U = triu(-A)+D; %令A = D-L-U,计算矩阵U
B = inv(D)*(L+U);
f = inv(D)*b;
%下面是迭代过程
while num <= Nmax
x = B*X + f; %用迭代格式进行迭代
if norm(x-X,'inf')<tol
disp('收敛!!!');
disp('迭代次数为');num
disp('方程组的解为');x
return;
%上面:达到精度要求就结束程序,输出迭代次数和方程组的解
end
X = x;
num = num+1;
end
%下面:如果达到最大迭代次数仍不收敛,输出警告语句及迭代的最终结果(并不是方程组的解)
disp('在最大迭代次数内不收敛!!!');
disp('最大迭代次数后的结果为');x
Jacobi迭代法解线性方程matlAB源码
需积分: 34 160 浏览量
2015-12-07
22:10:53
上传
评论
收藏 1KB ZIP 举报
一杯暖阳
- 粉丝: 0
- 资源: 8
最新资源
- NetOps-py通过sftp替换网络设备启动文件
- STM32单片机FPGA毕设电路原理论文报告任务驱动教学法在单片机课程教学中的应用
- STM32单片机FPGA毕设电路原理论文报告任务驱动法在单片机教学中的应用
- STM32单片机FPGA毕设电路原理论文报告人造金刚石压机智能化压力测控系统设计
- 以某列为依据匹配多项(Excel版)
- STM32单片机FPGA毕设电路原理论文报告人体短臂离心机实验台的显示控制系统
- STM32单片机FPGA毕设电路原理论文报告人工气候室监控系统的环境控制器研究
- STM32单片机FPGA毕设电路原理论文报告染整自动线张力控制系统的设计
- 数据挖掘与机器学习-实验
- 基于Linux系统Nginx的动态网站的LNMP环境源码包
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈