Hilbert矩阵病态线性代数方程组的求解教学文稿.docVIP

Hilbert矩阵病态线性代数方程组的求解 实验一 病态线性代数方程组的求解 1.估计Hilbert矩阵2-条件数与阶数的关系 运行tiaojianshu.m 输入m=10 可以得到如下表的结果 阶数 1 2 3 4 5 条件数 1 19.28 524.05 1.55e+4 4.76e+5 阶数 6 7 8 9 10 条件数 1.49e+7 4.75e+8 1.52e+10 4.93e+11 1.60e+13 2.选择不同维数，分别用Guass消去（LU分解），Jacobi迭代，GS迭代，SOR迭代求解，比较结果。 说明：Hx=b，H矩阵可以由matlab直接给出，为了设定参考解，我们先设x为分量全1的向量，求出b,然后将H和b作为已知量，求x，与设定的参考解对比。 对于Jacobi迭代，GS迭代，SOR迭代，取迭代初值x0为0向量，迭代精度eps=1.0e-6， 迭代次数100000, SOR迭代中w=1.2和0.8分别计算。 n=5 x分量 Gauss法 J迭代 GS迭代 SOR迭代 实际解x w=1.2 w=0.8 x(1) 1.0000 -8.8179 0.9998 0.9999 0.9998 1 x(2) 1.0000 -Inf 1.0028 1.0019 1.0031 1 x(3) 1.0000 -Inf 0.9879 0.9919 0.9863 1 x(4)