IIR和FIR两种滤波器有什么区别?

2024-06-11 08:52
文章标签 区别 两种 滤波器 fir iir

本文主要是介绍IIR和FIR两种滤波器有什么区别?,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

概念的区分

IIR(Infinite Impulse Response,无限脉冲响应)和FIR(Finite Impulse Response,有限脉冲响应)滤波器是两种常见的数字信号处理滤波器类型,它们在结构、性能和用途上有显著区别:

结构差异:

IIR滤波器具有反馈结构,即滤波器的输出不仅取决于当前和过去输入信号的值,还包括过去的输出信号。这意味着IIR滤波器的冲激响应永远不会完全衰减到零,理论上无穷期延续。
FIR滤波器则不包含反馈环路,其输出只依赖于当前和过去的输入信号。FIR滤波器的冲激响应在有限时间内会衰减至零。


稳定性:

FIR滤波器总是稳定的,因为没有反馈回路,不存在发散的可能性。
IIR滤波器的稳定性取决于其极点位置,如果所有的极点都在单位圆内,则滤波器是稳定的;否则,可能会导致不稳定。


相位响应:

FIR滤波器具有线性相位响应,即所有频率的信号通过滤波器后,其相位延迟与频率成线性关系,这对于保持信号的相位同步和时间对齐很重要。
IIR滤波器一般具有非线性相位响应,除非特别设计,否则相位延迟会随着频率变化,这在某些应用中可能会带来问题。


滤波特性:

IIR滤波器因其反馈结构可以实现较高的滤波效率,用较少的阶数就能达到所需的滤波特性,比如陡峭的滚降特性或接近理想滤波器的响应曲线。
FIR滤波器虽然阶数较高时才能达到类似的效果,但因其线性相位和无发散风险,更适合需要精确相位和稳态无误差的场合。


设计灵活性:

FIR滤波器的设计通常更为直观和灵活,易于实现严格的带宽、阻带衰减和线性相位响应。
IIR滤波器可以基于经典滤波器设计方法(如巴特沃兹、切比雪夫等)实现特定频率响应,但其非线性相位特性增加了设计难度。


在工程实践中,选择IIR还是FIR滤波器通常取决于具体的应用需求。如果对相位线性度要求较高,或担心稳定性问题,通常会选择FIR滤波器。而在对计算资源有限,需要以较低阶数实现较强滤波效果,且对相位非线性有一定的容忍度时,IIR滤波器则是一个合适的选择。在现代DSP应用中,由于计算能力的增强,FIR滤波器因其稳定性、线性相位和易于设计的特点而越来越受到青睐。

C语言实现

实现IIR(无限脉冲响应)和FIR(有限脉冲响应)带通滤波器通常涉及到计算每个新样本点时的滤波操作。以下是对这两种滤波器的基本实现框架的概述以及示例C代码片段:

IIR带通滤波器

IIR滤波器通常通过递归的方式实现,即输出不仅取决于当前的输入样本,还取决于过去若干个输入样本以及过去若干个输出样本。巴特沃斯、切比雪夫、椭圆滤波器等经典滤波器结构都可以被用来设计IIR带通滤波器。以下是一个二阶Butterworth IIR带通滤波器的例子:

#include <math.h>// 假设我们已经有了滤波器的系数
#define B0 0.5 * (1.0 + sqrt(2)) // 振幅系数B0
#define B1 0.0                  // 振幅系数B1
#define B2 (-1.0 * B0)          // 振幅系数B2
#define A1 (-1.0 * sqrt(2) / 2.0) // 阻尼系数A1
#define A2 (1.0 - sqrt(2) / 2.0)  // 阻尼系数A2// 初始化历史状态变量
float x_n_1 = 0.0; // 上一次输入样本
float x_n_2 = 0.0; // 上上次输入样本
float y_n_1 = 0.0; // 上一次输出样本
float y_n_2 = 0.0; // 上上次输出样本// IIR带通滤波器核心函数
float iir_bandpass_filter(float input)
{// 计算当前输出float output = B0 * input + B1 * x_n_1 + B2 * x_n_2 - A1 * y_n_1 - A2 * y_n_2;// 更新历史状态x_n_2 = x_n_1;x_n_1 = input;y_n_2 = y_n_1;y_n_1 = output;return output;
}// 使用滤波器
int main()
{float input_signal[signal_length];float output_signal[signal_length];// ... 填充input_signal...for (int i = 0; i < signal_length; ++i){output_signal[i] = iir_bandpass_filter(input_signal[i]);}// ... 处理output_signal...return 0;
}


注意,实际使用时你需要根据具体的设计参数(如截止频率、品质因数等)计算合适的滤波器系数B0、B1、B2、A1和A2。

FIR带通滤波器

FIR滤波器没有反馈部分,仅依赖于当前和过去若干个输入样本。带通FIR滤波器通常通过卷积实现,系数来自于窗函数法或其他设计方法。
 

#include <stdio.h>// 滤波器系数数组h[]
#define FILTER_ORDER 51 // FIR滤波器阶数
float h[FILTER_ORDER]; // 初始化滤波器系数
// ...省略了根据设计要求填充滤波器系数的过程...// FIR带通滤波器核心函数
void fir_bandpass_filter(float input, float* state, float* output)
{static const int STATE_SIZE = FILTER_ORDER - 1; // FIR滤波器的状态大小static float input_state[STATE_SIZE]; // 存储过去输入样本的状态// 移除最早的一个输入样本,添加最新的输入样本memmove(input_state + 1, input_state, sizeof(float) * (STATE_SIZE - 1));input_state[0] = input;// FIR滤波输出计算*output = 0.0;for (int i = 0; i < FILTER_ORDER; ++i){*output += h[i] * input_state[i];}
}// 使用滤波器
int main()
{float input_signal[signal_length];float output_signal[signal_length];float filter_state[FILTER_ORDER - 1] = {0}; // 初始化状态为零// ... 填充input_signal...for (int i = 0; i < signal_length; ++i){fir_bandpass_filter(input_signal[i], filter_state, &output_signal[i]);}// ... 处理output_signal...return 0;
}


同样,这里的FIR滤波器系数`h[]`需要根据目标频率响应曲线设计得到。在实际项目中,可以使用像MATLAB、Python的SciPy库或者其他DSP工具设计滤波器并获取这些系数。

请确保根据实际应用调整滤波器阶数、类型和系数,同时注意边界条件的处理,特别是对于实时流数据的第一段样本,可能需要预先填充足够的“零”值来初始化滤波器状态。

这篇关于IIR和FIR两种滤波器有什么区别?的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1050712

相关文章

Springboot @Autowired和@Resource的区别解析

《Springboot@Autowired和@Resource的区别解析》@Resource是JDK提供的注解,只是Spring在实现上提供了这个注解的功能支持,本文给大家介绍Springboot@... 目录【一】定义【1】@Autowired【2】@Resource【二】区别【1】包含的属性不同【2】@

Java中的String.valueOf()和toString()方法区别小结

《Java中的String.valueOf()和toString()方法区别小结》字符串操作是开发者日常编程任务中不可或缺的一部分,转换为字符串是一种常见需求,其中最常见的就是String.value... 目录String.valueOf()方法方法定义方法实现使用示例使用场景toString()方法方法

分辨率三兄弟LPI、DPI 和 PPI有什么区别? 搞清分辨率的那些事儿

《分辨率三兄弟LPI、DPI和PPI有什么区别?搞清分辨率的那些事儿》分辨率这个东西,真的是让人又爱又恨,为了搞清楚它,我可是翻阅了不少资料,最后发现“小7的背包”的解释最让我茅塞顿开,于是,我... 在谈到分辨率时,我们经常会遇到三个相似的缩写:PPI、DPI 和 LPI。虽然它们看起来差不多,但实际应用

Docker镜像pull失败两种解决办法小结

《Docker镜像pull失败两种解决办法小结》有时候我们在拉取Docker镜像的过程中会遇到一些问题,:本文主要介绍Docker镜像pull失败两种解决办法的相关资料,文中通过代码介绍的非常详细... 目录docker 镜像 pull 失败解决办法1DrQwWCocker 镜像 pull 失败解决方法2总

在C#中调用Python代码的两种实现方式

《在C#中调用Python代码的两种实现方式》:本文主要介绍在C#中调用Python代码的两种实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录C#调用python代码的方式1. 使用 Python.NET2. 使用外部进程调用 Python 脚本总结C#调

GORM中Model和Table的区别及使用

《GORM中Model和Table的区别及使用》Model和Table是两种与数据库表交互的核心方法,但它们的用途和行为存在著差异,本文主要介绍了GORM中Model和Table的区别及使用,具有一... 目录1. Model 的作用与特点1.1 核心用途1.2 行为特点1.3 示例China编程代码2. Tab

IDEA中Git版本回退的两种实现方案

《IDEA中Git版本回退的两种实现方案》作为开发者,代码版本回退是日常高频操作,IntelliJIDEA集成了强大的Git工具链,但面对reset和revert两种核心回退方案,许多开发者仍存在选择... 目录一、版本回退前置知识二、Reset方案:整体改写历史1、IDEA图形化操作(推荐)1.1、查看提

Nginx指令add_header和proxy_set_header的区别及说明

《Nginx指令add_header和proxy_set_header的区别及说明》:本文主要介绍Nginx指令add_header和proxy_set_header的区别及说明,具有很好的参考价... 目录Nginx指令add_header和proxy_set_header区别如何理解反向代理?proxy

Android自定义Scrollbar的两种实现方式

《Android自定义Scrollbar的两种实现方式》本文介绍两种实现自定义滚动条的方法,分别通过ItemDecoration方案和独立View方案实现滚动条定制化,文章通过代码示例讲解的非常详细,... 目录方案一:ItemDecoration实现(推荐用于RecyclerView)实现原理完整代码实现

Java中&和&&以及|和||的区别、应用场景和代码示例

《Java中&和&&以及|和||的区别、应用场景和代码示例》:本文主要介绍Java中的逻辑运算符&、&&、|和||的区别,包括它们在布尔和整数类型上的应用,文中通过代码介绍的非常详细,需要的朋友可... 目录前言1. & 和 &&代码示例2. | 和 ||代码示例3. 为什么要使用 & 和 | 而不是总是使