matalb实践（八）：SOR迭代解线性方程

最新推荐文章于 2024-04-15 22:52:56 发布

从零开始的奋豆

最新推荐文章于 2024-04-15 22:52:56 发布

阅读量237

点赞数 1

分类专栏： MATLAB 文章标签： matlab

本文链接：https://blog.csdn.net/m0_68926749/article/details/134431946

版权

MATLAB 专栏收录该内容

12 篇文章 9 订阅

订阅专栏

1.题目

2.SOR 迭代

$\mathrm{x_i^{(k+1)}=(1-o)x_i^{(k)}+\frac{o}{a_{ii}}~\left(b_i-\sum_{j=1}^{i-1}a_{ij}x_j^{(k+1)}-\sum_{j=i+1}^{n}a_{ij}x_j^{(k)}\right),\quad i=1,2,\cdots,n}$

3.代码实现

clc;clear;
a=[4 -1 0;-1 4 -1;0 -1 4];
b=-eye(3);
A=[a b;b a];

B=[2 1 2 2 1 2]';
x1=A^(-1)*B;
x=zeros(6,1);
w=1;
k1=1;
p=zeros(6,10);
while(max(abs(B-A*x))>=1e-5)
   for i=1:6
      x(i)=x(i)+w/4*(B(i)-sum(A(i,:)'.*x));
   end
   p(:,k1)=x;

   k1=k1+1;
end
disp(k1);
disp(x);
x2=zeros(6,1);
w=1.1;
k2=1;
while(sum(abs(B-A*x2))>=1e-5)
   for i=1:6
      x2(i)=x2(i)+w/4*(B(i)-sum(A(i,:)'.*x2));
   end
   k2=k2+1;
end
disp(k2);
disp(x);