线性方程数值解中jacobi迭代matlab实现
2007-11-24 12:09
579 查看
%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;
end
%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;
end
相关文章推荐
- 数值分析 jacobi迭代法求解线性方程组 MATLAB程序实现
- #公式与实现# Jacobi迭代 Gauss-Seidel迭代
- #公式与实现# Jacobi迭代与五点迭代
- Jacobi迭代算法matlab版
- Jacobi迭代和Gauss迭代 c语言实现
- Jacobi迭代和Gauss迭代 c语言实现
- 雅可比(Jacobi)迭代算法的C++实现
- 线性方程数值分析中高斯迭代matlab实现
- matlab实现不动点迭代、牛顿法、割线法
- #公式与实现# Jacobi迭代的并行实现
- #公式与实现# Jacobi迭代 Gauss-Seidel迭代
- #公式与实现# Jacobi迭代与五点迭代
- #公式与实现# Jacobi迭代的并行实现
- matlab实现ICP(3D迭代最近点算法)
- Matlab实现——严格对角占优三对角方程组求解(高斯赛尔德Gauss-Seidel迭代、超松弛)
- hog matlab代码1——实现hog+svm图像二分类
- Viola–Jones object detection framework--Rapid Object Detection using a Boosted Cascade of Simple Features中文翻译 及 matlab实现(见文末链接)
- 数据挖掘——线性神经网络的Matlab实现
- 动态星间链路分析及其STK/Matlab实现
- 最小二乘法曲线拟合以及Matlab实现