首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >c++中的FIR滤波器

c++中的FIR滤波器
EN

Stack Overflow用户
提问于 2014-06-20 07:26:21
回答 1查看 934关注 0票数 1

目前我正在尝试实现FIR低通滤波器。在MATLAB中计算了FIR系数。现在我需要在C++中实现FIR算法。

我将一个类定义为过滤器,并将FIR函数定义为:

代码语言:javascript
复制
double * Filter::FIR (double x[])
{

    unsigned int jj;
    unsigned int pp;
    double sum;


    pp = 0;

    // input values
    memcpy(&_x_sign, &x, sizeof(x));


    for (pp = 0; pp < order+1; pp++)

    {
        sum = 0 ;
        for (jj = 0; jj <  order+1; jj++)
        {
            sum += _b[jj] * _x_sign[pp - jj];
        }

        _y_sign[pp] = sum;

    }

    return (_y_sign);
}

当循环中(pp -jj)为负数时,_x_sign作为数组不起作用。如何修改FIR滤波器?

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2014-06-20 09:05:17

我看到你试图做一个卷积,如果你想通过在数组之外索引来避免分段错误,你需要改变表达式。例如,您可以使用而不是

代码语言:javascript
复制
_x_sign[pp -jj]
pp>jj? _x_sign[pp - jj] : 0

在本例中,通过使用三元操作符避免在内存外进行索引(如果jj大于pp,则不会计算_x_signpp - jj )。

您还可以更改

代码语言:javascript
复制
for (jj = 0; jj < order+1&& jj< pp+1; jj++)
    {
        sum += _b[jj] * _x_sign[pp - jj];
    }

第二种解决方案更好,因为当它无用时,您会避免做一个更大的循环。

票数 3
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/24322055

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档