北航数值分析大作业(二).docxVIP

《数值分析A》计算实习题目二姓名： 学号：学院： 2014年11月算法设计方案：算法结构概述：根据题目的设定，在整个程序中共需要四个函数：矩阵拟上三角化函数、矩阵QR分解函数、带双步位移QR分解求矩阵特征值函数（在这个函数中需要调用矩阵的QR分解函数）、列主元高斯消去法求矩阵特征值对应特征向量的函数，在设计带双步位移QR分解求矩阵特征值函数时，逻辑结构不采用教材上的结构，而选择老师在讲解时推荐的形式，在下一小节会详细介绍。程序总体逻辑结构框图：带双步位移QR分解求矩阵特征值函数的逻辑结构：将矩阵A拟上三角化，确定精度水平e，最大迭代次数L；，令k=1，m=N;若，则得到A的一个特征值，并置m=m-1，转（4）；否则转（5）；若m=1，则得到最后1或2个特征值，转（9）；否则，转（3）；若，则得到两个特征值，m=m-2，转（4）；否则转（6）；若k=L，则停止；否则转（7）；k=k+1；转（3）；结束程序源代码#includestdio.h#includemath.h#includestring.h#include stdlib.h#define N 10int i,j,k,l,m;double temp,c,h,tr,b,Tr,det;double A[N][N],A_orig[N][N],R[N][N],Q[N][N],RQ[N][N]; double temp1[N][N],t