用QR算法求矩阵的特征值(共7页).doc

一、实验名称：用QR算法求矩阵的特征值二、实验目的：1、通过实验进一步熟悉掌握求矩阵特征值的QR方法及原理。 2、理解QR方法的计算流程。 3、能够编程实现QR方法。三、实验内容：给定矩阵 , ，采用QR方法计算A和H矩阵的全部特征值。四、实验要求：（1） 根据QR算法原理编写程序求矩阵及矩阵的全部特征值（要求误差10）。（2） 直接用MATLAB的内部函数eig求矩阵及矩阵的全部特征值，并与（1）的结果比较。五、QR方法计算矩阵特征值的程序：function namda,time,data_na=qr_tz(A,tol)if nargin=1;tol=1e-5;endwucha=1;time=0;while (wuchatol)&(time500)q,r=qr(A);A1=r*q;tz0=diag(A1);tz1=diag(A);wucha=norm(tz0-tz1);A=A1;time=time+1;data_na(time,:)=