北航数值分析大作业 第一题 幂法与反幂法.docVIP

数 值 分 析（B） 大 作 业（一） 姓名： 学号： 电话： 1、算法设计： ①求、和的值： ：表示矩阵的按模最小特征值，为求得直接对待求矩阵A应用反幂法即可。 、：若矩阵A的特征值满足关系 且，要求、及时，可按如下方法求解： 对矩阵A用幂法，求得按模最大的特征值。 按平移量对矩阵A进行原点平移得矩阵，对矩阵B用反幂法求得B的按模最小特征值。 则：，即为所求。 ②求和A的与数最接近的特征值（k=0，1，…39）： 求矩阵A的特征值中与P最接近的特征值的大小，采用原点平移的方法： 先求矩阵 B=A-PI 对应的按模最小特征值，则+P即为矩阵A与P最接近的特征值。在本次计算实习中则是先求平移矩阵，对该矩阵应用反幂法求得，则与最接近的A的特征值为： 重复以上过程39次即可求得（k=0，1，…39）的值。 ③求A的（谱范数）条件数和行列式： 在（1）中用反幂法求矩阵A的按模最小特征值时，要用到Doolittle分解方法，在Doolittle分解完成后得到的两个矩阵分别为L和U，则A的行列式可由U阵求出，即：det(A)=det(U)。 求得det(A)不为0，因此A为非奇异的实对称矩阵，则： ，和分别为模最大特征值与模最小特征值。 2、程序源代码： #include Stdio.h #include Conio.h #include math.h //****************************************************************************// // 在存储带状矩阵时，下面的几个量在程序中反复用到，为方便编程故把它们定义成宏.// // M ：转换后的矩阵的行数，M=R+S+1。 // // N ：转换后的矩阵的列数，与原矩阵列数相等。 // // S ：上半带宽。 // // R ：下半带宽。 // // Epsilon：用来指定求解精度。 // //****************************************************************************// #define M 5 #define N 501 #define R 2 #define S 2 #define Epsilon 0.000000000001 void Creat_MatrixA(double array[M+1][N+1]); void Load_MatrixA(double arrayA[M+1][N+1],double arrayB[M+1][N+1]); void Load_vectoru(double u[N+1]); double Mifa(double u[N+1],double array[M+1][N+1]); double Fanmifa(double u[N+1],double arrayA[N+1][N+1]); void Solution_Yushu(double array[N+1][N+1],double u[N+1],double namda1,double namda501); void Doolittle_Dai(double array[M+1][N+1]); void Back_Doolittle_Dai(double array[M+1][N+1],double y[N+1],double x[N+1]); double max(double a,double b); double min(double a,double b); main(void) { //****************************************************************************// // 定义主程序中用到的变量 // //*****************************************************