最近在做一些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实现时会碰到这样那样的怪问题时,可以把相关问题发给我,我会跟你们一起来解决的。