iir数字滤波器_手把手教系列之一阶数字滤波器设计实现(附代码)
[导读] 前面分享了 IIR/FIR/mean/梳状数字滤波器的具体设计实现,这几种使用起来或许觉得计算量大,相对复杂。实际工程应用中通常有必要过滤来自传感器或音频流的数据,以抑制不必要的噪声。有的应用场景,可能只需要一个最简单的一阶滤波器即可。所以今天来分享一下怎么设计实现一阶数字滤波器。
一阶 RC 滤波?
小伙伴们一定都用过下面这个无源 RC 低通滤波电路:其拉普拉斯模型如下:由于
所以:
其幅频响应为:
由其传递函数可知,这是一个单极点系统,其阻带满足-20dB/10 倍频程斜率下降。其截止频率为:
如把 C/R 交换位置则变成了高通滤波器,其截止频率依然按上式进行计算。这里也分享一个可在线计算的网址给大家:
http://www.elecfans.com/tools/rclvboqijiezhipinlv.html
其通带增益为 0dB。为什么要先谈谈硬件的一阶滤波器呢? 因为这个是大家最为熟悉的东西,而且也一定学过对其进行幅频响应分析。
既然硬件很容易实现一阶低通或者高通滤波器,那么为什么还要讨论一阶数字滤波器呢?
- 硬件滤波器需要 RC 器件,R/C 的规格并不能随意选取,受厂家规格限制,其数值并不连续,特殊规格需要定制
- 数字滤波器非常灵活,一阶数字滤波器计算代价极低。随便一个单片机都可以玩的转。
- 在满足香农采样定理的前提下可灵活实现截止频率。
数字滤波器
这里直接把差分方程列出来,具体推导就不罗嗦了,有兴趣可以找书看看,比较容易:
其中
表示滤波时间常数,T 表示采样周期。
MATLAB 代码
clc; format compacts = tf('s');w = 50; % rad/sH = w/(s+w)T = 1/500;Hd = c2d(H,T,'zoh')opts = bodeoptions;opts.FreqUnits = 'rad/s';opts.XLim = [0.01, 10000];opts.Grid = 'on';bode(H,Hd, opts)
从其响应曲线看为一低通滤波器,相频响应不线性,从其差分方程也看出输出反馈参与运算了,所以其本质是 IIR 滤波器。
上代码
#include #include #include 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 { sim[i]=20*sin(2*PI*10*i/500)+rand()%5; }for(int i=0;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 滤波器,为啥要单列出来介绍一下呢?是因为其实现简单,实际使用时也不必进行复杂的仿真。
1.谨记:知者不惑,仁者不忧,勇者不惧
2.中芯国际透露:14nm或不能为某客户代工
3.数学之美:嵌入式编程凹凸性之妙用(附C代码)
4.MCU是如何从上电复位运行到main函数的?
5.不想在编译程序时浪费时间,这里有妙招!
6.STM32L5中如何关闭TrustZone?
免责声明:本文系网络转载,版权归原作者所有。如涉及作品版权问题,请与我们联系,我们将根据您提供的版权证明材料确认版权并支付稿酬或者删除内容。
iir数字滤波器_手把手教系列之一阶数字滤波器设计实现(附代码)相关推荐
- python代码在线回归中怎么运行_手把手教你用Python进行回归(附代码、学习资料)...
原标题:手把手教你用Python进行回归(附代码.学习资料) 作者: GURCHETAN SINGH翻译:张逸校对:丁楠雅 本文共5800字,建议阅读8分钟. 本文从线性回归.多项式回归出发,带你用P ...
- 怎么用python实现回归_手把手教你用Python进行回归(附代码、学习资料)-阿里云开发者社区...
我刚开始学习数据科学时,第一个接触到的算法就是线性回归.在把这个方法算法应用在到各种各样的数据集的过程中,我总结出了一些它的优点和不足. 首先,线性回归假设自变量和因变量之间存在线性关系,但实际情况却 ...
- 简单有趣的 NLP 教程:手把手教你用 PyTorch 辨别自然语言(附代码)
简单有趣的 NLP 教程:手把手教你用 PyTorch 辨别自然语言(附代码) 雷锋网(公众号:雷锋网)按:本文作者甄冉冉,原载于作者个人博客,雷锋网已获授权. 最近在学pyTorch的实际应用 ...
- python自助电影售票机_手把手教你用python抢票回家过年(代码简单)
首先看看如何快速查看剩余火车票? 当你想查询一下火车票信息的时候,你还在上12306官网吗?或是打开你手机里的APP?下面让我们来用Python写一个命令行版的火车票查看器, 只要在命令行敲一行命令就 ...
- 实战案例,手把手教你构建电商用户画像 | 附代码
导读:本文以真实案例,手把手教你搭建电商系统的用户画像. 先来看该电商用户画像用到的标签. 数据内容包括user_id(用户身份).item_id(商品).IDbehavior_type(用户行为类型 ...
- 英文期刊催稿信模板_手把手教你写投稿信,另附查尔斯沃思高质量模板
导语 本文是查尔斯沃思作者服务关于学术论文写作系列文章的最后一篇,我们邀请英国编辑团队资深成员,根据其自身丰富的撰稿经验,为中国作者呈现系统全面的写作指导建议,我们将其翻译成中文,方便大家理解.希望本 ...
- 手把手教系列之梳状滤波器设计实现
[导读]:前面一篇文章关于IIR/移动平均滤波器设计的文章.本文来聊一聊陷波滤波器,该滤波器在混入谐波干扰时非常有用,算法简单,实现代价低.本文来一探其在机理.应用场景. 注:尽量在每篇文章写写摘要, ...
- application实现网页计数_手把手教你利用爬虫爬网页(Python代码)
本文主要分为两个部分:一部分是网络爬虫的概述,帮助大家详细了解网络爬虫:另一部分是HTTP请求的Python实现,帮助大家了解Python中实现HTTP请求的各种方式,以便具备编写HTTP网络程序的能 ...
- python 接口设计_手把手教你在机器学习过程设计Python接口
前言 为了进行机器学习工程,首先要部署一个模型,在大多数情况下作为一个预测API.为了使此API在生产中工作,必须首先构建模型服务基础设施.这包括负载平衡.扩展.监视.更新等等. 乍一看,所有这些工作 ...
最新文章
- Visual Studio 2019设置回车代码补全
- linux和用户账户相关的系统文件
- HP, That's A Dream Or Not?
- tbase同步mysql_mysql主从同步
- 小程序promise封装post请求_Promise封装微信小程序的Request请求
- python解放二次开发_[转载]Python二次开发程序详解
- 【转载】(Git)用动图展示10大Git命令
- android chrome 不支持 audio/video的autoplay 属性
- [转载] Java static关键字详解
- unity C#修改脚本图标
- 28.java中Java8.0的新特性(附讲解和练习)
- “李记餐厅”微信点餐小程序的设计与实现
- 1050Ti解决csgo打不开、电脑无缘无故蓝屏的终极方法
- 200 一步步开始学习制作PCBPCBA--绘制原理图
- matlab音频信号导入,如何将语音信号导入simulink中
- SSH学习(个人笔记)
- 决策树后剪枝算法(四)最小错误剪枝MEP
- python利用以下公式求π的值_使用Python计算 π 值
- 给图像增加一种噪声并利用至少两种低通或高通滤波器实现频率域的滤波
- usb root hub的启用和设备的识别(1)
热门文章
- BeanShell使用json.jar包处理Json数据
- 国外物联网平台初探(四):Ayla Networks
- java,js,jstl,EL的简单交互
- leetcode 237: Delete Node in a Linked List
- MinGW-notepad++开发c/c++程序
- 忘记mysql 密码的取回方法
- asp页面和Asp.net页面传中文参数UrlEncode编码以及接收解码
- 模仿nginx修改进程名
- MFC系统托盘的实现
- 网页中插入VLC播放器播放rtsp视频流步骤