北航数值分析第二次大作业--QR分解.docxVIP

北航数值分析第二次大作业--QR分解 北航数值分析第二次大作业--QR分解 PAGE PAGEPAGE 16 北航数值分析第二次大作业--QR分解 PAGE 15 《数值分析A》 计算实习题目二 姓 名 学 号 联系方式 班 级 指导教师 2012年10月  算法设计方案 整个程序主要分为四个函数，主函数，拟上三角化函数，QR分解函数以及使用双步位移求解矩阵特征值、特征向量的函数。因为在最后一个函数中也存在QR分解，所以我没有采用参考书上把矩阵M进行的QR分解与矩阵Ak的迭代合并的方法，而是在该函数中调用了QR分解函数，这样增强了代码的复用性，减少了程序长度；但由于时间关系，对阵中方法的运算速度没有进行深入研究。 1.为了减少QR分解法应用时的迭代次数，首先对给定矩阵进行拟上三角化处理。 2.对经过拟上三角化处理的矩阵进行QR分解。 3.注意到计算特征值与特征向量的过程首先要应用前面两个函数，于是在拟上三角化矩阵的基础上对QR分解函数进行了调用。计算过程中，没有采用goto语句，而是根据流程图采用其他循环方式完成了设计，通过对迭代过程的合并，简化了程序的循环次数，最后在计算特征向量的时候采用了列主元高斯消去法。 二、源程序代码 #include #include #include int i,j,k,l,m; 2e ,A[i][j]); 2e ,R[i][j]); } printf( \n\n); } printf( Q矩阵为:\n\n); 2e ,Q[i][j]); } printf( \n\n); } printf( RQ矩阵为：\n\n); 2e ,RQ[i][j]); } printf( \n\n); } DoublestepsQR(A); 2e ,Re[j]); } printf(\n\n 特征值虚部依次为：\n\n ); 2e ,Im[j]); } 2e ,X[i][j]); } printf( \n\n); } getchar(); } //拟上三角化函数 void Quasiuppertriangular(double A[][10]) { for(j=0;j8;j++) { for(i=0;i10;i++) { U[i]=0; P[i]=0; T[i]=0; W[i]=0; } m=0; for(i=j+2;i10;i++) { if(A[i][j]!=0) { m=m+1; } } if(m==0) { continue; } d=0; for(i=j+1;i10;i++) { d=d+pow(A[i][j],2); } d=sqrt(d); c=-d; if(A[j+1][j]=0) { c=d; } h=c*(c-A[j+1][j]); U[j+1]=A[j+1][j]-c; for(i=j+2;i10;i++) { U[i]=A[i][j]; } for(i=0;i10;i++) { for(k=0;k10;k++) { P[i]=P[i]+U[k]*A[k][i]; } P[i]=P[i]/h; } t=0; for(i=0;i10;i++) { for(k=0;k10;k++) { T[i]=T[i]+U[k]*A[i][k]; } T[i]=T[i]/h; t=t+P[i]*U[i]; } t=t/h; for(i=0;i10;i++) { W[i]=T[i]-t*U[i]; for(k=0;k10;k++) { A[i][k]=A[i][k]-W[i]*U[k]-U[i]*P[k]; if(abs(A[i][k])1e-12) { A[i][k]=0; } } } } } //QR分解函数 void QRdecomposition(double A[][10]) { for(i=0;i10;i++) { for(j=0;