一种用于音频解码的定点修正离散余弦反变换IMDCT的方法

摘要

本发明公开了一种用于音频解码的定点修正离散余弦反变换的方法，包括以下步骤：步骤一，对Lbit的输入数据xr[N]做IMDCT处理；步骤二，利用Pbit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处理，最后输出Lbit的数据sample[N]和Lbit的更新的修正数据overlap_out[N]，其中L与P均为正整数，且L>P；步骤三，测得所述Lbit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K

权利要求

1.一种用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于，包括以下步骤： 步骤一，对L bit的输入数据xr[N]做IMDCT处理；

步骤二，利用P bit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处理，最后输出L bit的输出数据sample[N]和L bit的更新的修正数据overlap_out[N]，其中L与P均为正整数，且L>P； 步骤三，测得所述L bit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K

2.根据权利要求1所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于：所述方法还包括步骤五，将当前帧的更新的修正数据overlap_out[N]作为下一帧的修正数据overlap_in[N]，重复步骤一；所述K为当前帧的更新的修正数据overlap_out[N]中N个数里最大的溢出比特数；K的取值随帧的变化而变化。

3.根据权利要求2所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于：当K大于等于0小于等于M时，所述overlap_out[N]的定点存储格式在

1.(P-1)～（1+M）.（P-M-1）范围内动态变化，其中M

4.根据权利要求3所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于：所述L的值为32，P的值为16，K大于等于0小于等于3，所述overlap_out[N]的定点存储格式在1.15～4.12范围内动态变化。

5.根据权利要求3所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于：所述L的值为24，P的值为16，K大于等于0小于等于3，所述overlap_out[N]的定点存储格式在1.15～4.12范围内动态变化。

6.根据权利要求3所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于：所述L的值为16，P的值为8，K大于等于0小于等于3，所述overlap_out[N]的定点存储格式在1.7～4.4范围内动态变化。

7.根据权利要求3所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法，其特征在于：所述L的值为24，P的值为8，K大于等于0小于等于3，所述overlap_out[N]的定点存储格式在1.7～4.4范围内动态变化。

说明书

一种用于音频解码的定点修正离散余弦反变换IMDCT的方

法

技术领域

[0001] 本发明属于计算机技术领域，涉及一种用于音频解码的定点修正离散余弦反变换IMDCT的方法。

背景技术

[0002] 修正离散余弦反变换(IMDCT)在音视频信号编码中已得到广泛地应用，但其计算量很大，因此其快速算法在实时编码系统中显得尤为重要。在MP3,WMA,AAC（MP4）等主流音频解码方案中，修正离散余弦反变换（IMDCT）模块也是一种普遍使用的模块。对于其具体的软硬件实现，除了简单按照其原始公式来运算外，通常还有相应的快速算法。为了保证最终音频解码输出的PCM值的精度（通常为16位比特），IMDCT的输入输出及其中间结果一般都使用更高的比特数来存储。譬如将目前较常用的纯24位或者纯32位比特数作为数据的存储位宽，如图1和图2所示。

[0003] 图1中，输入数据xr[N]经过IMDCT处理后，再与修正数据overlap_in[N]经加窗及叠加处理，最后输出数据sample[N]和更新的修正数据overlap_out[N]。在整个运算过程中，所有的数据的存储位宽均为纯24位比特数。

[0004] 图2中，输入数据xr[N]经过IMDCT处理后，再与修正数据overlap_in[N]经加窗及叠加处理，最后输出数据sample[N]和更新的修正数据overlap_out[N]。在整个运算过程中，所有的数据的存储位宽均为纯32位比特数。

[0005] 事实上对于很多软件解码实现来说，由于主流CPU或者DSP的处理单元通常为16位或者32位而没有特别的24位处理单元，为了保证16位精度的PCM输出正常，IMDCT通常以32位进行计算和数据存储。典型的例子如著名的MP3开源软件LIBMAD，即如图2所示。无论是以24位还是32位来存储，现有IMDCT的计算和数据存储方法的一个特点是：其输入、输出以及中间结果都使用统一的比特数，即要么全部是24位，要么全部是32位。其存在的缺陷是：不能灵活利用存储资源，计算量大。

发明内容

[0006] 本发明所要解决的技术问题是：提供一种用于音频解码的定点修正离散余弦反变换IMDCT的方法，该方法能在确保最终PCM输出的高精度的同时降低存储资源的消耗。 [0007] 为解决上述技术问题，本发明采用如下技术方案。

[0008] 一种用于音频解码的定点修正离散余弦反变换IMDCT的方法，包括以下步骤： [0009] 步骤一，对L bit的输入数据xr[N]做IMDCT处理；

[0010] 步骤二，利用P bit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处理，最后输出L bit的输出数据sample[N]和L bit的更新的修正数据overlap_out[N]，其中L与P均为正整数，且L>P；

[0011] 步骤三，测得所述L bit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K

[0012] 步骤四，将L bitoverlap_out[N]按照K缩减为P位存储，定点存储格式为（1+K）.（P-K-1）；

[0013] 步骤五，将当前帧的更新的修正数据overlap_out[N]作为下一帧的修正数据overlap_in[N]，重复步骤一；所述K为当前帧的更新的修正数据overlap_out[N]中N个数里最大的溢出比特数；K的取值随帧的变化而变化。

[0014] 作为本发明的一种优选方案，当K大于等于0小于等于M时，所述定点存储格式在1.(P-1)～（1+M）.（P-M-1）范围内动态变化，其中M

[0015] 作为本发明的另一种优选方案，所述L的值为32，P的值为16，K大于等于0小于等于3，所述定点存储格式在1.15～4.12范围内动态变化。

[0016] 作为本发明的再一种优选方案，所述L的值为24，P的值为16，K大于等于0小于等于3，所述定点存储格式在1.15～4.12范围内动态变化。

[0017] 作为本发明的再一种优选方案，所述L的值为16，P的值为8，K大于等于0小于等于3，所述定点存储格式在1.7～4.4范围内动态变化。

[0018] 作为本发明的再一种优选方案，所述L的值为24，P的值为8，K大于等于0小于等于3，所述定点存储格式在1.7～4.4范围内动态变化。

[0019] 本发明的有益效果在于：本发明所述方法通过将overlap数据由高位宽改为低位宽存储，有效降低了存储要求；并通过检测overlap的最大溢出比特数，让低位宽存储的overlap能够动态使用，尽可能最大的定点化精度，保证了最终PCM输出的高精度。 附图说明

[0020] 图1为现有的纯24位IMDCT的计算和数据存储的方法流程图；

[0021] 图2为现有的纯32位IMDCT的计算和数据存储的方法流程图；

[0022] 图3为实施例一所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法流程图；

[0023] 图4为实施例二所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法流程图；

[0024] 图5为实施例三所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法流程图；

[0025] 图6为实施例四所述的用于音频解码的定点修正离散余弦反变换IMDCT的方法流程图。

具体实施方式

[0026] 在本发明中，overlap_in/overlap_out改由低位宽存储。由于overlap数组的长度跟xr一样，这样节省下来的存储用RAM将会相当可观。为了确保在降低存储要求的同时保证最终PCM输出的高精度，本发明将额外引入2比特的溢出标志K，以保证overlap存储的数据的动态范围得到合理表达。overlap_in和overlap_out在内存中共享同一块区域，即overlap区域既是IMDCT的输入也是IMDCT的输出。IMDCT在计算完成后将更新这一区域。以32位的存储为例，如果采用4.28格式作为定点化格式，即32位中的前4位表示整数部分(包括1比特符号位)，后28位表示小数部分。由于最终PCM输出的值域范围只能在-1到1之间，任何超出部分都只能被截断，所以只要前4位中出现有效值，就表示数据已经溢出。不过由于IMDCT的计算中的中间数据可能需要溢出以保证最终输出的正确性，那么该4位即用来作为溢出保护。

[0027] 当将IMDCT与加窗及叠加操作后的输出的overlap_out[N]由32位改为16位存储时，数据精度的损失将不可避免。其损失的大小，完全取决于16位的定点化格式。显而易见，1.15格式（1比特符号位+15位小数位）的损失最小，其表达的动态范围也最大（15位小数位），然而该格式将无任何比特数表达溢出情况。对于通常的音频编码系统来说，当前帧的overlap_out将作为下一帧的输入即下一帧的overlap_in。所以如果overlap_out中有值出现溢出的情况下，仍然用1.15格式表达的话，将在下一帧计算中出现严重失真。本发明充分考虑到这一情况，并使用K来记录当前帧的溢出情况，K将记录下overlap_out[N]中N个数里最大的溢出比特数。对于4.28格式的32位数据来说，其最大可能的溢出比特数为3，所 以K用2bits表达就已经足够。这样，16位的overlap_out的真正定点化格式是动态的，也即根据其溢出情况，使用尽可能最大化的精度来表达，其表达格式如下： [0028] (1+K).(15-K)

[0029] 这样overlap_out的定点化格式将在1.15～4.12范围内动态变化。下面结合附图对本发明的具体实施方式作进一步详细说明。

[0030] 实施例一

[0031] 本实施例提供一种用于音频解码的定点修正离散余弦反变换IMDCT的方法，如图3所示，包括以下步骤：

[0032] A1、对32bit的输入数据xr[N]做IMDCT处理；

[0033] A2、利用16bit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处理，最后输出32bit的数据sample[N]和32bit的更新的修正数据overlap_out[N]； [0034] A3、测得所述32bit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K<16；

[0035] A4、将32bitoverlap_out[N]按照K缩减为16位存储，定点存储格式为（1+K）.（15-K）；其中overlap_out[N]的定点化格式将在1.15～4.12范围内动态变化。 [0036] 举例来说，如果K值为0，也即N个overlap_out值都没有发生溢出，那么overlap_out[N]将使用1.15格式存储；如果K值为1，那么overlap_out[N]将使用2.14格式存储；如果K值为2，那么overlap_out[N]将使用3.13格式存储；如果K值为3，也即overlap_out里有值发生了最大可能的溢出，那么overlap_out[N]将使用4.12格式存储，以保证下一帧数据的计算正常。例如4.12格式的overlap[N]中发生最大溢出的那个数是"0101.010101010101",由于前4位中最大溢出位为3（加粗部分），那么K=3。 [0037] 实施例二

[0038] 本实施例与实施例一的区别在于，输入数据xr[N]和输出数据sample[N]均为纯24bit，如图4所示，包括以下步骤：

[0039] B1、对24bit的输入数据xr[N]做IMDCT处理；

[0040] B2、利用16bit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处 理，最后输出24bit的数据sample[N]和24bit的更新的修正数据overlap_out[N]；

[0041] B3、测得所述24bit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K<16；

[0042] B4、将24bitoverlap_out[N]按照K缩减为16位存储，定点存储格式为（1+K）.（15-K）；其中overlap_out[N]的定点化格式将在1.15～4.12范围内动态变化。 [0043] 实施例三

[0044] 本实施例与实施例一的区别在于，输入数据xr[N]和输出数据sample[N]均为纯16bit，修正数据overlap_in[N]为8bit，如图5所示，包括以下步骤：

[0045] C1、对16bit的输入数据xr[N]做IMDCT处理；

[0046] C2、利用8bit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处理，最后输出16bit的数据sample[N]和16bit的更新的修正数据overlap_out[N]； [0047] C3、测得所述16bit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K<8；

[0048] C4、将16bitoverlap_out[N]按照K缩减为8位存储，定点存储格式为（1+K）.（7-K）；其中overlap_out[N]的定点化格式将在1.7～4.4范围内动态变化。 [0049] 实施例四

[0050] 本实施例与实施例二的区别在于，修正数据overlap_in[N]位8bit，如图6所示，包括以下步骤：

[0051] D1、对24bit的输入数据xr[N]做IMDCT处理；

[0052] D2、利用8bit的修正数据overlap_in[N]对经IMDCT处理后的数据进行加窗及叠加处理，最后输出24bit的数据sample[N]和24bit的更新的修正数据overlap_out[N]； [0053] D3、测得所述24bit的更新的修正数据overlap_out[N]的N个值中的最大溢出位为K，其中K<8；

[0054] D4、将24bitoverlap_out[N]按照K缩减为8位存储，定点存储格式为（1+K）.（7-K）；其中overlap_out[N]的定点化格式将在1.7～4.4范围内动态变化。 [0055] 本发明针对音频解码中的IMDCT模块，提出了一种用于音频解码的定点修正离散余弦反变换IMDCT的方法，该方法能够降低存储要求且维持高精度输出，该方法具有如下优点：

[0056] 1、通过将overlap数据由高位宽（32位/24位）改为低位宽（16位）存储，有效降低了存储要求，节省了生产成本；

[0057] 2、通过检测overlap的最大溢出比特数，能让低位宽（16位）存储的overlap动态使用尽可能最大的定点化精度，保证了最终PCM输出的高精度要求；

[0058] 3、本发明与实现IMDCT的具体算法和长度（N）无关，也即无论IMDCT采用何种算法，何种长度，本发明都可以在有效降低存储要求的同时维持高精度输出； [0059] 4、本发明充分考虑到当前主流的CPU或者DSP的数据处理位宽，不需要使用非标准的比特位宽。

[0060] 本发明中建议的低存储位宽（16位）并非唯一，可以根据具体方案自由选择位宽，以取得输出精度和存储要求之间的最优平衡，如实施例二至四。

[0061] 这里本发明的描述和应用是说明性的，并非想将本发明的范围限制在上述实施例中。这里所披露的实施例的变形和改变是可能的，对于那些本领域的普通技术人员来说实施例的替换和等效的各种部件是公知的。本领域技术人员应该清楚的是，在不脱离本发明的精神或本质特征的情况下，本发明可以以其他形式、结构、布置、比例，以及用其他元件、材料和部件来实现。