0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

介绍一下FIR滤波器的使用

冬至子 来源:振动信号研究所 作者:一只小绿豆 2023-05-26 15:13 次阅读

1. 滤波的目的

信号实时处理中,滤波是十分必要的,因为信号中难免会由于各种原因混入噪声,干扰我们对信号进行分析。

这里强调了“实时”,是因为在一些场景中,我们可能需要对信号进行离线的、交互式的分析,此时仅仅是对原始信号进行采集即可,无需软件 进行滤波,只需要硬件在采集信号时做一下抗混叠滤波即可,这种场景不是本文所关心的。

在更多的场景下,比如我们的运动手表、手环在处理ppg信号计算心率时,是会进行实时计算的,那么就会进行实时滤波处理:因为大多数人的心率频率在0.7~3.6Hz范围内,至少在这之外的频率会在计算心率之前滤除掉。

2. 信号模拟

比如我们要处理一个信号,但是我们仅仅关心信号100Hz以下的频段,这时我们就需要一个低通滤波器了,此时我们先模拟出一个包含5Hz和3000Hz频率成分的信号,假设信号采样频率为8192,采样时间为1秒,共计8192个点。信号生成和展示的代码如下:

import numpy as np
from numpy import cos
import matplotlib.pyplot as plt

pi = np.pi
t = np.linspace(0, 1, 8192)

signal = 3 * cos(2 * pi * 5 * t + pi/3) + 7 * cos(2 * np.pi * 3000 * t + 3/8*pi)

plt.figure()
plt.plot(t, signal)
plt.show()

我们生成了这样的一个信号:

图片

生成的信号

3. FIR滤波器系数生成

这一步可以使用matlab进行辅助,本文仅仅是想要一个截止频率为10Hz的FIR低通滤波器,步骤如下:

  1. 打开matlab;
  2. 点击"APP";
  3. 找到滤波器设计工具,并点击;
  4. 选择响应类型、设计方法、阶数等。

图片

辅助设计界面

  1. 点击“文件” 、“导出”、“系数文件”,导出系数文件,我把导出的系数(FIR低通滤波系数仅有分子)画出来后如下图:

图片

系数波形图

3. FIR滤波原理

采用FIR进行滤波,从操作上看是进行卷积操作,对上述滤波器的系数进行FFT变换即可窥见一斑:

import numpy as np
import matplotlib.pyplot as plt
# b = [......] 101个系数组成的列表,此处省略
delta_f = 1
plt.plot([delta_f * i for i in range(4097)], abs(np.fft.rfft(b, 8192)))  # 单边FFT
plt.show()

图片

系数的傅里叶变换

这里得到的就是滤波器的幅频响应曲线,和滤波器辅助设计工具中所展示的幅频响应曲线是一致的。

4. 滤波计算代码与结果

把第二步生成的信号中高于100Hz的频率成分(即3000Hz的成分)滤除,得到的结果如下图所示,在结果的开头和结尾产生了失真,这是因为在卷积运算过程中,在刚开始和结尾计算时有效信息不足进行了补零操作导致的。

图片

滤波的结果

这一步的运算代码如下:

import numpy as np
from numpy import cos
import matplotlib.pyplot as plt
from filter_coeff import b
pi = np.pi

t = np.linspace(0, 1, 8192)
signal = 3 * cos(2 * pi * 5 * t + pi/3) + 7 * cos(2 * np.pi * 3000 * t + 3/8*pi)

fir_len = len(b)

res = []
for i in range(len(signal)):
    print(i)
    temp = 0
    if i < fir_len:
        data = (fir_len-i) * [0] + signal[0:i].tolist()
    elif i > len(signal) - fir_len:
        data = signal[i:].tolist() + (fir_len - len(signal) + i) * [0]
    else:
        data = signal[i:i+fir_len]
    for j in range(fir_len):
        temp += b[j] * data[-j-1]
    res.append(temp)


plt.figure()
plt.plot(t, res)
plt.show()
声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • 滤波器
    +关注

    关注

    161

    文章

    7804

    浏览量

    178058
  • 低通滤波器
    +关注

    关注

    14

    文章

    474

    浏览量

    47398
  • FIR
    FIR
    +关注

    关注

    4

    文章

    146

    浏览量

    33169
  • FFT变换
    +关注

    关注

    2

    文章

    10

    浏览量

    8755
收藏 人收藏

    评论

    相关推荐

    DSP in FPGA:FIR滤波器

    FIR 滤波器广泛应用于数字信号处理中,主要功能就是将不感兴趣的信号滤除,留下有用信号。##全并行FIR滤波器结构
    发表于 06-27 10:02 9510次阅读

    如何设计低通FIR滤波器

    此示例显示如何设计低通FIR滤波器。这里介绍的许多概念可以扩展到其他响应,如高通,带通等。FIR滤波器被广泛使用,因为它们具有强大的设计算法
    发表于 08-23 10:00

    并行FIR滤波器Verilog设计

    本文将简单介绍FIR滤波器的原理,详细介绍使用Verilog HDL设计并行FIR滤波器的流程和
    发表于 09-25 17:44

    什么是fir数字滤波器 什么叫FIR滤波器

    什么是fir数字滤波器 Part 1: Basics1.1 什么是FIR滤波器?FIR 滤波器
    发表于 01-16 09:42 1.7w次阅读

    高效FIR滤波器的设计与仿真-基于FPGA

    高效FIR滤波器的设计与仿真-基于FPGA 摘要:该文在介绍有限冲激响应(FIR)数字滤波器理论及常见实现方法的基础上,提出了
    发表于 01-16 09:56 1658次阅读
    高效<b class='flag-5'>FIR</b><b class='flag-5'>滤波器</b>的设计与仿真-基于FPGA

    用MATLAB设计FIR滤波器的方法

    用MATLAB设计FIR滤波器的方法 摘  要 介绍了利用MATLAB信号处理工具箱进行FIR滤波器设计的三种方法:程序设计法、FDATo
    发表于 01-16 18:12 1.5w次阅读
    用MATLAB设计<b class='flag-5'>FIR</b><b class='flag-5'>滤波器</b>的方法

    高阶FIR正交镜像滤波器的设计

    本文计论了高阶FIR型正交镜像滤波器的设计问题。根据FIR 正交镜像滤波器 设计的基本原理,将高阶正交镜像滤波器的设计问题转换为单变量的优化
    发表于 08-29 16:16 29次下载
    高阶<b class='flag-5'>FIR</b>正交镜像<b class='flag-5'>滤波器</b>的设计

    fir_滤波器sourc

    fir滤波器的有关资料 fir_滤波器sourc.rar
    发表于 12-14 14:12 24次下载

    基于MATLAB的FIR滤波器设计与滤波

    基于MATLAB的FIR滤波器设计与滤波
    发表于 12-14 22:08 64次下载

    详解FIR滤波器和IIR滤波器的区别

    数字滤波器广泛应用于硬件电路设计,般分为FIR滤波器和IIR滤波器。那么FIR
    发表于 05-03 11:36 20次下载

    FIR滤波器的FPGA设计与实现

    本文针对快速、准确选择参数符合项目要求的滤波器设计方法的目的,通过系统的介绍有限脉冲响应( Finite Impulse Response,FIR滤波器的原理、结构形式以及几种
    发表于 12-21 14:53 14次下载
    <b class='flag-5'>FIR</b><b class='flag-5'>滤波器</b>的FPGA设计与实现

    FPGA的FIR抽取滤波器设计详细教程

    介绍FIR抽取滤波器的工作原理,重点阐述了用XC2V1000实现FIR抽取滤波器的方法,并给出了仿真波形和设计特点。
    发表于 04-19 11:34 2314次阅读
    FPGA的<b class='flag-5'>FIR</b>抽取<b class='flag-5'>滤波器</b>设计详细教程

    FIR滤波器与IIR滤波器到底有什么区别

    你知道FIR滤波器与IIR滤波器的不同点吗?它有有什么特点?滤波器是工程师工作中必不可少的器件,滤波器分为很多种,本文详细
    发表于 08-09 14:15 3.3w次阅读
    <b class='flag-5'>FIR</b><b class='flag-5'>滤波器</b>与IIR<b class='flag-5'>滤波器</b>到底有什么区别

    FIR滤波器和IIR滤波器的区别与联系

    1.根据冲激响应的不同,将数字滤波器分为有限冲激响应(FIR滤波器和无限冲激响应(IIR)滤波器。对于FIR
    的头像 发表于 12-30 23:45 3867次阅读

    IIR滤波器FIR滤波器的区别

    数字滤波器是数字信号处理中最常用的种技术,可以对数字信号进行滤波、降噪、增强等处理,其中最常见的两种数字滤波器是IIR滤波器
    的头像 发表于 06-03 10:21 1.9w次阅读