解二次特征值问题的重启广义二阶双正交方式.pdf

…洲Ⅷ川Ⅲ㈣㈣洲0洲川 Y2375500 Contents AbstractinChinese………………………………………………． Abstract…………………………………………………………． iii 1 1 Introduction…．………………………………………………． 2 Second—order andthe second—order Krylovsubspace generalized ………………．．…………．．……… 5 biorthogonalizationprocess The second—order ．…………… 5 §2．1 generalized KryloVsubspace The second—order ．． 7 procedure §2．2 generalized biorthogonalization ofthe 12 §2．3 ImprovementMsOB(m)procedure……………………． 3 methodfor ．………………………… 14 obliqueprojection QEP 4 Restartedmodi6edsecond_order 16 biorthogonalizationprocedure 5 Numerical ………………………………………．． 18 experiments 6 Conclusionsanddiscussions．……………………………．…… 22 ……………………………………………………… 23 Bibliogr印hy ………………………………………………… 25 Acknowledgments 摘要 摘要 我们首先介绍了基于矩阵对A，B以及一对向量u1，乱2的广义二阶Krylov子 空间(‰(A，B；u1，札2))的相关知识，进而将其推广，引入了广义二阶左 间g。(A，B；u1，u2)是由以A，B为系数矩阵u1，u2为初始向量的二阶齐次线性递 推关系得到的向量张成的，它是对由矩阵对A，B以及一个初始向量乱所生成的二 阶Krylov子空间9。(A，B；u)的推广。 在给出了广义二阶左右Krylov子空间定义的基础上我们研究了其与标 准Krylov子空间在解二次特征值问题上的区别与联系，然后提出了一种基于 二阶双正交方法(SOB)f181的修正的二阶双正交的方法(MSOB)，并且对其进 行改进得到一种既节省存储又能减少运算量的修正二阶双正交方法。通过修正 的二阶双正交方法(MSOB)我们得到广义二阶左右Krylov子空问的一组双正交 基。 我们采用斜投影技术f12，171以及重启向量选择技巧[13]进一步提出可重启的修 优点有以下几点： 一、由于其直接应用于二次特征值问题，因此它保持了原有二次特征值问题 在结构及性质上的一些优点(对称性等)； 二、它可以同时求得QEP的左右特征向量； 三、重启方法的应用可以避免SOB方法随着迭代增加所带来的存储量的增加 以及双正交性不能保持所带来的误差。 最