资源描述:
《数值分析试验幂法及反幂法matlab》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、.-一、问题的描述及算法设计〔一〕问题的描述我所要做的课题是:对称矩阵的条件数的求解设计1、求矩阵A的二条件数问题A=2、设计容:1〕采用幂法求出A的.2〕采用反幂法求出A的.3〕计算A的条件数ⅡAⅡ2*ⅡA-1Ⅱ2=cond2〔A〕=/.〔精度要求为10-6〕3、设计要求1〕求出ⅡAⅡ2。2〕并进展一定的理论分析。〔二〕算法设计1、幂法算法〔1〕取初始向量u〔例如取u=(1,1,…1)〕,置精度要求,置k=1.〔2〕计算v=Au,m=max(v),u=v/m〔3〕假设
2、m=m
3、<,那么停顿计算〔m作为绝对值最大特征值,u作为相应的特征向量
4、〕否那么置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〕假设
5、m=m
6、<,那么停顿计算〔1/m作为绝对值最小特征值,u.可修编..-作为相应的特征向量〕;否那么置k=k+1,转〔3〕.二、算法的流程图〔一〕幂法算法的流程图noyes开场k=0;m1=0v=A*u[vmax,i]=max(abs(v))m=v(i);u=v/mabs(m-m1)<1e-6index=1;br
7、eak;输出:m,u,index完毕m1=m;k=k+1.可修编..-〔二〕反幂法算法的流程图noyes开场输入A;[m,u,index]=pow_inv(A,1e-6)k=0;m1=0v=invA*u[vmax,i]=max(abs(v))m=v(i);u=v/mabs(m-m1)<1e-6index=1;break;输出:m,u,index完毕m1=m;k=k+1输入A;[m,u,index]=pow(A,1e-6).可修编..-三、算法的理论依据及其推导〔一〕幂法算法的理论依据及推导幂法是用来确定矩阵的主特征值的一种迭代方法,也即,绝
8、对值最大的特征值。稍微修改该方法,也可以用来确定其他特征值。幂法的一个很有用的特性是它不仅可以生成特征值,而且可以生成相应的特征向量。实际上,幂法经常用来求通过其他方法确定的特征值的特征向量。1、幂法的迭代格式与收敛性质设n阶矩阵A的特征值,,…,是按绝对值大小编号的,x(i=1,2,…,n)为对应的特征向量,且为单根,即
9、
10、>
11、
12、≥…≥
13、
14、那么计算最大特征值与特征向量的迭代格式为v=Au,m=max(v),u=v/m〔1〕其中max(v)表示向量v绝对值的最大分量。2、对于幂法的定理按式〔1〕计算出m和u满足m=,u=〔二〕反幂法算法的理
15、论依据及推导反幂法是用来计算绝对值最小的特征值突然相应的特征向量的方法。是对幂法的修改,可以给出更快的收敛性。1、反幂法的迭代格式与收敛性质设A是非奇异矩阵,那么零不是特征值,并设特征值为
16、
17、≥
18、
19、≥…≥
20、
21、>
22、
23、那么按A的特征值绝对值的大小排序,有
24、
25、>
26、
27、≥…≥
28、
29、对A实行幂法,就可得A的绝对值最大的特征值1/和相应的特征向量,即A的绝对值最小的特征值和相应的特征向量。由于用A.可修编..-代替A作幂法计算,因此该方法称为反幂法,反幂法的迭代格式为v=Au,m=max(v),u=v/m〔2〕2、对于反幂法的定理按式〔2〕计算出的m和u满
30、足:m=,u=在式〔2〕中,需要用到A,这给计算带来很大的不方便,因此,把〔2〕式的第一式改为求解线性方程组Av=u〔3〕但由于在反幂法中,每一步迭代都需求解线性方程组〔3〕式,迭代做了大量的重复计算,为了节省工作量,可事先把矩阵A作LU分解,即A=LU所以线性方程组〔3〕改为Ly=u,Uv=y四、相关的数值结果〔一〕幂法程序的运行结果m=3.4142u=-0.7071index=11.0000-0.7071〔二〕反幂法程序的运行结果m=0.5858u=0.7071index=11.00000.7071〔三〕矩阵A的二条件数的结果ⅡAⅡ2*
31、ⅡA-1Ⅱ2=cond2〔A〕=m/m=3.4142/0.5858=5.828269五、数值计算结果的分析求n阶方阵A的特征值和特征向量,是实际计算中常常碰到的问题。对于n阶矩阵A,假设存在数和n维向量x满足Ax=x〔1〕.可修编..-那么称为矩阵A的特征值,x为相应的特征向量。由线性代数知识可知,特征值是代数方程
32、I-A
33、=+a+…+a+a=0〔2〕的根。从外表上看,矩阵特征值与特征向量的求解问题似乎很简单,只需求解方程〔2〕的根,就能得到特征值,再解齐次方程组〔I-A〕x=0〔3〕的解,就可得到相应的特征向量。上述方法对于n很小时是可以
34、的。但当n稍大时,计算工作量将以惊人的速度增大,并且由于计算带有误差,方程〔2〕未必是准确的特征方程,自然就不必说求解方程〔2〕与〔3〕的困难了。本次实验所用的幂法和反幂法分别是