特征值问题的预变换方法(I)_ 杨辉三角阵变换与二阶PDE 特征多项式.pdfVIP

中国科学: 数学 2011 年 第41 卷 第8 期: 701 724 特征值问题的预变换方法 杨辉三角阵变换 与二阶 特征多项式 孙家昶 中国科学院软件研究所并行计算实验室, 北京 100190 E-mail: sun@mail.rdcps.ac.cn 收稿日期: 2011-01-27; 接受日期: 2011-07-29 国家自然科学基金(批准号: 资助项目 摘要 本文提出一类求解特征值问题的下三角预变换方法, 目标是通过相似变换后矩阵下三角元素平 方和明显减少、且变换后的特征值及其特征向量较易求解, 使变换后的对角线可作为全体特征值很好 的一组初值, 其作用如同对于解方程组找到好的预条件子, 加速迭代收敛. 以二阶PDE 数值计算为例, 对于以Laplace 方程为代表的特征波向量组及正交多项式组有广泛的应用前景. 杨辉三角是我国古代数学家的一项重要成就. 本文引入杨辉三角矩阵作为预变换子, 给出一般矩 阵用杨辉三角矩阵作为左、右预变换子时变为上三角矩阵的充要条件, 给出了元素为行指标二次多项 式的两个矩阵类 (三对角线阵与五对角线阵) 中特征值何时保持二次多项式的充要条件, 并应用于构 造新的二元PDE 正交多项式. 关键词 特征问题预变换 二阶PDE 特征多项式 杨辉三角矩阵 主题分类 65F15, 65D99, 65N25, 15A18, 35P99 引言 众所周知, 方程组与特征值问题是科学计算研究与解法器的两大重要内容. 求解线性方程组的预 [1, 2] 条件子方法自上世纪七十年代末引入以来 , 已在算法效率与大规模实际计算中取得巨大的进展. 其 基本思想是把原始的病态方程 Au = f 变为另一个等价方程组 BAu = Bf 的求解, 以期改善条件数 [3] κ(BA) ≪ κ(A) 加速迭代. 预条件子的思想也正推广到某些稀疏非线性方程组求解 . 2010 年 11 月, 作者所在研究所的青年学者杨超在国产千万亿次计算机 “天河一号” 上测试 “地球外核热流动模拟” [4] 时主要用区域分解预条子已解出百亿未知数的方程组 . 如何处理特征值问题则要复杂得多. 特征值问题是一种特殊的非线性问题, 求解方式与相应的矩 阵或PDE 性质密切相关. 事实上, 基于第一原理, 不少物理学家想了很多办法构造特殊基底用以近似 求解 Schr¨odinger 方程特征值问题. 如我国物理学家王鼎盛等应用自洽和全电子的线性化缀加平面波 [5] 方法, 在密度泛函近似和间隙区完整势函数近似下, 计算了几十种元素晶体的结合能 , 并与作者合作 [6, 7] 在国产机上实现了特征值并行计算 . 国际上计算少量较小特征值有效的Jacobi-Davidson 算法首 [8] 先也是计算化学家提出的 . 英文引用格式 Sun J C. On pre-transformed methods for eigen-problems, I: Yanghui-triangle transform and 2nd order PDE eigen-problems (in Chinese). Sci Sin Math, 2011, 41(8): 701–724, doi: 10.1360/012011-76 孙家昶: 特征值问题的预变换方法 (I): 杨辉三角阵变换与二阶 PDE 特征多项式 在数学上, 类似于求解线性方程组的预条件子方法, 一个特征值问题 Au = λu, A = (a )