文件名称：幂法 反幂法 带双步位移
文件大小：10KB
文件格式：CPP
更新时间：2012-12-05 05:13:41
幂法 反幂法 带双步位移
void main(void)
{
cout.setf(ios::scientific); //输出值为科学计数法
cout.precision(12); //输出值精度为12位
char f;
int i;
double maxNorm,maxNorm1,lambda[N],lambda_s,miu[N],cond_A,det_A;
double A[UB+DB+1][N];
MatrixInit(A); //为矩阵A赋初值
maxNorm=PowerMethod(A); //求得模最大的特征根
maxNorm1=PMwithTr(A,maxNorm); //求得距离模最大特征根最远的根
lambda_s=AntiPowerMethod(A); //求得模最小特征根
if(maxNorm>maxNorm1) ////////////////////////////
{ //
lambda[0]=maxNorm1; // 判断maxNrom与maxNorm1
lambda[N-1]=maxNorm; // 的大小,大的即为最大特征值，
} // 小的即为最小特征值。
else //
lambda[0]=maxNorm; //
lambda[N-1]=maxNorm1; /////////////////////////////
for(i=1;i<=39;i++) //循环调用带原点平移的
{ //反幂法求最接近特征值
miu[i]=lambda[0]+i*(lambda[500]-lambda[0])/40;
lambda[i]=APMwithTr(A,miu[i]);
}
cond_A=fabs(maxNorm/lambda_s); //求条件数
det_A=Det_A(A); //求行列式
cout<<"lambda_1="<