开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > mpeg/mp3 > 查看源码
MP3_Decode.c
资源名称：mp3ruanjianjiema.rar　[点击查看]
上传用户：kingbiz
上传日期：2022-06-24
资源大小：2524k
文件大小：22k
源码类别：

mpeg/mp3

开发平台：

C/C++

MP3_Decode.c：源码内容
  4. #include "mp3_common.h"
  5. #include "mp3_decode.h"
  6.  
  7. /*============================================================================*/
  9. const struct BandIndex sfBandIndex[3]=
  10. {
  11. {{0,4,8,12,16,20,24,30,36,44,52,62,74,90,110,134,162,196,238,288,342,418,576},
  12. {0,4,8,12,16,22,30,40,52,66,84,106,136,192}},
  13. {{0,4,8,12,16,20,24,30,36,42,50,60,72,88,106,128,156,190,230,276,330,384,576},
  14. {0,4,8,12,16,22,28,38,50,64,80,100,126,192}},
  15. {{0,4,8,12,16,20,24,30,36,44,54,66,82,102,126,156,194,240,296,364,448,550,576},
  16. {0,4,8,12,16,22,30,42,58,78,104,138,180,192}}
  17. };
  19. /*============================================================================*/
  21. const static int bitrate[3][15] = {
  22. {0,32,64,96,128,160,192,224,256,288,320,352,384,416,448},
  23. {0,32,48,56,64,80,96,112,128,160,192,224,256,320,384},
  24. {0,32,40,48,56,64,80,96,112,128,160,192,224,256,320}
  25. };
  27. /*============================================================================*/
  28. /*============================================================================*/
  30. void decode_info()
  31. {
  32.     layer *hdr = fr_ps.header;
  33.     int x;
  35.     hdr->version = getbit(1);
  36.     hdr->lay = 4-getbit(2);
  37.     hdr->error_protection = !getbit(1); /* error protect. TRUE/FALSE */
  38.     hdr->bitrate_index = getbit(4);
  39.     hdr->sampling_frequency = getbit(2);
  40.     hdr->padding = getbit(1);
  41.     hdr->extension = getbit(1);
  42.     hdr->mode = getbit(2);
  43.     hdr->mode_ext = getbit(2);
  44.     hdr->copyright = getbit(1);
  45.     hdr->original = getbit(1);
  46.     hdr->emphasis = getbit(2);
  47. }
  49. /*============================================================================*/
  50. //攻取Side信息
  52. void III_get_side_info(III_side_info_t *si)
  53. {
  54. int ch, gr, i;
  55. int stereo = fr_ps.stereo;
  57. si->main_data_begin = getbit(9);  //SI 边信息结构
  58. if (stereo == 1)
  59. si->private_bits = getbit(5);
  60. else
  61. si->private_bits = getbit(3);
  63. for (ch=0; ch<stereo; ch++)         //第二颗粒是否传送信息
  64. for (i=0; i<4; i++)
  65. si->scfsi[ch][i] = getbit(1);
  67. for (gr=0; gr<2; gr++) {
  68. for (ch=0; ch<stereo; ch++) {
  69. grle[ch][gr].part2_3_length = getbit(12);
  70. grle[ch][gr].big_values = getbit(9);
  71. grle[ch][gr].global_gain = getbit(8);
  72. grle[ch][gr].scalefac_compress = getbit(4);
  73. grle[ch][gr].window_switching_flag = getbit(1);
  74. if (grle[ch][gr].window_switching_flag) {            //窗开关为1
  75. grle[ch][gr].block_type = getbit(2);
  76. grle[ch][gr].mixed_block_flag = getbit(1);
  77. for (i=0; i<2; i++)
  78. grle[ch][gr].table_select[i] = getbit(5);
  79. for (i=0; i<3; i++)
  80. grle[ch][gr].subblock_gain[i] = getbit(3);
  82. /* Set region_count parameters since they are implicit in this case. */
  84. if (grle[ch][gr].block_type == 0) {
  85. // printf("Side info bad: block_type == 0 in split block.n");
  86. // exit(0);
  87. }
  88. else if (grle[ch][gr].block_type == 2
  89. && grle[ch][gr].mixed_block_flag == 0)
  90. grle[ch][gr].region0_count = 8; /* MI 9; */     //region0和1的缺省配置
  91. else grle[ch][gr].region0_count = 7; /* MI 8; */
  92. grle[ch][gr].region1_count = 20 - grle[ch][gr].region0_count;
  93. }
  94. else {                                                       //窗开关为0
  95. for (i=0; i<3; i++)
  96. grle[ch][gr].table_select[i] = getbit(5);
  97. grle[ch][gr].region0_count = getbit(4);
  98. grle[ch][gr].region1_count = getbit(3);
  99. grle[ch][gr].block_type = 0;
  100. }
  101. grle[ch][gr].preflag = getbit(1);
  102. grle[ch][gr].scalefac_scale = getbit(1);
  103. grle[ch][gr].count1table_select = getbit(1);
  104.          }
  105. }
  106. }
  108. /*============================================================================*/
  109. /*============================================================================*/
  110. /*============================================================================*/
  111. //获取比例因子
  113. const struct
  114. {
  115. int l[5];
  116. int s[3];
  117. } sfbtable = {
  118. {0, 6, 11, 16, 21},
  119. {0, 6, 12}
  120. };
  122. const int slen[2][16]=
  123. {
  124. {0, 0, 0, 0, 3, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4},
  125. {0, 1, 2, 3, 0, 1, 2, 3, 1, 2, 3, 1, 2, 3, 2, 3}
  126. };
  128. void III_get_scale_factors(III_scalefac_t *scalefac, III_side_info_t *si, int gr, int ch)
  129. {
  130. int sfb, i, window;
  131. struct Granule *gr_info = &(grle[ch][gr]);
  133. if (gr_info->window_switching_flag && (gr_info->block_type == 2)) {
  134. if (gr_info->mixed_block_flag) { /* MIXED */ /* NEW - ag 11/25 */
  135. for (sfb = 0; sfb < 8; sfb++)
  136. (*scalefac)[ch].l[sfb] = hgetbits(
  137. slen[0][gr_info->scalefac_compress]);       //用 slen[0]
  138. for (sfb = 3; sfb < 6; sfb++)
  139. for (window=0; window<3; window++)
  140. (*scalefac)[ch].s[window][sfb] = hgetbits(
  141. slen[0][gr_info->scalefac_compress]);   //用 slen[0]
  142. for (sfb = 6; sfb < 12; sfb++)
  143. for (window=0; window<3; window++)
  144. (*scalefac)[ch].s[window][sfb] = hgetbits(
  145. slen[1][gr_info->scalefac_compress]);   //6-11用 slen[1]
  147. for (sfb=12,window=0; window<3; window++)
  148. (*scalefac)[ch].s[window][sfb] = 0;
  149. }
  150. else {  /* SHORT*/
  151. for (i=0; i<2; i++)
  152. for (sfb = sfbtable.s[i]; sfb < sfbtable.s[i+1]; sfb++)
  153. for (window=0; window<3; window++)
  154. (*scalefac)[ch].s[window][sfb] = hgetbits(
  155. slen[i][gr_info->scalefac_compress]);
  156. for (sfb=12,window=0; window<3; window++)
  157. (*scalefac)[ch].s[window][sfb] = 0;
  158. }
  159. }
  160. else {   /* LONG types 0,1,3 */
  161. for (i=0; i<4; i++) {
  162. if ((si->scfsi[ch][i] == 0) || (gr == 0))
  163. for (sfb = sfbtable.l[i]; sfb < sfbtable.l[i+1]; sfb++)
  164. (*scalefac)[ch].l[sfb] = hgetbits(
  165. slen[(i<2)?0:1][gr_info->scalefac_compress]);
  166. }
  167. (*scalefac)[ch].l[22] = 0;
  168. }
  169. }
  170.  
  171. /*============================================================================*/
  172. /*============================================================================*/
  173. //逆量化采样
  175. const int pretab[22] = {0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,2,2,3,3,3,2,0};
  176. static float is43_table[1024];
  178. void III_dequantize_sample(int is[SBLIMIT][SSLIMIT], float xr[SBLIMIT][SSLIMIT], III_scalefac_t *scalefac, struct Granule *gr_info, int ch)
  179. {
  180. int ss,sb,cb=0,sfreq=fr_ps.header->sampling_frequency;
  181. int next_cb_boundary, cb_begin, cb_width, sign;
  182. float temp;
  183. static int init=1;
  185. if(init)
  186. {
  187. for(init=0;init<1024;init++)
  188. {
  189.   is43_table[init]=1.33333333333*log(init);
  190. }
  191. init=0;
  192. }
  194. /* choose correct scalefactor band per block type, initalize boundary */
  196. if (gr_info->window_switching_flag && (gr_info->block_type == 2) )
  197. if (gr_info->mixed_block_flag)
  198. next_cb_boundary=sfBandIndex[sfreq].l[1];  /* LONG blocks: 0,1,3 */
  199. else
  200. {
  201. next_cb_boundary=sfBandIndex[sfreq].s[1]*3; /* pure SHORT block */
  202. cb_width = sfBandIndex[sfreq].s[1];
  203. cb_begin = 0;
  204. }
  205. else
  206. next_cb_boundary=sfBandIndex[sfreq].l[1];  /* LONG blocks: 0,1,3 */
  208. /* apply formula per block type */
  209. for (sb=0 ; sb < SBLIMIT ; sb++) {        //0-31
  210. for (ss=0 ; ss < SSLIMIT ; ss++) {    //0-17
  212. /*********              由sb和ss计算cb(比例因子带)和cb_width(比例因子带宽)   **************/
  213. if ( (sb*18)+ss == next_cb_boundary)
  214. { /* Adjust critical band boundary */
  215. if (gr_info->window_switching_flag && (gr_info->block_type == 2))
  216. {
  217. if (gr_info->mixed_block_flag)
  218. {
  219. if (((sb*18)+ss) == sfBandIndex[sfreq].l[8])
  220. {
  221. next_cb_boundary=sfBandIndex[sfreq].s[4]*3;
  222. cb = 3;
  223. cb_width = sfBandIndex[sfreq].s[cb+1] -
  224. sfBandIndex[sfreq].s[cb];
  225. cb_begin = sfBandIndex[sfreq].s[cb]*3;
  226. }
  227. else if (((sb*18)+ss) < sfBandIndex[sfreq].l[8])
  228. next_cb_boundary = sfBandIndex[sfreq].l[(++cb)+1];
  229. else {
  230. next_cb_boundary = sfBandIndex[sfreq].s[(++cb)+1]*3;
  231. cb_width = sfBandIndex[sfreq].s[cb+1] -
  232. sfBandIndex[sfreq].s[cb];
  233. cb_begin = sfBandIndex[sfreq].s[cb]*3;
  234. }
  235. }
  236. else
  237. {
  238. next_cb_boundary = sfBandIndex[sfreq].s[(++cb)+1]*3;
  239. cb_width = sfBandIndex[sfreq].s[cb+1] -
  240. sfBandIndex[sfreq].s[cb];
  241. cb_begin = sfBandIndex[sfreq].s[cb]*3;
  242. }
  243. }
  244.             else /* long blocks */
  245.            next_cb_boundary = sfBandIndex[sfreq].l[(++cb)+1];
  246. }
  247. /***************************************************************************************/
  248.   if(is[sb][ss]!=0)
  249.   {
  250. /* Compute overall (global) scaling. */
  251. xr[sb][ss] = (0.25 * (gr_info->global_gain - 210.0));
  254. /* Do long/short dependent scaling operations. */
  256. if (gr_info->window_switching_flag && (
  257. ((gr_info->block_type == 2) && (gr_info->mixed_block_flag == 0)) ||
  258. ((gr_info->block_type == 2) && gr_info->mixed_block_flag && (sb >= 2)) ))
  259. {
  261. xr[sb][ss] += (0.25 * -8.0 *gr_info->subblock_gain[(((sb*18)+ss) - cb_begin)/cb_width]);
  262. xr[sb][ss] += (0.25 * -2.0 * (1.0+gr_info->scalefac_scale)
  263. * (*scalefac)[ch].s[(((sb*18)+ss) - cb_begin)/cb_width][cb]);
  264. }
  265. else    /* LONG block types 0,1,3 & 1st 2 subbands of switched blocks */
  266. {
  267. xr[sb][ss] +=  (-0.5 * (1.0+gr_info->scalefac_scale)
  268. * ((*scalefac)[ch].l[cb]
  269. + gr_info->preflag * pretab[cb]));
  270. }
  272. /* Scale quantized value. */
  273.             //求出符号
  274. sign=0;
  275. if(is[sb][ss]<0)
  276. {
  277. sign=1;
  278.     is[sb][ss]=-is[sb][ss];
  279. }
  281. if(is[sb][ss]>1023)
  282. temp=is43_table[(unsigned)is[sb][ss]>>3]+2.7725887;// 4/3*ln8
  283. else
  284. temp=is43_table[is[sb][ss]];
  286. xr[sb][ss] *=0.69314318;        //log(2);
  287. xr[sb][ss] +=temp;
  288. xr[sb][ss] =exp(xr[sb][ss]);
  289. if (sign) xr[sb][ss] = -xr[sb][ss];
  290.   }
  291.   else
  292. xr[sb][ss] =0;
  294. }
  295. }
  296. }
  297.  
  298. /*============================================================================*/
  299. /*============================================================================*/
  300. //重新排序
  302. void III_reorder(SS xr ,  struct Granule *gr_info)
  303. {
  304.    int sfreq=fr_ps.header->sampling_frequency;
  305.    int sfb, sfb_start, sfb_lines;
  306.    int sb, ss, window, freq, src_line, des_line;
  307.    SS ro ;
  308.    ////
  309.    
  310.    if (gr_info->window_switching_flag && (gr_info->block_type == 2)) 
  311.    {
  312.       if (gr_info->mixed_block_flag) 
  313.   {
  314.          /* NO REORDER FOR LOW 2 SUBBANDS */
  315.          for (sb=0 ; sb < 2 ; sb++)
  316.             for (ss=0 ; ss < SSLIMIT ; ss++) 
  317.                ro[sb][ss] = xr[sb][ss];
  318.             
  319.          /* REORDERING FOR REST SWITCHED SHORT */
  320.          for(sfb=3,sfb_start=sfBandIndex[sfreq].s[3],
  321.             sfb_lines=sfBandIndex[sfreq].s[4] - sfb_start;
  322.             sfb < 13; sfb++,sfb_start=sfBandIndex[sfreq].s[sfb],
  323.             (sfb_lines=sfBandIndex[sfreq].s[sfb+1] - sfb_start))
  324.                for(window=0; window<3; window++)
  325.                   for(freq=0;freq<sfb_lines;freq++) {
  326.                      src_line = sfb_start*3 + window*sfb_lines + freq;
  327.                      des_line = (sfb_start*3) + window + (freq*3);
  328.                      ro[des_line/SSLIMIT][des_line%SSLIMIT] =
  329.                                     xr[src_line/SSLIMIT][src_line%SSLIMIT];
  330.                }
  331.       }
  332.       else {  /* pure short */
  333.          for(sfb=0,sfb_start=0,sfb_lines=sfBandIndex[sfreq].s[1];
  334.             sfb < 13; sfb++,sfb_start=sfBandIndex[sfreq].s[sfb],
  335.             (sfb_lines=sfBandIndex[sfreq].s[sfb+1] - sfb_start))
  336.                for(window=0; window<3; window++)
  337.                   for(freq=0;freq<sfb_lines;freq++) 
  338.   {
  339.                      src_line = sfb_start*3 + window*sfb_lines + freq;
  340.                      des_line = sfb_start*3 + window + (freq*3);
  341.                      ro[des_line/SSLIMIT][des_line%SSLIMIT] =
  342.                                     xr[src_line/SSLIMIT][src_line%SSLIMIT];
  343.   }
  344.       }
  345.       for (sb=0 ; sb < SBLIMIT ; sb++)
  346.          for (ss=0 ; ss < SSLIMIT ; ss++)
  347.             xr[sb][ss] = ro[sb][ss];
  349.    }
  350. }
  351.  
  352. /*============================================================================*/
  353. /*============================================================================*/
  354. /*============================================================================*/
  355. //立体声处理
  357. static int is_pos[576];
  358. static  float is_ratio[576];
  360. void III_stereo(float xr[2][SBLIMIT][SSLIMIT],III_scalefac_t *scalefac, struct Granule *gr_info)
  361. {
  362.    int sfreq = fr_ps.header->sampling_frequency;
  363.    int stereo = fr_ps.stereo;
  364.    int ms_stereo = (fr_ps.header->mode == MPG_MD_JOINT_STEREO) &&
  365.                    (fr_ps.header->mode_ext & 0x2);
  366.    int i_stereo = (fr_ps.header->mode == MPG_MD_JOINT_STEREO) &&
  367.                   (fr_ps.header->mode_ext & 0x1);
  368.    int sfb;
  369.    int i,j,sb,ss;
  370.    /////
  371.    
  372.    /* intialization */
  373.    for ( i=0; i<576; i++ )
  374.    {
  375.       is_pos[i] = 7;
  376.    }
  377.    ////
  379.    if ((stereo == 2) && i_stereo )
  380.    {  if (gr_info->window_switching_flag && (gr_info->block_type == 2))
  381.       {  if( gr_info->mixed_block_flag )
  382.          {  int max_sfb = 0;
  384.             for ( j=0; j<3; j++ )
  385.             {  int sfbcnt;
  386.                sfbcnt = 2;
  387.                for( sfb=12; sfb >=3; sfb-- )
  388.                {  int lines;
  389.                   lines = sfBandIndex[sfreq].s[sfb+1]-sfBandIndex[sfreq].s[sfb];
  390.                   i = 3*sfBandIndex[sfreq].s[sfb] + (j+1) * lines - 1;
  391.                   while ( lines > 0 )
  392.                   {  if ( xr[1][i/SSLIMIT][i%SSLIMIT] != 0.0 )
  393.                      {  sfbcnt = sfb;
  394.                         sfb = -10;
  395.                         lines = -10;
  396.                      }
  397.                      lines--;
  398.                      i--;
  399.                   }
  400.                }
  401.                sfb = sfbcnt + 1;
  403.                if ( sfb > max_sfb )
  404.                   max_sfb = sfb;
  406.                while( sfb<12 )
  407.                {  sb = sfBandIndex[sfreq].s[sfb+1]-sfBandIndex[sfreq].s[sfb];
  408.                   i = 3*sfBandIndex[sfreq].s[sfb] + j * sb;
  409.                   for ( ; sb > 0; sb--)
  410.                   {  is_pos[i] = (*scalefac)[1].s[j][sfb];
  411.                      if ( is_pos[i] != 7 )
  412.                         is_ratio[i] = tan( is_pos[i] * (PI / 12));
  413.                      i++;
  414.                   }
  415.                   sfb++;
  416.                }
  417.                sb = sfBandIndex[sfreq].s[11]-sfBandIndex[sfreq].s[10];
  418.                sfb = 3*sfBandIndex[sfreq].s[10] + j * sb;
  419.                sb = sfBandIndex[sfreq].s[12]-sfBandIndex[sfreq].s[11];
  420.                i = 3*sfBandIndex[sfreq].s[11] + j * sb;
  421.                for ( ; sb > 0; sb-- )
  422.                {  is_pos[i] = is_pos[sfb];
  423.                   is_ratio[i] = is_ratio[sfb];
  424.                   i++;
  425.                }
  426.              }
  427.              if ( max_sfb <= 3 )
  428.              {  i = 2;
  429.                 ss = 17;
  430.                 sb = -1;
  431.                 while ( i >= 0 )
  432.                 {  if ( xr[1][i][ss] != 0.0 )
  433.                    {  sb = i*18+ss;
  434.                       i = -1;
  435.                    } else
  436.                    {  ss--;
  437.                       if ( ss < 0 )
  438.                       {  i--;
  439.                          ss = 17;
  440.                       }
  441.                    }
  442.                 }
  443.                 i = 0;
  444.                 while ( sfBandIndex[sfreq].l[i] <= sb )
  445.                    i++;
  446.                 sfb = i;
  447.                 i = sfBandIndex[sfreq].l[i];
  448.                 for ( ; sfb<8; sfb++ )
  449.                 {  sb = sfBandIndex[sfreq].l[sfb+1]-sfBandIndex[sfreq].l[sfb];
  450.                    for ( ; sb > 0; sb--)
  451.                    {  is_pos[i] = (*scalefac)[1].l[sfb];
  452.                       if ( is_pos[i] != 7 )
  453.                          is_ratio[i] = tan( is_pos[i] * (PI / 12));
  454.                       i++;
  455.                    }
  456.                 }
  457.             }
  458.          } else
  459.          {  for ( j=0; j<3; j++ )
  460.             {  int sfbcnt;
  461.                sfbcnt = -1;
  462.                for( sfb=12; sfb >=0; sfb-- )
  463.                {  int lines;
  464.                   lines = sfBandIndex[sfreq].s[sfb+1]-sfBandIndex[sfreq].s[sfb];
  465.                   i = 3*sfBandIndex[sfreq].s[sfb] + (j+1) * lines - 1;
  466.                   while ( lines > 0 )
  467.                   {  if ( xr[1][i/SSLIMIT][i%SSLIMIT] != 0.0 )
  468.                      {  sfbcnt = sfb;
  469.                         sfb = -10;
  470.                         lines = -10;
  471.                      }
  472.                      lines--;
  473.                      i--;
  474.                   }
  475.                }
  476.                sfb = sfbcnt + 1;
  477.                while( sfb<12 )
  478.                {  sb = sfBandIndex[sfreq].s[sfb+1]-sfBandIndex[sfreq].s[sfb];
  479.                   i = 3*sfBandIndex[sfreq].s[sfb] + j * sb;
  480.                   for ( ; sb > 0; sb--)
  481.                   {  is_pos[i] = (*scalefac)[1].s[j][sfb];
  482.                      if ( is_pos[i] != 7 )
  483.                         is_ratio[i] = tan( is_pos[i] * (PI / 12));
  484.                      i++;
  485.                   }
  486.                   sfb++;
  487.                }
  489.                sb = sfBandIndex[sfreq].s[11]-sfBandIndex[sfreq].s[10];
  490.                sfb = 3*sfBandIndex[sfreq].s[10] + j * sb;
  491.                sb = sfBandIndex[sfreq].s[12]-sfBandIndex[sfreq].s[11];
  492.                i = 3*sfBandIndex[sfreq].s[11] + j * sb;
  493.                for ( ; sb > 0; sb-- )
  494.                {  is_pos[i] = is_pos[sfb];
  495.                   is_ratio[i] = is_ratio[sfb];
  496.                   i++;
  497.                }
  498.             }
  499.          }
  500.       } else
  501.       {  i = 31;
  502.          ss = 17;
  503.          sb = 0;
  504.          while ( i >= 0 )
  505.          {  if ( xr[1][i][ss] != 0.0 )
  506.             {  sb = i*18+ss;
  507.                i = -1;
  508.             } else
  509.             {  ss--;
  510.                if ( ss < 0 )
  511.                {  i--;
  512.                   ss = 17;
  513.                }
  514.             }
  515.          }
  516.          i = 0;
  517.          while ( sfBandIndex[sfreq].l[i] <= sb )
  518.             i++;
  519.          sfb = i;
  520.          i = sfBandIndex[sfreq].l[i];
  521.          for ( ; sfb<21; sfb++ )
  522.          {  sb = sfBandIndex[sfreq].l[sfb+1] - sfBandIndex[sfreq].l[sfb];
  523.             for ( ; sb > 0; sb--)
  524.             {  is_pos[i] = (*scalefac)[1].l[sfb];
  525.                if ( is_pos[i] != 7 )
  526.                   is_ratio[i] = tan( is_pos[i] * (PI / 12));
  527.                i++;
  528.             }
  529.          }
  530.          sfb = sfBandIndex[sfreq].l[20];
  531.          for ( sb = 576 - sfBandIndex[sfreq].l[21]; sb > 0; sb-- )
  532.          {  is_pos[i] = is_pos[sfb];
  533.             is_ratio[i] = is_ratio[sfb];
  534.             i++;
  535.          }
  536.       }
  537.    }
  540.    if (stereo==2)
  541.       for(sb=0;sb<SBLIMIT;sb++)
  542.          for(ss=0;ss<SSLIMIT;ss++)
  543.  {
  544.             float xr0,xr1;
  545. i = (sb*18)+ss;
  546. xr0=xr[0][sb][ss];
  547. xr1=xr[1][sb][ss];
  548.             if ( is_pos[i] == 7 ) {
  549.                if ( ms_stereo ) {
  550.                   xr[0][sb][ss] = (xr0+xr1)/1.41421356;
  551.                   xr[1][sb][ss] = (xr0-xr1)/1.41421356;
  552.                }
  553.                else {
  554.                   xr[0][sb][ss] = xr0;
  555.                   xr[1][sb][ss] = xr1;
  556.                }
  557.             }
  558.             else if (i_stereo ) {
  559.                xr[0][sb][ss] = xr0 * (is_ratio[i]/(1+is_ratio[i]));
  560.                xr[1][sb][ss] = xr0 * (1/(1+is_ratio[i]));
  561.             }
  562. //            else 
  563. //               printf("Error in streo processingn");
  564.             
  565.          }
  567. } 
  569. /*============================================================================*/
  570. /*============================================================================*/
  571. /*============================================================================*/
  573. float Ci[8]={-0.6,-0.535,-0.33,-0.185,-0.095,-0.041,-0.0142,-0.0037};
  574. void III_antialias(SS hybridIn, struct Granule *gr_info)
  575. {
  576.    static int    init = 1;
  577.    static float ca[8],cs[8];
  578.    float        bu,bd;  /* upper and lower butterfly inputs */
  579.    int           ss,sb,sblim;
  581.    if (init) {
  582.       int i;
  583.       float    sq;
  584.       for (i=0;i<8;i++) {
  585.          sq=sqrt(1.0+Ci[i]*Ci[i]);
  586.          cs[i] = 1.0/sq;
  587.          ca[i] = Ci[i]/sq;
  588.       }
  589.       init = 0;
  590.    }
  592.    /* clear all inputs */
  595.    if  (gr_info->window_switching_flag && (gr_info->block_type == 2) &&
  596.        !gr_info->mixed_block_flag ) return;
  598.    if ( gr_info->window_switching_flag && gr_info->mixed_block_flag &&
  599.      (gr_info->block_type == 2))
  600.       sblim = 1;
  601.    else
  602.       sblim = SBLIMIT-1;
  604.    /* 31 alias-reduction operations between each pair of sub-bands */
  605.    /* with 8 butterflies between each pair                         */
  607.    for(sb=0;sb<sblim;sb++)
  608.       for(ss=0;ss<8;ss++)
  609.   {
  610.          bu = hybridIn[sb][17-ss];
  611.          bd = hybridIn[sb+1][ss];
  612.          hybridIn[sb][17-ss] =(bu * cs[ss]) - (bd * ca[ss]);
  613.          hybridIn[sb+1][ss] = (bd * cs[ss]) + (bu * ca[ss]);
  614.       }
  615. }
  617. /*============================================================================*/
  618. /*============================================================================*/
  619. /*============================================================================*/
  622. /*
  624. void out_fifo(short  SAM[2][SSLIMIT][SBLIMIT], FILE *outFile)
  625. {
  626.     short *pcm=&SAM[0][0][0];
  627.     fwrite(pcm,2,SBLIMIT*SSLIMIT*2,outFile);
  628. }
  630. */
  631. void  buffer_CRC()
  632. {
  633.     getbit(16);
  634. }
  636. /*============================================================================*/
  637. /*============================================================================*/
  639. /* Return the number of slots for main data of current frame, */
  640. int main_data_slots()
  641. {
  642. unsigned int nSlots; 
  643.  
  644. nSlots = (unsigned int)((144 * bitrate[2][fr_ps.header->bitrate_index])
  645. / s_freq[fr_ps.header->sampling_frequency]);
  646. if (fr_ps.header->padding) nSlots++;
  647. nSlots -= 4;                         //减去4字节的头
  648. if (fr_ps.header->error_protection)  //如有CRC 则减2字节
  649. nSlots -= 2;
  650. if (fr_ps.stereo == 1)               //单声道-17字节的side 信息
  651. nSlots -= 17;
  652. else
  653. nSlots -=32;                     //单声道-32字节的side 信息
  654. return(nSlots);
  655. }
  657. /*============================================================================*/
  658. /*============================================================================*/
  659. /*============================================================================*/