期刊文献+

广义块Toeplitz特征值问题的基于sine变换的预处理子(英文)

Block Sine Transform Preconditioner for Generalized Block Toeplitz Eigenvalue Problem
下载PDF
导出
摘要 在求块Toeplitz矩阵束(A_(mn),B_(mn))特征值的Lanczos过程中,通过对移位块Toepltz矩阵A_(mn)-ρB_(mn)进行基于sine变换的块预处理,从而改进了位移块Toeplitz矩阵的谱分布,加速了Lanczos过程的收敛速度.该块预处理方法能通过快速算法有效快速执行.本文证明了预处理后Lanczos过程收敛迅速,并通过实验证明该算法求解大规模矩阵问题尤其有效. We employ the block sine transform-based preconditioner to precondition the shifted block Toeplitz matrix Amn - ρBmn involved in the Lanczes method to compute the minimum eigenvalue of the generalized block Toeplitz eigenvalue problem Anmx=λBmnx, where Amn, and Bmn are partitioned into m×m blocks with order n. The block sine transform-based preconditioner can improve the spectral distribution of the shifted block Toeplitz matrix and, hence, can speed up the convergence rate of the preconditioned Lanczos method. The block sine transform-based preconditioner can be implemented efficiently by the fast transform algorithm. A convergence analysis shows that the preconditioned Lanczos method converges sufficiently fast, and numerical results show that this method is highly effective for large matrix.
出处 《数学研究》 CSCD 2008年第3期240-250,共11页 Journal of Mathematical Study
基金 supported by National Natural Science Foundation of China Nos 10531080
关键词 分块Toeplitz矩阵 双对称 sine变换 预处理Lanczos方法 block Toeplitz matrix quadrantally symmetric sine transform preconditioned Lanczos method
  • 相关文献

参考文献8

  • 1Morgan R, Scott D. Preconditioning the Lanczos algorithms for sparse symmetric eigenvalue problems. SIAM J. Sci. Comput., 1993, 14: 585-593.
  • 2Chan R, Jin X. A family of block preconditioners or block systems. SIAM J. SCI. Comput., 1992, 13: 1218-1235.
  • 3Cybenko G, Van Loan C. Comptuting the minimum eigenvalue of a symmetric positive definte Toeplitz matrix. SIAM J. Sci. Stat. Comput., 1986, 7: 123-131.
  • 4Chan R, Ng M, Wong C. Sine transform based preconditioners for symmetric Toeplitz systems. Linear Algebra Appl., 1996, 232: 237-259.
  • 5Benedetto F. Preconditioning of block Toeplitz matrices bye sine transfoms. SIAM J. SCI. Comput., 1997, 18: 499-515.
  • 6Yip P, Rao K. Fast decimation-in-time algorithms for a family of discrete sine and cosine transforms. Circuits Systems Signal Process., 1984, 3: 387-408.
  • 7Yip P, Rao K. A fast computational algorithm for the discrete sine transform. IEEE Commun. Trans., 1980, 28: 304-307.
  • 8Golub G, Ye Q. An inverse free preconditioned krylov subspace method for symmetric generalized eigenvalue problems. SIAM J. Sci. Comput., 2002, 24: 312-334.

相关作者

内容加载中请稍等...

相关机构

内容加载中请稍等...

相关主题

内容加载中请稍等...

浏览历史

内容加载中请稍等...
;
使用帮助 返回顶部