手把手教系列之一阶数字滤波器设计实现(附代码)

嵌入式客栈 2020-06-12 00:00
注、星标 嵌入式客栈 ,干货及时送达

[导读] 前面分享了 IIR/FIR/mean/梳状数字滤波器的具体设计实现,这几种使用起来或许觉得计算量大,相对复杂。实际工程应用中通常有必要过滤来自传感器或音频流的数据,以抑制不必要的噪声。有的应用场景,可能只需要一个最简单的一阶滤波器即可。所以今天来分享一下怎么设计实现一阶数字滤波器。

一阶 RC 滤波?

小伙伴们一定都用过下面这个无源 RC 低通滤波电路:其拉普拉斯模型如下:由于

所以:

其幅频响应为:

由其传递函数可知,这是一个单极点系统,其阻带满足-20dB/10 倍频程斜率下降。其截止频率为:

如把 C/R 交换位置则变成了高通滤波器,其截止频率依然按上式进行计算。这里也分享一个可在线计算的网址给大家:

http://www.elecfans.com/tools/rclvboqijiezhipinlv.html

其通带增益为 0dB。为什么要先谈谈硬件的一阶滤波器呢? 因为这个是大家最为熟悉的东西,而且也一定学过对其进行幅频响应分析。

既然硬件很容易实现一阶低通或者高通滤波器,那么为什么还要讨论一阶数字滤波器呢?

  • 硬件滤波器需要 RC 器件,R/C 的规格并不能随意选取,受厂家规格限制,其数值并不连续,特殊规格需要定制
  • 数字滤波器非常灵活,一阶数字滤波器计算代价极低。随便一个单片机都可以玩的转。
  • 在满足香农采样定理的前提下可灵活实现截止频率。

数字滤波器

这里直接把差分方程列出来,具体推导就不罗嗦了,有兴趣可以找书看看,比较容易:

其中

表示滤波时间常数,T 表示采样周期。

MATLAB 代码

clc; format compact

s = tf('s');
w = 50% rad/s
H = w/(s+w)

T = 1/500;
Hd = c2d(H,T,'zoh')

opts = bodeoptions;
opts.FreqUnits = 'rad/s';
opts.XLim = [0.0110000];
opts.Grid = 'on';

bode(H,Hd, opts)

从其响应曲线看为一低通滤波器,相频响应不线性,从其差分方程也看出输出反馈参与运算了,所以其本质是 IIR 滤波器。

上代码

#include <stdio.h>
#include <math.h>
#include <string.h>
typedef struct _t_FSTO_FILTER
{

    float yn1;
    float a;
}t_FSTO_FILTER;

int init_first_order_lpf(t_FSTO_FILTER * pFilter,float Tf,float T)
{
    if(T<=0 || Tf<=0)
      return -1;

    pFilter->a   = Tf/(Tf+T);
    pFilter->yn1 = 0;
    return 0;
}
float fist_order_lpf(t_FSTO_FILTER * pFilter,float xn)
{
    float yn;
    yn = pFilter->a*pFilter->yn1+(1-pFilter->a)*xn;
    pFilter->yn1 = yn;
    return yn;
}

#define PI 3.1415f
#define SAMPLE_RATE 500.0f
#define SAMPLE_T (1/SAMPLE_RATE)
#define SAMPLE_SIZE (100)
int main()
{
    float sim[SAMPLE_SIZE];
    float out[SAMPLE_SIZE];
    t_FSTO_FILTER lpf;

    if(init_first_order_lpf(&lpf,0.005,SAMPLE_T)==-1)
      return -1;

    FILE *pFile=fopen("./simulationSin.csv","wt+");
    if(pFile==NULL)
    {
        printf("simulationSin.csv opened failed");
        return -1;
    }

    for(int i=0;i<SAMPLE_SIZE;i++)
    {
        sim[i]=20*sin(2*PI*10*i/500)+rand()%5;
    }
    for(int i=0;i<SAMPLE_SIZE;i++)
    {
        out[i]=fist_order_lpf(&lpf,sim[i]);
        fprintf(pFile,"%f,%f\n",sim[i],out[i]);
    }

    fclose(pFile);
    return 0;
}

取滤波时间常数为 0.005S,采样周期为 0.2S,为 40 倍关系,来看一下上述代码的滤波效果,波形未失真,效果棒棒哒~

如果将常数修改为 0.1S,看下效果:

由图可见,幅度已经衰减,波形已经失真,传递函数的幅频响应已进入衰减区。所以实际使用的时候,滤波器时间常数尽量取小于采样周期 10 倍为宜,具体可以仿真一下,或者类似上面测试程序测试一下为宜。

总结一下

一阶数字滤波计算简单,实现代价非常低。在滤除高频噪声时应用很广泛。其本质是 IIR 滤波器,为啥要单列出来介绍一下呢?是因为其实现简单,实际使用时也不必进行复杂的仿真。

本文辛苦原创分享,如果觉得有价值也请帮忙点赞转发支持,不胜感激!

END

往期精彩推荐,点击即可阅读




▲Linux内核中I2C总线及设备长啥样?  [推荐]
手把手教系列之状滤波器设计实现
手把手教系列之IIR数字滤波器设计实现

嵌入式客栈 欢迎关注嵌入式客栈,主要分享嵌入式Linux系统构建、嵌入式linux驱动开发、单片机技术、FPGA开发、信号处理、工业通讯等技术主题。欢迎关注,一起交流,一起进步!
评论
  • 在智能化技术快速发展当下,图像数据的采集与处理逐渐成为自动驾驶、工业等领域的一项关键技术。高质量的图像数据采集与算法集成测试都是确保系统性能和可靠性的关键。随着技术的不断进步,对于图像数据的采集、处理和分析的需求日益增长,这不仅要求我们拥有高性能的相机硬件,还要求我们能够高效地集成和测试各种算法。我们探索了一种多源相机数据采集与算法集成测试方案,能够满足不同应用场景下对图像采集和算法测试的多样化需求,确保数据的准确性和算法的有效性。一、相机组成相机一般由镜头(Lens),图像传感器(Image
    康谋 2024-12-12 09:45 98浏览
  • RK3506 是瑞芯微推出的MPU产品,芯片制程为22nm,定位于轻量级、低成本解决方案。该MPU具有低功耗、外设接口丰富、实时性高的特点,适合用多种工商业场景。本文将基于RK3506的设计特点,为大家分析其应用场景。RK3506核心板主要分为三个型号,各型号间的区别如下图:​图 1  RK3506核心板处理器型号场景1:显示HMIRK3506核心板显示接口支持RGB、MIPI、QSPI输出,且支持2D图形加速,轻松运行QT、LVGL等GUI,最快3S内开
    万象奥科 2024-12-11 15:42 91浏览
  • 一、SAE J1939协议概述SAE J1939协议是由美国汽车工程师协会(SAE,Society of Automotive Engineers)定义的一种用于重型车辆和工业设备中的通信协议,主要应用于车辆和设备之间的实时数据交换。J1939基于CAN(Controller Area Network)总线技术,使用29bit的扩展标识符和扩展数据帧,CAN通信速率为250Kbps,用于车载电子控制单元(ECU)之间的通信和控制。小北同学在之前也对J1939协议做过扫盲科普【科普系列】SAE J
    北汇信息 2024-12-11 15:45 119浏览
  • 本文介绍瑞芯微RK3588主板/开发板Android12系统下,APK签名文件生成方法。触觉智能EVB3588开发板演示,搭载了瑞芯微RK3588芯片,该开发板是核心板加底板设计,音视频接口、通信接口等各类接口一应俱全,可帮助企业提高产品开发效率,缩短上市时间,降低成本和设计风险。工具准备下载Keytool-ImportKeyPair工具在源码:build/target/product/security/系统初始签名文件目录中,将以下三个文件拷贝出来:platform.pem;platform.
    Industio_触觉智能 2024-12-12 10:27 90浏览
  • 全球智能电视时代来临这年头若是消费者想随意地从各个通路中选购电视时,不难发现目前市场上的产品都已是具有智能联网功能的智能电视了,可以宣告智能电视的普及时代已到临!Google从2021年开始大力推广Google TV(即原Android TV的升级版),其他各大品牌商也都跟进推出搭载Google TV操作系统的机种,除了Google TV外,LG、Samsung、Panasonic等大厂牌也开发出自家的智能电视平台,可以看出各家业者都一致地看好这块大饼。智能电视的Wi-Fi连线怎么消失了?智能电
    百佳泰测试实验室 2024-12-12 17:33 77浏览
  • 习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习笔记&记录学习习笔记&记学习学习笔记&记录学习学习笔记&记录学习习笔记&记录学习学习笔记&记录学习学习笔记记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记
    youyeye 2024-12-12 10:13 55浏览
  • 铁氧体芯片是一种基于铁氧体磁性材料制成的芯片,在通信、传感器、储能等领域有着广泛的应用。铁氧体磁性材料能够通过外加磁场调控其导电性质和反射性质,因此在信号处理和传感器技术方面有着独特的优势。以下是对半导体划片机在铁氧体划切领域应用的详细阐述: 一、半导体划片机的工作原理与特点半导体划片机是一种使用刀片或通过激光等方式高精度切割被加工物的装置,是半导体后道封测中晶圆切割和WLP切割环节的关键设备。它结合了水气电、空气静压高速主轴、精密机械传动、传感器及自动化控制等先进技术,具有高精度、高
    博捷芯划片机 2024-12-12 09:16 93浏览
  • 习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习笔记&记录学习习笔记&记学习学习笔记&记录学习学习笔记&记录学习习笔记&记录学习学习笔记&记录学习学习笔记记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记录学习学习笔记&记
    youyeye 2024-12-11 17:58 94浏览
  • 首先在gitee上打个广告:ad5d2f3b647444a88b6f7f9555fd681f.mp4 · 丙丁先生/香河英茂工作室中国 - Gitee.com丙丁先生 (mr-bingding) - Gitee.com2024年对我来说是充满挑战和机遇的一年。在这一年里,我不仅进行了多个开发板的测评,还尝试了多种不同的项目和技术。今天,我想分享一下这一年的故事,希望能给大家带来一些启发和乐趣。 年初的时候,我开始对各种开发板进行测评。从STM32WBA55CG到瑞萨、平头哥和平海的开发板,我都
    丙丁先生 2024-12-11 20:14 81浏览
  • 时源芯微——RE超标整机定位与解决详细流程一、 初步测量与问题确认使用专业的电磁辐射测量设备,对整机的辐射发射进行精确测量。确认是否存在RE超标问题,并记录超标频段和幅度。二、电缆检查与处理若存在信号电缆:步骤一:拔掉所有信号电缆,仅保留电源线,再次测量整机的辐射发射。若测量合格:判定问题出在信号电缆上,可能是电缆的共模电流导致。逐一连接信号电缆,每次连接后测量,定位具体哪根电缆或接口导致超标。对问题电缆进行处理,如加共模扼流圈、滤波器,或优化电缆布局和屏蔽。重新连接所有电缆,再次测量
    时源芯微 2024-12-11 17:11 123浏览
  • 应用环境与极具挑战性的测试需求在服务器制造领域里,系统整合测试(System Integration Test;SIT)是确保产品质量和性能的关键步骤。随着服务器系统的复杂性不断提升,包括:多种硬件组件、操作系统、虚拟化平台以及各种应用程序和服务的整合,服务器制造商面临着更有挑战性的测试需求。这些挑战主要体现在以下五个方面:1. 硬件和软件的高度整合:现代服务器通常包括多个处理器、内存模块、储存设备和网络接口。这些硬件组件必须与操作系统及应用软件无缝整合。SIT测试可以帮助制造商确保这些不同组件
    百佳泰测试实验室 2024-12-12 17:45 84浏览
我要评论
0
点击右上角,分享到朋友圈 我知道啦
请使用浏览器分享功能 我知道啦