毕业设计论文-- 关于方程Ax + d = λx的求解.docVIP

PAGE 3 - 本 科 生 毕 业 论 文（设计） 题目 关于方程Ax + d = λx的求解  目 录 摘要 Abstract 第1章 绪论 1.1 选题背景和意义 1.2 研究目的和思路 第2章 相关理论基础 2.1 线性方程理论 2.2 特征值和特征向量 第3章 根据A矩阵的情况对问题的讨论 3.1 一些特殊情况的解 3.2 比较一般的情况的解 第4章 复杂度分析和数值例子 4.1 对于确定算法的复杂度分析 4.2 数值例子 第5章 总结 参考文献 附录 摘 要 本文通过对矩阵不同形式的讨论，循序渐进的研究方程。在研究讨论这个方程的过程中，给出一系列的算法，并给出算法的复杂度分析和一些数值例子。在通过对不同类型的矩阵讨论这个方程的过程中，运用一些简单的恒等变换和分类讨论来简化问题，并且使用一些经典的数值方法来解决剩下的问题。 关键字：线性方程组； 特征值； 特征向量； 实对称； Jacobi方法 Abstract I study the equation step by step, through discussion of the different forms of matrix. In the research process of the equation, I give a series of algorithms, complexity analysis of these algorithms and some numerical examples. In the progress of discussing the equation by different types of matrix, I use some simple discussion of equation transformation and classification to simplify the problem, and use some classical numerical methods to solve the left problems. Keywords: Linear equation; eigenvalue; eigenvector; real symmetric; Jacobi Method 第1章 绪 论 1.1 选题背景和意义 本文主要研究方程Ax + d = λx的解，其中A是一个给定的n阶矩阵，且情况未知，d是一个给定的非零向量，所需要求解的未知量为向量x与一个数值λ，且向量x中满足x(n-1) = x(n)。如果把这个向量方程看做一个方程组并且把的x(n-1) = x(n)也看做一个方程的话，这个方程组一共有n+1个方程，n+1个未知量，但由于λ的存在，它并不是一个通常的线性方程。 1.2 研究目的和思路 通过对给定但是不确定具体形式的A矩阵进行讨论，根据A矩阵的性质来解原方程。由于原方程的形式非常接近线性方程组，如果先确定或者固定了λ的值，那么，来解这个方程剩余的x的值就很方便了。而我们也可以通过利用x(n-1) = x(n)和A矩阵的性质来确定λ的值或者给出λ的可行范围。这样我们就能够解决这个问题了。为了体现λ在方程中的意义，我将原方程变形为 (*) (λ - A)x = d 而这个形式就是我们在接下来的文章中讨论这个方程的解的主要形式。同时，这个形式很容易让我们想到最后的解x与λ的值，是否会与A矩阵的特征值和特征向量有关联。 第2章 相关理论基础 2.1 线性方程相关理论 2.1.1 线性方程组的Gauss消元法 所谓Gauss消元法就是利用初等变换把线性方程组化为对角矩阵。然后根据化出的形式判断方程是否有解，有解的话由底自上的解出每一个未知数。 2.1.2 线性方程组的直接解法 对于一般规模较小的矩阵我们都可以采用LU分解或者选主元的LU分解来进行。特殊的，如果矩阵是一个对称正定矩阵，那么就可以采用Cholesky分解，把原矩阵分解为LLT的形式。 2.1.3 线性方程组的迭代法 对于规模较大的矩阵，我们一般会采用迭代法的方法来解方程。 符号定义：Ax = b，A = D – L – U，D为对角元矩阵，-L为下三角（无对角元）矩阵，-U为上三角（无对角元）矩阵。记B = D-1(L+U)，M = (D - L)-1U。 （1）Jacobi迭代。 xk+1 = Bxk + D-1b （2）G-S迭代。 x k+1 = Mx k + (D - L) -1b。 （3）SOR迭代。 x k+1 = Lωx k + ω(D - ωL) -1b, 其中Lω = (D - ωL) -1