北航数值分析第二次大作业.docVIP

数值分析第二次大作业 姓名：李潇 学号：SY1303314 题目:使用带双步位移的QR分解法求矩阵的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知： (i,j=1,2,……,10) 1、算法设计： ①矩阵的拟上三角化： 对实矩阵A进行相似变换化为拟上三角矩阵，其变换矩阵采用Householder矩阵，变换过程如下： 若，则； 否则，， ， ， ， 。 当时，得，令又是对称正交矩阵，于是成立，因而与 相似。 ②矩阵的QR分解： 矩阵的QR分解过程与拟上三角化过程相似，在这里不再重复其原理。 ③求全部特征值 矩阵拟上三角化后利用带双步位移的QR方法，采用书本Page 63页具体算法实现。为了使编程方便，并体会goto语句使用的灵活性，程序中的跳转均使用goto Loop*实现。 ④求A的相应于实特征值的特征向量 求实特征值对应的特征向量，即是求解线性方程组，。因此，为得到全部实特征值对应的特征向量，解线性方程组的过程要循环n次（n为矩阵阶数）。线性方程组的求解采用列主元素Gauss消去法。 2、程序源代码： #include Stdio.h #include Conio.h #include math.h //****************************************************************************// // N ：待求矩阵的阶数 // // Epsilon：精度水平 // // MaxL ：限制求解的最大迭代次数 // //****************************************************************************// #define N 10 #define Epsilon 0.000000000001 #define MaxL 200 void Creat_MatrixA(double array[N+1][N+1]); void Load_MatrixA(double MatrixA[N+1][N+1],double arrayA[N+1][N+1]); void Clear_Elements(double array[N+1][N+1]); void Initial_to_I(double array[N+1][N+1]); void Nishang(double arrayA[N+1][N+1]); void QR_Fenjie(double arrayA[N+1][N+1],double Q[N+1][N+1],int m); void Solution_namda(double arrayA[N+1][N+1],double namda[N+1][2]); void Solution_A_Next(double arrayA[N+1][N+1],int m,double s,double t); void Solution_vector(double array[N+1][N+1],double namda[N+1][2],double x[N+1][N+1]); void GaussElemitation_select(int m,int flag,double array[N+1][N+1],double x[N+1][N+1]); void Print_test(double array[N+1][N+1]); void Print_testAB(double arrayA[N+1][N+1],double arrayB[N+1][N+1]); void Print_namda(double namda[N+1][2]); void Print_vector(double x[N+1][N+1],double namda[N+1][2]); int main(void) { //****************************************************************************// // 定义主程序中用到的变量 // //**********