滤波器的频率响应总是具有固有的“滚降”，因为您实际上无法实现完美矩形函数的通带。对于低通滤波器，频率响应幅度降至 -3dB 的点称为通带，超出的所有点称为阻带（从技术上讲，所有超出转角频率的点，但我们将采用转角频率为 -3dB 水平）。您的频率响应在通带之外衰减的速度取决于滤波器长度。

如果您查看滤波器的频率响应ht，您会发现它在 1000 Hz 时下降到 -6dB：

因此，滤波后功率下降 6dB 是有道理的，您在图中可以看到幅度减半。

如果您查看了fir1您使用的函数的文档，您也会意识到这一点（强调我的）：

B = fir1(N,Wn)设计一个N' 阶低通 FIR 数字滤波器，并返回长度N+1向量中的滤波器系数B。截止频率 Wn 必须在 之间0 < Wn < 1.0，1.0 对应于采样率的一半。滤波器B是真实的并且具有线性相位。 Wn 处滤波器的归一化增益为 -6 dB。

现在要创建更清晰的过滤器，其响应非常接近矩形，您必须使用 IIR 过滤器，它有自己的稳定性问题等问题，但绝对是一种选择。您可以在此处查看我的答案，了解有关实现离散形式 2、二阶部分 IIR 滤波器的一些想法，该滤波器给出了非常尖锐的角。该示例用于带通滤波器，但您可以阅读文档以了解其中使用的功能并自己实现低通版本。

滤波器需要时间来“决定”信号是刚好在滤波器转换内部还是外部。一种解决方案是将滤波器转换从任何感兴趣的信号移开，例如在两个测试信号之间的中间，在您的测试用例中，几乎没有信号会因不完美的决定而失真。

采用

ht = fir1(40,.5,'stop');

相反，F0 处的衰减非常低，F1 处的衰减非常高。