Python实现MATLAB filter函数：原理与代码详解

南乔

发布时间：2026-10-15| 5人浏览过

原创

本文详细解析MATLAB中filter函数的核心原理，并提供Python实现的完整代码。通过对比b=[8,-3.5,0.5]和a=[1,-1.5,0.25]系数的计算过程，展示如何手动实现IIR滤波算法。代码包含输入序列x的处理逻辑，结果与MATLAB原生函数完全一致，适用于需要跨平台验证或理解底层算法的技术场景。

MATLAB filter函数基础

在MATLAB中，filter函数用于实现数字滤波器，其标准调用格式如下所示。

y = filter(b,a,x)

Python实现方案

下面介绍如何在Python中实现与MATLAB filter函数等效的功能。

核心算法代码


def filter_matlab(b,a,x):
  y = []
  y.append(b[0] * x[0])
  for i in range(1,len(x)):
    y.append(0)
    for j in range(len(b)):
      if i >= j :
        y[i] = y[i] + b[j] * x[i - j ]
        j += 1
    for l in range(len(b)-1 ):
      if i >l:
        y[i] = (y[i] - a[l+1] * y[i -l-1])
        l += 1
    i += 1
  return y

测试用例与结果验证

具体示例如下：

设定滤波器系数b和a，以及输入信号x，调用自定义的filter_matlab函数进行计算。

b = [8,-3.5,0.5] a = [1,-1.5,0.25] x = [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] y = filter_matlab(b,a,x)

该Python实现函数的输出结果与MATLAB的filter函数计算结果完全一致，具体数值如下。


[8,
 24.5,
 52.25,
 94.75,
 156.5625,
 243.65625,
 363.84375,
 527.3515625,
 747.56640625,
 1042.01171875,
 1433.6259765625,
 1952.43603515625,
 2637.74755859375,
 3541.0123291015625,
 4729.581604003906,
 6291.619323730469,
 8342.533584594727,
 11033.395545959473,
 14561.959922790527,
 19187.090997695923]

以上为本文的全部内容，希望对读者有所帮助，欢迎继续支持本站。