以下是一个简单的C#语言实现带通滤波器算法的例子:
public class BandpassFilter
{
private double[] coefficients;
public BandpassFilter(int sampleRate, int centerFrequency, int bandwidth)
{
double w1 = 2.0 * Math.PI * (centerFrequency - bandwidth / 2) / sampleRate;
double w2 = 2.0 * Math.PI * (centerFrequency + bandwidth / 2) / sampleRate;
double m = 2.0 * Math.Cos((w2 + w1) / 2.0);
double b1 = -2.0 * Math.Cos((w2 - w1) / 2.0) / m;
double b2 = 1.0 / m;
coefficients = new double[] { b2, b1, 1.0, b1, b2 };
}
public double[] Filter(double[] input)
{
double[] output = new double[input.Length];
for (int i = 2; i < input.Length; i++)
{
output[i] = coefficients[0] * input[i - 2] + coefficients[1] * input[i - 1] + coefficients[2] * input[i] + coefficients[3] * output[i - 1] + coefficients[4] * output[i - 2];
}
return output;
}
}
该类包含了一个构造函数,用于初始化带通滤波器的参数,以及一个Filter方法,用于对输入数据进行滤波处理。使用时,可以按照以下步骤进行:
// 创建带通滤波器(采样率为44100Hz,中心频率为1000Hz,带宽为500Hz)
var filter = new BandpassFilter(44100, 1000, 500);
// 输入数据(长度为1024)
var input = new double[1024];
// TODO: 将数据按照需要进行赋值
// 进行滤波处理
var output = filter.Filter(input);
// TODO: 处理输出数据
注意,该实现仅适用于实时滤波,如果需要离线处理,请使用更加高效的算法。
特此记录
anlog
2023年3月22日