Implicit restart scheme for large scale Krylov subspace model reduction method

An implicitly restarted Krylov subspace model reduction method is presented. By this technique a stable, linear transfer function f(s) of order n can be approximated by an order m, where n≫m. The oblique projection of stable systems onto a Krylov subspace may generate unstable partial realizations along with undesirable modes. To ensure that the Krylov projected model for a stable system is also stable and contains useful modes only, further operations like stable projection and balanced truncation are then performed to compute a reduced order model f_r(s) of order r (where r\n\n\t\t