Abstract: A new algebraic approach to approximation of correlation matrices of discrete signals based on a special family of square matrices with sparse circulant alternation of elements is developed.