用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
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/21 21:14:18
用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
用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.比较计算结果
用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
雅克比迭代法
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