系列文章目录
提示:这里可以添加系列文章的所有文章的目录,目录需要自己手动添加
个人学习记录:matlab进行雅可比迭代计算
提示:写完文章后,目录可以自动生成,如何生成可参考右边的帮助文档
前言
提示:这里可以添加本文要记录的大概内容:
雅可比迭代法matlab编程
提示:以下是本篇文章正文内容,下面案例可供参考
一、雅可比迭代条件,
系数矩阵绝对对角占优。
二、代码
代码如下
A=[4 -1 1 ;4 -8 1 ;-2 1 5];
B=[7;-21;15];
tol=0.0000001;%容差设置
max1=500;%最大迭代次数
x0=[0;0;0];%正确计算结果给定
n=length(B);%B的长度,即B的列数
k=0;
x=zeros(n,1);%迭代初值给定
while k<max1
for i=1:n
x(i)=(B(i)-(A(i,[1:i-1,i+1:n])*x0([1:i-1,i+1:n],1)))/A(i,i);
end
if norm(x-x0,inf)<tol
break;
end
x0=x;
k=k+1;
end
总结
提示:这里对文章进行总结:
norm(A,p)函数用法,表示矩阵或向量的范数,p=inf表示2范数(行范数),返回最大行向量