原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
会计学;第1页/共41页;第2页/共41页;第3页/共41页;计算矩阵的主特征根及对应的特征向量; §1.1 乘幂法; 乘幂法的基本思想是取初始向量v(0)?Rn,作迭代
v(k+1) =Av(k) =Ak+1v(0) , k=0,1,2,…
产生迭代序列?v(k)?.; 1. 设|?1???2??????n? , 这时,(8.1)式可写成; 求矩阵A的按模最大的特征值; 对非零向量v,用max(v)表示v的按绝对值最大的分量,称向量u=v/max(v)为向量v的规范化向量.;所以;其收敛速度由比值|?2/?1|来确定,其值越小收敛越快.; 如用规范化乘幂法解例2,仍取u(0)=v(0)=(1,0)T,则有;k; 2. 设?1=?2=?=?r,且 |?1???r+1??????n? ,这时,(8.1)式可写成; v(2i)??12i(a1x1+a2x2) , v(2i+1)??12i+1(a1x1-a2x2); x1??k+1u(k+1)+?1u(k) , x2??k+1u(k+1)-?1u(k);K;§1.2 加速技术;可见,序列??k?线性收敛于?1 .;? k;则对B应用乘幂法可达到加速收敛的目的。;计算可得;这是因为|?2/?1|=1/2,而|m2/m1|=1/7,故对B应用乘幂法远比对A应用乘幂法收敛的快.;也可将上式改写成; 反幂法还可结合原点位移法应用.设已求得矩阵A的特征值??i的某个近似值;取初始向量u(0)=(1,0.714405,-0.249579)T,对B用反幂法计算可得:;? 反幂法 /* Inverse Power Method */; Jacobi方法是求实对称矩阵全部特征值和特征向量的一种矩阵变换方法。;R的第i个列向量恰为?i的特征向量;; Rpq(?)具有下列性质: ; 设实对称矩阵A=(aij)n?n ,记B=RpqT(?)ARpq(?)=(bij)n?n则它们元素之间有如下关系:;所以有;只需角?满足; 由式(8.7),令t=tan?,则t满足方程; ?是给定的精度要求,则A的特征值可取为?i?aii(k),i=1,2,…,n.;的全部特征值.;所以;;从而A的特征值可取为
?1?2.125825, ?2?8.388761, ?3?4.485401; Jacobi方法具有方法简单紧凑,精度高,收敛较快等优点, 是计算对称矩阵全部特征值和相应特征向量的有效方法,但计算量较大,一般适用于阶数不高的矩阵.
文档评论（0）