Iterative solution of nonlinear equations in several variables
Iterative solution of nonlinear equations in several variables
Hi-index | 22.14 |
An effective numerical computation of the steady-state Riccati matrix is based on the successive solutions of a Lyapunov equation using Newton's method. The requirements of this algorithm are an initial stabilizing matrix and the numerical solution of the associated Lyapunov equation. Computationally, the first requirement is the more influencing factor in solving the Riccati equation with reasonable accuracy and speed. In this paper an initial matrix, based on the parameter imbedded solution of the Riccati equation, is introduced for the Newton's algorithm. The imbedding Newton algorithm has been applied to a variety of system, both stable and unstable as well as high-dimensional, A matrices, one of which is reported here. The proposed modification has improved the required CPU time of previous initialization schemes by as much as a factor of 6 times for the same order of accuracy.