SSOR迭代法

SSOR迭代法 求解线性方程组 (MATLAB)

在这里插入图片描述

(数值分析作业自编)

代码如下:
*将代码分块,很容易理解。
左端是输出项,右端输入自定义项

function [R,k,x] = SSOR(A,b,error,w)
% SSOR迭代
% R为迭代矩阵B的谱半径;k为迭代步数;x为解向量
% A为系数矩阵;b为常数项;error为指定相邻两次迭代精度;w为松弛因子

%% w的输入控制
while (1) 
    if w>0 & w<=2
       % 用谱半径判断收敛性 不满足敛散性退出程序并报错
       B1=(diag(diag(A))-w*tril(-A,-1))\((1-w)*diag(diag(A))+w*triu(-A,1));
       B2=(diag(diag(A))-w*triu(-A,1))\((1-w)*diag(diag(A))+w*tril(-A,-1));
       R=max(sqrt(eig((B2*B1)'*(B2*B1))));
       if R<1
          break;
       else
          disp('错误,R>1迭代不收敛');
          disp('请尝试重新输入一个w:');
          w=input('请输入w:');
       end
    else 
       disp('请尝试重新输入一个w∈(0,2}:');
       w=input('再次输入w:');
    end  
end

%% 初始化
n=length(b);
x=zeros(n,1);x0=zeros(n,1);x1=zeros(n,1);
k=0;

%% 循环
while (1) 
    %x0上一次迭代的解向量 算法主体
     x0=x;
    for i=1:n
        x1(i)=(1-w)*x0(i)+w*(b(i)-A(i,1:i-1)*x1(1:i-1)-A(i,i+1:n)*x0(i+1:n))/A(i,i);
    end
    for i=n:-1:1
        x(i)=(1-w)*x1(i)+w*(b(i)-A(i,1:i-1)*x1(1:i-1)-A(i,i+1:n)*x(i+1:n))/A(i,i);
    end
    % 用相邻两次迭代解的无穷范数作为判敛原则
    if norm(x-x0,Inf)<error
        break;
    end
    k=k+1; 
end

end


输入:

A=[5 2 1; -1 4 2; 2 -3 10];
b=[-12 20 3]';
w=0.8;
error=1e-4;

调用

[R,k,x]=SSOR(A,b,error,w)

运行

>> test_SSOR
R =
    0.1553
k =
     6
x =
   -4.0000
    3.0000
    2.0000
>> 

在这里插入图片描述

  • 0
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值