Fast computational structures for an efficient implementation of the complete TDAC analysis/synthesis MDCT/MDST filter banks

A new fast computational structure identical both for the forward and backward modified discrete cosine/sine transform (MDCT/MDST) computation is described. It is the result of a systematic construction of a fast algorithm for an efficient implementation of the complete time domain aliasing cancellation (TDAC) analysis/synthesis MDCT/MDST filter banks. It is shown that the same computational structure can be used both for the encoder and the decoder, thus significantly reducing design time and resources. The corresponding generalized signal flow graph is regular and defines new sparse matrix factorizations of the discrete cosine transform of type IV (DCT-IV) and MDCT/MDST matrices. The identical fast MDCT computational structure provides an efficient implementation of the MDCT in MPEG layer III (MP3) audio coding and the Dolby Labs AC-3 codec. All steps to derive the computational structure are described in detail, and to put them into perspective a comprehensive list of references classified into categories is provided covering new research results achieved in the time period 1999-2008 in theoretical and practical developments of TDAC analysis/synthesis MDCT/MDST filter banks (general mathematical, symmetry and special properties, fast MDCT/MDST algorithms and efficient software/hardware implementations of the MDCT in MP3).

[1]  Po-sheng Wu,et al.  Efficient IMDCT core designs for audio signal processing , 2003, 2003 IEEE Workshop on Signal Processing Systems (IEEE Cat. No.03TH8682).

[2]  Gerhard Fettweis,et al.  New Recursive Algorithms for the Unified Forward and Inverse MDCT/MDST , 2003, J. VLSI Signal Process..

[3]  John Princen,et al.  Subband/Transform coding using filter bank designs based on time domain aliasing cancellation , 1987, ICASSP '87. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[4]  Weibei Dou,et al.  An improved audio encoding architecture based on 16-bit fixed-point DSP , 2002, IEEE 2002 International Conference on Communications, Circuits and Systems and West Sino Expositions.

[5]  Wonyong Sung,et al.  Speaking partner: an ARM7-based multimedia handheld device , 2002, IEEE Workshop on Signal Processing Systems.

[6]  Henrique S. Malvar,et al.  Fast algorithms for orthogonal and biorthogonal modulated lapped transforms , 1998, 1998 IEEE Symposium on Advances in Digital Filtering and Signal Processing. Symposium Proceedings (Cat. No.98EX185).

[7]  Peng Liu,et al.  Embedded software optimization for MP3 decoder implemented on RISC core , 2004, IEEE Transactions on Consumer Electronics.

[8]  Suleyman Malki,et al.  A hardware implementation of an MP3 decoder , 2003 .

[9]  Szu-Wei Lee.  Improved algorithm for efficient computation of the forward and backward MDCT in MPEG audio coder , 2001 .

[10]  B. Hunt,et al.  The discreteW transform , 1985 .

[11]  Wonyong Sung,et al.  Software optimization of MPEG audio layer-III for a 32 bit RISC processor , 2002, Asia-Pacific Conference on Circuits and Systems.

[12]  Jar-Ferr Yang,et al.  Recursive architectures for the forward and inverse modified discrete cosine transforms , 2000, 2000 IEEE Workshop on SiGNAL PROCESSING SYSTEMS. SiPS 2000. Design and Implementation (Cat. No.00TH8528).

[13]  Huazhong Shu,et al.  Radix-3 Algorithm for the Fast Computation of Forward and Inverse MDCT , 2007, IEEE Signal Processing Letters.

[14]  Prabindh Sundareson.  A Re-evaluation of Fundamental Transform Structures for Efficient Implementation on Semi-parallel DSP Architectures , 2002 .

[15]  Gerald Schuller,et al.  Improved integer transforms using multi-dimensional lifting [audio coding examples] , 2004, 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[16]  Che-Hong Chen,et al.  Recursive architectures for realizing modified discrete cosine transform and its inverse , 2003, IEEE Trans. Circuits Syst. II Express Briefs.

[17]  Young-Cheol Park,et al.  Design optimization of a dual MP3/AAC decoder , 2002, 2002 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[18]  Vladimir Britanak.  The fast DCT-IV/DST-IV computation via the MDCT , 2003, Signal Process..

[19]  G. Fettweis,et al.  Improved implementation of MDCT in MP3 audio coding , 2004, APCC/MDMC '04. The 2004 Joint Conference of the 10th Asia-Pacific Conference on Communications and the 5th International Symposium on Multi-Dimensional Mobile Communications Proceeding.

[20]  G. Plonka,et al.  Fast and numerically stable algorithms for discrete cosine transforms , 2005 .

[21]  Y. Wang,et al.  Some peculiar properties of the MDCT , 2000, WCC 2000 - ICSP 2000. 2000 5th International Conference on Signal Processing Proceedings. 16th World Computer Congress 2000.

[22]  Heng-Ming Tai,et al.  Design and implementation of a fast algorithm for modulated lapped transform , 2002 .

[23]  Mu-Huo Cheng,et al.  Fast IMDCT and MDCT algorithms - a matrix approach , 2003, IEEE Trans. Signal Process..

[24]  Trieu-Kien Truong,et al.  Fast algorithm for computing the forward and inverse MDCT in MPEG audio coding , 2006, Signal Process..

[25]  K.R. Rao,et al.  An efficient implementation of the forward and inverse MDCT in MPEG audio coding , 2001, IEEE Signal Processing Letters.

[26]  H. Tai,et al.  Fast algorithm for computing modulated lapped transform , 2001 .

[27]  Yin-Tsung Hwang,et al.  A novel MDCT/IMDCT computing kernel design , 2005, IEEE Workshop on Signal Processing Systems Design and Implementation, 2005..

[28]  Fuhuei Lin,et al.  A unified computing kernel for MDCT/IMDCT in modern audio coding standards , 2007, 2007 International Symposium on Communications and Information Technologies.

[29]  Jar-Ferr Yang,et al.  Unified selectable fixed-coefficient recursive structures for computing DCT, IMDCT and subband synthesis filtering , 2004, 2004 IEEE International Symposium on Circuits and Systems (IEEE Cat. No.04CH37512).

[30]  Young-Cheol Park,et al.  Design optimization of main-profile MPEG-2 AAC decoder , 2001, 2001 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (Cat. No.01CH37221).

[31]  Ralf Geiger,et al.  Audio Coding based on Integer Transforms , 2001 .

[32]  Steven G. Johnson,et al.  Type-IV DCT, DST, and MDCT algorithms with reduced numbers of arithmetic operations , 2007, Signal Process..

[33]  Gerald Schuller,et al.  Improved integer transforms for lossless audio coding , 2003, The Thrity-Seventh Asilomar Conference on Signals, Systems & Computers, 2003.

[34]  Y.-C. Yang,et al.  Low power and cost effective VLSI design for an MP3 audio decoder using an optimised synthesis-subband approach , 2004 .

[35]  Miikka Vilermo,et al.  Modified Discrete Cosine Transform: Its Implications for Audio Coding and Error Concealment , 2003 .

[36]  A. Spanias,et al.  Perceptual coding of digital audio , 2000, Proceedings of the IEEE.

[37]  Jürgen Herre,et al.  MPEG-4 High-Efficiency AAC Coding , 2007 .

[38]  Henrique S. Malvar,et al.  Fast algorithm for the modulated complex lapped transform , 2003, IEEE Signal Processing Letters.

[39]  Peng Liu,et al.  A hardware/software co-optimization approach for embedded software of MP3 decoder , 2007 .

[40]  Vladimir Britanak.  An efficient computing of oddly stacked MDCT/MDST via evenly stacked MDCT/MDST and vice versa , 2005, Signal Process..

[41]  Chi-Min Liu,et al.  A Unified Fast Algorithm for Cosine Modulated Filter Banks in Current Audio Coding Standards , 1999 .

[42]  M. Dietz,et al.  MPEG-4 high-efficiency AAC coding [Standards in a Nutshell] , 2008 .

[43]  Vladimir Britanak.  A Note on the MDCT/MDST and Pseudoinverse Matrix , 2004, Comput. Artif. Intell..

[44]  Henrique S. Malvar,et al.  Signal processing with lapped transforms , 1992 .

[45]  Yuriy A. Reznik,et al.  Efficient implementation of a class of MDCT/IMDCT filterbanks for speech and audio coding applications , 2008, 2008 IEEE International Conference on Acoustics, Speech and Signal Processing.

[46]  A. W. Johnson,et al.  Adaptive transform coding incorporating Time Domain Aliasing Cancellation , 1987, Speech Commun..

[47]  Gerhard Fettweis,et al.  Computation of forward and inverse MDCT using Clenshaw's recurrence formula , 2003, IEEE Trans. Signal Process..

[48]  K. Rao,et al.  Discrete Cosine and Sine Transforms: General Properties, Fast Algorithms and Integer Approximations , 2006 .

[49]  Xinjian Chen,et al.  New algorithm for modulated complex lapped transform with symmetrical window function , 2004, IEEE Signal Processing Letters.

[50]  N. Rama Murthy,et al.  A parallel/pipelined algorithm for the computation of MDCT and IMDCT , 2003, Proceedings of the 2003 International Symposium on Circuits and Systems, 2003. ISCAS '03..

[51]  Vladimir Britanak,et al.  A new fast algorithm for the unified forward and inverse MDCT/MDST computation , 2002, Signal Process..

[52]  V. Gurkhe.  Optimization of an MP3 decoder on the ARM processor , 2003, TENCON 2003. Conference on Convergent Technologies for Asia-Pacific Region.

[53]  T.W. Fox,et al.  Goertzel implementations of the forward and inverse modified discrete cosine transform , 2004, Canadian Conference on Electrical and Computer Engineering 2004 (IEEE Cat. No.04CH37513).

[54]  J. Yang,et al.  Regular implementation algorithms of time domain aliasing cancellation , 1996 .

[55]  Heng-Ming Tai,et al.  Design and Efficient Implementation of a Modulated Complex Lapped Transform Processor Using Pipelining Technique , 2001 .

[56]  Xin Yang,et al.  Tradeoffs in modified discrete cosine transform implementations , 2001, ASICON 2001. 2001 4th International Conference on ASIC Proceedings (Cat. No.01TH8549).

[57]  Hi-Seok Kim,et al.  A new optimized algorithm for computation of MDCT and its inverse transform , 2004, Proceedings of 2004 International Symposium on Intelligent Signal Processing and Communication Systems, 2004. ISPACS 2004..

[58]  Ajay Kumar Bangla,et al.  A new efficient implementation of TDAC synthesis filterbank based on Radix-2 FFT , 2006, 2006 14th European Signal Processing Conference.

[59]  Gerhard Fettweis,et al.  New recursive algorithms for the forward and inverse MDCT , 2001, 2001 IEEE Workshop on Signal Processing Systems. SiPS 2001. Design and Implementation (Cat. No.01TH8578).

[60]  Hi-Seok Kim,et al.  An optimized algorithm for computing the modified discrete cosine transform and its inverse transform , 2004, 2004 IEEE Region 10 Conference TENCON 2004..

[61]  Henrique S. Malvar.  A modulated complex lapped transform and its applications to audio processing , 1999, 1999 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings. ICASSP99 (Cat. No.99CH36258).

[62]  F. R. Gantmakher.  The Theory of Matrices , 1984 .

[63]  Lotfi Senhadji,et al.  A fast algorithm for the computation of 2-D forward and inverse MDCT , 2008, Signal Process..

[64]  Ted Painter,et al.  Audio Signal Processing and Coding , 2007 .

[65]  Marina Bosi,et al.  Introduction to Digital Audio Coding and Standards , 2004, J. Electronic Imaging.

[66]  Gene H. Golub,et al.  Matrix computations , 1983 .

[67]  L. Yaroslavsky,et al.  On the relationship between MDCT, SDPT and DFT , 2000, WCC 2000 - ICSP 2000. 2000 5th International Conference on Signal Processing Proceedings. 16th World Computer Congress 2000.

[68]  Xinjian Chen,et al.  A Novel DCT-Based Algorithm for Computing the Modulated Complex Lapped Transform , 2006, IEEE Transactions on Signal Processing.

[69]  Peter Lancaster,et al.  The theory of matrices , 1969 .