数值分析试验幂法与反幂法matlab(共10页).doc

精选优质文档-倾情为你奉上一、问题的描述及算法设计（一）问题的描述我所要做的课题是：对称矩阵的条件数的求解设计1、求矩阵A的二条件数问题 A=2、设计内容：1）采用幂法求出A的 .2）采用反幂法求出A的.3）计算A的条件数 A2* A-12=cond2（A）=/.（精度要求为10-6）3、设计要求1）求出A2。2）并进行一定的理论分析。（二）算法设计1、幂法算法（1）取初始向量u（例如取u=(1,1,1)）,置精度要求，置k=1. （2）计算v=Au,m=max(v), u= v/ m（3）若| m= m|，则停止计算（m作为绝对值最大特征值，u作为相应的特征向量）否则置k=k+1，转（2）2、反幂法算法（1）取初始向量u（例如取u=(1,1,1)）,置精度要求，置k=1. （2）对A作LU分解，即A=LU（3）解线性方程组 Ly=u,Uv=y（4）计算 m=max(v), u= v/ m（5）若|m=m|，则停止计算（1