北航数值分析大作业第二题精解.doc

《北航数值分析大作业第二题精解.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、word格式目标：使用带双步位移的QR分解法求矩阵的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知：(i,j=1,2,……,10)开始输入矩阵A用函数“nishangsanjiadiv”将矩阵A拟上三角化即为A（n-1）用函数“characteristic”求解矩阵A（n-1）即A的所有特征值用函数“qrdiv”将A（n-1）QR分解用函数“characteristicvector”求解矩阵A（n-1）即A的所有实特征值对应的特征向量输出A的所有特征值λ、A的所有实特征值对应的特征向量、拟上三角

2、矩阵A（n-1）、及其Q、R和R*Q算法：以上是程序运作的逻辑，其中具体的函数的算法，大部分都是数值分析课本上的逻辑，在这里特别写出矩阵A的实特征值对应的一个特征向量的求法：由于，因此在经过选主元的高斯消元以后，即的最后一行必然为零，左上方变为n-1阶单位矩阵，右上方变为n-1阶向量，然后令，则。....word格式这样即求出所有A所有实特征值对应的一个特征向量。#include#include#include#defineN10#defineE1.0e-

3、12#defineMAX10000//以下是符号函数doublesgn(doublea){doublez;if(a>E)z=1;elsez=-1;returnz;}//以下是矩阵的拟三角分解voidnishangsanjiaodiv(doubleA[N][N]){inti,j,k;intm=0;doubled,c,h,t;doubleu[N],p[N],q[N],w[N];for(i=0;i

4、i))continue;for(j=i+1,d=0;j

5、}for(j=0;j

6、doubleA[N][N],doubleQ[N][N],doubleR[N][N]){inti,j,k;//intm=0;doubled,c,h;doubleu[N],w[N],p[N];for(i=0;i

7、j][i]<=E)m=m+1;//if(m==(N-1-i))continue;for(j=i,d=0;j

8、];}for(j=0;j