MPEG音频实验

一.实验原理

• MPEG-1 Audio LayerII编码器原理

• 多相滤波器组
  将PCM样本变换到32个子带的频域信号
  
• 心理声学模型
  1、将样本变换到频域
  2、确定声压级别
  

3、考虑安静时阈值
4、将音频信号分解成“乐音(tones)” 和“非乐音/噪声”部分：因为两种信号的掩蔽能力不同

5、音调和非音调掩蔽成分的消除
6、单个掩蔽阈值的计算
7、全局掩蔽阈值的计算

还要考虑别的临界频带的影响。一个掩蔽信号会对其它频带上的信号产生掩蔽效应。这种掩蔽效应称为掩蔽扩散。

8、每个子带的掩蔽阈值
9、计算每个子带信号掩蔽比(signal-to-mask ratio, SMR)，并将SMR传递给编码单元。（SMR = 信号能量 / 掩蔽阈值）

二.实验要求

1.输出音频的采样率和目标码率
2.选择三个不同特性的音频文件
噪声（持续噪声、突发噪声）
音乐
混合
3. 某个数据帧，输出
该帧所分配的比特数
该帧的比例因子
该帧的比特分配结果

三. 代码实现

1.输出音频的采样率和目标码率

fprintf (stderr, "--------------------------------------------\n");
  fprintf (stderr, "Input File : '%s'   %.1f kHz\n",
	   (strcmp (inPath, "-") ? inPath : "stdin"),
	   s_freq[header->version][header->sampling_frequency]);
  fprintf (stderr, "Output File: '%s'\n",
	   (strcmp (outPath, "-") ? outPath : "stdout"));
  fprintf (stderr, "%d kbps ", bitrate[header->version][header->bitrate_index]);
  fprintf (stderr, "%s ", version_names[header->version]);
  if (header->mode != MPG_MD_JOINT_STEREO)
    fprintf (stderr, "Layer II %s Psycho model=%d  (Mode_Extension=%d)\n",
	     mode_names[header->mode], *psy, header->mode_ext);
  else
    fprintf (stderr, "Layer II %s Psy model %d \n", mode_names[header->mode],
	     *psy);

输出：
采样率：22.1khz
目标码率：96 kbps

2.选择三个不同特性的音频文件某个数据帧，输出
添加代码：

if (frameNum == 7) {
  printf("采样率：%lfkhz\n", s_freq[header.version][header.sampling_frequency]);
  printf("目标码率：%dMbps\n", bitrate[header.version][header.bitrate_index]);
  printf("该帧所分配的比特数：%d\n", lg_frame * 8);
  for (i = 0; i < nch; i++)
  {
   printf("声道%d\n", i);
   for (j = 0; j < SBLIMIT; j++) {
    printf("子带%d\t", j + 1);
    printf("比例因子：");
    switch (scfsi[i][j]) {
    case 0:
     for (k = 0; k < 3; k++)
      printf("%3d ", scalar[i][k][j]);
     break;
    case 1:
     printf("%3d ", scalar[i][0][j]);
     printf("%3d ", scalar[i][0][j]);
     printf("%3d ", scalar[i][2][j]);
     break;
    case 2:
     printf("%3d ", scalar[i][0][j]);
     printf("%3d ", scalar[i][0][j]);
     printf("%3d ", scalar[i][0][j]);
     break;
    case 3:
     printf("%3d ", scalar[i][0][j]);
     printf("%3d ", scalar[i][2][j]);
     printf("%3d ", scalar[i][2][j]);
     break;
    }
    printf(" 比特分配结果: ");
    printf("%3d\n", bit_alloc[i][j]);
   }
  }
 }

噪声：

输出：

音乐：

输出：

音乐加噪声：

输出：

四.结论

频率越高，分配到的比特数越少