clc
clear
n = input('请输入矩阵阶数:\n');
for i = 1:n
fprintf('请输入矩阵第%d行\n',i);
A(i,:) = input('');
end
A
B(:,1) = input('请输入B向量:\n');
B
for i = 1:n
x(i) = 0;
x2(i) = 0;
end
for i =1:4
for j = 1:n
for k = 1:n
if j ~= k
x(j) = x(j)-A(j,k)*x2(k);
end
end
x(j)=x(j) + B(j);
x(j)=x(j)/A(j,j);
end
x2=x;
for j = 1:n
x(j)=0;
end
fprintf('第%d次迭代结果:\n',i);
x2
end
【雅可比(jacobi)迭代法 matlab实现】
推荐阅读
- 数值分析|有限元法,有限差分法和有限体积法的区别
- 数学|matlab三维山峰/山脉/山地曲面数据图
- 数值分析|MATLAB unifrnd 与 rand函数的区别
- 数值分析|Matlab函数linspace吃亏记录
- 判断并将矩阵转化为严格对角占优矩阵