数据压缩试验六：MPEG音频编码

一、实验原理

1、感知编码的基本思路

（1）心理声学模型：
听觉系统中存在一个听觉阈值电平，低于这个电平的声音信号就听不到。
听觉阈值的大小随声音频率的改变而改变。
一个人是否听到声音取决于声音的频率，以及声音的幅度是否高于这种频率下的听觉阈值。
听觉掩蔽特性：即听觉阈值电平是自适应的，会随听到的不同频率声音而发生变化。

声音压缩算法可以确立这种特性的模型来取消更多的冗余数据。

（2）人耳听觉特性
人类听觉系统大致等效于一个在0Hz到20KHz频率范围内由25个重叠的带通滤波器组成的滤波器组。
人耳不能区分同一频带内同时发生的不同声音；
人耳频带被称为临界频带（critical band）；
临界频带是心理声学模型分析的基本单位。
掩蔽效应在一定频率范围内不随带宽增大而改变，直至超过某个频率值。

2.程序设计的整体框架
MPEG-1 Audio Layer II编码器原理框图如下图所示：

二、代码及分析

#ifdef NEWENCODE
    // 计算比例因子
    scalefactor_calc_new(*sb_sample, scalar, nch, frame.sblimit);
    find_sf_max (scalar, &frame, max_sc);
    if (frame.actual_mode == MPG_MD_JOINT_STEREO) {
      /* this way we calculate more mono than we need */
      /* but it is cheap */
      combine_LR_new (*sb_sample, *j_sample, frame.sblimit);
      scalefactor_calc_new (j_sample, &j_scale, 1, frame.sblimit);
    }
#else
    scale_factor_calc (*sb_sample, scalar, nch, frame.sblimit);
    pick_scale (scalar, &frame, max_sc);
	

    if (frame.actual_mode == MPG_MD_JOINT_STEREO) {
      /* this way we calculate more mono than we need */
      /* but it is cheap */
      combine_LR (*sb_sample, *j_sample, frame.sblimit);
      scale_factor_calc (j_sample, &j_scale, 1, frame.sblimit);
    }
#endif
	/* 添加 */
	FILE *output = fopen("output.txt", "wb");
	if (frameNum == 20)
	{
		fprintf(output,"输出音频的采样率为：%.1fkhz \r\n", s_freq[header.version][header.sampling_frequency]);
		fprintf(output, "目标码率为：%dMbps \r\n", bitrate[header.version][header.bitrate_index]);
		fprintf(output, "选择的数据帧为：%d \r\n", frameNum);
		fprintf(output, "该帧所分配的比特数为：%dbits\r\n", adb);
		for (int k = 0; k < nch; k++)//nch即num channels: 1 for mono, 2 for stereo
		{
			fprintf(output, "声道[%d] \r\n", k + 1);
			for (int j = 0; j < frame.sblimit; j++)//子带总数
			{
				fprintf(output, "子带[%d]:    ", j + 1);
				for (int m = 0; m < 3; m++)//每个子带有3个比例因子
				{
					fprintf(output, "%d\t", scalar[k][m][j]);
				}
				fprintf(output, "\r\n");
			}
		}
	}
#ifdef NEWENCODE
    sf_transmission_pattern (scalar, scfsi, &frame);
    main_bit_allocation_new (smr, scfsi, bit_alloc, &adb, &frame, &glopts);
    //main_bit_allocation (smr, scfsi, bit_alloc, &adb, &frame, &glopts);

    if (error_protection)
      CRC_calc (&frame, bit_alloc, scfsi, &crc);

    write_header (&frame, &bs);
    //encode_info (&frame, &bs);
    if (error_protection)
      putbits (&bs, crc, 16);
    write_bit_alloc (bit_alloc, &frame, &bs);
    //encode_bit_alloc (bit_alloc, &frame, &bs);
    write_scalefactors(bit_alloc, scfsi, scalar, &frame, &bs);
    //encode_scale (bit_alloc, scfsi, scalar, &frame, &bs);
    subband_quantization_new (scalar, *sb_sample, j_scale, *j_sample, bit_alloc,
    			  *subband, &frame);
    //subband_quantization (scalar, *sb_sample, j_scale, *j_sample, bit_alloc,
    //	  *subband, &frame);
    write_samples_new(*subband, bit_alloc, &frame, &bs);
    //sample_encoding (*subband, bit_alloc, &frame, &bs);
#else
transmission_pattern (scalar, scfsi, &frame);
//根据掩蔽值为子带进行迭代比特分配让掩噪比最低之后输出每个子带的比特分配信息
    main_bit_allocation (smr, scfsi, bit_alloc, &adb, &frame, &glopts);
	/* 添加 */
	if (frameNum == 20)
	{
		fprintf(output, "\n比特分配:\n");		
		for (int k = 0; k < nch; k++)
		{
			fprintf(output, "channel[%d] \r\n", k + 1);
			for (int j = 0; j < frame.sblimit; j++)
			{
				fprintf(output, "subband[%d]:%d\r\n", j, bit_alloc[k][j]);
			}
		}
	}

三、实验结果及分析

1、实验结果

2、结果分析
实验输出结果表明，同一子带的三个比例因子相差较小，所以可以用比例因子选择的方法压缩文件；从帧比特分配来看，高频子带分配比特数较少，低频子带分批比特数较多。