目前我正在尝试实现FIR低通滤波器。在MATLAB中计算了FIR系数。现在我需要在C++中实现FIR算法。
我将一个类定义为过滤器,并将FIR函数定义为:
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滤波器?
发布于 2014-06-20 09:05:17
我看到你试图做一个卷积,如果你想通过在数组之外索引来避免分段错误,你需要改变表达式。例如,您可以使用而不是
_x_sign[pp -jj]
pp>jj? _x_sign[pp - jj] : 0在本例中,通过使用三元操作符避免在内存外进行索引(如果jj大于pp,则不会计算_x_signpp - jj )。
您还可以更改
for (jj = 0; jj < order+1&& jj< pp+1; jj++)
{
sum += _b[jj] * _x_sign[pp - jj];
}第二种解决方案更好,因为当它无用时,您会避免做一个更大的循环。
https://stackoverflow.com/questions/24322055
复制相似问题