幂法求特征值实验报告4.docVIP

本科实验报告 课程名称： 计算机数值方法B 实验项目： 幂法求特征值 实验地点： ZSA401 专业班级： 软件19班 学号： 学生姓名： \ 指导教师： 李志 2012年 4月 16日 幂法求特征值 一、实验目的和要求 (1) 培养理论联系实际的思路，训练综合运用已经学过的理论和实际知识去分析和解决实际问题的能力。 (2) 帮助学生全面消化已学的相关课程内容，深刻理解计算数值方法课程的内涵，培养使用电子计算机进行科学计算和解决问题的能力。 (3) 进行基本技能训练和巩固。使学生得到选择算法、编写程序、分析数值结果、写数值实验报告、课堂讨论等环节的综合训练。 二、实验内容和原理 使用幂法求A模为最大的特征值及其相应的特征向量。 三、主要仪器设备 使用的计算机：HP ProBook 6450b，软件环境：Win-tc 四、操作方法与实验步骤 #include stdio.h #include math.h main() { int i,j,k; float x[3],y[3],temp[3],m; int a[3][3]; printf(please input the array\n); for(i=0;i3;i++) for(j=0;j3;j++) scanf(%d,a[i][j]); printf(please input the starting values\n); for(i=0;i3;i++) scanf(%f,x[i]); for(k=0;k7;k++) { y[0]=x[0]*a[0][0]+x[1]*a[0][1]+x[2]*a[0][2]; y[1]=x[0]*a[1][0]+x[1]*a[1][1]+x[2]*a[1][2]; y[2]=x[0]*a[2][0]+x[1]*a[2][1]+x[2]*a[2][2]; for(i=0;i3;i++) printf(y[%d]=%f ,i,y[i]); printf(\n); for(i=0;i3;i++) temp[i]=y[i]; for(i=0;i2;i++) if(fabs(temp[i])fabs(temp[i+1])) temp[i+1]=temp[i]; m=temp[2]; printf(m=%f\n,m); for(i=0;i3;i++) x[i]=y[i]/m; for(i=0;i3;i++) printf(x[%d]=%f ,i,x[i]); printf(\n); printf(\n); } getch(); } 五、实验数据记录和处理 实验结果： 六、实验结果与分析 优点算法简便，便于理解计算，实现了课程设计的算法要求及功能。缺点是有很多问题，由于时间问题该算法做得比较粗糙，还不能很好的处理的问题。 七、讨论、心得 刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。 这段时间的实验课提高了我的分析问题，解决问题的能力，特别提高了对一个程序的整体操作能力，对程序的细微之处有了明显的提高认识的地方，追求最实用的程序，弥补学习上的不足，同时认识到还应深入理解课本上的知识，学过的东西要知道理论与实践相结合，增加动手能力。