北航数值第一次大作业课程.doc

一、算法的设计方案： （一）各所求值得计算方法 1、最大特征值λ501，最小特征值λ1，按模最小特征值λs的计算方法 首先使用一次幂法运算可以得到矩阵的按模最大的特征值λ，λ必为矩阵A的最大或最小特征值，先不做判断。对原矩阵A进行一次移项，即（A-λI），在进行一次幂法运算，可以得到另一个按模最大特征值λ0。比较λ和λ0的大小，较大的即为λ501，较小的即为λ1。 对矩阵A进行一次反幂法运算，即可得到按模最小特征值λs。 2、A与μk值最接近的特征值λik的计算方法 首先计算出k所对应的μk值，对原矩阵A进行一次移项，即（A-μkI），得到一个新的矩阵，对新矩阵进行一次反幂法运算，即可得到一个按模最小特征值λi。则原矩阵A与μk值最接近的特征值 λik=λi+μk。 3、A的（谱范数）条件数cond（A）2的计算方法 其中 QUOTE 矩阵A的按模最大和按模最小特征值。 （二）程序编写思路。 由于算法要求A的零元素不存储，矩阵A本身为带状矩阵，所以本题的赋值，LU分解，反幂法运算过程中，均应采用Doolittle分解法求解带状线性方程组的算法思路。 幂法、反幂法和LU分解均是多次使用，应编写子程序进行反复调用。 二、源程序： #includestdio.h #includeiostream #includestdlib.h