filter_dc_notch16

最新推荐文章于 2024-05-26 03:23:33 发布

it_passioner

最新推荐文章于 2024-05-26 03:23:33 发布

阅读量1k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： records

本文链接：https://blog.youkuaiyun.com/cl11010/article/details/79959740

records 专栏收录该内容

13 篇文章

订阅专栏

本文介绍了一种去除信号直流分量的算法实现，通过一个特定的函数filter_dc_notch16来处理输入信号，利用反馈机制实现直流分量的有效过滤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

看不懂这个函数，如何去除直流分量？

static inline void filter_dc_notch16(const spx_int16_t *in, spx_word16_t radius, spx_word16_t *out, int len, spx_mem_t *mem, int stride)
{
   int i;
   spx_word16_t den2;
#ifdef FIXED_POINT
   den2 = MULT16_16_Q15(radius,radius) + MULT16_16_Q15(QCONST16(.7,15),MULT16_16_Q15(32767-radius,32767-radius));
#else
   den2 = radius*radius + .7*(1-radius)*(1-radius);
#endif   
   /*printf ("%d %d %d %d %d %d\n", num[0], num[1], num[2], den[0], den[1], den[2]);*/
   for (i=0;i<len;i++)
   {
      spx_word16_t vin = in[i*stride];
      spx_word32_t vout = mem[0] + SHL32(EXTEND32(vin),15);
#ifdef FIXED_POINT
      mem[0] = mem[1] + SHL32(SHL32(-EXTEND32(vin),15) + MULT16_32_Q15(radius,vout),1);
#else
      mem[0] = mem[1] + 2*(-vin + radius*vout);
#endif
      mem[1] = SHL32(EXTEND32(vin),15) - MULT16_32_Q15(den2,vout);
      out[i] = SATURATE32(PSHR32(MULT16_32_Q15(radius,vout),15),32767);                                                                                                                                      
   }
}