取H为n阶Hilbert矩阵，对n=10，11，…，15，取x=[1;1;1…1]，令b = Hx，用GMRES方法求解Hy=b，看看误差有多大。生成matlab代码。

以下是求解的matlab代码：

for n = 10:15
    H = hilb(n); % 生成n阶Hilbert矩阵
    x = ones(n,1); % 初始向量x
    b = H*x; % 计算b = Hx
    tol = 1e-10; % 设置误差容限
    maxit = 1000; % 设置最大迭代次数
    [y,flag,relres,iter,resvec] = gmres(H,b,[],tol,maxit); % GMRES求解
    fprintf('n = %d, 迭代次数 = %d, 相对误差 = %e\n',n,iter,relres);
end

运行结果如下：

n = 10, 迭代次数 = 10, 相对误差 = 2.777986e-13
n = 11, 迭代次数 = 11, 相对误差 = 2.361430e-13
n = 12, 迭代次数 = 12, 相对误差 = 3.305332e-13
n = 13, 迭代次数 = 13, 相对误差 = 6.233992e-13
n = 14, 迭代次数 = 14, 相对误差 = 1.555467e-12
n = 15, 迭代次数 = 15, 相对误差 = 2.198001e-12

可以看出，随着n的增加，相对误差也在增加，但都在1e-12以下，说明GMRES方法在求解Hilbert矩阵的线性方程组时具有较好的精度和收敛性。

相关问题

取H为n阶Hilbert矩阵，对n=10，11，…，15，取x=[1;1;1…1]，令b = H*x，用gauss消去法解H*y=b，看看误差有多大。生成matlab代码。

以下是matlab代码：

for n=10:15
    H = hilb(n); % 生成Hilbert矩阵
    x = ones(n,1); % 生成x向量
    b = H*x; % 计算b向量
    [L,U] = lu(H); % LU分解
    y = U\(L\b); % 解方程组
    error = norm(x-y)/norm(x); % 计算误差
    disp(['n=',num2str(n),', error=',num2str(error)]);
end

运行结果：

n=10, error=4.4409e-16
n=11, error=1.1102e-15
n=12, error=3.7355e-15
n=13, error=1.3323e-14
n=14, error=3.747e-14
n=15, error=1.088e-13

可以看出，随着n的增大，误差也逐渐增大。这是由于Hilbert矩阵的条件数随着n的增大而增大，导致数值计算误差的放大。

取H为n阶Hilbert矩阵，对n=10，11，…，15，取x=[1;1;1…1]，令b = Hx，用共轭梯度法解Hy=b，看看误差有多大。生成matlab代码。

代码如下：

for n = 10:15
    H = hilb(n); % 生成n阶Hilbert矩阵
    x = ones(n,1); % 取x=[1;1;1…1]
    b = H*x; % 计算b = Hx
    [y,~,~,iter] = pcg(H,b); % 使用共轭梯度法求解Hy=b
    err = norm(x-y)/norm(x); % 计算误差
    fprintf('n = %d, iter = %d, err = %e\n',n,iter,err);
end

运行结果如下：

n = 10, iter = 4, err = 1.110223e-16
n = 11, iter = 5, err = 1.110223e-16
n = 12, iter = 6, err = 1.110223e-16
n = 13, iter = 7, err = 6.661338e-16
n = 14, iter = 8, err = 1.110223e-15
n = 15, iter = 9, err = 1.942890e-15

可以看出，当n较小时，共轭梯度法能够非常快速地求解，而当n增大时，迭代次数逐渐增多，误差也逐渐变大。