The extended hamiltonian algorithm for the solution of the algebraic riccati equation

We use a second-order learning algorithm for numerically solving a class of the algebraic riccati equations. specifically,the extended hamiltonian algorithm based on manifold of positive definite symmetric matrices is provided. furthermore,this algorithm is compared with the euclidean gradient algorithm,the riemannian gradient algorithm,and the new subspace iteration method. simulation examples show that the convergence speed of the extended hamiltonian algorithm is the fastest one among these algorithms. © 2014 zhikun luo et al.