IIR 滤波器的实现(C++)

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
下面是一个简单的C++代码,用于实现二阶IIR滤波器的设计和滤波: ```cpp #include <vector> #include <cmath> class IIRFilter { private: double b0, b1, b2, a1, a2; double x1, x2, y1, y2; public: IIRFilter(double fc, double fs, double Q, int type) { double w0 = 2 * M_PI * fc / fs; double alpha = sin(w0) / (2 * Q); double a = exp(-alpha * w0); double beta = sqrt(1 - pow(a, 2)); switch (type) { case 0: // Low-pass filter b0 = (1 - a) / 2; b1 = 1 - a; b2 = (1 - a) / 2; a1 = -2 * a; a2 = pow(a, 2); break; case 1: // High-pass filter b0 = (1 + a) / 2; b1 = -(1 + a); b2 = (1 + a) / 2; a1 = -2 * a; a2 = pow(a, 2); break; case 2: // Band-pass filter b0 = alpha; b1 = 0; b2 = -alpha; a1 = -2 * beta * cos(w0); a2 = pow(beta, 2); break; case 3: // Band-stop filter b0 = 1; b1 = -2 * cos(w0); b2 = 1; a1 = -2 * beta * cos(w0); a2 = pow(beta, 2); break; default: break; } } double filter(double x) { double y = b0 * x + b1 * x1 + b2 * x2 - a1 * y1 - a2 * y2; x2 = x1; x1 = x; y2 = y1; y1 = y; return y; } }; ``` 在这个类中,我们使用了一个向量来存储滤波器的前向和反馈系数,还定义了一个滤波函数来实现滤波过程。在构造函数中,我们根据滤波器类型和参数计算出前向和反馈系数,并将它们存储到类的私有变量中。在滤波函数中,我们通过更新输入和输出历史值来实现滤波过程。 使用这个类的示例代码如下: ```cpp #include <iostream> int main() { double fc = 1000; // Cutoff frequency double fs = 44100; // Sampling frequency double Q = 1 / sqrt(2); // Quality factor int type = 0; // Low-pass filter IIRFilter filter(fc, fs, Q, type); // Generate a test signal (sine wave) double f = 440; // Frequency double duration = 5; // Duration (seconds) int samples = duration * fs; // Number of samples std::vector<double> x(samples); for (int i = 0; i < samples; i++) { x[i] = sin(2 * M_PI * f * i / fs); } // Apply the filter to the test signal std::vector<double> y(samples); for (int i = 0; i < samples; i++) { y[i] = filter.filter(x[i]); } // Output the filtered signal for (int i = 0; i < samples; i++) { std::cout << y[i] << std::endl; } return 0; } ``` 在这个示例中,我们首先创建了一个IIRFilter对象,并设置了滤波器的类型、截止频率、采样频率和品质因数。然后,我们生成一个测试信号(正弦波),并将其送入滤波器中进行滤波。最后,我们将滤波后的信号输出到控制台。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值