Volume 50, Issue 4
Optimal Parameters for Doubling Algorithms

Tsung-Ming Huang, Ren-Cang Li, Wen-Wei Lin & Linzhang Lu

J. Math. Study, 50 (2017), pp. 339-357.

Published online: 2018-04

Export citation
  • Abstract

In using the structure-preserving algorithm (SDA) [Linear Algebra Appl., 2005, vol. 396, pp. 55–80] to solve a continuous-time algebraic Riccati equation, a parameter-dependent linear fractional transformation $z→(z−\gamma)/(z+\gamma)$ is first performed in order to bring all the eigenvalues of the associated Hamiltonian matrix in the open left half-plane into the open unit disk. The closer the eigenvalues are brought to the origin by the transformation via judiciously selected parameter $\gamma,$ the faster the convergence of the doubling iteration will be later on. As the first goal of this paper, we consider several common regions that contain the eigenvalues of interest and derive the best $\gamma$ so that the images of the regions under the transform are closest to the origin. For our second goal, we investigate the same problem arising in solving an $M$-matrix algebraic Riccati equation by the alternating-directional doubling algorithm (ADDA) [SIAM J. Matrix Anal. Appl., 2012, vol. 33, pp. 170–194] which uses the product of two linear fractional transformations $(z_1,z_2)→[(z_2−\gamma_2)/(z_2+\gamma_1)][(z_1−\gamma_1)/$$(z_1+\gamma_2)]$ that involves two parameters. Illustrative examples are presented to demonstrate the efficiency of our parameter selection strategies.

  • AMS Subject Headings

15A24, 65F30, 65H10

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address

min@ntnu.edu.tw (Tsung-Ming Huang)

rcli@uta.edu (Ren-Cang Li)

wwlin@am.nctu.edu.tw (Wen-Wei Lin)

lzlu@xmu.edu.cn (Linzhang Lu)

  • BibTex
  • RIS
  • TXT
@Article{JMS-50-339, author = {Huang , Tsung-MingLi , Ren-CangLin , Wen-Wei and Lu , Linzhang}, title = {Optimal Parameters for Doubling Algorithms}, journal = {Journal of Mathematical Study}, year = {2018}, volume = {50}, number = {4}, pages = {339--357}, abstract = {

In using the structure-preserving algorithm (SDA) [Linear Algebra Appl., 2005, vol. 396, pp. 55–80] to solve a continuous-time algebraic Riccati equation, a parameter-dependent linear fractional transformation $z→(z−\gamma)/(z+\gamma)$ is first performed in order to bring all the eigenvalues of the associated Hamiltonian matrix in the open left half-plane into the open unit disk. The closer the eigenvalues are brought to the origin by the transformation via judiciously selected parameter $\gamma,$ the faster the convergence of the doubling iteration will be later on. As the first goal of this paper, we consider several common regions that contain the eigenvalues of interest and derive the best $\gamma$ so that the images of the regions under the transform are closest to the origin. For our second goal, we investigate the same problem arising in solving an $M$-matrix algebraic Riccati equation by the alternating-directional doubling algorithm (ADDA) [SIAM J. Matrix Anal. Appl., 2012, vol. 33, pp. 170–194] which uses the product of two linear fractional transformations $(z_1,z_2)→[(z_2−\gamma_2)/(z_2+\gamma_1)][(z_1−\gamma_1)/$$(z_1+\gamma_2)]$ that involves two parameters. Illustrative examples are presented to demonstrate the efficiency of our parameter selection strategies.

}, issn = {2617-8702}, doi = {https://doi.org/10.4208/jms.v50n4.17.04}, url = {http://global-sci.org/intro/article_detail/jms/11322.html} }
TY - JOUR T1 - Optimal Parameters for Doubling Algorithms AU - Huang , Tsung-Ming AU - Li , Ren-Cang AU - Lin , Wen-Wei AU - Lu , Linzhang JO - Journal of Mathematical Study VL - 4 SP - 339 EP - 357 PY - 2018 DA - 2018/04 SN - 50 DO - http://doi.org/10.4208/jms.v50n4.17.04 UR - https://global-sci.org/intro/article_detail/jms/11322.html KW - Doubling algorithm, SDA, ADDA, optimal parameter, algebraic Riccati equation. AB -

In using the structure-preserving algorithm (SDA) [Linear Algebra Appl., 2005, vol. 396, pp. 55–80] to solve a continuous-time algebraic Riccati equation, a parameter-dependent linear fractional transformation $z→(z−\gamma)/(z+\gamma)$ is first performed in order to bring all the eigenvalues of the associated Hamiltonian matrix in the open left half-plane into the open unit disk. The closer the eigenvalues are brought to the origin by the transformation via judiciously selected parameter $\gamma,$ the faster the convergence of the doubling iteration will be later on. As the first goal of this paper, we consider several common regions that contain the eigenvalues of interest and derive the best $\gamma$ so that the images of the regions under the transform are closest to the origin. For our second goal, we investigate the same problem arising in solving an $M$-matrix algebraic Riccati equation by the alternating-directional doubling algorithm (ADDA) [SIAM J. Matrix Anal. Appl., 2012, vol. 33, pp. 170–194] which uses the product of two linear fractional transformations $(z_1,z_2)→[(z_2−\gamma_2)/(z_2+\gamma_1)][(z_1−\gamma_1)/$$(z_1+\gamma_2)]$ that involves two parameters. Illustrative examples are presented to demonstrate the efficiency of our parameter selection strategies.

Tsung-Ming Huang, Ren-Cang Li, Wen-Wei Lin & Linzhang Lu. (2019). Optimal Parameters for Doubling Algorithms. Journal of Mathematical Study. 50 (4). 339-357. doi:10.4208/jms.v50n4.17.04
Copy to clipboard
The citation has been copied to your clipboard