(C++ 实现matlab filter()函数)
C++ 实现matlab filter()函数
笔者在做信号处理的过程中,用到了一个带通滤波器,通过matlab计算设计好参数之后,直接调用filter(b,a,x)即可得到滤波结果,在将matlab代码转换为C++实现时,对于滤波器的结构参数a,b可以通过matlab计算得到,因而只需要实现filter()这个函数。笔者实现的filter()函数如下:
实现的filter()函数如下:
//a,b由matlab设计滤波器得到的滤波器参数,x 输入,y输出,n滤波器阶数,s_n输入数据长度
void filter(doublea,doubleb,doublex,doubley,int n,int s_n) {
for (int i = 0; i if (i < n){
for (int k = 0; k <= i; k++){
y[i] += b[k] * x[i - k];}
for (int k = 1; k <= i; k++){
y[i]-= a[k] * y[i - k];}
}
else{
for (int m = 0; m < n; m++){
y[i] += b[m] * x[i - m];}
for (int m = 1; m < n; m++){
y[i] -= a[m] * y[i - m];}
}
}
}
一个简单的demo:
#includeusing namespace std;
void filter(double*a,double*b,double*x,double*y,int n,int s_n) {
for (int i = 0; i
对应的matlab代码:
y=filter(b,a,x)
测试结果我就不贴了,自己去运行
参考:https://ww2.mathworks.cn/help/matlab/ref/filter.html