作者:
Zhang, JuanLi, ShifengXiangtan Univ
Dept Math & Computat Sci Hunan Key Lab Computat & Simulat Sci & Engn Xiangtan 411105 Hunan Peoples R China
In this paper, applying special properties of doubling transformation, a structure-preserving doubling algorithm is developed for computing the positive definite solutions for a nonlinear matrix equation. Further, by ...
详细信息
In this paper, applying special properties of doubling transformation, a structure-preserving doubling algorithm is developed for computing the positive definite solutions for a nonlinear matrix equation. Further, by mathematical induction, we establish the convergence theory of the structure-preserving doubling algorithm. Finally, we offer corresponding numerical examples to illustrate the effectiveness of the derived algorithm. (C) 2020 Elsevier Ltd. All rights reserved.
The discretized Bethe-Salpeter eigenvalue problem arises in the Green's function evaluation in many body physics and quantum chemistry. Discretization leads to a matrix eigenvalue problem for H is an element of C-...
详细信息
The discretized Bethe-Salpeter eigenvalue problem arises in the Green's function evaluation in many body physics and quantum chemistry. Discretization leads to a matrix eigenvalue problem for H is an element of C-2n(x2n) with a Hamiltonian-like structure. After an appropriate transformation of H to a standard symplectic form, the structure-preserving doubling algorithm, originally for algebraic Riccati equations, is extended for the discretized Bethe-Salpeter eigenvalue problem. Potential breakdowns of the algorithm, due to the ill condition or singularity of certain matrices, can be avoided with a double-Cayley transform or a three-recursion remedy. A detailed convergence analysis is conducted for the proposed algorithm, especially on the benign effects of the double-Cayley transform. Numerical results are presented to demonstrate the efficiency and the structure-preserving nature of the algorithm.
A new iterative doubling algorithm for the solution of the discrete time Riccati equation is proposed. The algorithm is based on the Cyclic Reduction Method (CRM). The proposed doubling algorithm does not require non-...
详细信息
ISBN:
(纸本)9781728166957
A new iterative doubling algorithm for the solution of the discrete time Riccati equation is proposed. The algorithm is based on the Cyclic Reduction Method (CRM). The proposed doubling algorithm does not require non-singularity of the transition matrix and is faster than the classical doubling algorithm. The method can be applied to infinite measurement noise case, where the Riccati equation takes the form of the Lyapunov equation. In this case, the classical doubling algorithm is faster.
A highly accurate doubling algorithm to solve the most fundamental quadratic matrix equation in the quasi-birthand-death (QBD) process is developed. It follows from the general framework of the doubling algorithm for ...
详细信息
A highly accurate doubling algorithm to solve the most fundamental quadratic matrix equation in the quasi-birthand-death (QBD) process is developed. It follows from the general framework of the doubling algorithm for the first standard form (SF1) but can be implemented to compute the minimal nonnegative solution with high entrywise relative accuracy for all entries, large or tiny. The algorithm is globally and quadratically convergent, except for QBD equations in the critical case where convergence is linear with the linear rate 1/2. Numerical examples are presented to demonstrate and confirm our claims. The development here parallels the recent work of Xue and Li (2017) [3] on the M-matrix algebraic Riccati equation. (C) 2019 Elsevier Inc. All rights reserved.
We consider computing the minimal nonnegative solution of the nonsymmetric algebraic Riccati equation with *** is well known that such equations can be efficiently solved via the structure-preserving doubling algorith...
详细信息
We consider computing the minimal nonnegative solution of the nonsymmetric algebraic Riccati equation with *** 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 *** 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 ***,the doubling algorithm based on the proposed generalized transformation is presented and shown to be ***,the convergence result and the comparison theorem on convergent rate are *** 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.
作者:
Zhang, JuanLi, ShifengXiangtan Univ
Dept Math & Computat Sci Hunan Key Lab Computat & Simulat Sci & Engn Xiangtan 411105 Hunan Peoples R China
In the paper, we apply a structure-preserving doubling algorithm to solve the continuous coupled algebraic Riccati equation (CCARE). Using the existence and uniqueness of the CCARE, we show that the iteration solution...
详细信息
In the paper, we apply a structure-preserving doubling algorithm to solve the continuous coupled algebraic Riccati equation (CCARE). Using the existence and uniqueness of the CCARE, we show that the iteration solution of the CCARE are positive semi-definite, symmetric, and unique. Further, we discuss the convergent analysis of the structure-preserving doubling algorithm. Moreover, we present two modified structure-preserving doubling algorithms. Finally, we offer corresponding numerical examples to illustrate the effectiveness of the derived numerical algorithms.
We consider the numerical solution of large-scale discrete-time algebraic Riccati equations. The doubling algorithm is adapted, with the iterates for A not computed explicitly but recursively. The resulting algorithm ...
详细信息
We consider the numerical solution of large-scale discrete-time algebraic Riccati equations. The doubling algorithm is adapted, with the iterates for A not computed explicitly but recursively. The resulting algorithm is efficient, with computational complexity and memory requirement proportional to the size of the problem, and essentially converges quadratically. An error analysis, on the truncation of iterates, and some numerical results are presented. (C) 2014 Elsevier B.V. All rights reserved.
Recently, Guo and Lin [SIAM J. Matrix Anal. Appl., 31 (2010), 2784-2801] proposed an efficient numerical method to solve the palindromic quadratic eigenvalue problem (PQEP) ((2)A(T)+Q + A)z = 0 arising from the vibrat...
详细信息
Recently, Guo and Lin [SIAM J. Matrix Anal. Appl., 31 (2010), 2784-2801] proposed an efficient numerical method to solve the palindromic quadratic eigenvalue problem (PQEP) ((2)A(T)+Q + A)z = 0 arising from the vibration analysis of high speed trains, where A, Q is an element of C-nxn have special structures: both Q and A are, among others, m x m block matrices with each block being k x k (thus, n = mk), and moreover, Q is block tridiagonal, and A has only one nonzero block in the (1,m)th block position. The key intermediate step of the method is the computation of the so-called stabilizing solution to the n x n nonlinear matrix equation X + A(T)X(-1)A = Q via the doubling algorithm. The aim of this article is to propose an improvement to this key step through solving a new nonlinear matrix equation having the same form but of only k x k in size. This new and much smaller matrix equation can also be solved by the doubling algorithm. For the same accuracy, it takes the same number of doubling iterations to solve both the larger and the new smaller matrix equations, but each doubling iterative step on the larger equation takes about 4.8 as many flops than the step on the smaller equation. Replacing Guo's and Lin's key intermediate step by our modified one leads to an alternative method for the PQEP. This alternative method is faster, but the improvement in speed is not as dramatic as just for solving the respective nonlinear matrix equations and levels off as m increases. Numerical examples are presented to show the effectiveness of the new method. Copyright (c) 2014 John Wiley & Sons, Ltd.
In this paper, we propose and discuss a new class of complex nonsymmetric algebraic Riccati equations (NAREs) whose four coefficient matrices form a matrix with its omega-comparison matrix being an irreducible singula...
详细信息
In this paper, we propose and discuss a new class of complex nonsymmetric algebraic Riccati equations (NAREs) whose four coefficient matrices form a matrix with its omega-comparison matrix being an irreducible singular M-matrix. We also prove that the extremal solutions of the NAREs exist uniquely in the noncritical case and exist in the critical case. Some good properties of the solutions are also shown. Besides, some classical numerical methods, including the Schur methods, Newton's method, the fixed-point iterative methods and the doubling algorithms, are also applied to solve the NAREs, and the convergence analysis of these methods are given in details. For the doubling algorithms, we also give out the concrete parameter selection strategies. The numerical results show that our methods are efficient for solving the NAREs.
We consider the solution of the large-scale nonsymmetric algebraic Riccati equation XCX - XD - AX + B = 0, with M equivalent to [D, -C;-B, A] is an element of R(n1+n2)x(n1+n2) being a nonsingular M-matrix, and A, D be...
详细信息
We consider the solution of the large-scale nonsymmetric algebraic Riccati equation XCX - XD - AX + B = 0, with M equivalent to [D, -C;-B, A] is an element of R(n1+n2)x(n1+n2) being a nonsingular M-matrix, and A, D being sparse-like, with the products A(-1)u, A(-T)u, D(-1)v and D(-T)v computable in O(n(1)) or O(n(2)) complexity, for some vectors u and v. In the nonsymmetric algebraic Riccati equation arose from a two-dimensional transport model, B, C are low-ranked corrections of some invertible diagonal matrices. The structure-preserving doubling algorithm by Guo, Lin and Xu (2006) is adapted, with the appropriate applications of the Sherman-Morrison-Woodbury formula and the sparse plus-low-rank representations of various iterates. The resulting large-scale doubling algorithm has an O(n) computational complexity and memory requirement per iteration (with n = max{n(1), n(2)}) and converges essentially quadratically, as illustrated by the numerical examples. (c) 2021 Elsevier B.V. All rights reserved.
暂无评论