实验6反幂法求矩阵按模最小特征值

西华大学计算机系学生上机实践报告 西华数学与计算机学院上机实践报告课程名称：计算方法A年级：上机实践成绩：指导教师：严常龙姓名： 上机实践名称：反幂法求矩阵按模最小特征值学号：上机实践日期：上机实践编号：1上机实践时间：一、目的1通过本实验加深对反幂法的构造过程的理解；2能对反幂法提出正确的算法描述编程实现，得到计算结果。二、内容与设计思想自选方阵，用反幂法求解其按模最小特征值。可使用实例：三、使用环境操作系统：Windows XP软件环境：Microsoft Visual C+四、核心代码及调试过程#include#include#define MAX_N 20 /矩阵最大维数#define MAXREPT 100#define epsilon 0.00001 /求解精度int main()int n;int i,j,k;double xmax,oxmax;static double aMAX_NMAX_N;static double lMAX_NMAX_N,uMAX_NMAX_N;static double xMAX_N,nxMAX_N;printf(n请输入矩阵阶数n:); /输入矩阵维数scanf(%d,&n);if(nMAX_N) printf(the input n is larger than MAX_N,please redefine the MAX_N.n); return 1;if(n=0) printf(please input a number between 1 and %d.n,MAX_N); return 1; /输入A矩阵printf(请输入矩阵的值aij i,j=0.%d;n,n-1);for(i=0;in;i+) for(j=0;jn;j+) scanf(%lf,&aij);for(i=0;in;i+) xi=1;oxmax=0;for(i=0;iMAXREPT;i+) for(j=0;jn;j+) /幂乘 nxj=0; for(k=0;kn;k+) nxj+=ajk*xk; xmax=0.0; for(j=0;jxmax) xmax=fabs(nxj); for(j=0;jn;j+) nxj/=xmax; for(j=0;jn;j+) xj=nxj; if(fabs(xmax-oxmax)epsilon) printf(solve.max lamda=%lfn,xmax); /输出 printf(the vector is:n); for(i=0;in;i+) printf(%fn,nxi); break; /return 0; oxmax=xmax;/printf(after %d repeat ,max no result .n,MAXREPT);for(i=0;in;i+) uii=1; /U矩阵对角元为for(k=0;kn;k+) for(i=k;in;i+) /计算L矩阵 lik=aik; for(j=0;j=k-1;j+) lik-=(lij*ujk); for(j=k+1;jn;j+) /计算U矩阵 ukj=akj; for(i=0;i=k-1;i+) ukj-=(lki*uij); ukj/=lkk; for(i=0;in;i+) xi=1;for(i=0;iMAXREPT;i+) for(j=0;jn;j+) /反幂乘 nxj=xj; for(k=0;k=0;j-) xj=nxj; for(k=j+1;kn;k+) xj-=ujk*xk; xmax=0.0; for(j=0;jxmax) xmax=fabs(xj); for(j=0;jn;j+) xj/=xmax; if(fabs(xmax-oxmax)epsilon) printf(solve. min lamda=%lfn,1/xmax); /输出 printf(the vector is:n); for(i=0;in;i+) printf(%fn,xi); break; /return 0; oxmax=xmax;return 1;五、总结本次试验利用C语言实现了用反幂法求解其按模最小特