Jacobi迭代与SOR迭代求解希尔伯特矩阵

最新推荐文章于 2023-07-17 15:22:42 发布

qqlu_did_lq

最新推荐文章于 2023-07-17 15:22:42 发布

阅读量1.2w

点赞数 18

分类专栏：矩阵文章标签：迭代希尔伯特

本文链接：https://blog.csdn.net/qqlu_did/article/details/45300819

版权

本文通过MATLAB程序演示了Jacobi迭代和SOR迭代法求解希尔伯特矩阵线性方程组的过程。实验表明，Jacobi迭代法导致解发散，而SOR迭代法（ω在1到1.5之间变化）得到的解是收敛的，且收敛速度与矩阵维度和松弛因子ω有关。

摘要由CSDN通过智能技术生成

给出线性方程组 Hn*x = b，其中系数矩阵Hn为希尔伯特矩阵：

假设 x ∗ =(1, 1, . . . , 1)T，b = Hnx ∗。若取 n = 6,8, 10，分别用 Jacobi

迭代法及 SOR迭代（ω = 1, 1:25,1:5）求解，比较计算结果。

MATLAB源码如下，运行Demo_Jacobi_SOR.m文件，附件包含另外两个函数文件，分别为：Jacobi.m与SOR.m。

Demo_Jacobi_SOR.m

clear all
clc
n=[6 8 10];
for i=1:length(n)
H{i}=hilb(n(i));
size_H{i}=size(H{i},1);
x_true{i}=ones(1,size_H{i});
b{i}=x_true{i}*H{i};
x_ini{i}=zeros(1,size_H{i});
%accuracy=5*(10^-6);
accuracy=0.01;
end
%Jacobi
disp('Jacobi');
for i=1:length(n)
    fprintf('The dimension is %d\n',n(i))
    [x{i},k{i}]=Jacobi(H{i},b{i},x_ini{i},accuracy,x_true{i});
end
%SOR
disp('SOR');
w=1;
disp('SOR w=1');
for i=