作业帮 > 综合 > 作业

用Matlab 或C语言 都可以 最好能再加点分析过程

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/26 04:37:26
用Matlab 或C语言 都可以 最好能再加点分析过程
给出线性方程组HnX=b,其中系数矩阵Hn为希尔伯特矩阵:
Hn=(hij),hij=1/(i+j-1),i,j=1,2,……n.
假设X*=(1,1,…,1)T,b=HnX*.若取n=6,8,10,分别用雅可比迭代法及SOR迭代(w=1,1.25,1.比较计算结果
雅克比迭代法
function [x,n]=jacobi(A,b,x0,eps,varargin)
if nargin==3
eps= 1.0e-6;
M = 200;
elseif nargin=eps
x0=x;
x=B*x0+f;
n=n+1;
if(n>=M)
disp('Warning:迭代次数太多,可能不收敛!');
return;
end
end
sor
function [x,n]=SOR(A,b,x0,w,eps,M)
if nargin==4
eps= 1.0e-6;
M = 200;
elseif nargin=eps
x0=x;
x =B*x0+f;
n=n+1;
if(n>=M)
disp('Warning:迭代次数太多,可能不收敛!');
return;
end
end