Rapid update of odd DCT and DST for real-time signal processing

Abstract

When processing a signal or an image using the Discrete Cosine Transform (DCT) or Discrete Sine Transform (DST), a typical approach is to extract a portion of the signal by windowing and then form the DCT or DST of the window contents. By shifting the window point by point over the signal, the entire signal may be processed. DCTs and DSTs are defined where the denominator in the transform kernel is either an odd or an even integer, resulting in transforms known as the even DCT (EDCT), even DST (EDST), odd DCT (ODCT) and odd DST (ODST). Each is available in types I to IV, for a total of 16 different transforms. The widely used transform commonly called the "DCT" is actually the EDCT-II. In this paper we extend our previous work using the EDCT-II and EDST-II, and show that a similar approach yields algorithms for the ODCT-II and ODST-II. We develop algorithms to "update" the ODCT-II and ODST-II simultaneously to reflect the modified window contents using less computation than directly evaluating the modified transform via standard Fast Transform algorithms. These algorithms are able to handle arbitrary step sizes up to the length of the transform, i.e. the algorithm simultaneously updates the ODCT-II and ODST-II to reflect inclusion of r, where 1 ≤ r ≤ N-1, additional data points and removal of r old points from the signal. Examples of applications where this algorithm would be useful include target recognition where time constraints may not permit the immediate processing of every incoming data point, adaptive system identification, etc.