qr方法求矩陣全部特征值(編輯修改稿)
2024-09-17 13:00 本頁面
【文章內(nèi)容簡介】 =0) return 0。 else{ printf( 矩陣A[][]的行列式為:%.5f\n,det*L[n1][n1])。 return det*L[n1][n1]。 }}int Non_singularMatrix(double L[n][n]){//判斷是否是非奇異矩陣 printf( 判斷矩陣A[][]的行列式是否為0?\n)。 if(Det(L)!=0) return 1。 return 0。}void Disp(double H[n][n]){//輸出矩陣 for(int i=0。in。i++){ for(int j=0。jn。j++) printf(%.6f ,H[i][j])。 printf(\n)。 }}int IsZero(double a[],int j){//判斷數(shù)組是否全為0 for(int i=0。ij。i++) if(a[i]!=0) return 0。 return 1。}int sgn(double y){//符號函數(shù) if(y0) return 1。 else if(y==0) return 0。 else return 1。}void Hessenberg(double A[n][n]){//將矩陣化為上Hessenberg矩陣 printf(原矩陣為:\n)。Disp(A)。//輸出原矩陣 double T[n][n],B[n][n],C[n][n],c[n1],v[n1],u[n1],R[n1][n1],I[n1][n1],t,w,s。 int i,j,k,m。 for(k=0。kn2。k++){ for(i=0。ink1。i++) for(j=0。jnk1。j++) if(i==j) I[i][j]=1。 else I[i][j]=0。//定義單位陣I[][] double max=fabs(A[k+1][k])。//求最大值 for(i=0。ink1。i++){ if(maxfabs(A[i+k+1][k])) max=fabs(A[i+k+1][k])。 } for(i=0。ink1。i++) c[i]=A[i+k+1][k]/max。//標準化數(shù)組 if(IsZero(c,nk1))//判斷數(shù)組是否全為0 continue。//數(shù)組為0，則這一步不需要約化 for(i=0,t=。ink1。i++) t+=c[i]*c[i]。 v[k]=sgn(A[k+1][k])*sqrt(t)。// u[0]=c[0]+v[k]。 for(j=1。jnk1。j++) u[j]=c[j]。//求矩陣u[][] w=v[k]*(c[0]+v[k])。 for(i=0。ink1。i++) for(j=0。jnk1。j++) R[i][j]=I[i][j]u[i]*u[j]/w。 for(i=0。in。i++) for(j=0。jn。j++){ if(i==j) T[i][j]=1。 else T[i][j]=0。//定義矩陣T[][]為單位陣 } for(i=0。ink1。i++) for(j=0。jnk1。j++) T[i+k+1][j+k+1]=R[i][j]。//初等反射陣T[][] for(i=
醫(yī)療健康相關推薦
鄂ICP備17016276號-1