matlab filter 函数,C++ 实现matlab filter()函数

Matlab中的filter函数是一个数字滤波器，可以将一个输入信号通过一个滤波器进行滤波处理，得到一个输出信号。C++中可以使用各种数字滤波算法来实现类似的功能。

以下是一个示例代码，使用C++实现了类似于Matlab filter函数的数字滤波器：

#include<iostream>
#include<math.h>

using namespace std;

double filter(double *b, int nb, double *a, int na, double *x, double *y, int n)
{
    int i,j;
    double sum;
    for(i=0;i<n;i++)
    {
        sum=b[0]*x[i];
        for(j=1;j<=nb;j++)
        {
            if(i-j>=0)
                sum+=b[j]*x[i-j];
        }
        for(j=1;j<=na;j++)
        {
            if(i-j>=0)
                sum-=a[j]*y[i-j];
        }
        y[i]=sum;
    }
    return y[n-1];
}

int main()
{
    double b[3]={0.1,0.2,0.1};
    double a[3]={1,-0.2,0.1};
    double x[10]={1,2,3,4,5,6,7,8,9,10};
    double y[10]={0};
    int nb=2,na=2,n=10;
    double output=filter(b,nb,a,na,x,y,n);
    cout<<"Output: "<<output<<endl;
    return 0;
}

此处，我们使用差分方程的方法实现数字滤波器，其中b和a分别是分子和分母多项式的系数，x和y分别是输入和输出信号的数组，n是信号的长度。

在主函数中，我们定义了一组滤波器的系数和输入信号，然后调用filter函数进行滤波处理，并输出结果。

需要注意的是，这只是一个简单的实现示例，实际应用中可能需要考虑更多的因素，比如滤波器的稳定性、滤波器的阶数等等。