ADC常用的十大滤波算法(C语言)

news2024/10/7 10:25:32

 

目录

​编辑

​一、限幅滤波法

二、中位值滤波法

三、算术平均滤波法

四、递推平均滤波法

五、中位值平均滤波法

六、限幅平均滤波法

七、一阶滞后滤波法

八、加权递推平均滤波法

九、消抖滤波法

十、限幅消抖滤波法


​一、限幅滤波法

1、方法:

    • 根据经验判断两次采样允许的最大偏差值(设为A)

    • 每次检测到新值时判断:

a. 如果本次值与上次值之差<=A,则本次值有效

b. 如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值

2、优点:

    • 能有效克服因偶然因素引起的脉冲干扰

3、缺点

    • 无法抑制那种周期性的干扰

    • 平滑度差

 

/* A值根据实际调,Value有效值,new_Value当前采样值,程序返回有效的实际值 */
#define A 10
char Value;
char filter()
{
    char new_Value;
    new_Value = get_ad(); // 获取采样值
    if( abs(new_Value - Value) > A)   
        return Value;     // abs()取绝对值函数
    return new_Value;
}

二、中位值滤波法

1、方法:

    • 连续采样N次(N取奇数),把N次采样值按大小排列

    • 取中间值为本次有效值

2、优点:

    • 能有效克服因偶然因素引起的波动干扰

    • 对温度、液位的变化缓慢的被测参数有良好的滤波效果

3、缺点:

    • 对流量、速度等快速变化的参数不宜

 

#define N 11
char filter()
{
    char value_buf[N];
    char count, i, j, temp;
    for(count = 0; count < N; count ++) //获取采样值
    {
        value_buf[count] = get_ad();
        delay();
    }
    for(j = 0; j < (N-1); j++)
    {
        for(i = 0; i < (n-j); i++)
        {
            if(value_buf[i] > value_buf[i+1])
            {
                temp = value_buf[i];
                value_buf[i] = value_buf[i+1];
                value_buf[i+1] = temp;
            }
        }
    }
    return value_buf[(N-1)/2];
}

三、算术平均滤波法

1、方法:

    • 连续取N个采样值进行算术平均运算

    • N值较大时:信号平滑度较高,但灵敏度较低

    • N值较小时:信号平滑度较低,但灵敏度较高

    • N值的选取:一般流量,N=12;压力:N=4

2、优点:

    • 适用于对一般具有随机干扰的信号进行滤波

    • 这样信号的特点是有一个平均值,信号在某一数值范围附近上下波动

3、缺点:

    • 对于测量速度较慢或要求数据计算速度较快的实时控制不适用

    • 比较浪费RAM

 

#define N 12
char filter()
{
    int sum = 0;
    for(count = 0; count < N; count++)
    {
        sum += get_ad();
    } 
    return (char)(sum/N);
}

四、递推平均滤波法

1、方法:

    • 把连续取N个采样值看成一个队列

    • 队列的长度固定为N

    • 每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则)

    • 把队列中的N个数据进行算术平均运算,就可获得新的滤波结果

    • N值的选取:流量,N=12;压力:N=4;液面,N=4 ~ 12;温度,N=1 ~ 4

2、优点:

    • 对周期性干扰有良好的抑制作用,平滑度高

    • 适用于高频振荡的系统

3、缺点:

    • 灵敏度低

    • 对偶然出现的脉冲性干扰的抑制作用较差

    • 不易消除由于脉冲干扰所引起的采样值偏差

    • 不适用于脉冲干扰比较严重的场合

    • 比较浪费RAM

 

#define N 10

u16 value_buf[N];

u16 sum=0;

u16 curNum=0;

u16 moveAverageFilter()

{

    if(curNum < N)

    {

        value_buf[curNum] = getValue();

        sum += value_buf[curNum];

        curNum++;

        return sum/curNum;

    }

    else

    {

        sum -= sum/N;

        sum += getValue();

        return sum/N;

    }

}

五、中位值平均滤波法

1、方法:

    • 相当于“中位值滤波法”+“算术平均滤波法”

    • 连续采样N个数据,去掉一个最大值和一个最小值

    • 然后计算N-2个数据的算术平均值

    • N值的选取:3~14

2、优点:

    • 融合了两种滤波法的优点

    • 对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差

3、缺点:

    • 测量速度较慢,和算术平均滤波法一样

    • 比较浪费RAM

 

char filter()
{
    char count, i, j;
    char Value_buf[N];
    int sum = 0;
    for(count = 0; count < N; count++)
    {
        Value_buf[count] = get_ad();
    } 
    for(j = 0; j < (N-1); j++)
    {
        for(i = 0; i < (N-j); i++)
        {
            if(Value_buf[i] > Value_buf[i+1])
            {
                temp = Value_buf[i];
                Value_buf[i] = Value_buf[i+1];
                Value_buf[i+1] = temp;
            }
        }  
    }    
    for(count = 1; count < N-1; count ++)
    {
        sum += Value_buf[count];
    }
    return (char)(sum/(N-2));
}

六、限幅平均滤波法

1、方法:

    • 相当于“限幅滤波法”+“递推平均滤波法”

    • 每次采样到的新数据先进行限幅处理,

    • 再送入队列进行递推平均滤波处理

2、优点:

    • 融合了两种滤波法的优点

    • 对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差

3、缺点:

    • 比较浪费RAM

 

#define A 10
#define N 12
char value, i = 0;
char value_buf[N];
char filter()
{
    char new_value, sum = 0;
    new_value = get_ad();
    if(Abs(new_value - value) < A)  
        value_buf[i++] = new_value;
    if(i==N)  
        i=0;
    for(count = 0; count < N; count++)
    {
        sum += value_buf[count];
    }
    return (char)(sum/N);
}

七、一阶滞后滤波法

1、方法:

    • 取a=0~1

    • 本次滤波结果=(1-a)本次采样值+a上次滤波结果

2、优点:

    • 对周期性干扰具有良好的抑制作用

    • 适用于波动频率较高的场合

3、缺点:

    • 相位滞后,灵敏度低

    • 滞后程度取决于a值大小

    • 不能消除滤波频率高于采样频率的1/2的干扰信号

 

/*为加快程序处理速度,取a=0~100*/
#define a 30
char value;
char filter()
{
    char new_value;
    new_value = get_ad();
    return ((100-a)*value + a*new_value);
}

八、加权递推平均滤波法

1、方法:

    • 是对递推平均滤波法的改进,即不同时刻的数据加以不同的权

    • 通常是,越接近现时刻的数据,权取得越大。

    • 给予新采样值的权系数越大,则灵敏度越高,但信号平滑度越低

2、优点:

    • 适用于有较大纯滞后时间常数的对象

    • 和采样周期较短的系统

3、缺点:

    • 对于纯滞后时间常数较小,采样周期较长,变化缓慢的信号

    • 不能迅速反应交易系统当前所受干扰的严重程度,滤波效果差

 

/* coe数组为加权系数表 */
#define N 12
char code coe[N] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
char code sum_coe = {1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9 + 10 + 11 + 12};
char filter()
{
    char count;
    char value_buf[N];
    int sum = 0;
    for(count = 0; count < N; count++)
    {
        value_buf[count] = get_ad();
    }
    for(count = 0; count < N; count++)
    {
        sum += value_buf[count] * coe[count];
    } 
    return (char)(sum/sum_coe);
}

九、消抖滤波法

1、方法:

    • 设置一个滤波计数器

    • 将每次采样值与当前有效值比较:

    • 如果采样值=当前有效值,则计数器清零

    • 如果采样值>或<当前有效值,则计数器+1,并判断计数器是否>=上限N(溢出)

    • 如果计数器溢出,则将本次值替换当前有效值,并清计数器

2、优点:

    • 对于变化缓慢的被测参数有较好的滤波效果,

    • 可避免在临界值附近控制器的反复开/关跳动或显示器上数值抖动

3、缺点:

    • 对于快速变化的参数不宜

    • 如果在计数器溢出的那一次采样到的值恰好是干扰值,则会将干扰值当作有效值导入交易系统

 

#define N 12
char filter()
{
    char count = 0, new_value;
    new_value = get_ad();
    while(value != new_value)
    {
        count++;
        if(count >= N) 
            return new_value;
        new_value = get_ad();
    }
    return value;
}

十、限幅消抖滤波法

1、方法:

    • 相当于“限幅滤波法”+“消抖滤波法”

    • 先限幅,后消抖

2、优点:

    • 继承了“限幅”和“消抖”的优点

    • 改进了“消抖滤波法”中的某些缺陷,避免将干扰值导入系统

3、缺点:

    • 对于快速变化的参数不宜

 

#define A 10
#define N 12
char value;
char filter()
{
    char new_value, count = 0;
    new_value = get_ad();
    while(value != new_value)
    {
        if(Abs(value - new_value) < A)
        {
            count++;
            if(count >= N) 
                return new_value;
            new_value = get_ad();
        }
        return value;
    }
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/569292.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

系列1——geotools简介

geotools简介 官网 https://docs.geotools.org/latest/userguide/index.html 架构图 常用功能列表 ModulePurposegt-renderImplements of Java2D rendering engine to draw a mapgt-jdbcImplements for accessing spatial databasegt-mainImplements for accessing spatial…

搭建一个Fabric网络

生成相关证书文件 首先Fabric提供了一些工具用于生成所需的文件&#xff0c;cryptogen模块会根据提供的配置文件生成后续模块运行过程中需要的证书和数据文件。 cryptogen常用命令&#xff1a; generate用于根据配置文件生成证书。 showtemplate用于提供一个模板。 具体操作如…

ROS学习——rviz打开bag文件

一、首先可以在这个网站中下载.bag文件用于学习。 二、显示.bag文件信息 在终端1中启动&#xff1a; roscore 在终端2中输入&#xff1a; rosbag info <bag_file_name> 这里要把<bag_file_name>换成你自己的bag文件名字&#xff0c;之后在终端就会显示出bag文件…

Lidar AI Solution环境配置

目录 Lidar AI Solution环境配置前言1. Lidar AI Solution1.1 Pipeline overview1.2 GetStart 2. CUDA-BEVFusion2.1 3D目标检测(nuScenes验证集)2.2 演示2.3 模型和数据2.4 前置条件2.5 快速开始推理2.5.1 下载模型和数据到CUDA-BEVFusion文件夹2.5.2 配置environment.sh2.5.3…

力扣贪心算法专题(三)力扣题 452、435、763、56、738、968、714 思路及C++实现

文章目录 贪心算法452. 用最少数量的箭引爆气球435. 无重叠区间做法1 右边界排序 不重叠区间做法2 右边界排序 不重叠区间做法3 左边界排序 重叠区间 763.划分字母区间做法1做法2 56. 合并区间738.单调递增的数字暴力解法贪心算法 968.监控二叉树714.买卖股票的最佳时机含手续费…

CSDN周赛55期 - 简单分析

上期周赛被判抄袭了。。。无语。我从第一期周赛一步步走来&#xff0c;所有题我都做过&#xff0c;我还需要抄袭&#xff1f; 虽然第一时间去申诉&#xff0c;但貌似并没什么用。算了&#xff0c;C站的审核也就这样了&#xff0c;失望寒心。 本期还是《计算之魂》主题赛——不得…

JavaScript 基础 DOM (四)

正则表达式正则表达式 正则基本使用 定义规则 const reg /表达式/其中/ /是正则表达式字面量正则表达式也是对象 使用正则 test()方法 用来查看正则表达式与指定的字符串是否匹配 如果正则表达式与指定的字符串匹配 &#xff0c;返回true&#xff0c;否则false reg.test(…

eclipse环境配置

eclipse环境配置 0. 前言1. 下载eclipse2. 下载sdk3. 配置sdk环境变量4. 验证sdk其他 0. 前言 本节记录如何配置eclipse的sdk配置 操作系统&#xff1a;Windows10 专业版 开发环境&#xff1a;eclipse 1. 下载eclipse 来这边下载一下eclipse 菜鸟教程 完事之后打开eclips…

回归测试概念和4种回归测试策略——你想知道的都在这里啦!

前言&#xff1a; 回归测试是指修改了旧代码后&#xff0c;重新进行测试以确认修改没有引入新的错误或导致其他代码产生错误。自动回归测试将大幅降低系统测试、维护升级等阶段的成本。回归测试作为软件生命周期的一个组成部分&#xff0c;在整个软件测试过程中占有很大的工作…

Tomcat源码:Acceptor与Poller、PollerEvent

参考资料&#xff1a; 《Tomcat源码解析系列&#xff08;十一&#xff09;ProtocolHandler》 《Tomcat源码解析系列&#xff08;十二&#xff09;NioEndpoint》 前文&#xff1a; 《Tomcat源码&#xff1a;启动类Bootstrap与Catalina的加载》 《Tomcat源码&#xff1a;容器…

绝对不能错过的7个零基础免费的ChatGPT镜像网站

还在为打不开openai官网烦心&#xff1f;本文帮你实现ChatGPTMidJourney自由(&#xffe3;∇&#xffe3;)/ &#x1f4d2;收集了一些截至目前(2023年5月25日午12:00)可以免费访问&#xff0c;并且零基础也能正常使用的镜像网站&#xff0c;后续将持续维护更新(&#xff61;&a…

(2)lordPE脱壳

1.寻找OEP 1.1第一种情况 1.修改OD选项&#xff0c;调试设置"事件"为系统断点&#xff0c;直接打开"查看"->"内存",设置00400000下F2断点&#xff0c;单步F8找到0040****开头的OEP例如&#xff1a;00401528 1.2第二种情况 进入od后如果直…

Github疯传!200本经典计算机书籍!

好书在精不在多&#xff0c;每一本经典书籍都值得反复翻阅&#xff0c;温故而知新&#xff01; 下面分享几本计算机经典书籍&#xff0c;都是我自己看过的。 重构 改善既有代码的设计 就像豆瓣评论所说的&#xff0c;看后有种醍醐灌顶、欲罢不能的感觉。无论你是初学者&#…

数据结构与算法02:数组和链表

目录 【数组】 为什么Go语言的切片是成倍扩容&#xff1f; 【链表】 单链表 循环链表 双向链表 双向循环链表 数组和链表如何选择&#xff1f; 如何使用链表实现 LRU 缓存淘汰算法&#xff1f; 链表的一些操作 【每日一练】 【数组】 数组&#xff08;Array&#…

JMeter参数化四种实现方式

1 参数化释义 什么是参数化&#xff1f;从字面上去理解的话&#xff0c;就是事先准备好数据&#xff08;广义上来说&#xff0c;可以是具体的数据值&#xff0c;也可以是数据生成规则&#xff09;&#xff0c;而非在脚本中写死&#xff0c;脚本执行时从准备好的数据中取值。 参…

Sentinel热点key

1.基本介绍 官方文档 何为热点&#xff1f;热点即经常访问的数据。很多时候我们希望统计某个热点数据中访问频次最高的 Top K 数据&#xff0c;并对其访问进行限制。比如&#xff1a; 商品 ID 为参数&#xff0c;统计一段时间内最常购买的商品 ID 并进行限制用户 ID 为参数&…

HyDE、UDAPDR(LLM大模型用于信息检索)

本篇博文继续整理LLM在搜索推荐领域的应用&#xff0c;往期文章请往博主主页查看更多。 Precise Zero-Shot Dense Retrieval without Relevance Labels 这篇文章主要做zero-shot场景下的稠密检索&#xff0c;通过借助LLM的力量不需要Relevance Labels&#xff0c;开箱即用。作…

【数据结构】---几分钟简单几步学会手撕链式二叉树(上)

文章目录 前言&#x1f31f;一、二叉树链式结构的实现&#x1f30f;1.1 前置说明&#x1f4ab;快速创建一棵简单的二叉树 &#x1f30f;1.2 二叉树的遍历的时间、空间复杂度&#x1f30f;1.3 二叉树的遍历&#x1f4ab;1.3.1 前序、中序以及后序遍历&#xff1a;&#x1f4ab;1…

深入理解hashmap底层实现原理

目录 总体介绍 HashMap元素的存储 在hashmap中添加元素 HashMap的扩容机制 HashMap的线程安全性 1.添加和删除元素时存在不安全性 2.进行扩容操作时存在不安全性 3.哈希冲突存在不安全性 4.线程之间的不可见性导致安全问题 总体介绍 HashMap是我们用于元素映射使用频率最…

MySQL——子查询

来一篇 MySQL-子查询 记录一下这个美好的时光,学习记录篇,下文中所有SQL 语句 均可在 MySQL DB 学习Demo 此处下载SQL语句执行,有相关DB 与 表。 1. 需求分析与问题解决 1.1 实际问题 现有解决方式一: SELECT salary FROM employees WHERE last_name = Abel SELECT last…