The Rayleigh-Ritz method, refinement and Arnoldi process for periodic matrix pairs

  • Authors:
  • Eric King-Wah Chu;Hung-Yuan Fan;Zhongxiao Jia;Tiexiang Li;Wen-Wei Lin

  • Affiliations:
  • School of Mathematical Sciences, Building 28, Monash University, VIC 3800, Australia;Department of Mathematics, National Taiwan Normal University, Taipei 116, Taiwan;Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China;Department of Mathematics, Southeast University, Nanjing 211189, China;CMMSC and NCTS, Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, Taiwan

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2011

Quantified Score

Hi-index 7.29

Visualization

Abstract

We extend the Rayleigh-Ritz method to the eigen-problem of periodic matrix pairs. Assuming that the deviations of the desired periodic eigenvectors from the corresponding periodic subspaces tend to zero, we show that there exist periodic Ritz values that converge to the desired periodic eigenvalues unconditionally, yet the periodic Ritz vectors may fail to converge. To overcome this potential problem, we minimize residuals formed with periodic Ritz values to produce the refined periodic Ritz vectors, which converge under the same assumption. These results generalize the corresponding well-known ones for Rayleigh-Ritz approximations and their refinement for non-periodic eigen-problems. In addition, we consider a periodic Arnoldi process which is particularly efficient when coupled with the Rayleigh-Ritz method with refinement. The numerical results illustrate that the refinement procedure produces excellent approximations to the original periodic eigenvectors.