古典Jacobi方法.doc

#includeiostream #includemath.h using namespace std; const int n=3; double e=0.0005,a1,a2,a3,a4,max; int i,j,k,d; double a[n][n]={{1,1,0.5},{1,1,0.25},{0.5,0.25,2}},s[n][n],t[n][n],s1[n][n],t1[n][n],Q[n][n]; double feiduimax(double A[][n]); void selectarray(int p,int q); void main() { double sum; feiduimax(a); a1=(a[j][j]-a[i][i])/(2*a[i][j]); if(a1=0) a2=1/(fabs(a1)+sqrt(1+pow(a1,2))); else a2=-1/(fabs(a1)+sqrt(1+pow(a1,2))); a3=1/sqrt(1+pow(a2,2)); a4=a2*a3; selectarray(i,j); printf(s=\n); for(i=0;in;i++) for(j=0;jn;j++) { printf( %f ,s[i][j]); if(j==n-1) printf(\n); } printf(\n\n); for(i=0;in;i++) { for(j=0;jn;j++) { Q[i][j]=s[i][j]; } } while(max=e) { for(i=0;in;i++) { for(j=0;jn;j++) { s1[i][j]=s[j][i]; } } for(i=0;in;i++) { for(j=0;jn;j++) { sum=0; for(k=0;kn;k++) sum=sum+s1[i][k]*a[k][j]; t1[i][j]=sum; } } for(i=0;in;i++) { for(j=0;jn;j++) { sum=0; for(k=0;kn;k++) sum=sum+t1[i][k]*s[k][j]; a[i][j]=sum; } } printf(a=\n); for(i=0;in;i++) { for(j=0;jn;j++) { printf( %f ,a[i][j]); if(j==n-1) printf(\n); } } printf(\n\n); feiduimax(a); a1=(a[j][j]-a[i][i])/(2*a[i][j]); if(a1=0) a2=1/(fabs(a1)+sqrt(1+pow(a1,2))); else a2=-1/(fabs(a1)+sqrt(1+pow(a1,2))); a3=1/sqrt(1+pow(a2,2)); a4=a2*a3; selectarray(i,j); printf(s=\n); for(i=0;in;i++) { for(j=0;jn;j++) { printf( %f ,s[i][j]); if(j==n-1) printf(\n); } } printf(\n\n); for(i=0;in;i++) { for(j=0;jn;j++) { sum=0; for(k=0;kn;k++) sum=sum+Q[i][k]*s[k][j]; Q[i][j]=sum; } } } printf(Q=\n); for(i=0;in;i++) { for(j=0;jn;j++) { printf( %f ,Q[i][j]); if(j==n-1) printf(\n); } } printf(\n\n); for(i=0;in;i++) { printf( 特征值r%d=%f\n,i,a[i][i]); printf( 特征向量x%d=\n,i); for(j=0;jn;j++) printf( %f\n,Q[j][i]); printf(\n\n); } } double feiduimax(double A[][3]) { int m,k1; max=fabs(A[0][1]