病态希尔伯特方程组解法-高等数值分析大作业.pdfVIP

病态方程组的求解 理论分析表明，数值求解病态线性方程组很困难，考虑求解如下的线 性方程组的求解： 1 h = , , = 1,2 … . . Hx=b，其中H 是Hilbert 矩阵，H=(hij)n×n，ij i+j−1 1.估计矩阵的2-条件数和矩阵阶数的关系； 2.对不同的n，取x = (1,1,1 …,1) ∈Rn，利用Hx=b 得出b.然后分别用 Guass 消去，Jacobi 迭代，Gauss-seidel 迭代，SOR 迭代求解方程： Ax=b ； 把近似解x (k) 与准确解x 对比，计算结果； 3.综合计算结果，讨论病态方程组的求解。 解：1.用 MATLAB 计算 1-9 阶 Hilbert 矩阵的 2-条件数得到的结果如下表所示 （cond2 （H）表示n 阶Hilbert 矩阵的2-条件数）： 1-9 阶hilbert 矩阵的2-条件数 n 1 2 3 4 5 6 7 8 9 cond2(H ) 1 19.2815 524.0568 1.55E+04 4.77E+05 1.49E+07 4.75E+08 1.53E+10 4.93E+11 对表中数据进行分析可发现，随矩阵的阶数的增加，2-条件数越来越大，后项 与前项之比约为30，2-条件数呈近似指数式增长。对cond2 （H）取以10 为底 的对数，可发现在前几项log10(cond2(Hn))随n 的增加约呈线性增加关系。 n 为20 时，log10(cond2(Hn))与n 的关系与下图所示： 此时可发现当n ≤13时，log10(cond2(Hn))与 n 接近于线性关系，n 13后， log10(cond2(Hn))变化十分缓慢。为了对图像的线性进行比较，对n 为1-13 部分 进行线性拟合，结果如下图所示： 从上图可看出，n 较小时，log10(cond2(Hn))与n 成近似线性关系。 继续增大n，n=200 时，log10(cond2(Hn))与n 如下图所示： 由图可看出，n 为200 时，log10(cond2(Hn))与n 的关系与n=20 时相近，在 n 较小时 log10(cond2(Hn))与 n 的关系接近线性关系，当 n 大于一定的值后 log10(cond2(Hn))变化十分缓慢，但有上升的趋势。 因此hilbert 矩阵的 2-条件数与阶数的关系是：hilbert 矩阵的 2-条件数本身随 矩阵的阶数的增加会以近似于指数形式增长，随阶数的增加，其2-条件数会变得 很大，但对其做一定处理，如取hilbert 矩阵的2-条件数以10 为底的对数后，可 发现，2-条件数以 10 为底的对数整体很小，在阶数较小时其值与阶数的关系接 近线性关系，当阶数大于一定的值后，其值还是增加得十分缓慢，变化较小。 2.使用Guass 消去，Jacobi 迭代，Gauss-seidel 迭代，SOR 迭代求解 近似解： T 定迭代初值为（0.99,0.99,0.99,0.99 ），迭代次数为20000 ，判断迭代法是否收 || （）− （−1）||2 敛的标准都是 ≤1−06 ，即前后两项二范数的