A symplectic orthogonal method for single input or single ouput discrete time optimal quadratic control problems

A new, numerically stable, structure preserving method for the discrete linear quadratic control problem with single input or single output is introduced, which is similar to Byers’ method in the continuous case and faster than the general $QZ$-algorithm approach of Pappas, Laub, and Sandell.