在 0.667 和 3.833 Hz 范围内快速实现 IIR 滤波器

Question

我需要去除 0.667 - 3.833 Hz 范围之外的频率(心率范围:40 至 230 bpm)。我查看了 Apple 文档，但只看到了 AVFoundation Butterworth 过滤器。显然，“值的有效范围是 20 Hz 到 (SampleRate/2)。”该滤波器来自 AVAudioUnitEQ，所以这是可以理解的……我似乎找不到任何在我正在使用的频率范围内工作的滤波器。

我是否必须实现我自己的 Butterworth 过滤器，或者我是否遗漏了什么？

编辑:所以我深入研究了文档，发现了 IIR filter和 FIR filter ，但这并不是我所期望的...... IIR 滤波器是一个“级联双二阶 IIR 滤波器”，而不是 Butterworth filter (这是我更喜欢的)。我做了一个little reading我不知道这个过滤器应该如何工作。如果我使用这些系数，是否可以实现与 Butterworth 滤波器类似的功能？如果是这样，我该怎么做？目标是拥有一个没有波纹的平坦通带。

我大致关注的博文:http://www.ignaciomellado.es/blog/Measuring-heart-rate-with-a-smartphone-camera

包含一些双二阶滤波器建议的博文:http://www.ignaciomellado.es/blog/Measuring-heart-rate-with-a-smartphone-camera

Answer 1

IIR 滤波器的实现在任何语言中都非常简单。 IIR 的设计既是挑战又是乐趣所在。可以做出很多权衡。双二阶 IIR 滤波器是非常好的二阶滤波器(可以使用巴特沃斯系数进行配置)，对于许多应用来说，一个就足够了，但是您也可以级联它们以制作更高阶的滤波器，如果您需要。

首先要弄清楚你真正需要什么，这通常意味着你需要了解你想要过滤的噪声在频域中是什么样子，然后设计一个抑制噪声的滤波器。

有很多工具可以让您了解给定的一组 IIR 系数的频率响应是什么样的。

这里是 one许多来自快速谷歌搜索。

一个关键参数是您对数据源进行采样的频率(采样率)

一旦你有了系数，实现就非常简单，它只是当前输入样本和一些先前输入和输出的函数。

你可以找到一些细节here , 但这是一个函数。

y[n] = b0* x[n] + b1 * x[n-1] + b2 * x[n-2] - a1 * y[n-1] - a2 * y[n-2]

其中 x[n] 是您当前的样本，x[n-1] 是之前的样本，依此类推，y[n] 是您当前的输出值，而 y[n-1]、y[n- 2] 是以前的输出。

如果您真的想了解过滤器设计的细节，我找到的最好的资源之一是使用 SciPy ，同样也有很多关于如何执行此操作的教程。

我没看过，但这是一个 youtube video描述如何用 Python 设计巴特沃斯滤波器。

正如我之前所说，一旦你有了系数，任何语言的实现都非常简单。例如，这是我用 c++ 做的一个非常简单的 Biquad，用于过滤 Arduino 项目中的数据。这可以很容易地转换为 swift 或任何其他语言。

class BiQuad {

public:
  BiQuad(double b0, double b1, double b2, double a1, double a2)
      : b0_(b0), b1_(b1), b2_(b2), a1_(a1), a2_(a2) {}


  double processSample(double xn) {
    xn_ = xn;
    yn_ = b0_ * xn_ + b1_ * xnm1_ + b2_ * xnm2_ - a1_ * ynm1_ - a2_ * ynm2_;
    ynm2_ = ynm1_;
    ynm1_ = yn_;
    xnm2_ = xnm1_;
    xnm1_ = xn_;
    return yn_;
  }

  double getCurrentY() { return yn_; }

  void prime(double value) {
    xn_ = xnm1_ = xnm2_ = yn_ = ynm1_ = ynm2_ = value;
  }

private:
  double b0_;
  double b1_;
  double b2_;
  double a1_;
  double a2_;
  double xn_;
  double xnm1_ = 0.0;
  double xnm2_ = 0.0;
  double yn_ = 0.0;
  double ynm1_ = 0.0;
  double ynm2_ = 0.0;
};

希望这对您有所帮助。