资源描述：
《数值分析幂法课程设计---用幂法求矩阵的最大特征值》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、课程设计报告课程名称计算方法课程设计课题名称用幂法求矩阵的最大特征值专业班级学号姓名指导教师13目录1.设计内容与设计要求………………………………………………………32.问题的描述及算法设计……………………………………………………43.算法的流程图………………………………………………………………44.算法的理论依据及其推导…………………………………………………45.相关的数值结果……………………………………………………………86.数值计算结果的分析………………………………………………………97.附件…………………………………………………………………………1
2、08.总结…………………………………………………………………………1213一、设计内容与设计要求1．设计内容：对课程《计算方法》中的常见算法进行综合设计或应用（具体课题题目见后面的供选题目）。2．设计要求：l课程设计报告正文内容a.问题的描述及算法设计；b.算法的流程图（要求画出模块图）；c.算法的理论依据及其推导；d.相关的数值结果（通过程序调试)，；e.数值计算结果的分析；f.附件（所有程序的原代码，要求对程序写出必要的注释）。l书写格式a．要求用A4纸打印成册b．正文格式：一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。c．
3、正文的内容:正文总字数要求在3000字左右（不含程序原代码）。d．封面格式如下页。l考核方式指导老师负责验收程序的运行结果，并结合学生的工作态度、实际动手能力、创新精神和设计报告等进行综合考评，并按优秀、良好、中等、及格和不及格五个等级给出每位同学的课程设计成绩。具体考核标准包含以下几个部分：a．平时出勤（占10%）b．系统需求分析、功能设计、数据结构设计及程序总体结构合理与否（占10%）c．程序能否完整、准确地运行，个人能否独立、熟练地调试程序（占40%）d．设计报告（占30%）注意：不得抄袭他人的报告（或给他人抄袭），一旦发现，成绩为零分。13e．独立
4、完成情况（占10%）。l课程验收要求a．判定算法设计的合理性，运行相关程序，获得正确的数值结果。b．回答有关问题。c．提交课程设计报告。d．提交软盘（源程序、设计报告文档）。e．依内容的创新程度，完善程序情况及对程序讲解情况打分。3、进度安排1、班级：信息与计算科学：0501、0502、05032、主讲教师：聂存云3、辅导教师：聂存云、赵钍焱4、时间安排：第16周星期一8时：00分——16时：30分星期三8时：00分——16时：30分星期五8时：00分——16时：30分星期天8时：00分——11时：00分13一、问题的描述及算法设计对于给定矩阵A用幂法求解
5、矩阵A的最大特征值；本算法用于求矩阵按模最大的特征值及其相应特征向量。程序执行后，先通过键盘输入矩阵、迭代初值向量、精度控制和迭代允许最大次数，程序即可给出每次迭代的次数和对应的迭代特征值、特征向量及误差序列，它们都按10位有效数输出。其中最后输出的结果即为所求的特征值和特征向量序列。如果迭代超出次还没有求出满足精度的根则输出迭代超限提示，此时可以根据输出序列判别收敛情况。二、算法的描述；（1）输入矩阵、初始向量，误差；（2）；（3）计算；（4）；（5）；（6）如果,则显示特征值和对应的特征向量),终止；（7）,转（3）注：如上算法中的符号表示取向量中绝对
6、值最大的分量。本算法使用了数据规范化处理技术以防止计算过程中出现益出错误。三、算法的理论依据及其推导；幂法是通过求矩阵特征向量来求出特征值的一种迭代法.其基本思想是:若我们求某个n阶方阵A的特征值和特征向量,先任取一个初始向量X(0),构造如下序列:13X(0),X(1)=AX(0),X(2)=AX(1),…,X(K)=AX(K+1),…⑴当k增大时,序列的收敛情况与绝对值最大的特征值有密切关系,分析这一序列的极限,即可求出按模最大的特征值和特征向量.假定矩阵A有n个线性无关的特征向量.n个特征值按模由大到小排列:│λ1│>=│λ2│>=…>=│λn│⑵其
7、相应的特征向量为:V1,V2,…,Vn⑶它们构成n维空间的一组基.任取的初始向量X(0)由它们的线性组合给出X(0)=a1V1+a2V2+…+anVn⑷由此知,构造的向量序列有X(k)=AX(k-1)=A2X(k-2)=…=AkX(0)=a1λ1kV1+a2λ2kV2+…+anλnkVn⑸下面按模最大特征值λ1是单根的情况讨论:由此公式(5)可写成X(k)=λ1k(a1V1+a2(λ2/λ1)kV2+…+an(λn/λ1)kVn)⑹若a1≠0，由于
8、λi/λ1
9、<1(i≥2),故k充分大时，X(k)=λ1k(a1V1+εk)其中εk为一可以忽略的小量，这说
10、明X(k)与特征向量V1相差一个常数因子，即使a1=0，由于计算过
