11种软件滤波方法的示例程序

假定从8位AD中读取数据(如果是更高位的AD可定义数据类型为int),子程序为get_ad();

1、限副滤波

/*  A值可根据实际情况调整

value为有效值，new_value为当前采样值

滤波程序返回有效的实际值  */

#define A 10

char value;

char filter()

{

char  new_value;

new_value = get_ad();

if ( ( new_value - value > A ) || ( value - new_value > A )

return value;

return new_value;

}

2、中位值滤波法

/*  N值可根据实际情况调整

排序采用冒泡法*/

#define N  11

char filter()

{

char value_buf[N];

char count,i,j,temp;

for ( count=0;count

{

value_buf[count] = get_ad();

delay();

}

for (j=0;j

{

for (i=0;i

{

if ( value_buf>value_buf[i+1] )

{

temp = value_buf;

value_buf = value_buf[i+1];

value_buf[i+1] = temp;

}

}

}

return value_buf[(N-1)/2];

}

3、算术平均滤波法

/*

*/

#define N 12

char filter()

{

int  sum = 0;

for ( count=0;count

{

sum + = get_ad();

delay();

}

return (char)(sum/N);

}

4、递推平均滤波法(又称滑动平均滤波法)

/*

*/

#define N 12

char 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

sum = value_buf[count];

return (char)(sum/N);

}

5、中位值平均滤波法(又称防脉冲干扰平均滤波法)

/*

*/

#define N 12

char filter()

{

char count,i,j;

char value_buf[N];

int  sum=0;

for  (count=0;count

{

value_buf[count] = get_ad();

delay();

}

for (j=0;j

{

for (i=0;i

{

if ( value_buf>value_buf[i+1] )

{

temp = value_buf;

value_buf = value_buf[i+1];

value_buf[i+1] = temp;

}

}

}

for(count=1;count

sum += value[count];

return (char)(sum/(N-2));

}

6、限幅平均滤波法

/*

*/

略 参考子程序1、3

7、一阶滞后滤波法

/* 为加快程序处理速度假定基数为100，a=0~100 */

#define a 50

char value;

char filter()

{

char  new_value;

new_value = get_ad();

return (100-a)*value + a*new_value;

}

8、加权递推平均滤波法

/* 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

{

value_buf[count] = get_ad();

delay();

}

for (count=0,count

sum += value_buf[count]*coe[count];

return (char)(sum/sum_coe);

}

9、消抖滤波法

#define N 12

char filter()

{

char count=0;

char new_value;

new_value = get_ad();

while (value !=new_value);

{

count++;

if (count>=N)   return new_value;

delay();

new_value = get_ad();

}

return value;

}

10、限幅消抖滤波法

/*

*/

略 参考子程序1、9

11、IIR滤波例子

int  BandpassFilter4(int InputAD4)

{

int  ReturnValue;

int  ii;

RESLO=0;

RESHI=0;

MACS=*PdelIn;

OP2=1068; //FilterCoeff4[4];

MACS=*(PdelIn+1);

OP2=8;    //FilterCoeff4[3];

MACS=*(PdelIn+2);

OP2=-2001;//FilterCoeff4[2];

MACS=*(PdelIn+3);

OP2=8;    //FilterCoeff4[1];

MACS=InputAD4;

OP2=1068; //FilterCoeff4[0];

MACS=*PdelOu;

OP2=-7190;//FilterCoeff4[8];

MACS=*(PdelOu+1);

OP2=-1973; //FilterCoeff4[7];

MACS=*(PdelOu+2);

OP2=-19578;//FilterCoeff4[6];

MACS=*(PdelOu+3);

OP2=-3047; //FilterCoeff4[5];

*p=RESLO;

*(p+1)=RESHI;

mytestmul<<=2;

ReturnValue=*(p+1);

for  (ii=0;ii<3;ii++)

{

DelayInput[ii]=DelayInput[ii+1];

DelayOutput[ii]=DelayOutput[ii+1];

}

DelayInput[3]=InputAD4;

DelayOutput[3]=ReturnValue;

//  if (ReturnValue<0)

//  {

//  ReturnValue=-ReturnValue;

//  }

return ReturnValue;

}

关于中值滤波算法，以及C语言实现(转)

源:关于中值滤波算法,以及C语言实现 1.什么是中值滤波? 中值滤波是对一个滑动窗口内的诸像素灰度值排序,用其中值代替窗口中心象素的原来灰度值,它是一种非线性的图像平滑法,它对脉冲干扰级椒盐噪声的抑制 ...

SSE图像算法优化系列二十二：优化龚元浩博士的曲率滤波算法，达到约1000 MPixels&sol;Sec的单次迭代速度

2015年龚博士的曲率滤波算法刚出来的时候,在图像处理界也曾引起不小的轰动,特别是其所说的算法的简洁性,以及算法的效果.执行效率等方面较其他算法均有一定的优势,我在该算法刚出来时也曾经有关注,不过 ...

基于Vivado HLS在zedboard中的Sobel滤波算法实现

 基于Vivado HLS在zedboard中的Sobel滤波算法实现 平台:zedboard  + Webcam 工具:g++4.6  + VIVADO HLS  + XILINX EDK + ...

OPENCV基本滤波算法

图像滤波的主要目的是为了在保留图像细节的情况下尽量的对图像的噪声进行消除,从而是后来的图像处理变得更加的方便. 图像的滤波效果要满足两个条件:1.不能损坏图像的轮廓和边缘这些重要的特征信息.2.图像的 ...

基于FPGA的均值滤波算法实现

我们为了实现动态图像的滤波算法,用串口发送图像数据到FPGA开发板,经FPGA进行图像处理算法后,动态显示到VGA显示屏上,前面我们把硬件平台已经搭建完成了,后面我们将利用这个硬件基础平台上来实现基于 ...

基于FPGA的中值滤波算法实现

在这一篇开篇之前,我需要解决一个问题,上一篇我们实现了基于FPGA的均值滤波算法的实现,最后的显示效果图上发现有一些黑白色的斑点,我以为是椒盐噪声,然后在做基于FPGA的中值滤波算法的实验时,我发现黑 ...

linux内核中的C语言常规算法&lpar;前提&colon;你的编译器要支持typeof和type&rpar;

学过C语言的伙伴都知道,曾经比较两个数,输出最大或最小的一个,或者是比较三个数,输出最大或者最小的那个,又或是两个数交换,又或是绝对值等等,其实这些算法在linux内核中通通都有实现,以下的代码是我从 ...

C语言经典算法 - 多维矩阵转一维矩阵的代码

下边内容内容是关于C语言经典算法 - 多维矩阵转一维矩阵的内容,应该能对码农也有好处. #include #include int mai ...

一张有料的图片！！！附文件-图片合成器C语言实现算法

凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 1.打开F盘,先找两个文件,一个是压缩文件(2.rar),另一个是一张图片(1.jpg),将其放到F盘目录下, ...

随机推荐

Apriori——python3实现

最近看了关联算法中的Apriori没看懂,这次看了一些论文总算看懂了,不过还是没能够自己实现.在github搜到一些代码看,看的不很懂,这里先贴上(当中有自己加的注释),有时间再补充研究. # -*- ...

Hibernate jar包详解

Hibernate一共包括了23个jar包,令人眼花缭乱.本文将详细讲解Hibernate每个jar包的作用,便于你在应用中根据自己的需要进行取舍. 下载Hibernate,例如2.0.3稳定版本,解 ...

lib 合并

方法1: 命令如下:lib.exe lib1.lib lib2.lib ... /out:hw.lib 另一种合并的方法: 思路:lib.exe 可以把lib文件里的obj文件取出来,然后把这些obj ...

Atom编辑器入门到精通&lpar;四&rpar; Atom使用进阶

在本节中将介绍Atom提供的更高级的使用技巧,通过这些技巧将会进一步提高你的代码编写效率 代码片段(Snippets) Snippets是一种在代码中快捷插入代码块的方式,下面是维基百科中对Snipp ...

Matplotlib下载地址

http://www.lfd.uci.edu/~gohlke/pythonlibs/#matplotlib

Java中的变量与变量的作用域

关于Java中的变量及变量的作用域 关于Java中的变量及变量的作用域 0. 变量的概念 在程序运行期间,系统可以为程序分配一块内存单元,用来存储各种类型的数据.系统分配的内存单元要使用一个标记符来标 ...

Android四大组件之 --- Service入门

在前面了解了什么是多线程和异步消息处理机制后,我们来看一下到底什么是Service. 1. 什么是service(service的作用)一个服务是一个具有较长生命周期但没有用户界面的程序,例如:一个正 ...

C&num;简单线程

一.实例1 static void Main(string [] args) { Console.WriteLine("开始线程"); startFunc(); Console.W ...

SWUST OJ&lpar;952&rpar;

单链表的插入操作实现 #include #include typedef struct LinkList { int data; st ...

java线程安全— synchronized和volatile

java线程安全— synchronized和volatile package threadsafe; public class TranditionalThreadSynchronized { pu ...