A vectorized systolic array for block constrained RLS

A new vectorized systolic array for computing the weights of the block recursive least squares (RLS) with p linear constraints is proposed. Using the idea of M-invariant matrices, an upper triangular array for updating R_n of the modified QR decomposition of the data matrix is first derived where the first p rows corresponding to the constraints are in a frozen mode. Next, a lower triangular array for updating R_n^-T is derived where the first p rows are in no operation mode. Connecting these two arrays, the overall near rhombic array is obtained for computing the weights. Also, the numerical stability of the algorithm is examined by simulations