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

《实验6反幂法求矩阵按模最小特征值》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、西华大学计算机系学生上机实践报告西华数学与计算机学院上机实践报告课程名称：计算方法A年级：上机实践成绩：指导教师：严常龙姓名：上机实践名称：反幂法求矩阵按模最小特征值学号：上机实践日期：上机实践编号：1上机实践时间：一、目的1．通过本实验加深对反幂法的构造过程的理解；2．能对反幂法提出正确的算法描述编程实现，得到计算结果。二、内容与设计思想自选方阵，用反幂法求解其按模最小特征值。可使用实例：三、使用环境操作系统：WindowsXP软件环境：MicrosoftVisualC++四、核心代码及调试过程#include#include#defin

2、eMAX_N20//矩阵最大维数#defineMAXREPT100#defineepsilon0.00001//求解精度intmain(){intn;inti,j,k;doublexmax,oxmax;staticdoublea[MAX_N][MAX_N];staticdoublel[MAX_N][MAX_N],u[MAX_N][MAX_N];staticdoublex[MAX_N],nx[MAX_N];printf("请输入矩阵阶数n:");//输入矩阵维数scanf("%d",&n);if(n>MAX_N){第5页共5页西华大学计算机系学生上机实践报告printf(

3、"theinputnislargerthanMAX_N,pleaseredefinetheMAX_N.");return1;}if(n<=0){printf("pleaseinputanumberbetween1and%d.",MAX_N);return1;}//输入A矩阵printf("请输入矩阵的值a[i][j]i,j=0...%d;",n-1);for(i=0;i

4、{for(j=0;jxmax)xmax=fabs(nx[j]);for(j=0;j

5、=0;i

6、)//计算U矩阵{u[k][j]=a[k][j];for(i=0;i<=k-1;i++)u[k][j]-=(l[k][i]*u[i][j]);u[k][j]/=l[k][k];}}for(i=0;i=0;j--){x[j]=nx[j];for(k=j+1;k

7、k];}xmax=0.0;第5页共5页西华大学计算机系学生上机实践报告for(j=0;jxmax)xmax=fabs(x[j]);for(j=0;j