C++ 实现matlab filter()函数

(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：

#include
using 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