音频通过滤波器后产生Noise analysis

       最近在做一些digital filter research, 进而用这些滤波器来实现一些音频功能。碰到这样一个问题，算法在matlab下用.m代码进行仿真时，效果不错，达到设计的预想，但写成C 代码时会出现一些garbage bits, 在听觉表现来上就是出现Noise,在频域来看会出现一些频域上的impluse.

      分析下来原因是不同数据类型转换时做得不够好，其实就是一些简单的细节问题，但往往细节就决定了成败。

      例如：audio 常见的是16bits short数据，但在C中的运算往往会引入float 或double型数据做为中间运算的存储值，这就有可能在中间或最后结果碰到将float or double型数据转成 short型，在这中间就会saturate operate. 

     double result;

     short  output;

      if(result > 32768.0 )

         result  = 32768.0;

      else if(result < -32768.0)

         result = -32768.0;

 

     output = (short) result;

 

     这过程中隐含重大错误，因为0x8000 = 32768. 但对short 来说: ox8000是个负数，上面的操作会导致把原来最大的数据变成一负数。

     直接后果会出现Noise. 应该改成下面形式：

     double result;

     short  output;

      if(result > =32768.0 )

         result  = 32767.0;

      else if(result < -32768.0)

         result = -32768.0;

 

     output = (short) result;

 

     小问题往往是难倒Hero, 在这记录一下，这就是经验积累。 如有同行，在做算法设计时，碰到matlab算法仿真是好的，但C实现时会碰到这样那样的怪问题时，可以把相关问题发给我，我会跟你们一起来解决的。