一、问题描述
利用迭代方法求解线性方程组。
二、实验目的
掌握Jacobi 方法和Gauss-Seidel 方法的原理，能够编写代码实现两种迭代方法；能够利用代码分析线性方程组求解中的误差情况。
三、实验内容及要求
-
用代码实现：对下列方程中重新组织，并用Jacobi 方法和Gauss-Seidel 方法求解。（注意重新组织的过程也需要利用代码实现，每次迭代都输出产生的近似解）
u − 8v − 2w = 1
u + v + 5w = 4
3u − v + w =− 2 -
对于元素为 Aij = 5/(i + 2j − 1) 的n × n 矩阵， 令x* = [1, . . . , 1]^T， b = Ax* . 使用MATLAB 中的反斜线命令计算双精度解Xc。代码计算 Ax=b 的误差放大因子并与A 的条件数进行比较，观察Xc与误差放大因子之间的关系。
(a) n = 6 (b) n = 10
四、算法原理
-
给出Jacobi 方法和Gauss-Seidel 方法的原理:
当使用迭代方法求解线性方程组时，Jacobi方法和Gauss-Seidel方法是两种常见的迭代方法。
Jacobi 方法：
Jacobi方法是一种迭代法，用于求解线性方程组 Ax = b，其中A是系数矩阵，x是未知向量，b是右侧常数向量。迭代过程如下：
-
将系数矩阵A分解为三个部分：A = D + L + U，其中D是对角矩阵，L是严格下三角矩阵（主对角线以下元素为零），U是严格上三角矩阵（主对角线以上元素为零）。
-
将线性方程组写成迭代形式：x^(k+1) = D^(-1) * (b - (L + U) * xk)，其中k是迭代步数，x(k)是第k次迭代的近似解。
-
重复迭代步骤2，直到满足预定的收敛准则（例如，达到最大迭代次数或误差小于某个阈值）。
Gauss-Seidel 方法：
Gauss-Seidel方法是Jacobi方法的改进版本，它也是一种迭代法用于求解线性方程组 Ax = b。迭代过程如下：
-
将系数矩阵A分解为三个部分：A = D + L + U，其中D是对角矩阵，L是严格下三角矩阵，U是严格上三角矩阵。
-
将线性方程组写成迭代形式：x^(k+1) = (D - L)^(-1) * U * x^k + (D - L)^(-1) * b，其中k是迭代步数，x^(k)是第k次迭代的近似解。
-
重复迭代步骤2，直到满足预定的收敛准则。
Gauss-Seidel方法相较于Jacobi方法的优势在于它使用了当前迭代步中已经更新的近似解，而不是等待整个迭代步完成后再更新。这使得Gauss-Seidel通常比Jacobi方法更快地收敛。然而，对于某些情况，Jacobi方法可能更适用。选择使用哪种方法通常取决于问题的性质和收敛的要求。
- 给出求解线性方程组中误差分析涉及的概念及计算方法。
当进行线性方程组的误差分析时，我们关注一系列关键概念以评估数值解的准确性和稳定性。这些概念涉及不同方面的误差，从残差到相对误差，以及与矩阵条件数相关的内容。
残差是近似解对应方程组的左侧与右侧的差异。通过计算残差，我们可以评估近似解的逼近程度，即 r = b − A x r = b - Ax r=b−Ax。
在误差分析中，我们关注的另一个概念是误差向量，表示真解与近似解之间的差异。误差向量 e = x ∗ − x e = x* - x e=x∗−x提供了关于数值解准确性的信息。
为了量化误差的大小，我们使用误差范数，它是误差向量的范数。不同的范数（如欧几里得范数和无穷范数）提供了不同的视角，用于评估误差的大小。
相对误差是一种衡量真解与近似解之间相对差异的指标。通过计算相对误差，我们可以了解数值解相对于真实解的准确度。
与误差相关的另一个关键概念是条件数，它描述了线性方程组对输入数据的敏感性。矩阵条件数是评估问题的数值稳定性的指标。
最后，误差放大因子用于衡量误差在迭代求解过程中如何被放大。这有助于我们了解数值方法的收敛性和数值解的稳定性。
通过综合考虑这些概念，我们能够进行全面的误差分析，评估数值解的质量，并了解在数值计算中可能引入的误差。
五、测试数据及结果
-
贴图给出迭代情况展示两种迭代方法的收敛性。
-
给出 n = 6 和 n = 10 的近似值及对应的误差放大因子和条件数。
六、总结与思考
在完成实验中，我涉及了线性方程组的数值解法，具体而言是Jacobi方法和Gauss-Seidel方法的实现。通过这个实验，我得到了一些总结与思考。
首先，在实现Jacobi和Gauss-Seidel方法时，我深刻理解了这两种迭代方法的原理。Jacobi方法与Gauss-Seidel方法相比，更为简单，每次迭代都是基于上一次迭代的近似解。而Gauss-Seidel方法则更为迅速，因为它在每次迭代中使用了当前步骤已更新的近似解。这让我对这两种方法的迭代过程和计算方式有了更深入的认识。
其次，通过实验，我对于线性方程组的误差分析有了更全面的了解。我学到了如何计算残差、误差向量以及使用误差范数和相对误差来评估数值解的准确性。此外，我了解了条件数的概念，它是衡量线性方程组数值稳定性的关键指标。
在编写MATLAB代码的过程中，我进一步熟悉了MATLAB语言的语法和一些常用函数，如矩阵运算和迭代。在解决代码错误和逻辑问题的过程中，我提高了调试和排错的能力。
总体而言，这个实验为我提供了一个深入学习线性代数、数值计算和MATLAB编程的机会。通过亲自实践，我更好地理解了这些概念和方法，同时也提高了在数值计算中解决实际问题的能力。这不仅对我的学术研究有益，还增强了我在编程和算法方面的技能。
