SOR方法：输入一个方阵，将其分解为对角矩阵、下三角矩阵和上三角矩阵 - MATLAB开发
5星
- 浏览量: 0
- 大小:None
- 文件类型：None
简介：
本MATLAB项目实现SOR（Successive Over-Relaxation）方法，用于将给定的方阵分解成对角矩阵、下三角矩阵和上三角矩阵，适用于线性代数问题求解。
函数[x] = SOR_HW(A,b,x_0,omega)
% 输入方阵A、向量b以及初始x值和松弛因子omega
N = 1000; % 迭代次数上限
n = length(A); % 矩阵维度
tol = 0.0001; % 收敛容许误差
x = zeros(n, 1);
% 将方阵A分解为三个矩阵：对角矩阵(D)、严格下三角矩阵(L)和严格上三角矩阵(U)
D = diag(diag(A));
L = -tril(A,-1);
U = -triu(A,1);
a = (D-omega*L);
for i=1:N
x = a\(((1-omega)*D + omega*U)*x_0) + omega*(a\b);
if norm(x-x_0)
全部评论 (0)
还没有任何评论哟~
