解码一帧Layer3第4步：逆量化和重排序 -- requantizer方法
1.逆量化 逆量化的功能是把哈夫曼解码得到的值还原成576个频谱值。长块、短块分别用不同的公式：
其中ISi为第i个哈夫曼码值，XRi为第i个逆量化值。
长块 576个频谱值分为21个增益因子（scale factor）频带（band），由intWidthLong[]保存这21个频带的子带（subband）宽度。intWidthLong[]的初始化见《（六）用JAVA编写MP3解码器——帧数据结构》Layer3.java第203--204行。
短块 576个频谱值分为13个增益因子频带，由intWidthShort[]保存这13个频带的子带宽度。每一个频带包含3个窗。intWidthShort[]的初始化见《（六）用JAVA编写MP3解码器——帧数据结构》Layer3.java第205--206行。
混合块 576个频谱值分为17个增益因子频带。前8个频带为长块，各用一个增益因子逆量化，这8个子带的总和为36；后9个频带为短块。
子带可以理解为一个增益因子逆量化频率线的条数（Subband: Subdivision of the audio frequency band）。
2.重排序 经逆量化后得到的频谱值不是按时间顺序排列的，在编码的MDCT过程中，对于长块产生的频谱值先子带后频率排列；对于短块，产生的频谱值按子带、窗、频率的顺序排列。为了提供哈夫曼编码效率，短窗中的数据被重新排列，按照子带、频率、窗的顺序排列。重排序就是将频率值按时间顺序重新排列，重排序按窗、子带、频带顺序进行排序。
将逆量化和重排序放一块儿同时进行可以减少函数调用，更重要的是可以减少循环，使程序运行效率更高一点。这样带来的不好的一面是代码看起来很复杂，让人理不清头绪。
从逆量化公式可以看出这一模块要进行大量的浮点运算，为减少浮点运算提高运行效率，逆量化采用查表法实现，用到的floatPow2和floatPowIS两张表在构造方法中已经在Layer3的构造方法内初始化过了，公式中用到的其它输入变量的值已经在解码帧边信息这一步初始化过了，以输入变量的值作为索引直接查表，就得到逆量化的结果了；逆量化得到的576个频谱值保存到形参的xr[]中。
class Layer3的requantizer方法源码如下：
//4.
//>>>>REQUANTIZATION & REORDER=============================================
private static float[][][] xr; // [2][32][18];
private static float[] floatPow2; // [256 + 118 + 4]
private static float[] floatPowIS; // [8207]
private static int[] intWidthLong; // [22] 长块的增益因子频带(用一个增益因子逆量化频率线的条数)
private static int[] intWidthShort; // [13] 短块的增益因子频带
private static int rzero_bandL;
private static int[] rzero_bandS = new int[3];
/*
* ISO/IEC 11172-3 ANNEX B,Table 3-B.6. Layer III Preemphasis
*/
private final int[] pretab = {0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,2,2,3,3,3,2,0};
/*
* requantizer方法：逆量化并对短块(纯短块和混合块中的短块)重排序
*
* 在 requantizer 方法内赋值的变量:
* rzero_bandL -- 长块非零哈夫曼值的频带数,用于强度立体声(intensity stereo)处理
* rzero_bandS -- 短块非零哈夫曼值的频带数,用于强度立体声(intensity stereo)处理
* nozero_index -- 非零哈夫曼值的"子带"数
*
* Layer3 逆量化公式(ISO/IEC 11172-3, 2.4.3.4.7.1)
* LONG blocks:
* xr[i] = sign(is[i]) * abs(is[i])^(4/3) * 2^((1/4) * (global_gain - 210)) *
* 2^-((scalefac_scale + 1) / 2 * (scalefac_l[sfb] + preflag * pretab[sfb]))
* SHORT blocks:
* xr[i] = sign(is[i]) * abs(is[i])^(4/3) * 2^((1/4) * (global_gain - 210 - 8 *
* subblock_gain[w])) * 2^-((scalefac_scale + 1) / 2 * scalefac_s[sfb][w])
*/
private void requantizer(int ch, int gr, float[][] xr) {
int bi, wi, band, width, pre, v, is_idx = 0, sb = 0, ss = 0;
float requ_v;
int sb_start, ss_start; // 用于计算短块重排序后的下标
final GRInfo gr_info = objSI.ch[ch].gr[gr];
final int preflag = gr_info.preflag;
final int shift = 1 + gr_info.scalefac_scale;
final int l[] = scfL[ch];
final int s[][] = scfS[ch];
final int[] sub_gain = gr_info.subblock_gain;
sub_gain[0] <<= 3;
sub_gain[1] <<= 3;
sub_gain[2] <<= 3;
int pow2_idx = 256 - gr_info.global_gain;
if( objHeader.isMSStereo())
pow2_idx += 2; // ms_stereo,相当于除以根2
// pure SHORT blocks:
// window_switching_flag = 1, block_type = 2, mixed_block_flag=0
if ((gr_info.window_switching_flag == 1) && (gr_info.block_type == 2)) {
if (gr_info.mixed_block_flag == 1) {
//-------------------------------------------------------------
// 混合块
//-------------------------------------------------------------
rzero_bandS[0] = rzero_bandS[1] = rzero_bandS[2] = 2;
rzero_bandL = -1;
/*
* 混合块的前8个频带是长块。 前8块各用一个增益因子逆量化，这8个增益因子
* 的频带总和为36， 这36条频率线用长块公式逆量化。
*/
for (band = 0; band < 8; band++) {
pre = (preflag == 0) ? 0 : pretab[band];
requ_v = floatPow2[pow2_idx + ((l[band] + pre) << shift)];
width = intWidthLong[band];
for (bi = 0; bi < width; bi++) {
v = is[is_idx];
if (v < 0) {
xr[sb][ss] = -requ_v * floatPowIS[-v];
rzero_bandL = band;
}
else if (v > 0) {
xr[sb][ss] = requ_v * floatPowIS[v];
rzero_bandL = band;
}
else
xr[sb][ss] = 0;
is_idx++;
if (++ss == 18) {
ss = 0;
sb++;
}
}
}
sb_start = 2; ss_start = 0; // 为短块重排序准备好下标
/*
* 混合块的后9个频带是被加窗的短块，其每一块内3个值的增益因子频带相同。
* 后9块增益因子对应的频率子带值为intWidthShort[3..11]
*/
for (band = 3; band < 12; band++) {
width = intWidthShort[band];
for (wi = 0; wi < 3; wi++) {
requ_v = floatPow2[pow2_idx + sub_gain[wi] + (s[wi][band] << shift)];
/*
* 计算当前窗wi内重排序后的下标
*/
sb = sb_start;
ss = ss_start + wi;
if (ss >= 18) {
ss -= 18;
sb++;
}
for (bi = 0; bi < width; bi++) {
v = is[is_idx];
if (v < 0) {
xr[sb][ss] = -requ_v * floatPowIS[-v];
rzero_bandS[wi] = band;
}
else if (v > 0) {
xr[sb][ss] = requ_v * floatPowIS[v];
rzero_bandS[wi] = band;
}
else
xr[sb][ss] = 0;
is_idx++;
ss += 3;
if (ss >= 18) {
ss -= 18;
sb++;
}
}
}
/*
* 调整下标,准备开始在下一增益因子频带内重排序
*/
ss -= 2; //wi=2时已执行过ss += 3
if (ss < 0) {
ss = 0;
sb--;
}
sb_start = sb; ss_start = ss;
}
} else {
//-------------------------------------------------------------
// 纯短块
//-------------------------------------------------------------
rzero_bandS[0] = rzero_bandS[1] = rzero_bandS[2] = rzero_bandL = -1;
sb_start = ss_start = 0;
for (band = 0; band < 12; band++) {
width = intWidthShort[band];
for (wi = 0; wi < 3; wi++) {
requ_v = floatPow2[pow2_idx + sub_gain[wi] + (s[wi][band] << shift)];
sb = sb_start;
ss = ss_start + wi;
if (ss >= 18) {
ss -= 18;
sb++;
}
for (bi = 0; bi < width; bi++) {
v = is[is_idx];
if (v < 0) {
xr[sb][ss] = -requ_v * floatPowIS[-v];
rzero_bandS[wi] = band;
}
else if (v > 0) {
xr[sb][ss] = requ_v * floatPowIS[v];
rzero_bandS[wi] = band;
}
else
xr[sb][ss] = 0;
is_idx++;
ss += 3;
if (ss >= 18) {
ss -= 18;
sb++;
}
}
}
ss -= 2;
if (ss < 0) {
ss = 0;
sb--;
}
sb_start = sb; ss_start = ss;
}
}
rzero_bandS[0]++; rzero_bandS[1]++; rzero_bandS[2]++; rzero_bandL++;
} else {
//-------------------------------------------------------------
// 长块
//-------------------------------------------------------------
rzero_bandL=-1;
for (band = 0; band < 21; band++) {
pre = (preflag == 0) ? 0 : pretab[band];
requ_v = floatPow2[pow2_idx + ((l[band] + pre) << shift)];
width = intWidthLong[band];
for (bi = 0; bi < width; bi++) {
v = is[is_idx];
if (v < 0) {
xr[sb][ss] = -requ_v * floatPowIS[-v];
rzero_bandL = band;
}
else if (v > 0) {
xr[sb][ss] = requ_v * floatPowIS[v];
rzero_bandL = band;
}
else
xr[sb][ss] = 0;
is_idx++;
if (++ss == 18) {
ss = 0;
sb++;
}
}
}
rzero_bandL++;
}
/*
* is_index之后xr[][]清零
*
* iso11172-3上指出应该用"缺省值"逆量化剩余部分,"缺省值"是0还是pow2[pow2_idx]?
* 取pow2[pow2_idx]时有的歌曲有爆音,应该取0才对?
*/
for (; sb < 32; sb++) {
for (; ss < 18; ss++)
xr[sb][ss] = 0;
ss = 0;
}
}
//<<<<REQUANTIZATION & REORDER=============================================
官方文档中没有实现重排序方法的描述，如何进行重排序手册是对我们没有帮助的。我在调试运行时打印出纯短块的576条频率线重排序前后下标的对应关系，想了解重排序细节请先看附件里的pure_SHORT_reorder.txt，琢磨出重排序前后下标的变化规律，再回过来看源码就容易理解其流程了。逆量化究竟如何进行，我这儿就不作讲解了，相信你能从这个.txt文件中琢磨出来。
相关推荐
纯java实现的MP3解码器，装逼使用，面试使用，很好的实现。
java的mp3解码器(非JMF、控制台播放音乐) 用于在你自己的程序中引入,播放MP3文件 运行前先编译---build.bat/build.sh
纯java写的解码器 不含任何jar包
基于通用可编程GPU的视频编解码器——架构、算法与实现
jmp123是一个解码速度很快的纯JAVA的开源的MP3解码器。在2核2.9G CPU上跑，解码速度达到每秒15000帧以上；在4核 CPU上会跑得更快。
C++实现的Mp3解码器
大体流程是从SD卡读取MP3格式的文件，然后提取里边音频数据，经过哈夫曼解码、哈夫曼信息解码、尺度因子解码、然后逆量化、重排序，立体声解码，频率翻转、合成多项滤波，最终生成PCM码流，通过IIS将数据丢给DAC放大...
好不容易找到的，虽然使用matlab实现，但对理解MP3格式很有帮助！
1、好用的mp3播放器，可以试试 2、mp3解码器
H.264的三大编解码器之一JM最新版，为VC下C工程源码，可直接使用或二次开发。 文件中包含说明文档
Huffman编解码问题——讲解.docx
Netty 框架学习 —— 预置的 ChannelHandler 和编解码器（csdn）————程序
JCodec 是一个实现一组流行的视频和音频编解码器的库。目前 JCodec 支持 MP4 视频格式。JCodec 是在 FreeBSD 许可证下分发的自由软件。支持的功能：从视频文件中解码特定的图像帧。性能/质量注意事项 由于 JCodec 是...
许多刚接触mina的朋友，对于mina的编解码器的编写很迷惑.希望这个文档可以帮助朋友们少走弯路。 资源中是一个比较典型的编解码器写法。生成了可执行文件。并对编解码器的代码有详细注释。
Netty 框架学习 —— 编解码器框架（csdn）————程序
用vc6.0编写的MP3播放器和MP3解码器。 单独将MP3解码器编译成静态链接库lib文件， 然后MP3播放器中再调用lib。 很有参考价值。。。
mp3解码库,千千静听,kugou等mp3播放器使用的解码库
java的TLV包的编码和解码，能够对Tag Length Value对应的值进行编码和解码。包内包含jar包和源码包，还有一个java写的带有详细注释的Demo
jmp123 JAVA MP3解码,jmp123 JAVA MP3解码