数值分析实验报告基本QR算法求全部特征值.docVIP

PAGE  PAGE 3 数值分析实验报告 专业 信息与计算科学 班级 信计101 姓名 学号 协作队员 实验日期 2013 年 1 月 5 日 星期 六 成绩评定 教师签名 批改日期 题目 问题提出 给定矩阵, （1）用Matlab函数“eig”求矩阵全部特征值。 （2）用幂法求A的主特征值及对应的特征向量。 （3）用基本QR算法求全部特征值（可用Matlab函数“qr”实现矩阵的QR分解）。 模型建立 用幂法求A的主特征值及对应的特征向量的模型:选取v0=(1,1,1,1,1)T,按照下列公式构造向量序列{uk}{vk} v0=u0≠0vk=Auk-1μk=maxvkuk=vk/μk 则有 limk→∞uk=x1maxx1 limk→∞μk=λ1 循环足够多次后，可以近似得出， λ1=μk x1=uk 求解方法 (1) A=[2 3 4 5 6;4 4 5 6 7;0 3 6 7 8;0 0 2 8 9;0 0 0 1 0]; a=eig(A) (2) pmethod.m function [l,v,