收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

分布式存储环境下矩阵广义特征值问题的并行计算

魏立峰  
【摘要】: 矩阵广义特征值问题是当前迅速发展的计算机科学和数值代数中的一个非常活跃的研究课题。它在很多应用中扮演非常重要的角色,从数学角度来看,矩阵特征值问题的应用大多来自数学物理方程、差分方程、Markov过程等。目的是为了计算固体、流体、电磁、微观粒子、系统控制等重大问题。在实际的科学研究与工程应用中,比如在建筑工程、航空航天研究、生物科学、计算物理以及石油勘探中,都要涉及到大规模矩阵广义特征值问题的计算。 大规模并行处理系统(MPP)和PC机群为并行求解矩阵广义特征值问题提供了分布式存储环境。本文从理论和实验两方面深入研究了分布式环境下实矩阵广义特征值问题Ax=λBx的并行计算,提出了一些新的算法。本文的主要研究成果概括如下: (1)关于实对称定三对角矩阵广义特征值问题 ⅰ) 提出了一种基于等规模矩阵划分策略的分治算法,子问题的规模之和等于原问题的规模。该算法利用割线法迭代计算广义特征值,具有计算量少,计算速度快的优点,并且算法具有内在并行性。理论分析和数值实验均证明该算法比K.Y.Li算法快,在矩阵规模较大时,计算时间约可以减少10%-20%。当n→∞时,加速比S_p→p,在分布式环境下可以获得很好的线性加速比和稳定的效率。 ⅱ) 提出了一种具有全新矩阵划分策略的分治算法。在这个划分策略中,子问题的规模之和等于原问题的规模减1,文中给出了特征值分割定理及其证明。该算法所提供的特征区间要优于第一种算法,因而在采用相同的迭代方式时,算法具有计算量少、计算时间短的优点。在矩阵规模较大时,计算时间约可以减少30%甚至更多。同时,算法也具有良好的内在并行性,当n→∞时,加速比S_p→p,在分布式环境下可以获得很好的线性加速比和稳定的效率。 (2)关于实对称带状矩阵广义特征值问题 ⅰ) 提出了一种结合多分法的并行分治算法,给出了特征值分割定理及其证明。在该算法中,子问题的规模之和等于原问题的规模。算法利用二分法结合广义Rayleigh商迭代,计算包含在特征区间内的每个特征对。在计算全部特征值和特征向量的情况下,算法的计算复杂性为O(n~2r~2),当r<<n时,优于 坛一一一~一续攀遨粼望廷二生掣夔 LApACK的计算复杂性O(n’)。在并行实现时,算法结合并行性好的多分法, 在分布式环境下获得了很高的并行加速比和效率。 ii)提出了另一种结合多分法的并行分治算法,给出了特征值分割定理及其证明。 在该算法中,子问题的规模之和等于原问题的规模减;,其中;为半带宽。在 计算全部特征值和特征向量的情况下,算法的计算复杂性为。(。’:’),并且该 算法的计算时间绝大多数情况都少于前一种分治算法。结合多分法,算法具 有很好的并行性,在分布式环境下,获得了很高的并行加速比和效率。 iii)提出了一种基于分治策略的并行同伦连续方法。通过限制步长、抛弃特征值 束的方法,提高算法效率。在矩阵规模较大时、算法的计算时间明显比 LAPACK少。在算法的并行实现时,充分利用计算与通信的重叠,在分布式 环境下获得了很高的并行加速比和稳定的效率。 (3)提出了一种计算实对称定矩阵广义特征值问题的并行分治算法。该算法是在计 算矩阵B的谱分解和MDR约化的基础上,将问题转化为对称三对角一正定对角 矩阵广义特征值问题,然后利用分治算法结合割线法迭代进行计算。根据扰动 理论,给出了每个特征值的包含区间。数值实验表明,在分布式环境下,该算 法可以获得很高的并行加速比和效率。 (4)提出了一种计算非对称矩阵广义特征值问题的分治算法。该算法在分治策略的 基础上,利用Laguerre迭代计算广义特征值,并采用各种方法尽可能减少特征 值的遗漏。采用收缩方法计算遗漏的特征值。数值实验表明,该算法计算速度 快,并行效果好。 (5)在分布式环境下实现了本文提出的所有算法,并形成了一个符合规范的软件包。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 刘涛;周明全;武仲科;;针对颅面点云数据的Deluney三角网构建方法[J];计算机工程与应用;2011年22期
2 高振国;蔡绍滨;李香;赵蕴龙;王慧强;;非完全无线报文重传分析研究及其分治策略[J];电子学报;2011年07期
3 黄勇;;稀疏保留判决分析在人脸表情识别中的应用[J];计算机工程;2011年14期
4 ;[J];;年期
5 ;[J];;年期
6 ;[J];;年期
7 ;[J];;年期
8 ;[J];;年期
9 ;[J];;年期
10 ;[J];;年期
11 ;[J];;年期
12 ;[J];;年期
13 ;[J];;年期
14 ;[J];;年期
15 ;[J];;年期
16 ;[J];;年期
17 ;[J];;年期
18 ;[J];;年期
19 ;[J];;年期
20 ;[J];;年期
中国重要会议论文全文数据库 前10条
1 忻欣;冯纯伯;;奇异H~∞控制问题二次矩阵不等式的可解判据的研究[A];1994中国控制与决策学术年会论文集[C];1994年
2 高鹰;谢胜利;姚振坚;;一种新的基于广义特征值分解的盲信号分离算法[A];第十二届全国信号处理学术年会(CCSP-2005)论文集[C];2005年
3 高尚;;基于分治算法的选择问题研究[A];2005年全国理论计算机科学学术年会论文集[C];2005年
4 官玉才;陈璞;;子空间迭代法的一种新改进方案[A];计算爆炸力学进展[C];2006年
5 樊久铭;何晓峰;;基于模态区间方法的动力特征值算法研究[A];2010振动与噪声测试峰会论文集[C];2010年
6 钟慧湘;庞云阶;冯月萍;;估计基本矩阵的一个新方法[A];中国仪器仪表学会第六届青年学术会议论文集[C];2004年
7 刘东卓;;改进的小波神经网络TLS-ESPRIT算法[A];2008通信理论与技术新发展——第十三届全国青年通信学术会议论文集(下)[C];2008年
8 谭茂金;张庚骥;石耀霖;;阵列电法测井的垂直模式匹配理论研究[A];中国地球物理学会第二十四届年会论文集[C];2008年
9 陈青青;常卓;陆歌皓;;基于改进Delaunay三角的Robocup机器人策略站位[A];2009年中国智能自动化会议论文集(第二分册)[C];2009年
10 姜源源;易丽清;魏泳涛;冯国英;;正方形多芯光子晶体光纤的有限元分析[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
中国博士学位论文全文数据库 前10条
1 魏立峰;分布式存储环境下矩阵广义特征值问题的并行计算[D];中国人民解放军国防科学技术大学;2002年
2 王顺绪;特征值问题的并行计算[D];南京航空航天大学;2008年
3 王元媛;求解Toeplitz矩阵束广义特征值问题的预处理方法[D];厦门大学;2008年
4 钟慧湘;基本矩阵计算方法的研究[D];吉林大学;2005年
5 孙立双;矿体三维建模及储量计算关键问题研究[D];东北大学;2008年
6 黄永辉;基于实时计算方法的结构动力快速分析方法及其应用研究[D];湖南大学;2008年
7 徐忠海;结构特征灵敏度分析若干问题研究[D];吉林大学;2008年
8 杨锋;自适应广义特征分解及在通信中的应用[D];中国科学技术大学;2007年
9 管子玉;基于图学习的Web信息检索技术研究[D];浙江大学;2010年
10 杨其雨;基于非标航片的三维数字地图生成技术研究[D];上海大学;2008年
中国硕士学位论文全文数据库 前10条
1 聂永明;对称广义特征值问题的数值求解算法的研究[D];南昌大学;2012年
2 陆昊辉;广义特征值问题中的预处理方法[D];吉林大学;2005年
3 刘丹;关于Wilkinson定理的一些推广[D];中国海洋大学;2008年
4 张峰超;一种基于LOBPCG的改进算法[D];复旦大学;2010年
5 邢五丰;基于归纳学习的结构损伤识别方法研究[D];武汉理工大学;2006年
6 王正国;计算几个低阶特征对的动力缩聚法的分析和探讨[D];南京航空航天大学;2006年
7 刘红;同伦方法求解广义特征值问题[D];新疆大学;2008年
8 郝玉珏;非对称广义特征值问题及其在可坍塌槽道流研究中的应用[D];天津大学;2007年
9 安方;固定界面模态综合法的理论分析[D];南京航空航天大学;2008年
10 成斐鸣;基于ICA的扫描图像去噪方法研究[D];西安理工大学;2008年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978