关于方Ax+d=λx的求解 毕设答辩.ppt

* 浙江大学 毕业论文 关于方程Ax + d = λx的求解 答辩人： 导师： 专业：数学系 信息与计算科学 论文框架 第1章 绪论 第2章 相关理论基础 第3章 根据A矩阵的情况对问题的讨论 第4章 复杂度分析和数值例子 第5章 总结 参考文献 附录 第1章 绪论 本文主要研究方程Ax + d = λx的解，其中A是一个给定的n阶矩阵，且情况未知，d是一个给定的非零向量，所需要求解的未知量为向量x与一个数值λ，且向量x中满足x(n-1) = x(n)。如果把这个向量方程看做一个方程组并且把的x(n-1) = x(n)也看做一个方程的话，这个方程组一共有n+1个方程，n+1个未知量，但由于λ的存在，它并不是一个通常的线性方程。 通过对给定但是不确定具体形式的A矩阵进行讨论，根据A矩阵的性质来解原方程。由于原方程的形式非常接近线性方程组，如果先确定或者固定了λ的值，那么，来解这个方程剩余的x的值就很方便了。 为了体现λ在方程中的意义，我将原方程变形为 (*) (λ - A)x = d 这个形式很容易让我们想到最后的解x与λ的值，是否会与A矩阵的特征值和特征向量有关联。 第2章 相关理论基础 线性方程相关理论 线性方程组的解法（直接法、迭代法） 特征值和特征向量相关理论 矩阵的相似与相似对角化 矩阵特征值的计算方法 牛顿法等数值算法 第3章 根据A矩阵的情况对问题的讨论 当A矩阵是对角阵的情况 若A矩阵是一个对角阵，那么我们不妨设对角线上的元素为a1,a2…an-1,an。那么(*)就可以转化为n个非常简单的方程：(λ – ak)xk = dk, k = 1,2,…,n-1,n。然后利用x(n-1) = x(n)这个条件，我们先分析(λ – an-1)xn-1 = dn-1和(λ – an)xn = dn这最后的两个式子，不妨令x = xn-1 = xn，那么这两个方程可以重写作： (λ – an-1)x = dn-1 (λ - an)x = dn 这是个二元方程组，而这个方程里面λ和x的解的情况（是否有解，有解的话是唯一解还是无穷多解）和系数an-1,an,dn-1,dn有相当大的关系。我们只要分情况讨论就可以了。分类情况步骤： （1）dn-1≠0，dn≠0。如果dn-1和dn不相等，那么λ = (an-1 * dn - an * dn-1) / (dn - dn-1)。如果dn-1和dn相等，那就看一下an-1和an的情况：如果an-1和an不相等，那么就一定无解了；如果an-1和an相等，那么就相当于这两个方程其实是一个方程，我们只要找出和前n-1个ak不相等的λ即可。 第3章 根据A矩阵的情况对问题的讨论 当A矩阵是对角阵的情况 （2）dn-1=0，dn≠0。由于dn≠0，所以x一定不能等于0，而dn-1=0，则只能让λ取an-1的值了。 （3）dn-1≠0，dn=0。类似于（2）的情况只能取λ=an。 （4）dn-1=0，dn=0。那么就找出和前n-2个ak不相等的λ即可。 这样，只要求出了λ的值（或者已经判断无解），那么求x的值就非常方便了，只需要把求出的λ值分别代入方程(λ – ak)xk = dk，就可以解出{xk},k= 1,...,n了。而且由于λ是利用x(n-1) = x(n)的性质求出来的，所以这也就保证了解出来的 {xk},k= 1,...,n里面一定是满足x(n-1) = x(n)的。这样，对于A矩阵是对角阵的情况我们就完美解决了。 第3章 根据A矩阵的情况对问题的讨论 当A矩阵是上三角矩阵的情况 我们假定A = {ai,j}i,j = 1,…,n，且当ij时，ai,j = 0。让我们写出方程组的一般表达形式为：λxi – Σai,jxj = di, i = 1,…,n。然后我们看当i = n-1和i = n时的两项，分别为： λxn-1 – an-1,n-1xn-1 – an-1,nxn = dn-1 λxn – an,nxn = dn 令xn-1 = xn = x，并且重新整理上述两项可以得到： (λ – an-1,n-1 – an-1,n)x = dn-1 (λ – an,n)x = dn 那么，我们很容易的得到了和我们在讨论A矩阵是对角阵的时候产生的非常相似的两个方程，那么类似的，我们利用相同的分类讨论方法就可以求出λ或者判断方程是否是无解的了，这里就不重复叙述了。 第3章 根据A矩阵的情况对问题的讨论 当A矩阵是实对称矩阵的情况 当A是实对称矩阵的时候，我们再去写出具体方程组然后考察性质似乎遇到了一点麻烦。那么我们可以利用A可以相似对角化成一个实对角阵的性质来帮助我们解决问题。即