数值分析幂法和反幂法_第1页
数值分析幂法和反幂法_第2页
数值分析幂法和反幂法_第3页
数值分析幂法和反幂法_第4页
数值分析幂法和反幂法_第5页
已阅读5页,还剩15页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

幂法和反幂法求矩阵特征值随机产生一对称矩阵,对不同的原点位移和初值(至少取3个)分别使用幂法求计算矩阵的主特征值及主特征向量,用反幂法求计算矩阵的按模最小特征值及特征向量,并比较不同的原点位移和初值说明收敛。认真读题,了解问题的数学原形;选择合适问题求解的数值计算方法;设计程序并进行计算;对结果进行解释说明;采采用方法及结果说明对于幂法和反幂法求解矩阵特征值和特征向量的问题将从问题分析,算法设计和流程图,理论依据,程序及结果进行阐述该问题。一.问题的分析:求n阶方阵A的特征值和特征向量,是实际计算中常常碰到的问题,如:机械、结构或电磁振动中的固有值问题等。对于n阶矩阵A,若存在数人和n维向量x满足TOC\o"1-5"\h\zAx=人x(1)则称人为矩阵A的特征值,x为相应的特征向量。由高等代数知识可知,特征值是代数方程|人I-A|=人n+a人n-i+—+a人+a=0(2)的根。从表面上看,矩阵特征值与特征向量的求解问1题似乎很简单,只需求解方程(2)的根,就能得到特征值人,再解齐次方程组(人I-A)x=0(3)的解,就可得到相应的特征向量。上述方法对于n很小时是可以的。但当n稍大时,计算工作量将以惊人的速度增大,并且由于计算带有误差,方程(2)未必是精确的特征方程,自然就不必说求解方程(2)与(3)的困难了。幂法是一种计算矩阵主特征值(矩阵按模最大的特征值)及对应特征向量的迭代方法,特别是用于大型稀疏矩阵。反幂法是计算海森伯格阵或三角阵的对应一个给定近似特征值的特征向量的有效方法之一。二.算法设计及流程图1、幂法算法(1)取初始向量u(o)(例如取u(o)=(1,1,...1)Q,置精度要求e,置k=1.(2)计算v(k)=Au(k-i),m=max(v(k)),u(k)=v(k)/m(3)若|m广m.1|<e,则停止计算(m.作为绝对值最大特征值人1,u俄)作为相应的特征向量)否则置k=k+1,转(2)2、反幂法算法(1)取初始向量U(0)(例如取U(0)^1,1,…1)T),置精度要求E,置k=1.(2)对A作LU分解,即A=LU(3)解线性方程组Ly(k)=u(k-i),Uv(k)=y(k)(4)计算m=max(v(k)),u(k)=v(k)/m(5)若|m广m;1|<e,则停止计算(1/m「作为绝对值最小特征值牛u(k)作为相应的特征向量);否则置k=k+1,转(3)."幂法流程图:开始m1=m;k=k+1结束反幂法流程图输入入;[m,u,index]=pow(A,1e-6)k=0;m1=rv=A*ur[vmax,i]=max(abs(v))1rm=v(i);u=v/m一abs(m-m1)Z1e-e'''''''''''vindex=1;break;输出:m,u,indexabs(m-m1)<1e-6开始k=0;m1=0v=invA*um1=m;k=k+1输入A;结束(一)幂法算法的理论依据及推导index=1;break;输出:m=v(i);u=v/mm,abs(m-m1)<1e-6开始k=0;m1=0v=invA*um1=m;k=k+1输入A;结束(一)幂法算法的理论依据及推导index=1;break;输出:m=v(i);u=v/mm,u,index三、算法的理论依据及其推导[vmax,i]=max(abs(v))输入A;[m,u,index]=pow_inv(A,1e-6)幂法是用来确定矩阵的主特征值的一种迭代方法,也即,绝对值最大的特征值。稍微修改该方法,也可以用来确定其他特征值。幂法的一个很有用的特性是它不仅可以生成特征值,而且可以生成相应的特征向量。实际上,幂法经常用来求通过其他方法确定的特征值的特征向量。1、幂法的迭代格式与收敛性质设n阶矩阵A的特征值%,人2,…,人〃是按绝对值大小编号的,x‘(iF,2,…,n)为对应人.的特征向量,且七为单根,即IX|>|人|2,・・2|人|则计算最大特征值与特征向量的迭代格式为v(k)=Au(k-i),m=max(v(k)),u(k)=v(k)/m⑴上'其中max(v(k))表示向量v(k)绝对值的最大分量。2、对于幂法的定理按式(1)计算出mk和u(k)满足limm=X,limu(k)=——ik—Xk—X(二)反幂法算法的理论依据及推导反幂法是用来计算绝对值最小的特征值忽然相应的特征向量的方法。是对幂法的修改,可以给出更快的收敛性。1、反幂法的迭代格式与收敛性质设A是非奇异矩阵,则零不是特征值,并设特征值为|X|n|x|n,・・n|x|>|X|TOC\o"1-5"\h\z则按A-i的特征值绝对'值的大小排序,有“^|:|>|《|>W|1|nn一11对A-i实行幂法,就可得A-i的绝对值最大的特征值1/Xn和相应的特征向量,即A的绝对值最小的特征值和相应的特征向量。n由于用A-i代替A作幂法计算,因此该方法称为反幂法,反幂法的迭代格式为v(k)=A-iu(k-i),m=max(v(k)),u(k)=v(k)/m(2)2、对于反幂法的定理*按式(2)计算出的mk和u(k)满足:1xk-*k厂'k-*max(x)nn在式(2)中,需要用到A-i,这给计算带来很大的不方便,因此,把(2)式的第一式改为求解线性方程组Av(k)=u(k-i)(3)但由于在反幂法中,每一步迭代都需求解线性方程组(3)式,迭代做了大量的重复计算,为了节省工作量,可事先把矩阵A作LU分解,即A=LU所以线性方程组(3)改为Ly(k)=u(k-i),Uv(k)=y(k)四、算法程序设计代码幂法程序,在matlab中建立一个M文件并保存。%pow.mfunction[m,u,index,k]=pow(A,u,ep,it_max)ifnargin<4it_max=1000;endifnargin<3ep=1e-5;endn=length(A);index=0;k=0;m1=0;m0=0;I=eye(n);T=A-m0*I;whilek<=it_maxv=T*u;[vmax,i]=max(abs(v));m=v(i);u=v/m;ifabs(m-m1)<ep;index=1;break;endm=m+m0;m1=m;k=k+1;end在matlab输入面板,输入A=rand(4);%产生一个4维随机矩阵B=A+A’;u=[1111]’;%设立初始向量[m,u,index,k]=pow(B,u,ep,it_max)%最多可省略2个参数程序结束。在M文件中可以通过改变m0的值改变原点位移,从而达到原点位移加速。反幂法程序设计代码:

在matlab中建立一个M文件并保存。%pow_inv.mfunction[m,u,index,k]=pow_inv(A,u,ep,it_max)ifnargin<4it_max=1000;endifnargin<3ep=1e-5;endn=length(A);index=0;k=0;m1=0;m0=0;I=eye(n);T=A-m0*I;invT=inv(T);whilek<=it_maxv=invT*u;[vmax,i]=max(abs(v));m=v(i);u=v/m;ifabs(m-m1)<epindex=1;break;endm1=m;k=k+1;endm=1/m;m=m+m0;在matlab输入面板,输入A=rand(4);%产生一个4维随机矩阵B=A+A’;u=[1111]’;%设立初始向量[m,u,index,k]=pow_inv(B,u,ep,it_max)%最多可省略2个参数[m,u,index,k]=pow_inv(B,u,ep,it_max)%最多可省略2个参数程序结束。在M文件中可以通过改变m0的值改变原点位移,从而达到原点位移加速。【结果显示】%在M0=1e-4>>B=rand(4);>>A=B+B’A=0.26750.57760.63441.31300.57761.15030.76410.13670.63440.76410.02570.41931.31300.13670.41931.2248>>u=[1111]';>>[m,u,index,k]=pow(A,u)m=2.6813u=0.85760.69340.56231.0000index1=k=49修改M0=1e-3m=2.6814u=0.85760.6934

0.56231.0000二k=1001修改M0=0%此时为幂法m=2.6815u=0.85760.69350.56231.0000inTk=10修改U=[1234]修改M0=1e-4m=2.6813u=

0.85760.69340.56231.0000k=9修改M0=1e-3m=2.6805u=0.85760.69340.56221.0000inTk=修改M0=0m=2.6814u=0.85760.6934

0.56231.0000index=k=9修改U=[3567]修改M0=1e-4m=2.6819u=0.85770.69370.56241.0000ind:k=修改M0=1e-3m=2.6814=xapui0000'I死9如0ZG69'0ZZG80=n0源9・乙=in0=0/现割1001=习0=xapui0000'I£乙9如0州69'09药80

Um0muindexk[0.85760.69340.5623[1111]0.00012.68131.0000][0.58760.69340.56231490.0012.68141.0000][0.85760.69350.56230100102.68151.0000][0.85760.69340.5623110[1234]0.00012.68131.0000][0.85760.69340.5622190.0012.68051.0000][0.85760.69340.56231702.68141.0000][0.85770.69370.562419[3567]0.00012.68191.0000][0.85760.69340.5623170.0012.69141.0000][0.85770.69370.56240100102.692反幂法结果显示:在m0为0时1.0000]17

>>A=rand(4)A=0.81800.34120.8385□.54660.66020.53410.56810.44490.34200.72710.3704□.63460.28970.30930.T027□.6213»E=A+A'B=1.63591.00141.18050.83631.00141.Ot5S21.Z9520.75421.18051.29520.MOB1.33730.83630.TC421.39F31.242G>>u=[111I]311—1111>>.JiijUjiiideK,k]=pow_inv(Bjle—4)m=0.3037u=-0.54751.CICICICI0.12S6-0.5040ind.es=1k二16M0=0.001U=[1111]1.63591.00141.1805□.83631.00141.06S21.2952□.5M21.18051.2952□.fias1.53730.83630.57421.39731.2426»u=[l11iru=1111>>[n,UjindeKjk]=povr_inv(BjLe-4)n=0.3847u=-0.8996i.oooa0.272S-0.2364index=1k=1.63591.00141.1805□.83631.00141.06821.2952□.5M21.1S051.29520.74081.53730.S3630.5M21.33^31.2426»n=[l11ir1111[iRjUjLndeXjk]=pov_mv(Bj1e-4)n=0.3847u=-0.89951.00000.2725-0.2364irides=1B=1.C3591.00141.18050.83631.C014i.oess1.29520.57421.18051.29520.74081.39730.83630.57421.39731.2426»u.=[l357]nu=1357>>[ikjv.sindeZjk]=pow_inv(瓦u)m=0.384?u=-0.89961.00000.2726-0.2354iniex-1B=1.63591.00141.1S050.03631.00141.06S21.29520.57421.10051.29520.74081.39730.83630.57421.39731.2426»u=[l35r],u一1357>>[iRj乌indeXjk]=pow_inv(BJu)m=0.384?u=-0.89951.00000.2726-0.2364index=1B=1.63591.00141.10050.S3631.00141.06821.29520.57421.1S051.29520.740S1.39730.83630.57421.39731.2426»u.=[l357]nu=1357»[iRjUjindeKjk]=pmv(Bju)in=0.3347-0.09951.00000.2F26-0.2364indeK=1B=1.63591.OOM1.18050.83631.00141.06821.29520.57421.18051.29520.74081.39730.83C30.57421.39731.2426»u=[2345],u-234-EiiLjUjinde^jk]=pov_inv(Bju)0.3847-0.89961.00000.2726-0.23C4indeK=1k二1.63591.00141.18050.S3631.00141.06321.29520.57421.18051.29E20.74081.39730.03630.57421.39731.2426»34S]‘u=2345>>[mjUjindsKjk]=pcw_inT7(BJn=0.3047u=-0.B9951.00000.2726-0.2364indeK=1k=1.63591.00141.1S050.03631.00141.06821.29520.57421.18051.29520.74081.39730.83630.574

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论