QR.rar_特征向量

#include "stdio.h" #include "5sstq.c" #include "5strq.c" main() { int i,j,k,l=60; static double q[5][5],b[5],c[5]; static double a[5][5]={ {10.0,1.0,2.0,3.0,4.0}, {1.0,9.0,-1.0,2.0,-3.0},{2.0,-1.0,7.0,3.0,-5.0}, {3.0,2.0,3.0,12.0,-1.0},{4.0,-3.0,-5.0,-1.0,15.0}}; double eps=0.000001; strq(a,5,q,b,c); k=sstq(5,b,c,q,eps,l); printf("MAT A IS:\n"); for (i=0; i<=4; i++) { for (j=0; j<=4; j++) printf("%13.7e ",a[i][j]); printf("\n"); } printf("\n"); if (k>0) { printf("MAT B IS:\n"); for (i=0; i<=4; i++) printf("%13.7e ",b[i]); printf("\n\n"); printf("MAT Q IS:\n"); for (i=0; i<=4; i++) { for (j=0; j<=4; j++) printf("%13.7e ",q[i][j]); printf("\n"); } printf("\n"); } } int sstq(n,b,c,q,eps,l) int n,l; double b[],c[],q[],eps; { int i,j,k,m,it,u,v; double d,f,h,g,p,r,e,s; c[n-1]=0.0; d=0.0; f=0.0; for (j=0; j<=n-1; j++) { it=0; h=eps*(fabs(b[j])+fabs(c[j])); if (h>d) d=h; m=j; while ((m<=n-1)&&(fabs(c[m])>d)) m=m+1; if (m!=j) { do { if (it==l) { printf("fail\n"); return(-1); } it=it+1; g=b[j]; p=(b[j+1]-g)/(2.0*c[j]); r=sqrt(p*p+1.0); if (p>=0.0) b[j]=c[j]/(p+r); else b[j]=c[j]/(p-r); h=g-b[j]; for (i=j+1; i<=n-1; i++) b[i]=b[i]-h; f=f+h; p=b[m]; e=1.0; s=0.0; for (i=m-1; i>=j; i--) { g=e*c[i]; h=e*p; if (fabs(p)>=fabs(c[i])) { e=c[i]/p; r=sqrt(e*e+1.0); c[i+1]=s*p*r; s=e/r; e=1.0/r; } else { e=p/c[i]; r=sqrt(e*e+1.0); c[i+1]=s*c[i]*r; s=1.0/r; e=e/r; } p=e*b[i]-s*g; b[i+1]=h+s*(e*g+s*b[i]); for (k=0; k<=n-1; k++) { u=k*n+i+1; v=u-1; h=q[u]; q[u]=s*q[v]+e*h; q[v]=e*q[v]-s*h; } } c[j]=s*p; b[j]=e*p; } while (fabs(c[j])>d); } b[j]=b[j]+f; }