码迷,mamicode.com
首页 > 其他好文 > 详细

fir滤波器

时间:2018-05-03 10:29:18      阅读:141      评论:0      收藏:0      [点我收藏+]

标签:fir

 // reference implementation of an FIR     template <typename SampleType, typename NumericType>     static void reference (const NumericType* firCoefficients, size_t numCoefficients,                            const SampleType* input, SampleType* output, size_t n) noexcept     {         if (numCoefficients == 0)         {//输出清零             zeromem (output, sizeof (SampleType) * n);             return;         } //分配内存         HeapBlock<SampleType> scratchBuffer (numCoefficients                                             #if JUCE_USE_SIMD                                              + (SIMDRegister<NumericType>::SIMDRegisterSize / sizeof (SampleType))                                             #endif                                              );        #if JUCE_USE_SIMD         SampleType* buffer = reinterpret_cast<SampleType*>  (SIMDRegister<NumericType>::getNextSIMDAlignedPtr  (reinterpret_cast<NumericType*> (scratchBuffer.getData())));        #else         SampleType* buffer = scratchBuffer.getData();        #endif //内存清零         zeromem (buffer, sizeof (SampleType) * numCoefficients);         for (size_t i = 0; i < n; ++i)         { //滑动,将数据向左移动一个位置             for (size_t j = (numCoefficients - 1); j >= 1; --j)                 buffer[j] = buffer[j-1];             buffer[0] = input[i];//最新数据             SampleType sum (0);//累加和             for (size_t j = 0; j < numCoefficients; ++j)                 sum += buffer[j] * firCoefficients[j];             output[i] = sum;         }     }


fir滤波器

标签:fir

原文地址:http://blog.51cto.com/9233403/2111978

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!