没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
JACOBI迭代法:
X = [0,0,0,0,0,0,0,0,0,0]';
A=[4,2,-3,-1,2,1,0,0,0,0;8,6,-5,-3,6,5,0,1,0,0;4,2,-2,-1,3,2,-1,0,3,1;0,-2,1,5,-1,3,-1,1,9,4;-4,2,6,-1,6,7,-3,3,2,3;8,6,-8,5,7,17,2,6,-3,5;0,2,-1,3,-4,2,5,3,0,1;16,10,-11,-9,17,34,2,-1,2,2;4,6,2,-7,13,9,2,0,12,4;0,0,-1,8,-3,-24,-8,6,3,-1];
b = [5,12,3,2,3,46,13,38,19,-21]';
X1 = A\b;
t = [];
e = 1;
f = 0;
while e
X0=X;
for j = 1 : 10
X(j) = (b(j)-A(j,[1:j-1,j+1:10]) * X0([1:j-1,j+1:10]))/A(j,j);
end
t = [t,X];f = f+1;
for j = 1 : 10
if abs(X(j)-X0(j))>e break; end
end
if j == 10 break;end
end
A
disp('方程组精确解:')
X1
disp('方程组达到精度的解:')
X
disp('方程组每次迭代的解:')
t
disp('方程组达到精度所需迭代次数:')
f
X = [0,0,0,0,0,0,0,0,0,0]';
A=[4,2,-3,-1,2,1,0,0,0,0;8,6,-5,-3,6,5,0,1,0,0;4,2,-2,-1,3,2,-1,0,3,1;0,-2,1,5,-1,3,-1,1,9,4;-4,2,6,-1,6,7,-3,3,2,3;8,6,-8,5,7,17,2,6,-3,5;0,2,-1,3,-4,2,5,3,0,1;16,10,-11,-9,17,34,2,-1,2,2;4,6,2,-7,13,9,2,0,12,4;0,0,-1,8,-3,-24,-8,6,3,-1];
b = [5,12,3,2,3,46,13,38,19,-21]';
X1 = A\b;
t = [];
e = 1;
f = 0;
while e
X0=X;
for j = 1 : 10
X(j) = (b(j)-A(j,[1:j-1,j+1:10]) * X0([1:j-1,j+1:10]))/A(j,j);
end
t = [t,X];f = f+1;
for j = 1 : 10
if abs(X(j)-X0(j))>e break; end
end
if j == 10 break;end
end
A
disp('方程组精确解:')
X1
disp('方程组达到精度的解:')
X
disp('方程组每次迭代的解:')
t
disp('方程组达到精度所需迭代次数:')
f
资源评论
- PlyTools2017-01-06是Matla代码,三种方法挺全,但是上来试了个Jacobi发现是结果错的。。。
littleloris
- 粉丝: 0
- 资源: 2
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功