In this paper, we develop a matrix framework to solve the problem of finding orthonormal rational function vectors with prescribed poles yk,l ∈ C, with respect to a certain discrete inner product that is defined by a set of data points zi,j ∈ C and corresponding weight vectors ~ wi,j . Our algorithm for solving the problem is recursive, and it is of complexity O(n). If all data points are real ...