北航数值分析大作业一

;Р?} Р?printf("λ1=%.12e\n",d1);Р?printf("λ501=%.12e\n",d501);Р//用反幂法求λsР?chuzhi(u);Р?ds=fmifa(u,MatrixC,0);Р?printf("λs=%.12e\n",ds);Р//用带原点平移的反幂法求λikР?for(i=1;i<=39;i++)Р?{Р??a=d1+(i*(d501-d1))/40;Р??assignment(MatrixC);Р??chuzhi(u);Р??d=fmifa(u,MatrixC,a);Р??printf("与μ%02d=%+.12e最接近的特征值λi%02d=%+.12e\n",i,a,i,d);Р?}Р//求A的条件数Р?d=fabs((d1/ds));Р?printf("A的（谱范数）条件数cond<A>2=%.12e\n",d);Р//求detAР?assignment(MatrixC);Р?LU(MatrixC);Р?a=1;Р?for(i=0;i<=500;i++)Р?{Р??a*=MatrixC[2][i];Р?}Р?printf("行列式detA=%.12e\n",a);РР//测试不同迭代初始向量对λ1计算结果的影响。РР?printf("改变迭代初始向量对λmax计算结果的测试如下：\n");Р assignment(MatrixC);Р?for(i=0;i<=500;i++)Р?{Р chuzhi2(u,i);Р?d1=mifa(u,MatrixC,0);Р printf("u%03d,λmax=%+e ",i,d1);Р?if(((i+1)%3)==0)Р?printf("\n");Р?}Р?printf("Press any key to continue\n");Р?getchar();Р}Р 程序结果：