outf.cpp

用QR法求矩阵全部特征值,qr算法是幕法的推广和变形

#include <stdio.h>
#include  <malloc.h>
void inarry(float x[], int *n);
void outf(float r[],int n)
{
 float *x,*x0;
 x=(float*)calloc(n*n,4);
 x0=(float*)calloc(n*n,4);
 inarry(x,&n);	
 float f;
 float eps=0.5;
 int i,j,k;
for(k=0;k<n*2;k++)
 { 
   for(i=0;i<n;i++)
   {for(j=0;j<n;j++)
    {if(i==j)
     x0[i*n+j]=x[i*n+j]-r[k];
     else
     x0[i*n+j]=x[i*n+j];
     
   }
  }
  
 f=x0[0]*x0[4]*x0[8]+x0[1]*x0[5]*x0[6]+x0[2]*x0[3]*x0[7]-x0[2]*x0[4]*x0[6]- x0[1]*x0[3]*x0[8]-x0[0]*x0[5]*x0[7];
 printf("验证求的的R是否为真解f=%f",f);
 if(f<eps)
 printf("   是真实解true\n");
 else 
 printf("   不是真实解false\n");
 k++;

}
free(x0);
}

//f=x0[0]*x0[5]*x0[10]*x0[15]+x0[1]*x0[6]*x0[11]*x0[12]+x0[2]*x0[7]*x0[8]*x0[13]+x0[3]*x0[4]*x0[9]*x0[14]-x0[3]*x0[6]*x0[9]*x0[12]-x0[2]*x0[5]*x0[8]*x0[15]-x0[1]*x0[4]*x0[11]*x0[14]-x0[0]*x0[7]*x0[10]*x0[13];