DiscreteRiccatiSolve

DiscreteRiccatiSolve[{a, b}, {q, r}]

DiscreteRiccatiSolve[{a, b}, {q, r, p}]