一种ADC采样算法，中位值平均滤波+递推平均滤波

已于 2024-03-21 07:58:56 修改
阅读量1k 收藏 20
点赞数 1
分类专栏： STM32 外设学习+项目实践 文章标签： 算法
于 2023-11-11 11:11:43 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45817947/article/details/134345497
版权

前言

在实际AD采集场景中，会出现周期性变化和偶然脉冲波动干扰对AD采集的影响

这里使用中位值平均滤波+递推平均滤波的结合

参考前人写好的代码框架，也参考博主GuYH_下面这篇博客，在此基础上稍作修改，写出这篇博客，能应用于实际项目。常用ADC采样数字滤波算法最全汇总！！！【❤️建议收藏❤️】_adc滤波算法_GuYH_的博客-CSDN博客

以后有机会，搞个上位机对比几种算法的优劣。

有纰漏请指出，转载请说明。

学习交流请发邮件 1280253714@qq.com

算法讲解

中位值滤波：采样10次，去掉最大和最小的值，求和再除8，对脉冲波有抑制作用。

递推平均滤波：把连续取N个采样值看成一个队列，队列的长度固定为N，每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据(先进先出原则)， 把队列中的N个数据进行算术平均运算,就可获得新的滤波结果。对周期性干扰有良好的抑制作用，平滑度高；适用于高频振荡的系统。

优点：适用于周期性变化和偶然脉冲波动干扰的AD采集。

缺点：灵敏度低；对偶然出现的脉冲性干扰的抑制作用较差，不适于脉冲干扰较严重的场合。

代码

因为项目用的是比51还low的单片机，所以不能进行太高频的AD采样，否则会大量占用CPU；采样次数也不能设置太高，否则占用大量RAM空间。

每10ms进行一次中位值滤波

每100ms进行一次递推平均滤波

adc.h

#ifndef __ADC_H
#define	__ADC_H

#include "sys.h"


#define AD_Channel_Num  4 	// ADC的通道数
#define AD_Sample_Num   4 	// 采样次数

void AdcInit(void);
u16 ADC_Sample(u8 adch);
void AdcLoopTask(void);

typedef struct {
    u16 V_Chip_0V6;
    u16 V_Bat;
    u16 V_In;
    u16 I_Cur;
    u8  sampleCnt;                      //采样计数
    u8  bSampleInit;					// 去除前几次的采样值的标志位，0为去除，1为开始滤波
    u8  u8ArraryIndex;					// 当前采样的索引			
	u16 u16SampleValue[AD_Channel_Num];	// 规则转换时DMA搬运的目标数组
	u16 u16ValueSum[AD_Channel_Num];	// ADC采样值的总和
	u16	u16Value[AD_Channel_Num];		// 平均滤波后的AD值	
} ADC_S;

extern ADC_S stAdc;

#endif	//__ADC_H

adc.c

#include "sys.h"

ADC_S stAdc;
/********************************************
*	@函数名	AdcInit
*	@描述	ADC循环任务
*	@参数	无
*	@返回值	无
*	@注意	无
********************************************/
void AdcInit(void)
{
	TRISA0 = 1;		//将RA0设置为输入
	TRISA2 = 1;		//将RA2设置为输入
	TRISA5 = 1;		//将RA5设置为输入
}

/********************************************
*	@函数名	AD_Sample
*	@描述	AD采样
*	@参数	adch - 检测通道
*	@返回值	ad_result - 8次AD平均值 
*	@注意	采样通道需自行设置为模拟口，采样10次,取中间八次的平均值为采样结果存于adresult中
********************************************/
u16 ADC_Sample(u8 adch)
{
	u32 adsum = 0;
	u16 admin = 4096,admax = 0;
	u8 adtimes = 0;
	u16 ad_temp,adresult;
	u8 j;
	ADCON1 = 0B00000101;			//左对齐,选用VDD 2.0V做AD基准
	ADCON0 = 0X81 | (adch << 2);	//ADCclk = Fosc/32
	asm("nop");
	asm("nop");

	for(j=0;j<10;j++)
	{
		GODONE = 1;						//开始转换
		u8 i = 64;
		
		while(GODONE)
		{
			__delay_us(100);	//延时100us(编译器内置函数)
			if(0 == (--i))		//延时6.4ms仍没有AD转换结束，跳出程序
				return;
		}
		
		ad_temp=(ADRESH<<4)+(ADRESL>>4);	//计算12位AD值
		adresult = ad_temp;
		
		if(ad_temp > admax)
			admax = ad_temp;				//AD采样最大值
		else if(ad_temp < admin)
			admin = ad_temp;				//AD采样最小值
		
		adsum += ad_temp;
	}
	adsum = adsum - admax - admin;
	adresult = adsum >> 3;		//8次平均值作为最终结果
	return adresult;
}

/********************************************
*	@函数名	AdcSampleAllChanel
*	@描述	对所有ADC通道都进行采样
*	@参数	无
*	@返回值	无
*	@注意	无
********************************************/
static void AdcSampleAllChanel(void)
{
	stAdc.u16SampleValue[0] 	= ADC_Sample(0x0F);		//芯片内部0.6V电压
	IO_Vbat_Sample_Ctrl = 1;	
	stAdc.u16SampleValue[1] 	= ADC_Sample(0x00);		//电池电压
	IO_Vbat_Sample_Ctrl = 0;
	stAdc.u16SampleValue[2] 	= ADC_Sample(0x05);		//输入电压
	stAdc.u16SampleValue[3] 	= ADC_Sample(0x02);		//电流采样电压
}

/********************************************
*	@函数名	AdcMoveAverageFilter
*	@描述	ADC递推平均滤波法
			把连续取N个采样值看成一个队列，队列的长度固定为N
			每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据(先进先出原则)
			把队列中的N个数据进行算术平均运算,就可获得新的滤波结果
			优点：对周期性干扰有良好的抑制作用，平滑度高；试用于高频振荡的系统
			缺点：灵敏度低；对偶然出现的脉冲性干扰的抑制作用较差，不适于脉冲干扰较严重的场合
*	@参数	无
*	@返回值	无
*	@注意	无
********************************************/
static void AdcMoveAverageFilter(void)
{
	u8 i = 0;
/******************************* u16ValueSum加上最新的值并剔除最旧的值 *******************************/
	for (i=0; i<AD_Channel_Num; i++)
	{
		stAdc.u16SampleValue[i] &= 0x0FFF;
		stAdc.u16ValueSum[i] += stAdc.u16SampleValue[i];
		stAdc.u16ValueSum[i] -= stAdc.u16ValueArray[i][stAdc.u8ArraryIndex];
	}

/******************************* 最近采样的AD_Sample_Num个数取平均 *******************************/	
	if ( 0 == stAdc.bSampleInit )
	{
		for (i=0; i<AD_Channel_Num; i++)
		{
			stAdc.u16Value[i] = stAdc.u16SampleValue[i];
		}
	} 
	else
	{
		for (i=0; i<AD_Channel_Num; i++)
		{
			stAdc.u16Value[i] = stAdc.u16ValueSum[i]>>2;	//在均值滤波的基础上进行递推均值滤波
		}
	}

/******************************* u16ValueArray添加最新的值 *******************************/	
	for (i=0; i<AD_Channel_Num; i++)
	{
		stAdc.u16ValueArray[i][stAdc.u8ArraryIndex] = stAdc.u16SampleValue[i];
	}	

/******************************* AD游标更新 *******************************/	
    stAdc.u8ArraryIndex++;
    if( stAdc.u8ArraryIndex >= AD_Sample_Num ) {
        stAdc.u8ArraryIndex = 0;
        stAdc.bSampleInit = 1;
    }	
}

typedef struct {
	u16 V_Chip_0V6;
    u16 V_Bat;
    u16 V_In;
    u16 I_Cur;
} ADC_TEST_S;
ADC_TEST_S stAdc1;
/********************************************
*	@函数名	AdcLoopTask
*	@描述	ADC循环任务
*	@参数	无
*	@返回值	无
*	@注意	无
********************************************/
void AdcLoopTask(void)
{
	
	AdcSampleAllChanel();
	stAdc.sampleCnt++;
	if (stAdc.sampleCnt > 10)
	{
		stAdc.sampleCnt = 0;
		AdcMoveAverageFilter();
		stAdc.V_Chip_0V6 = stAdc.u16Value[0];
		stAdc.V_Bat		 = stAdc.u16Value[1];
		stAdc.V_In		 = stAdc.u16Value[2];
		stAdc.I_Cur		 = stAdc.u16Value[3];
	}	
}

改进递推平均滤波占用大量RAM空间的方法

减去的不是队首的值，而是上一次得到的平均值

如果有10个通道，递推平均滤波个数为10个，那么，改进后，可以省去10*10*2=200byte个字节

static void AdcMoveAverageFilter(void)
{
	u8 i = 0;
/******************************* u16ValueSum加上最新的值并剔除最旧的值 *******************************/
	for (i=0; i<AD_Channel_Num; i++)
	{
		stAdc.u16SampleValue[i] &= 0x0FFF;
		stAdc.u16ValueSum[i] += stAdc.u16SampleValue[i];
		stAdc.u16ValueSum[i] -= stAdc.u16Value[i];
	}

/******************************* 最近采样的AD_Sample_Num个数取平均 *******************************/	
	if ( 0 == stAdc.bSampleInit )
	{
		for (i=0; i<AD_Channel_Num; i++)
		{
			stAdc.u16ValueSum[i] 	= stAdc.u16SampleValue[i] * AD_Sample_Num;
			stAdc.u16Value[i] 		= stAdc.u16SampleValue[i];
		}
	} 
	else
	{
		for (i=0; i<AD_Channel_Num; i++)
		{
			stAdc.u16Value[i] = stAdc.u16ValueSum[i]>>3;	//在均值滤波的基础上进行递推均值滤波
		}
	}

    if (stAdc.bSampleInit == 0) {
		stAdc.bSampleInit = 1;
	}	
}

这里我采集了4个通道，递推平均滤波个数为4个（8个在改进前RAM溢出了 哭笑 ）

改进前

改进后：

省去了4*4*2=32byte个字节的空间

递推平均滤波的直观体现

随机噪声

我用Excel模拟了一个滤波前的随机噪声，滤波后的值为当前值与其前七个值的平均数，权重为1/8。

规则序列

这里我选取了36组数据，生成sin从0到360度的数据，可以看到，滤波后的数据出现了相位滞后。且递推平均滤波个数越多，相位滞后越明显。

规则序列出现噪声后的滤波效果

sin（180°）时本来是0，这里我改成0.5，可以看出滤波后有一定的抑制作用，但是180°到250°的数据整体被抬高了一点。

改进型递推平均滤波

TianYaKe-天涯客
关注
  • 1
    点赞
  • 踩
  • 20
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
单片机上常用的10个滤波算法
m0_56770931的博客
01-21 1191
一、限幅滤波 1）方法 根据经验判断两次采样允许的最大偏差值A。 每次采新值时判断： 若本次值与上次值之差<=A，则本次有效； 若本次值与上次值之差>A，本次无效，用上次值代替本次。 2）优缺点 克服脉冲干扰，无法抑制周期性干扰，平滑度差。 3）代码 /* A值根据实际调，Value有效值，new_Value当前采样值，程序返回有效的实际值*/#defineA10charValue;char filter(){charnew...
通用软件滤波算法-递推平均滤波法
Hsy的博客
07-04 5879
一、处理思想 回顾一下算术平均滤波法的处理思想：采集N组数据求算术平均值作为单次数据的有效值。 假设N为10，采集某组数据时，若第2次采集时环境已发生较大变化，而我们需要再采集8次数据，才能将数据变化反映到有效值上。可以看出改算法的实时性较差。 而当N取不同值时，大致有一下规律： N↑ 实时性↓ 滤波平滑度↑ 为了解决实时性问题，有一种改进算法：递推平均滤波法 其处理...
AD处理-递推中位值平均一阶滤波
02-17
单路或者n路AD读取，m次递推中位值平均一阶滤波，n、m值可设置。
ADC常用的十大滤波算法(C语言）
最新发布
willhu2008的专栏
06-20 770
根据经验判断两次采样允许的最大偏差值（设为A）. 如果本次值与上次值之差<=A，则本次值有效. 如果本次值与上次值之差>A，则本次值无效，放弃本次值，用上次值代替本次值能有效克服因偶然因素引起的脉冲干扰无法抑制那种周期性的干扰平滑度差。
递推中位平均值滤波算法(C语言)
qq_41861711的博客
05-12 499
介绍了递推中位平均值滤波算法
C语言十大滤波算法
wzk456的专栏
03-15 3425
1、限幅滤波法 方法：根据经验判断，确定两次采样允许的最大偏差值（设为A），每次检测到新值时判断：如果本次值与上次值之差＜=A,则本次值有效。如果本次值与上次值之差＞A,则本次值无效,放弃本次值,用上次值代替本次值 *函数名称：AmplitudeLimiterFilter()-限幅滤波法 *优点：能有效克服因偶然因素引起的脉冲干扰 *缺点：无法抑制那种周期性的干扰，且平滑度差 *说明： 1、调用函数 GetAD(),该函数用来取得当前值 2、变量说明 Value:最近一次有效...
中位值平均滤波法
weixin_47225948的博客
09-21 1206
中位值平均滤波法是信号处理中的一种滤波方法，其原理是将信号中每个采样点周围一定窗口大小内的值取中位数，然后用中位数替代原来采样点的值。这样可以有效地去除随机噪声、粗糙噪声和冲击噪声。中位值平均滤波法的原理可以简单概括如下：将信号中每个采样点周围的一定窗口大小内的值排序；取排序后的中间值作为该采样点的新值；重复以上步骤，将所有采样点都经过中位值平均滤波。中位值平均滤波法的优点是能够有效地去除噪声，且不会造成太大的延迟。
卡尔曼滤波及其改进算法、中位值平均滤波(防脉冲干扰平均滤波法)matlab程序.zip
06-26
卡尔曼滤波，扩展卡尔曼滤波，无迹卡尔曼滤波程序matlab，注释详细！适合新手！中位值平均滤波(防脉冲干扰平均滤波法)纯属看心情附赠哈哈，单独需要的可以看我的博文https://blog.csdn.net/weixin_42859360/article/details/105695648，里面有代码直接copy！
滤波算法 | 递推平均滤波法
Systemcall名屋
06-11 6951
基本方法：采用队列作为测量数据存储器 , 设队列的长度为 N ,每进行一次测量 ,把测量结果放于队尾 ,而扔掉原来队首的一个数据 ,这样在队列中始终就有 N 个 “最新” 的数据。当计算平均值时 ,只要把队列中的 N 个数据进行算数平均 ,就可得到新的算数平均值。这样每进行一次测量 ,就可得到一个新的算术平均值。 其程序如下所示： #include <stdio.h> // 定...
递推平均滤波法
weixin_47225948的博客
09-21 1341
递推平均滤波法是一种数字信号处理方法，其原理是通过输入信号的平均值来减少噪声的影响。该方法通过维护一个固定长度的缓冲区，每次新来一个数据时，将最老的数据出队并将新数据入队，再计算当前缓冲区内所有数据的平均值并输出。这个输出值就是经过平滑处理后的信号输出。递推平均滤波法的优点是占用的内存较小，且能够实现实时平滑处理，但缺点是滤波效果受到缓冲区长度的限制，当长度太短时，噪声减少的效果不太显著；当长度太长时，信号的变化被平滑处理的过度，可能导致一些细节信息的丢失。因此，选择合适的缓冲区长度是关键。
【量亿数据-金融数据】几种经典的滤波算法
weixin_37997007的博客
07-25 1349
1、限幅滤波法（又称程序判断滤波法） A、方法： 根据经验判断，确定两次采样允许的最大偏差值（设为A） 每次检测到新值时判断： 如果本次值与上次值之差 如果本次值与上次值之差>A，则本次值无效，放弃本次值，用上次值代替本次值 B、优点： 能有效克服因偶然因素引起的脉冲干扰 C、缺点 无法抑制那种周期性的干扰 平滑度差 2、中位值滤波法 A、方法： 连续采样N次（N取奇
基于STM32单片机的ADC采样滤波程序限幅滤波法中位值滤波法 算术平均滤波法等十大滤波算法源代码.zip
06-18
基于STM32单片机的ADC采样滤波程序限幅滤波法中位值滤波法 算术平均滤波法等十大滤波算法源代码： 方法一：限幅滤波法 方法：根据经验判断，确定两次采样允许的最大偏差值（设为A），每次检测到新值时判断： 如果...
ADC滤波算法.rar
05-03
一、该模块包含滤波算法有：中位值滤波、中位值平均滤波、递推平均滤波、一阶滞后滤波。用户可根据项目不同情况选用不同的滤波算法。 1.1、中位值滤波：连续采样N次(N取奇数)，把N次采样值按大小排列，取中间值为...
11种AD采样的软件滤波方法
04-20
中位值平均滤波法是一种融合了中位值滤波法和算术平均滤波法的滤波方法，该方法连续采样N个数据，去掉一个最大值和一个最小值，然后计算N-2个数据的算术平均值。该方法的优点是融合了两种滤波法的优点，对于偶然出现...
10种软件滤波法
08-03
2. 中位值滤波法：连续采样N次后，将这些值排序，选取中间值作为有效值。这种滤波方式对偶发性波动干扰有良好效果，适用于慢变参数，如温度和液位，但不适用于快速变化的参数。 3. 算术平均滤波法：连续N次采样的...
ADC 转换_软件滤波
10-27
本文将介绍三种常见的软件滤波技术：限幅滤波法、中位值滤波法和算术平均滤波法。 1. **限幅滤波法**： 限幅滤波法基于设定的偏差阈值A来判断采样值的有效性。当连续两次采样的差值不超过A时，新值被认为是有效的...
滤波算法（二）—— 中位值滤波算法
热门推荐
bigmarshal的博客
03-01 1万+
1、算法介绍 中位值滤波算法的实现方法是采集N个周期的数据，去掉N个周期数据中的最大值和最小值，取剩下的数据的平均值。中位值滤波算法特别适用于会偶然出现异常值的系统。中位值滤波算法应用比较广泛，比如用于一些比赛的评分，经常是去掉一个最高分去掉一个最低分，将其他评分取平均值作为选手的最终得分。 优点：相比于平均值滤波算法，中位值滤波算法能够有效滤除偶然的脉冲干扰。 缺点：与平均...
加权递推平均滤波法
weixin_47225948的博客
09-21 1797
加权递推平均滤波法是一种数字滤波方法，它的原理是将当前样本值与前几个样本值加权平均得到滤波结果。具体来说，假设当前样本值为x(n)，前k个样本值为x(n-1), x(n-2),…,x(n-k)，则加权递推平均滤波法的计算公式为：其中，y(n)为滤波后的结果，a是滤波系数，通常取值为0.8~0.99之间。根据计算公式可以看出，滤波结果是当前样本值和前一次滤波结果的加权平均。由于加权递推平均滤波法考虑了多个样本值的加权平均，因此对于周期性噪声有较好的抑制效果。
10种常用滤波算法特点及例程
罗伯特祥的博客
11-16 8515
文章目录1. 限幅滤波法（1）原理（2）特点（3）例程2. 中位值滤波法（1）原理（2）特点（3）例程3. 算术平均滤波法（1）原理（2）说明（3）特点（4）例程4. 一阶滞后滤波法（1）原理（2）特点（3）例程5. 限幅消抖滤波法（1）原理（2）特点（3）例程 1. 限幅滤波法 （1）原理 通过给定的最大偏差值A来实现，如果本次值与上次差值<=A，则本次有效；否则无效，使用上次值代替。 （2）特点 优点：能有效克服偶然因素导致的脉冲干扰 缺点：无法抑制周期性干扰，平滑性差 （3）例程 int F
单片机ADC采样算法----递推平均值采样法程序
06-02
下面是使用递推平均值采样法实现单片机ADC采样的示例程序： ``` #define N 10 // 采样的个数 unsigned int adc_buf[N]; // 保存采样值的数组 unsigned int adc_average(void) { unsigned int sum = 0; unsigned char i; // 获取N个采样值并累加 for(i = 0; i < N; i++) { adc_buf[i] = ADC_read(); // 获取采样值 sum += adc_buf[i]; // 累加 } // 计算平均值 for(i = 0; i < N - 1; i++) { adc_buf[i] = adc_buf[i+1]; // 向前移位 } adc_buf[N-1] = ADC_read(); // 获取最新的采样值 sum += adc_buf[N-1]; // 累加 return sum / N; // 计算平均值 } ``` 其中，ADC_read()为获取ADC采样值的函数，需要根据具体的单片机型号和硬件接口进行实现。该程序中，通过数组保存之前N个采样值，并通过向前移位的方式实现了递推平均值的计算。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值