流式DSP滤波：121抽头FIR滤波器状态保存疑问

首先咱们得把核心逻辑掰明白：FIR滤波器的状态本质就是输入信号的延迟线缓存，长度正好是抽头数减1（你的121抽头对应120长度的状态）。这个缓存的作用是记住上一批数据里，需要和下一批新数据卷积的前序输入样本——只有保存好它，流式处理的结果才能和一次性处理完全一致。

一、Matlab里的状态对应逻辑 #

Matlab的filter函数其实已经帮你封装好了状态传递，你可能没注意到它的第二个返回值：

% 一次性处理全量数据
[y_full, ~] = filter(b, 1, x_full); % b是121抽头的FIR系数，a=1表示FIR结构

% 分批次流式处理
[y_batch1, zf1] = filter(b, 1, x_batch1); % zf1就是处理完第一批后的状态
[y_batch2, zf2] = filter(b, 1, x_batch2, zf1); % 传入上一批的状态zf1继续处理

这里的zf1就是延迟线的当前状态：如果x_batch1长度≥120，那zf1就是x_batch1的最后120个样本；如果x_batch1长度<120，zf1会是x_batch1加上初始零填充凑够120长度——这就是你之前看到状态和data_in不匹配的原因，短批次处理时初始零会混入状态里。

二、手动实现时的状态保存方法 #

看你写的C代码片段，手动做卷积的话，需要自己维护一个长度为120的状态数组（比如叫delay_line），步骤如下：

1. 初始化状态：一开始delay_line全填0，对应Matlab的初始零状态。
2. 处理每个输入样本：
   假设你的FIR系数是filter[0]到filter[120]，输出计算逻辑如下：

   data_out[out] = filter[0] * data_in[current_sample];
   for (int j = 1; j < 121; j++) {
       data_out[out] += filter[j] * delay_line[j-1];
   }
   

   计算完输出后，更新延迟线：把当前输入样本加入缓存，移除最旧的样本（用循环队列效率更高，这里先给个简单的右移实现）：

   // 右移延迟线，把新样本放到最前面
   for (int j = 119; j >= 0; j--) {
       delay_line[j+1] = delay_line[j];
   }
   delay_line[0] = data_in[current_sample];
   

3. 保存/加载状态：每批数据处理完后，当前的delay_line就是要保存的状态；下一批处理前，把保存的状态赋值给delay_line即可。

三、验证状态正确性的小技巧 #

你可以做个简单测试确认状态是否正确：

• 生成一段测试输入，比如长度130（120+10）的随机数组x。
• 一次性处理得到y_full。
• 分两次处理：第一次处理前10个样本，保存状态；第二次处理后120个样本，加载保存的状态。
• 把两次的输出拼接，对比和y_full是否完全一致（前10个输出是初始零状态下的结果，后120个输出应该和y_full的对应部分完全相同）。

这样就能确保你的流式处理结果和一次性处理完全对齐了。

内容的提问来源于stack exchange，提问作者BeautyQueen