qr方法求矩陣全部特征值(文件)
【正文】 f(不是上Hessenberg矩陣!將其化為上Hessenberg矩陣.\n)。 QRAlgorithm(A)。}MATLAB程序為：并且，與MATLAB軟件求解的結(jié)果進行比較，兩者的結(jié)果是一致的，這也證明了程序的正確性。 程序中的標識符的說明（類型、含義等）見程序。用此組數(shù)據(jù)時define n 5*/ double A1[3][3]={6,2,1,2,3,1,1,1,1}。//將矩陣A轉(zhuǎn)化為上Hessenberg矩陣 } else printf(是上Hessenberg矩陣!不需要將其再化為上Hessenberg矩陣!\n)。 } else printf(矩陣A[][]是非奇異矩陣!符合分解條件!\n)。j++) Z[i][j]=A[i][j]。 for(int i=0。 } printf(\n)。 for(i=0。Disp(Q)。j++) if(fabs(A[i][j])max) max=fabs(A[i][j])。//求A矩陣下三角的絕對值最大值 for(i=1。k++) t+=R[i][k]*Q[k][j]。jn。j++) Q[i][j]=V[i][j]。 for(i=0。in。jn。 F[i][j]=y。kn。jn。p[m+1][m+1]=c。 s=A[m+1][m]/sqrt(A[m][m]*A[m][m]+A[m+1][m]*A[m+1][m])。jn。mn1。jn。 double p[n][n],Q[n][n],R[n][n],F[n][n],V[n][n],c,s,t,y,max。j++) if(E[i][j]!=0) return 0。}int IsHessenberg(double E[n][n]){//判斷是否是上Hessenberg矩陣 for(int i=2。jn。 C[i][j]=s。jn。 B[i][j]=s。jn。j++) T[i+k+1][j+k+1]=R[i][j]。//定義矩陣T[][]為單位陣 } for(i=0。i++) for(j=0。jnk1。//求矩陣u[][] w=v[k]*(c[0]+v[k])。// u[0]=c[0]+v[k]。//數(shù)組為0，則這一步不需要約化 for(i=0,t=。 } for(i=0。//定義單位陣I[][] double max=fabs(A[k+1][k])。i++) for(j=0。 for(k=0。}void Hessenberg(double A[n][n]){//將矩陣化為上Hessenberg矩陣 printf(原矩陣為:\n)。 return 1。 printf(\n)。in。 }}int Non_singularMatrix(double L[n][n]){//判斷是否是非奇異矩陣 printf( 判斷矩陣A[][]的行列式是否為0?\n)。 } det*=L[k][k]。 L[i][k]=t。 } det*=1。jn。j++) if(maxL[j][k]){ max=L[j][k]。k++){ double max=L[k][k]。 數(shù)據(jù)結(jié)構(gòu)描述主要數(shù)據(jù)成員說明double A[n][n] 存放矩陣Adouble Q[n][n],存放QR分解式的正交矩陣Qdo
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1