Jacobi矩阵特征值并行算法.pdf

第24卷第 1期 纺 织 高 校 基 础 科 学 学 报 Vo1．24，No．1 2011年3月 BASICSCIENCESJOURNALOFTEXTILEUNIVERSITEIS March，2011 文章编号：1006-8341(2011)01-0021-05 Jacobi矩阵特征值的并行算法 刘艳红，吕全义 (西北工业大学 应用数学系，陕西 西安710129) 摘要 ：提 出了并行求解实三对角矩阵特征值方法，该方法主要针对Jacobi矩 阵．应用求多项式根 的Sturm法，将矩阵特征多项式的求根区间隔离成单根区间；对已隔离出的单根区间先用二分法 求解，达到一定精度后再用牛顿法精确求解．考虑到处理机负载平衡问题，将求根 区间分成若干 等分，然后按区间循环地将其分给各个处理机．各处理机并行地进行求根计算，它们之间无通信． 通过此方法实现了处理机负载平衡，算法并行效率达0．85以上．数值算例表明了此并行算法的 高效性． 关键词：Jacobi矩阵；Sturm法；牛顿法；并行算法；并行效率 中图分类号 ：0246 文献标识码 ：A 矩阵特征值问题不仅可直接解决诸如非线性规划、优化、常微分方程等各类数学计算问题 ，而且在结 构力学、工程设计、计算物理和量子力学中具有广泛的应用．由于它具有重要的理论和实际意义，所 以矩阵 特征值问题成为当前国内外高性能计算机的主要计算任务之一．国内外许多学者对矩阵特征值问题进行 了深入研究，其大多可归结于计算三对角矩阵的特征值．目前，三对角矩阵特征值的求解方法主要有Qn 法…、分而治之法 ]、分裂一合并法 3【]、Sturm序列法L4等．文献[1]介绍了带Wilkinson位移的QR法；文献 [2—3]基于分治思想将矩阵做降阶分块处理，然后并行求解 ；文献[4]基于二分法提出了并行求解算法，但 负载平衡问题没有很好解决。关于大型三对角矩阵特征值的并行求解 问题，实现负载平衡仍然面临很大挑 战．本文提出并行求解Jacobi矩阵全部特征值的方法，首先将求根区间分成幻(为正整数，P为处理机台 数)等分 ，然后按区间循环地将其分给各个处理机．各处理机在明确 自己的所有求根区问后，进行并行求 根计算．计算过程中处理机之问不需通信，此方法有效地实现了负载平衡，且达到最小耗时效果． 1 预备知识 定义 1[5 实系数多项式序列f0() () ()，… ()，如果满足条件 (1) ()在(a，6)内没有实根 ； (2)相邻的2个多项式 一() () (1≤k≤m)在 (a，6)内没有公共根 ； (3)如果 是序列中某多项式 () (1≤k≤m一1)的根，则A一()和 +。()异号；则称该多项 式序列为 )在区间(a，b)的Sturm序列． 定义26【 设三对角矩阵 收稿 日期 ：2010-09-25 基金项 目：陕西省 自然科学基金资助项 目(2009JM1008) 通讯作者：吕全义 (1963一)，女，辽宁省沈阳市人 ，西北工业大学教授．E—mail；luquan@nwpu。cdu．cn 纺 织 高 校 基 础 科 学 学 报 第24l卷 61 c1 口l bz 白 A ： … ’一 C l 一 I 如果 ，，都是实数，且 ciD，则称A为Jaeobi矩阵． 设 Jacobi矩阵A的特征多项式为 (入)=I八f—A1．记 (^)为4的前i阶主子矩阵的特征多项式， 则有递推关系 (A)=det(AJ—AI)= (J)L一