首页 | 本学科首页   官方微博 | 高级检索  
     


A new two‐phase structure‐preserving doubling algorithm for critically singular M‐matrix algebraic Riccati equations
Authors:Tsung‐Ming Huang  Wei‐Qiang Huang  Ren‐Cang Li  Wen‐Wei Lin
Affiliation:1. Department of Mathematics, National Taiwan Normal University, Taipei, Taiwan;2. Department of Applied Mathematics, National Chiao Tung University, Hsinchu, Taiwan;3. Department of Mathematics, University of Texas at Arlington, Arlington, USA;4. School of Mathematical Science, Xiamen University, Xiamen, Fujian, China
Abstract:Among numerous iterative methods for solving the minimal nonnegative solution of an M‐matrix algebraic Riccati equation, the structure‐preserving doubling algorithm (SDA) stands out owing to its overall efficiency as well as accuracy. SDA is globally convergent and its convergence is quadratic, except for the critical case for which it converges linearly with the linear rate 1/2. In this paper, we first undertake a delineatory convergence analysis that reveals that the approximations by SDA can be decomposed into two components: the stable component that converges quadratically and the rank‐one component that converges linearly with the linear rate 1/2. Our analysis also shows that as soon as the stable component is fully converged, the rank‐one component can be accurately recovered. We then propose an efficient hybrid method, called the two‐phase SDA, for which the SDA iteration is stopped as soon as it is determined that the stable component is fully converged. Therefore, this two‐phase SDA saves those SDA iterative steps that previously have to have for the rank‐one component to be computed accurately, and thus essentially, it can be regarded as a quadratically convergent method. Numerical results confirm our analysis and demonstrate the efficiency of the new two‐phase SDA. Copyright © 2015 John Wiley & Sons, Ltd.
Keywords:M‐matrix algebraic Riccati equation  minimal nonnegative solution  critical case  two‐phase structure‐preserving doubling algorithm  M‐matrix
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号