函数[x] = SOR_HW(A,b,x_0,omega)% 输入方阵A,b,初始x和omega的值格式长; N = 1000; %迭代次数n = 长度(A); tol = 0.0001; x =零(n,1); %将方阵A分解为三个矩阵:对角矩阵(D); 严格下三角矩阵(L); 严格上三角矩阵(U) D = 诊断(诊断(A)); L =-tril(A,-1); U = -triu(A,1); a = (D-欧米茄*L); 对于 i=1:N x = a\(((1-omega)*D + omega*U)*x_0) + omega*(a\b); 如果范数(x-x_0)<tol 休息; 结尾x_0=x; 结尾结尾
2021-11-07 13:32:29
1KB
matlab
1