用SSOR迭代法和预优共轭梯度法求解Ax=b方程组(附MATLAB代码)

最新推荐文章于 2025-06-04 09:02:56 发布

原创

最新推荐文章于 2025-06-04 09:02:56 发布 · 4.6k 阅读

29 ·

CC 4.0 BY-SA版权

一、 SSOR 迭代法

SSOR迭代法的迭代格式为

$(D-\omega C_L)x_{m-\frac{1}{2}}=\{\omega C_U+(1-\omega)D\}x_{m-1}+\omega b,\\ (D-\omega C_U)x_m= \{\omega C_L+(1-\omega)D\}x_{m-\frac{1}{2}}+\omega b.$

对应的MATLAB代码为

function [x,k,error,time]=myssor(A,b,x,w,tol,max_it)
%% SSOR迭代算法
%% To solve the equation Ax=b;
%% Iterative formula : (D-wL)x_(k+1)=((1-w)D+wU)x_k+wb; (D-wU)x_(k+1)=((1-w)D+wL)x_k+wb;
% A: input matrix;  b:right vector 
% x: initial vector;
% max_it : the maximum number of iterations;
% tol : Accuracy;
% w : relaxation factor;
% k : number of iterations at termination;
% error : the