一种有效的MDCT 快速算法.docVIP

一种有效的MDCT快速算法 1)祝平平刘建国多谱信息处理技术国防科技重点实验室图像信息处理与智能控制518129） 摘 要: 本文给出了一种有效的改进型离散余弦变换(MDCT)的快速算法。本算法在Britanak[1]提出的用离散余弦变换(DCT)和离散正弦变换(DCT)来计算改进型离散余弦变换(MDCT)的算法思想上，进一步改进算法流图，改变了扩展因子的位子，减少了乘法运算的次数，从而提高了算法的计算度。本算法能很好的计算N≠2n的MDCT,并在G.711宽带扩展层的编解码器中得到运用。 关键词: 改进型离散余弦变换(MDCT)；离散余弦变换(DCT)；离散正弦变换(DST) Efficient Fast Algorithm of the Forward MDCT 1)ZHU Ping-ping，1)PAN Chao，1)LIU Zhen-bing，2)Hu Chen 1)（Institute forPattern Recognition and Artificial Intelligence, Huazhong University of Science and Technology, State Key Lab for Multi-spectral Information Processing Technologies, Key Lab of State Education Ministry for Image Processing and Intelligent Control, Wuhan 430074） 2) （Huawei Technologies CO.,LTD，Shenzhen 518129） Abstract: An efficient fast algorithm of the forward MDCT is presented. We improve the flow graph of Britanak’s algorithm which employs discrete cosine transform of type II (DCT-II) and discrete sine transform of type II (DST-II) to compute the forward MDCT. In order to obtain more precision, we change the location of the scaling factors and reduce the amount of multiplication operators. This algorithm is suitable for N≠2 MDCT, and is used in CODEC in G.711 wide-brand scaling lay. Keywords: Modified Discrete Cosine Transform (MDCT)；； 1 引言 MDCT算法的推导 2.1 MDCT算法的定义和分析 对于G.711 宽带扩展层中MDCT的定义： (1) 其中 ， 是满足时域混叠消除的窗函数。在 G.711 宽带扩展层中要求 。 主要是计算核： (2) 本文将只讨论对(2)式的实现。 [1]介绍了一种用DCT和DST来计算2N点MDCT的算法。先将2N点MDCT转化成N点DCT-IV计算，再引入扩展因子将DCT-IV转化成N/2点的DCT-II和N/2点的DST-II计算。但是，由于[1]算法引入的扩展因子组成了蝶形计算单元，如果计算N=40的MDCT的情况蝶形计算单元就会过于复杂，而且这么多的乘法运算将很大程度上提高算法的计算复杂度并降低算法的计算精度。[2]也是将MDCT转化成DCT来计算。MDCT系数在经过扩展因子处理后可用DCT-II计算。算法[1][2]的共同特点是数据在进行DCT计算前都乘以了扩展系数，而由于乘法运算会损失一定的末尾数据信息，使得后续的DCT和DST的计算精度明显下降。 为此，我们提出的算法将乘法运算移到DST变换之后，这一点主要在图3中得到体现。由于我们的扩展因子较少的影响了运算，本算法会得到更高的计算精度。 2.2 MDCT算法的实现 图1给出了整个算法的计算思想流图，我们给出了将2N点MDCT转换为N/2个点的DST-III和N/2个点的DST-IV的算法。而根据这种算法思想可以不断分解直至少数个奇数点的DST-III和DST-IV计算。下面具体给出算法的推导。 图1. 整个算法思想计算流图 2N个点的MDCT预处理得到N个点的DCT-IV。： (3) 那么： (4) 上式就是个点的DCT-IV。可以用流图表示为： 图2.MDCT计