Question

求助FIR滤波器算法

我现在手头上有个问题，需要用FIR滤波器解决。大概是这样的：
首先我要求解时域上的一个微分方程，得到一组f(t)序列，同时需要对f(t)利用FIR滤波器进行处理，但是现在有个问题，因为我的时域序列是按时间步一点一点的求出来的。比如现在进行到第M步，那M+1及其之后的f(t)数值都不知道。那根据FIR滤波器的实现方法，这一点经过滤波之后的值应该怎么求？如果按照下面的公式（a是滤波器系数，N是阶数，系数在n=0点偶对称分布）：

那我M点之后的值是不知道的，是不是应该设置为0？
还是直接先把所有f(t)全部求出来之后再进行滤波？
现在挺迷惑，希望大神指教！！

Answer 1

现在进行到第M步，只能得出前面M-N时刻的滤波结果。
补0肯定是不行的。

滤波结果比采样结果延迟N点。
直接先把所有f(t)全部求出来之后再进行滤波当然更方便了。

写成

y(n)=求和a（n）*f(M-N)
更合理

Answer 2

你先去看一下卷积运算是怎么定义的，就应该知道怎么FIR滤波怎么算了。

追问

我这个不太一样，因为数据不是已知的，而是一边滤波一边求。新的数据是根据之前的数据求出来的。上面说的全部求出来再进行滤波的方法不行。