C语言十大滤波算法

合集下载

十大滤波算法

十大滤波算法
9.D、缺点:
10.?? ?对流量、速度等快速变化的参数不宜。
11.E、整理:shenhaiyu 2013-11-01
12.*/
13.?
14.intFilter_Value;
15.?
16.voidsetup(){
17.??Serial.begin(9600);? ?? ?//初始化串口通信
18.??randomSeed(analogRead(0));//产生随机种子
24.}
25.?
26.voidloop(){
27.??Filter_Value = Filter();? ?? ?//获得滤波器输出值
28.??Serial.println(Filter_Value);//串口输出
29.??delay(50);
30.}
31.?
32.//用于随机产生一个300左右的当前值
43.??filter_buf[FILTER_N]= Get_AD();
44.??for(i =0; i < FILTER_N; i++){
45.? ? filter_buf[i]= filter_buf[i +1];//所有数据左移,低位仍掉
46.? ? filter_sum += filter_buf[i];
19.}
20.?
21.voidloop(){
22.??Filter_Value = Filter();? ?? ?//获得滤波器输出值
23.??Serial.println(Filter_Value);//串口输出
24.??delay(50);
25.}
26.?
27.//用于随机产生一个300左右的当前值

C语言十大滤波算法

C语言十大滤波算法

C语言十大滤波算法Company Document number:WTUT-WT88Y-W8BBGB-BWYTT-19998十大滤波算法程序大全(精心整理版)(转自网络)11、限幅滤波法*************************************************** *函数名称:AmplitudeLimiterFilter()-限幅滤波法*优点:能有效克服因偶然因素引起的脉冲干扰*缺点:无法抑制那种周期性的干扰,且平滑度差*说明:1、调用函数GetAD(),该函数用来取得当前值2、变量说明Value:最近一次有效采样的值,该变量为全局变量NewValue:当前采样的值ReturnValue:返回值3、常量说明A:两次采样的最大误差值,该值需要使用者根据实际情况设置*入口:Value,上一次有效的采样值,在主程序里赋值*出口:ReturnValue,返回值,本次滤波结果****************************************************/#defineA10unsignedcharValueunsignedcharAmplitudeLimiterFilter(){unsignedcharNewValue;unsignedcharReturnValue;NewValue=GatAD();if(((NewValue-Value)>A))||((Value-NewValue)>A)))ReturnValue=Value;elseReturnValue=NewValue;return(ReturnValue);}2、中位值滤波法/*****************************************************函数名称:MiddlevalueFilter()-中位值滤波法*优点:能有效克服因偶然因素引起的波动干扰;对温度、液位等变化缓慢的被测参数有良好的滤波效果*缺点:对流量,速度等快速变化的参数不宜*说明:1、调用函数GetAD(),该函数用来取得当前值Delay(),基本延时函数2、变量说明ArrDataBuffer[N]:用来存放一次性采集的N组数据Temp:完成冒泡法试用的临时寄存器i,j,k:循环试用的参数值3、常量说明N:数组长度*入口:*出口:value_buf[(N-1)/2],返回值,本次滤波结果*****************************************************/ #defineN11unsignedcharMiddlevalueFilter(){unsignedcharvalue_buf[N];unsignedchari,j,k,temp;for(i=0;i<N;i++){value_buf[i]=get_ad();delay();}for(j=0;j<N-1;j++){for(k=0;k<N-j;k++){if(value_buf[k]>value_buf[k+1]){temp=value_buf[k];value_buf[k]=value_buf[k+1];value_buf[k+1]=temp;}}}returnvalue_buf[(N-1)/2];}3、算术平均滤波法/*********************************************************说明:连续取N个采样值进行算术平均运算优点:试用于对一般具有随机干扰的信号进行滤波。

十种软件滤波的算法

十种软件滤波的算法

十种软件滤波的算法软件滤波在嵌入式的数据采集和处理中有着很重要的作用,这10种方法各有优劣,根据自己的需要选择。

同时提供了C语言的参考代码,希望对各位能有帮助。

1、限幅滤波法(又称程序判断滤波法)A、方法:根据经验判断,确定两次采样允许的最大偏差值(设为A)每次检测到新值时判断:如果本次值与上次值之差<=A,则本次值有效如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值B、优点:能有效克服因偶然因素引起的脉冲干扰C、缺点无法抑制那种周期性的干扰平滑度差/* A值可根据实际情况调整value为有效值,new_value为当前采样值滤波程序返回有效的实际值*/#define A 10char value;char filter(){char new_value;new_value =get_ad();if ( ( new_value - value > A ) || ( value- new_value > A )return value;returnnew_value;}2、中位值滤波法A、方法:连续采样N次(N取奇数)把N次采样值按大小排列取中间值为本次有效值B、优点:能有效克服因偶然因素引起的波动干扰对温度、液位的变化缓慢的被测参数有良好的滤波效果C、缺点:对流量、速度等快速变化的参数不宜/* N值可根据实际情况调整排序采用冒泡法*/#define N 11char filter(){charvalue_buf[N];charcount,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];}3、算术平均滤波法A、方法:连续取N个采样值进行算术平均运算N值较大时:信号平滑度较高,但灵敏度较低N值较小时:信号平滑度较低,但灵敏度较高N值的选取:一般流量,N=12;压力:N=4B、优点:适用于对一般具有随机干扰的信号进行滤波这样信号的特点是有一个平均值,信号在某一数值范围附近上下波动C、缺点:对于测量速度较慢或要求数据计算速度较快的实时控制不适用比较浪费RAM#define N 12char filter(){int sum = 0;for (count=0;count<N;count++){sum + =get_ad();delay();}return(char)(sum/N);}4、递推平均滤波法(又称滑动平均滤波法)A、方法:把连续取N个采样值看成一个队列队列的长度固定为N每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则)把队列中的N个数据进行算术平均运算,就可获得新的滤波结果N值的选取:流量,N=12;压力:N=4;液面,N=4~12;温度,N=1~4B、优点:对周期性干扰有良好的抑制作用,平滑度高适用于高频振荡的系统C、缺点:灵敏度低对偶然出现的脉冲性干扰的抑制作用较差不易消除由于脉冲干扰所引起的采样值偏差不适用于脉冲干扰比较严重的场合比较浪费RAM#define N 12char value_buf[N];char i=0;char filter(){char count;int sum=0;value_buf[i++] =get_ad();if ( i == N) i = 0;for (count=0;count<N,count++)sum =value_buf[count];return(char)(sum/N);}5、中位值平均滤波法(又称防脉冲干扰平均滤波法)A、方法:相当于“中位值滤波法”+“算术平均滤波法”连续采样N个数据,去掉一个最大值和一个最小值然后计算N-2个数据的算术平均值N值的选取:3~14B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C、缺点:测量速度较慢,和算术平均滤波法一样比较浪费RAM#define N 12char filter(){char count,i,j;charvalue_buf[N];int sum=0;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;}}}for(count=1;count<N-1;count++)sum +=value[count];return(char)(sum/(N-2));}6、限幅平均滤波法A、方法:相当于“限幅滤波法”+“递推平均滤波法”每次采样到的新数据先进行限幅处理,再送入队列进行递推平均滤波处理B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C、缺点:比较浪费RAM/**/ 略参考子程序1、37、一阶滞后滤波法A、方法:取a=0~1本次滤波结果=(1-a)*本次采样值+a*上次滤波结果B、优点:对周期性干扰具有良好的抑制作用适用于波动频率较高的场合C、缺点:相位滞后,灵敏度低滞后程度取决于a值大小不能消除滤波频率高于采样频率的1/2的干扰信号/* 为加快程序处理速度假定基数为100,a=0~100 */#define a 50char value;char filter(){char new_value;new_value =get_ad();return(100-a)*value + a*new_value;}8、加权递推平均滤波法A、方法:是对递推平均滤波法的改进,即不同时刻的数据加以不同的权通常是,越接近现时刻的数据,权取得越大。

C语言十大滤波算法

C语言十大滤波算法

C语言十大滤波算法C语言是一种广泛应用于嵌入式系统、图形界面、游戏开发等领域的编程语言。

在信号处理和图像处理等领域,滤波算法是一种重要的处理方式。

滤波算法可以对信号进行去噪、平滑、边缘检测等操作,从而提高信号的质量和准确度。

在C语言中,有许多优秀的滤波算法被广泛应用。

下面将介绍C语言中的十大滤波算法,并讨论它们的原理和应用领域。

1.均值滤波算法:均值滤波是一种简单有效的滤波算法,通过计算像素周围若干个邻域像素的平均值作为滤波结果。

均值滤波适用于去除高频噪声,但会造成图像细节的模糊。

2.中值滤波算法:中值滤波算法通过计算像素周围若干个邻域像素的中值作为滤波结果。

中值滤波可以有效去除椒盐噪声,但不能处理高斯噪声。

3.高斯滤波算法:高斯滤波算法利用高斯函数对图像进行滤波,以平滑图像并去除噪声。

高斯滤波在保持图像边缘信息的同时,能够有效降低噪声。

4.自适应中值滤波算法:自适应中值滤波算法根据像素邻域内像素的不同情况选择中值滤波器的大小,对不同噪声情况进行适应性处理。

5.双边滤波算法:双边滤波算法是一种非线性滤波算法,通过同时考虑空间信息和灰度差异信息,可在去噪的同时保持图像的边缘信息。

6.快速傅里叶变换(FFT)滤波算法:FFT滤波是一种频域滤波算法,通过将信号从时域转换到频域,对频谱进行滤波后再进行逆变换,能够有效去除周期性噪声。

7.小波变换滤波算法:小波变换是一种时频联合分析方法,将信号分解为不同频率的子带,通过阈值处理可以实现去噪。

8.自适应滤波算法:自适应滤波算法根据图像中的纹理复杂度自动选择合适的滤波器,能够在保持图像细节的同时去除噪声。

9.协同滤波算法:协同滤波算法是一种基于用户行为数据的推荐算法,通过分析用户的历史数据和相似用户群体的数据,对用户进行个性化推荐。

10.卡尔曼滤波算法:卡尔曼滤波算法是一种利用动态模型对状态进行推断的滤波算法,适用于系统状态估计、信号恢复等应用。

以上是C语言中的十大滤波算法,它们在不同领域的应用有所差异,但都能够有效地处理信号和数据,提高数据质量和准确度。

十一种滤波方法及C语言程序

十一种滤波方法及C语言程序

一.十一种通用滤波算法(转)1、限幅滤波法(又称程序判断滤波法)A、方法:根据经验判断,确定两次采样允许的最大偏差值(设为A)每次检测到新值时判断:如果本次值与上次值之差<=A,则本次值有效如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值B、优点:能有效克服因偶然因素引起的脉冲干扰C、缺点无法抑制那种周期性的干扰平滑度差2、中位值滤波法A、方法:连续采样N次(N取奇数)把N次采样值按大小排列取中间值为本次有效值B、优点:能有效克服因偶然因素引起的波动干扰对温度、液位的变化缓慢的被测参数有良好的滤波效果C、缺点:对流量、速度等快速变化的参数不宜3、算术平均滤波法A、方法:连续取N个采样值进行算术平均运算N值较大时:信号平滑度较高,但灵敏度较低N值较小时:信号平滑度较低,但灵敏度较高N值的选取:一般流量,N=12;压力:N=4B、优点:适用于对一般具有随机干扰的信号进行滤波这样信号的特点是有一个平均值,信号在某一数值范围附近上下波动C、缺点:对于测量速度较慢或要求数据计算速度较快的实时控制不适用比较浪费RAM4、递推平均滤波法(又称滑动平均滤波法)A、方法:把连续取N个采样值看成一个队列队列的长度固定为N每次采样到一个新数据放入队尾,并扔掉原来队首的一次数据.(先进先出原则) 把队列中的N个数据进行算术平均运算,就可获得新的滤波结果N值的选取:流量,N=12;压力:N=4;液面,N=4~12;温度,N=1~4 B、优点:对周期性干扰有良好的抑制作用,平滑度高适用于高频振荡的系统C、缺点:灵敏度低对偶然出现的脉冲性干扰的抑制作用较差不易消除由于脉冲干扰所引起的采样值偏差不适用于脉冲干扰比较严重的场合比较浪费RAM5、中位值平均滤波法(又称防脉冲干扰平均滤波法)A、方法:相当于“中位值滤波法”+“算术平均滤波法”连续采样N个数据,去掉一个最大值和一个最小值然后计算N-2个数据的算术平均值N值的选取:3~14B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C、缺点:测量速度较慢,和算术平均滤波法一样比较浪费RAM6、限幅平均滤波法A、方法:相当于“限幅滤波法”+“递推平均滤波法”每次采样到的新数据先进行限幅处理,再送入队列进行递推平均滤波处理B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C、缺点:比较浪费RAM7、一阶滞后滤波法A、方法:取a=0~1本次滤波结果=(1-a)*本次采样值+a*上次滤波结果B、优点:对周期性干扰具有良好的抑制作用适用于波动频率较高的场合C、缺点:相位滞后,灵敏度低滞后程度取决于a值大小不能消除滤波频率高于采样频率的1/2的干扰信号8、加权递推平均滤波法A、方法:是对递推平均滤波法的改进,即不同时刻的数据加以不同的权通常是,越接近现时刻的数据,权取得越大。

10种简单的数字滤波算法(C++源程序)

10种简单的数字滤波算法(C++源程序)

10种简单的数字滤波算法(C++源程序)以下是10种简单的数字滤波算法C++实现示例:1. 均值滤波均值滤波是数字滤波算法的一种常见形式,它可以通过计算一定范围内像素值的平均值来平滑图像。

其C++实现如下:#include <iostream>#include <opencv2/opencv.hpp>using namespace std;using namespace cv;// Function to implement mean filtervoid meanBlur(Mat& img, Mat& result, int k_size){int img_rows = img.rows;int img_cols = img.cols;// Create a same sized blank imageresult.create(img_rows, img_cols, img.type());for(int r=0; r<img_rows; r++){for(int c=0; c<img_cols; c++){// Define the window of radius k_sizeint r_min = max(0, r-k_size/2);int r_max = min(img_rows-1, r+k_size/2);int c_min = max(0, c-k_size/2);int c_max = min(img_cols-1, c+k_size/2);// Calculate the mean valueint sum = 0;int count = 0;for (int i=r_min; i<=r_max; i++){for (int j=c_min; j<=c_max; j++){sum += img.at<uchar>(i,j);count++;}}result.at<uchar>(r,c) = (uchar) (sum/count);}}}int main(int argc, char** argv){// Load the imageMat img = imread("image.jpg", 0);// Check if image is loaded properlyif(!img.data){cout << "Failed to load image" << endl;return -1;}// Define the kernel sizeint k_size = 3;// Apply mean filterMat result;meanBlur(img, result, k_size);// Display the resultnamedWindow("Original Image", WINDOW_NORMAL);namedWindow("Mean Filtered Image", WINDOW_NORMAL);imshow("Original Image", img);imshow("Mean Filtered Image", result);waitKey(0);return 0;}在上述代码中,`meanBlur()` 函数接收一个输入图像`img` 和一个输出图像`result`,以及一个整数参数`k_size`,该参数指定滤波器的大小,即窗口的半径。

十大滤波算法

十大滤波算法
5,305);
35.}
36.?
37.//限幅滤波法(又称程序判断滤波法)
38.#define FILTER_A1
39.intFilter(){
40.??intNewValue;
41.??NewValue = Get_AD();
42.??if(((NewValue - Value)> FILTER_A)||((Value - NewValue)> FILTER_A))
48.? ?? ???filter_buf[i +1]= filter_temp;
49.? ?? ?}
50.? ?}
51.??}
52.??returnfilter_buf[(FILTER_N -1)/2];
53.}
3、算术平均滤波法
ARDUINO代码
1./*
2.A、名称:算术平均滤波法
3.B、方法:
十大滤波算法程序大全
1、限幅滤波法(又称程序判断滤波法)
ARDUINO代码
1./*
2.A、名称:限幅滤波法(又称程序判断滤波法)
3.B、方法:
4.?? ?根据经验判断,确定两次采样允许的最大偏差值(设为A),
5.?? ?每次检测到新值时判断:
6.?? ?如果本次值与上次值之差<=A,则本次值有效,
7.?? ?如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值。
36.??inti, j;
37.??intfilter_temp;
38.??for(i =0; i < FILTER_N; i++){
39.? ? filter_buf[i]= Get_AD();

十一种滤波方法及C语言程序

十一种滤波方法及C语言程序

一. 十一种通用滤波算法(转)1、限幅滤波法(又称程序判断滤波法)A、方法:根据经验判断,确定两次采样允许的最大偏差值(设为A)每次检测到新值时判断:如果本次值与上次值之差<=A, 则本次值有效如果本次值与上次值之差>A, 则本次值无效, 放弃本次值, 用上次值代替本次值B、优点:能有效克服因偶然因素引起的脉冲干扰C缺点无法抑制那种周期性的干扰平滑度差2、中位值滤波法A、方法:连续采样N次(N取奇数)把N 次采样值按大小排列取中间值为本次有效值B、优点:能有效克服因偶然因素引起的波动干扰对温度、液位的变化缓慢的被测参数有良好的滤波效果C缺点:对流量、速度等快速变化的参数不宜3、算术平均滤波法A、方法:连续取N 个采样值进行算术平均运算N 值较大时:信号平滑度较高,但灵敏度较低N 值较小时:信号平滑度较低,但灵敏度较高N值的选取:一般流量,N=12;压力:N=4B、优点:适用于对一般具有随机干扰的信号进行滤波这样信号的特点是有一个平均值,信号在某一数值范围附近上下波动C缺点:对于测量速度较慢或要求数据计算速度较快的实时控制不适用比较浪费RAM4、递推平均滤波法(又称滑动平均滤波法)A、方法:把连续取N个采样值看成一个队列队列的长度固定为N 每次采样到一个新数据放入队尾, 并扔掉原来队首的一次数据.(先进先出原则)把队列中的N个数据进行算术平均运算,就可获得新的滤波结果N值的选取:流量,N=12;压力:N=4;液面,N=4〜12;温度,N=1〜4 B、优点:对周期性干扰有良好的抑制作用,平滑度高适用于高频振荡的系统C缺点:灵敏度低对偶然出现的脉冲性干扰的抑制作用较差不易消除由于脉冲干扰所引起的采样值偏差不适用于脉冲干扰比较严重的场合比较浪费RAM5、中位值平均滤波法(又称防脉冲干扰平均滤波法)A、方法:相当于“中位值滤波法” +“算术平均滤波法”连续采样N个数据,去掉一个最大值和一个最小值然后计算N-2 个数据的算术平均值N 值的选取:3〜14B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C缺点:测量速度较慢,和算术平均滤波法一样比较浪费RAM6、限幅平均滤波法A、方法:相当于“限幅滤波法” +“递推平均滤波法” 每次采样到的新数据先进行限幅处理,再送入队列进行递推平均滤波处理B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C缺点:比较浪费RAM7、一阶滞后滤波法A、方法:取a=0〜1本次滤波结果=(1-a)*本次采样值+a*上次滤波结果B、优点:对周期性干扰具有良好的抑制作用适用于波动频率较高的场合C缺点:相位滞后,灵敏度低滞后程度取决于 a 值大小不能消除滤波频率高于采样频率的1/2 的干扰信号8、加权递推平均滤波法A、方法:是对递推平均滤波法的改进,即不同时刻的数据加以不同的权通常是,越接近现时刻的数据,权取得越大。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

十大滤波算法程序大全(精心整理版)(转自网络)1推荐1、限幅滤波法****************************************************函数名称:AmplitudeLimiterFilter()-限幅滤波法*优点:能有效克服因偶然因素引起的脉冲干扰*缺点:无法抑制那种周期性的干扰,且平滑度差*说明:?1、调用函数????GetAD(),该函数用来取得当前值?2、变量说明????Value:最近一次有效采样的值,该变量为全局变量????NewValue:当前采样的值????ReturnValue:返回值?3、常量说明????A:两次采样的最大误差值,该值需要使用者根据实际情况设置*入口:Value,上一次有效的采样值,在主程序里赋值*出口:ReturnValue,返回值,本次滤波结果****************************************************/#define?A??10unsignedcharValueunsignedcharAmplitudeLimiterFilter(){??unsignedcharNewValue;??unsignedcharReturnValue;??NewValue=GatAD();??if(((NewValue-Value)>A))||((Value-NewValue)>A)))??ReturnValue=Value;??elseReturnValue=NewValue;??return(ReturnValue);}2、中位值滤波法/*****************************************************函数名称:MiddlevalueFilter()-中位值滤波法*优点:能有效克服因偶然因素引起的波动干扰;对温度、液??????位等变化缓慢的被测参数有良好的滤波效果*缺点:对流量,速度等快速变化的参数不宜*说明:?1、调用函数????GetAD(),该函数用来取得当前值????Delay(),基本延时函数?2、变量说明????ArrDataBuffer[N]:用来存放一次性采集的N组数据????Temp:完成冒泡法试用的临时寄存器????i,j,k:循环试用的参数值?3、常量说明????N:数组长度*入口:*出口:value_buf[(N-1)/2],返回值,本次滤波结果*****************************************************/ #defineN11unsignedcharMiddlevalueFilter(){?unsignedcharvalue_buf[N];?unsignedchari,j,k,temp;?for(i=0;i<N;i++)?{???value_buf[i]=get_ad();???delay();?}?for(j=0;j<N-1;j++)?{??for(k=0;k<N-j;k++)??{???if(value_buf[k]>value_buf[k+1])????{??????temp=value_buf[k];??????value_buf[k]=value_buf[k+1];??????value_buf[k+1]=temp;????}??}?}?returnvalue_buf[(N-1)/2];}3、算术平均滤波法/*********************************************************说明:连续取N个采样值进行算术平均运算优点:试用于对一般具有随机干扰的信号进行滤波。

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

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

**********************************************************/#defineN12charfilter(){?unsignedintsum=0;?unsignedchari;???sum+=get_ad();???delay();?}?return(char)(sum/N);}4、递推平均滤波法(又称滑动平均滤波法)/***************************************************说明:把连续N个采样值看成一个队列,队列长度固定为N。

?????每次采样到一个新数据放入队尾,并扔掉队首的一?????次数据。

把队列中的N各数据进行平均运算,既获得?????新的滤波结果。

优点:对周期性干扰有良好的抑制作用,平滑度高;试用于高频振荡的系统缺点:灵敏度低;对偶然出现的脉冲性干扰的抑制作用较差,不适于脉冲干?????扰较严重的场合****************************************************/#defineN12unsignedcharvalue_buf[N];unsignedcharfilter(){?unsignedchari;?unsignedcharvalue;?intsum=0;?value_buf[i++]=get_ad();??????//采集到的数据放入最高位???value_buf[i]=value_buf[i+1];??//所有数据左移,低位扔掉???sum+=value_buf[i];?}?value=sum/N;?return(value);}5、中位值平均滤波法(又称防脉冲干扰平均滤波法)/********************************************?说明:采一组队列去掉最大值和最小值?优点:融合了两种滤波的优点。

对于偶然出现的脉冲性干扰,可消??????除有其引起的采样值偏差。

对周期干扰有良好的抑制作用,??????平滑度高,适于高频振荡的系统。

?缺点:测量速度慢*********************************************/#defineN12ucharfilter(){?unsignedchari,j,k,l;?unsignedchartemp,sum=0,value;?unsignedcharvalue_buf[N],;?for(i=0;i<N;i++)?{???value_buf[i]=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;?????}???}?}?for(i=1;i<N-1;i++)?sum+=value_buf[i];?value=sum/(N-2);?return(value);}6、递推中位值滤波法/************************************************?优点:对于偶然出现的脉冲性干扰,可消除由其引起的采样值偏差。

??????对周期性干扰有良好的抑制作用,平滑度高;试用于高频振荡??????的系统?缺点:测量速度慢*************************************************/ charfilter(charnew_data,charqueue[],charn){?charmax,min;?charsum;?chari;?queue[0]=new_data;?max=queue[0];?min=queue[0];?sum=queue[0];?for(i=n-1;i>0;i--)?{???if(queue[i]>max)???max=queue[i];???elseif(queue[i]<min)???min=queue[i];???sum=sum+queue[i];???queue[i]=queue[i-1];?}?i=n-2;?sum=sum-max-min+i/2;????//说明:+i/2的目的是为了四舍五入?sum=sum/i;?return(sum);}7、限幅平均滤波法/************************************************?优点:对于偶然出现的脉冲性干扰,可消除有其引起的采样值偏差。

*************************************************/#defineA10#defineN12unsignedchardata[];unsignedcharfilter(data[]){?unsignedchari;?unsignedcharvalue,sum;?data[N]=GetAD();?if(((data[N]-data[N-1])>A||((data[N-1]-data[N])>A))?data[N]=data[N-1];?//elsedata[N]=NewValue;?for(i=0;i<N;i++)?{???data[i]=data[i+1];???sum+=data[i];?}?value=sum/N;?return(value);}8、一阶滞后滤波法/**************************************************** *函数名称:filter()-一阶滞后滤波法*说明:?1、调用函数????GetAD(),该函数用来取得当前值????Delay(),基本延时函数?2、变量说明????Or_data[N]:采集的数据????Dr0_flag、Dr1_flag:前一次比较与当前比较的方向位????coeff:滤波系数????F_count:滤波计数器?3、常量说明????N:数组长度????Thre_value:比较门槛值*入口:*出口:*****************************************************/ #defineThre_value?10#define?N??50floatOr_data[N];unsignedcharDr0_flag=0,Dr1_flag=0;voidabs(floatfirst,floatsecond){?floatabs;?if(first>second)?{??abs=first-second;??Dr1_flag=0;?}?else?{??abs=second-first;??Dr1_flag=1;?}?return(abs);}?voidfilter(void){?uchari=0,F_count=0,coeff=0;?floatAbs=0.00;?//确定一阶滤波系数?for(i=1;i<N;i++)???{?????Abs=abs(Or_data[i-1],Or_data[i]);?????if(!(Dr1_flag^Dr0_flag))???????????????????//前后数据变化方向一致?????{????????F_count++;???????if(Abs>=Thre_value)???????{?????????F_count++;?????????F_count++;???????}???????if(F_count>=12)???????F_count=12;???????coeff=20*F_count;???????}?????else???????????????????????????????????????//去抖动?????coeff=5;?????//一阶滤波算法?????if(Dr1_flag==0)????????????????????????????//当前值小于前一个值?????Or_data[i]=Or_data[i-1]-coeff*(Or_data[i-1]-Or_data[i])/256; ?????else?????Or_data[i]=Or_data[i-1]+coeff*(Or_data[i]-Or_data[i-1])/256;??? ???????????F_count=0;?????????????????????????????????//滤波计数器清零?????Dr0_flag=Dr1_flag;???}}9、加权递推平均滤波法/************************************************************?coe:数组为加权系数表,存在程序存储区。

相关文档
最新文档