北航数值分析Jacobi法计算矩阵特征值.docVIP

北航数值分析Jacobi法计算矩阵特征值.doc

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
准备工作 ? 算法设计 矩阵特征值的求法有幂法、Jacobi法、QR法等,其中幂法可求得矩阵按模最大的特征值(反幂法可求得按模最小特征值),Jacobi法则可以求得对称阵的所有特征值。 分析一:由题目中所给条件λ1≤λ2≤…≤λn,可得出λ1、λn按模并不一定严格小于或大于其他特征值,且即使按模严格小于或大于其他特征值,也极有可能出现|λs|λ1||λn |或|λs|λn||λ1 |的情况,导致按幂法和反幂法无法求解λ1或λn二者中的一者; 分析二:题目要求求解与数μk =λ1+k(λn-λ1)/40最接近的特征值λik(k=1,2,3…39),这个问题其实可以转换为求A-μk 按模最小的特征值的问题,但因为在第一个问题中无法确定能肯定的求得λ1和λn,所以第二个问题暂先搁浅; 分析三:cond(A) = ||A|| * ||A|| =|λ|max * |λ|min,这可以用幂法和反幂法求得,det(A) =λ1 *λ2 * … *-12 λn,这需要求得矩阵A的所有特征值。 由以上分析可知,用幂法和反幂法无法完成所有问题的求解,而用Jacobi法求得矩阵所有特征值后可以求解题目中所给的各个问题。所以该题可以用Jacobi法求解。 ? 模块设计 由Jacobi法的经典4步骤,设计以下模块: 矩阵初始化模块 void init(double *m) 迭代模块 查找非对角线最大模 void find_pq(double *m) 计算旋转角度 void calCosSin(double *m) 更新矩阵 void refreshMextrix(double *m) 计算最终结果模块 void calFinal(double *m) ? 数据结构设计 由于矩阵是对称阵,上下带宽均为2,所以可以考虑用二维数组压缩存储矩阵上半带或下半带。但由于Jacobi法在迭代过程中会破坏矩阵的形态,所以原来为零的元素可能会变为非零,这就导致原来的二维数组无法存储迭代后的矩阵。基于此的考虑,决定采用一维数组存储整个下三角阵,以此保证迭代的正确进行。 完整代码如下(编译环境windows10 + visual studio2010): word 编辑版. 完整代码 // math.cpp : 定义控制台应用程序的入口点。 // 椣据畬敤尠瑳慤硦栮 #includestdio.h #includemath.h #includetime.h #define N 501 #define V (N+1)*N/2+1 #define e 2.718281828459045235360287471352662497757247093699959574966967627724076630353 #define a(i) (1.64 - 0.024 * (i)) * sin(0.2 * (i)) - 0.64 * pow(e , 0.1 / (i)) #define b 0.16 #define c -0.064 #define eps pow((double)10.0,-12) 搣晥湩??楢獴尠ㄥ?昵尠 #define PFrols 5 #define PFe %.11e #define FK 39 int p; int q; double cosz; double sinz; double MAX; int kk; //#define PTS pts word 编辑版. #ifdef PTS void PTS(double *m) { printf(-----------------------------------------------------------------------\n); printf( 迭代第%d次\n,kk); for(int i = 1 ; i = PFrols ; i++) { for( int j = (i-1)*i/2+1 ; j = (i+1)*i/2 ; j++) { printf(PFbits,m[j]); } putchar(10); } for(int i = 1 ; i = PFrols+1 ; i++) { printf( ... ); } putchar(10); printf( . .\n); . .\n); printf( printf( . .\n); fo

文档评论(0)

jinzhuang + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档