We consider computing the minimal nonnegative solution of the nonsymmetric algebraic Riccati equation with M-matrix.It is well known that such equations can be efficiently solved via the structure-preserving doubling ...We consider computing the minimal nonnegative solution of the nonsymmetric algebraic Riccati equation with M-matrix.It is well known that such equations can be efficiently solved via the structure-preserving doubling algorithm(SDA)with the shift-and-shrink transformation or the generalized Cayley transformation.In this paper,we propose a more generalized transformation of which the shift-and-shrink transformation and the generalized Cayley transformation could be viewed as two special cases.Meanwhile,the doubling algorithm based on the proposed generalized transformation is presented and shown to be well-defined.Moreover,the convergence result and the comparison theorem on convergent rate are established.Preliminary numerical experiments show that the doubling algorithm with the generalized transformation is efficient to derive the minimal nonnegative solution of nonsymmetric algebraic Riccati equation with M-matrix.展开更多
基金The work of B.Tang was supported partly by Hunan Provincial Innovation Foundation for Postgraduate(No.CX2016B249)Hunan Provincial Natural Science Foundation of China(No.2018JJ3019)+1 种基金The work of N.Dong was supported partly by the Hunan Provincial Natural Science Foundation of China(Nos.14JJ2114,2017JJ2071)the Excellent Youth Foundation and General Foundation of Hunan Educational Department(Nos.17B071,17C0466).
文摘We consider computing the minimal nonnegative solution of the nonsymmetric algebraic Riccati equation with M-matrix.It is well known that such equations can be efficiently solved via the structure-preserving doubling algorithm(SDA)with the shift-and-shrink transformation or the generalized Cayley transformation.In this paper,we propose a more generalized transformation of which the shift-and-shrink transformation and the generalized Cayley transformation could be viewed as two special cases.Meanwhile,the doubling algorithm based on the proposed generalized transformation is presented and shown to be well-defined.Moreover,the convergence result and the comparison theorem on convergent rate are established.Preliminary numerical experiments show that the doubling algorithm with the generalized transformation is efficient to derive the minimal nonnegative solution of nonsymmetric algebraic Riccati equation with M-matrix.