带通FIR滤波器

bandpass FIR filter

本文关键字：滤波器 FIR 带通      更新时间：2023-10-16

我需要制作一个简单的带通音频滤波器。现在我使用了这个简单的C++类：http://www.cardinalpeak.com/blog/a-c-class-to-implement-low-pass-high-pass-and-band-pass-filters

它运行良好并切断了所需的条带。但是当我尝试以小步长更改上限或下限时，在某些极限值上，我听到错误的结果 - 频率衰减或偏移（不对应于电流限制）声音。

脉冲响应计算功能：

void Filter::designBPF()
{
    int n;
    float mm;
    for(n = 0; n < m_num_taps; n++){
        mm = n - (m_num_taps - 1.0) / 2.0;
        if( mm == 0.0 ) m_taps[n] = (m_phi - m_lambda) / M_PI;
        else m_taps[n] = (   sin( mm * m_phi ) -
                             sin( mm * m_lambda )   ) / (mm * M_PI);
    }
    return;
}

哪里 

m_lambda = M_PI * Fl / (Fs/2);
m_phi = M_PI * Fu / (Fs/2);

Fs - 采样率 （44.100）Fl - 下限福 - 上限

和简单的过滤功能：

float Filter::do_sample(float data_sample)
{
    int i;
    float result;
    if( m_error_flag != 0 ) return(0);
    for(i = m_num_taps - 1; i >= 1; i--){
        m_sr[i] = m_sr[i-1];
    }   
    m_sr[0] = data_sample;
    result = 0;
    for(i = 0; i < m_num_taps; i++) result += m_sr[i] * m_taps[i];
    return result;
}

我需要使用任何窗口函数（布莱克曼等）吗？如果是，我该怎么做？我试图增加我对布莱克曼窗口的冲动反应：

m_taps[n] *= 0.42 - 0.5 * cos(2.0 * M_PI * n / double(N - 1)) +
                0.08 * cos(4.0 * M_PI * n / double(N - 1));

但结果是错误的。我需要规范化水龙头吗？

我发现了FIR滤波器的一个很好的免费实现：http://www.iowahills.com/A7ExampleCodePage.html

。此窗口FIR滤波器C代码分为两部分，第一部分是 计算矩形窗口的脉冲响应（低 通过、高通、带通或陷波）。然后是一扇窗户（凯撒、汉宁、 等）应用于脉冲响应。有几个窗口 选择...

y[i] = 波形[i] × （0.42659071 – 0.49656062cos（w） + 0.07684867cos（2w））

其中w = （2）i/n，n是波形中的元素数

试试这个，我从以下位置获取代码：http://zone.ni.com/reference/en-XX/help/370592P-01/digitizers/blackman_window/

我希望这有所帮助。

相关文章：