幂法_c语言

编写了幂法的子函数，单位化时采用二范数，相关参数可以根据需求改动较为灵活

int main() 
{
	double a[3][3]={{6,-12,6},{-21,-3,24},{-12,-12,51}};
        A_lamada1(a);
	return 0;
}
double A_lamada1(double a[3][3])       //根据幂法计算矩阵A的按模最大的特征值lamada1 
{
    int n,i,j;
    double u[3]={1};//向量初始化[1 0 0]
    double y[3];  //单位化矩阵 
    double l;     //模长
    double temp=0;
    double lamada=0;
    double lamada1=0;
    double e;
	for(n=0;n<=100;n++)
	{
	    for(i=0;i<3;i++)                //计算模长 
	    {
	     temp=temp+u[i]*u[i];
	    }
        l=sqrt(temp);
        temp=0;
        for(i=0;i<3;i++)
        {
          y[i]=u[i]/l;
          u[i]=0;
	    }
	    for(i=0;i<3;i++)              //计算u（k+1） 
	    {
		  for(j=0;j<3;j++)
		  {
		    u[i]=u[i]+a[i][j]*y[j];	
		  }
		} 
	    for (i=0;i<3;i++)                //计算lamada1 
	    {
		 lamada1=lamada1+y[i]*u[i];
	    }
	    e=fabs(lamada1-lamada)/fabs(lamada1);         //精度控制 
	    if(e<eps)                             
	    {
		printf(" %d  lamada1=%10.6e",n,lamada1);
		return lamada1;
		break;
	    }
	    else
	    lamada=lamada1;
	    lamada1=0;
    }
   if(n>100)
    {
  	    printf("error");
    }	
}

100为最大计算次数，计算次数用n记录，例子中的矩阵的精确计算结果为45，eps为计算精度，本例中为1e-6，计算结果如下：