下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上1、限幅滤波法(又称程序判断滤波法)A、方法:根据经验判断,确定两次采样允许的最大偏差值(设为A)每次检测到新值时判断:如果本次值与上次值之差<=A,则本次值有效如果本次值与上次值之差>A,则本次值无效,放弃本次值,用上次值代替本次值B、优点:能有效克服因偶然因素引起的脉冲干扰C、缺点无法抑制那种周期性的干扰平滑度差#define A 10 char value;char filter() char new_value; new_value = get_ad(); if ( ( new_value - value > A ) | ( value -
2、new_value > A ) return value; return new_value; 2、中位值滤波法A、方法:连续采样N次(N取奇数)把N次采样值按大小排列取中间值为本次有效值B、优点:能有效克服因偶然因素引起的波动干扰对温度、液位的变化缓慢的被测参数有良好的滤波效果C、缺点:对流量、速度等快速变化的参数不宜#define N 11 char filter() char value_bufN; char count,i,j,temp; for ( count=0;count<N;count+) value_bufcount = get_ad(); delay(); fo
3、r (j=0;j<N-1;j+) for (i=0;i<N-j;i+) if ( value_buf>value_bufi+1 ) temp = value_buf; value_buf = value_bufi+1; value_bufi+1 = temp; return value_buf(N-1)/2; 3、算术平均滤波法A、方法:连续取N个采样值进行算术平均运算N值较大时:信号平滑度较高,但灵敏度较低N值较小时:信号平滑度较低,但灵敏度较高N值的选取:一般流量,N=12;压力:N=4B、优点:适用于对一般具有随机干扰的信号进行滤波这样信号的特点是有一个平均值,信号在某
4、一数值范围附近上下波动C、缺点:对于测量速度较慢或要求数据计算速度较快的实时控制不适用比较浪费RAM#define N 12 char 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;压
5、力:N=4;液面,N=412;温度,N=14B、优点:对周期性干扰有良好的抑制作用,平滑度高适用于高频振荡的系统C、缺点:灵敏度低对偶然出现的脉冲性干扰的抑制作用较差不易消除由于脉冲干扰所引起的采样值偏差不适用于脉冲干扰比较严重的场合比较浪费RAM#define N 12 char value_bufN;char i=0; char filter() char count; int sum=0; value_bufi+ = get_ad(); if ( i = N ) i = 0; for ( count=0;count<N,count+) sum = value_bufcount; r
6、eturn (char)(sum/N); 5、中位值平均滤波法(又称防脉冲干扰平均滤波法)A、方法:相当于“中位值滤波法”+“算术平均滤波法”连续采样N个数据,去掉一个最大值和一个最小值然后计算N-2个数据的算术平均值N值的选取:314B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C、缺点:测量速度较慢,和算术平均滤波法一样比较浪费RAM#define N 12 char filter() char count,i,j; char value_bufN; int sum=0; for (count=0;count<N;count+) val
7、ue_bufcount = get_ad(); delay(); for (j=0;j<N-1;j+) for (i=0;i<N-j;i+) if ( value_buf>value_bufi+1 ) temp = value_buf; value_buf = value_bufi+1; value_bufi+1 = temp; for(count=1;count<N-1;count+) sum += valuecount; return (char)(sum/(N-2); 6、限幅平均滤波法A、方法:相当于“限幅滤波法”+“递推平均滤波法”每次采样到的新数据先进行限幅
8、处理,再送入队列进行递推平均滤波处理B、优点:融合了两种滤波法的优点对于偶然出现的脉冲性干扰,可消除由于脉冲干扰所引起的采样值偏差C、缺点:比较浪费RAM略 参考子程序1、37、一阶滞后滤波法(低通滤波)A、方法:取a=01本次滤波结果=(1-a)*本次采样值+a*上次滤波结果B、优点:对周期性干扰具有良好的抑制作用适用于波动频率较高的场合C、缺点:相位滞后,灵敏度低滞后程度取决于a值大小不能消除滤波频率高于采样频率的1/2的干扰信号#define a 50 char value; char filter() char new_value; new_value = get_ad(); retu
9、rn (100-a)*value + a*new_value; 8、加权递推平均滤波法A、方法:是对递推平均滤波法的改进,即不同时刻的数据加以不同的权通常是,越接近现时刻的数据,权取得越大。给予新采样值的权系数越大,则灵敏度越高,但信号平滑度越低B、优点:适用于有较大纯滞后时间常数的对象和采样周期较短的系统C、缺点:对于纯滞后时间常数较小,采样周期较长,变化缓慢的信号不能迅速反应系统当前所受干扰的严重程度,滤波效果差#define N 12 char code coeN = 1,2,3,4,5,6,7,8,9,10,11,12;char code sum_coe = 1+2+3+4+5+6+7
10、+8+9+10+11+12; char filter() char count; char value_bufN; int sum=0; for (count=0,count<N;count+) value_bufcount = get_ad(); delay(); for (count=0,count<N;count+) sum += value_bufcount*coecount; return (char)(sum/sum_coe); 9、消抖滤波法A、方法:设置一个滤波计数器将每次采样值与当前有效值比较:如果采样值当前有效值,则计数器清零如果采样值<>当前有效值
11、,则计数器+1,并判断计数器是否>=上限N(溢出)如果计数器溢出,则将本次值替换当前有效值,并清计数器B、优点:对于变化缓慢的被测参数有较好的滤波效果,可避免在临界值附近控制器的反复开/关跳动或显示器上数值抖动C、缺点:对于快速变化的参数不宜如果在计数器溢出的那一次采样到的值恰好是干扰值,则会将干扰值当作有效值导入系统#define N 12 char filter() char count=0; char new_value; new_value = get_ad(); while (value !=new_value); count+; if (count>=N) return
12、 new_value; delay(); new_value = get_ad(); return value; 10、限幅消抖滤波法A、方法:相当于“限幅滤波法”+“消抖滤波法”先限幅,后消抖B、优点:继承了“限幅”和“消抖”的优点改进了“消抖滤波法”中的某些缺陷,避免将干扰值导入系统C、缺点:对于快速变化的参数不宜第11种方法:IIR 数字滤波器A. 方法:确定信号带宽, 滤之。Y(n) = a1*Y(n-1) + a2*Y(n-2) + . + ak*Y(n-k) + b0*X(n) + b1*X(n-1) + b2*X(n-2) + . + bk*X(n-k)B. 优点:高通,低通,带
13、通,带阻任意。设计简单(用matlab)C. 缺点:运算量大。int BandpassFilter4(int InputAD4) int ReturnValue; int ii; RESLO=0; RESHI=0; MACS=*PdelIn; OP2=1068; /FilterCoeff44; MACS=*(PdelIn+1); OP2=8; /FilterCoeff43; MACS=*(PdelIn+2); OP2=-2001;/FilterCoeff42; MACS=*(PdelIn+3); OP2=8; /FilterCoeff41; MACS=InputAD4; OP2=1068; /
14、FilterCoeff40; MACS=*PdelOu; OP2=-7190;/FilterCoeff48; MACS=*(PdelOu+1); OP2=-1973; /FilterCoeff47; MACS=*(PdelOu+2); OP2=-19578;/FilterCoeff46; MACS=*(PdelOu+3); OP2=-3047; /FilterCoeff45; *p=RESLO; *(p+1)=RESHI; mytestmul<<=2; ReturnValue=*(p+1); for (ii=0;ii<3;ii+) DelayInputii=DelayInpu
15、tii+1; DelayOutputii=DelayOutputii+1; DelayInput3=InputAD4; DelayOutput3=ReturnValue; / if (ReturnValue<0) / / ReturnValue=-ReturnValue; / return ReturnValue; 第12种方法:RC滤波的一种实现.RcDigital(double & X, double & Y) static int MidFlag; static double Yn_1,Xn_1; double MyGetX=0,MyGetY=0; double A
16、lfa; Alfa=0.7; if(X=0|Y=0) MidFlag=0; Xn_1=0; Yn_1=0; MyGetX=0; MyGetY=0; if(X>0&&Y>0) if(MidFlag=1) MyGetY = (1 - Alfa) * Y + Alfa * Yn_1; MyGetX = (1 - Alfa) * X + Alfa * Xn_1; Xn_1 = MyGetX; Yn_1 = MyGetY; else MidFlag=1; MyGetX = X; MyGetY = Y; Xn_1 = X; Yn_1 = Y; X = MyGetX; Y = M
17、yGetY; 实际上的应用并没有单一的应用,例如:实际的AD上存在一定的冲击干扰,这时就中位值平均滤波,又由于效率的问题对其添加滑动平均滤波法,提高其动态性。sint32 slAD =0; void filter() static sint32 slaTemp10 = 0; static uint8 k =0 ; static sint32 slFirst = 0; static uint8 ucFirstFlag =0; uint8 queue10=0,1,2,3,4,5,6,7,8,9; uint8 ucTemp =0; uint8 i,j =0; sint32 slTemp =0; slaTempk = get_ad()/2; k+; for(j=0;j<10;j+) for(i=j+1;i&
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 电信多路复用器市场发展预测和趋势分析
- 纱笼市场发展现状调查及供需格局分析预测报告
- 竹篮产业规划专项研究报告
- 《品牌价值超越竞争》课件
- 沙狐球饼产业深度调研及未来发展现状趋势
- 《循环系统影像诊断》课件
- 电磁铁市场洞察报告
- 笔筒产品入市调查研究报告
- 白垩粉产业运行及前景预测报告
- 照相用回光灯市场发展现状调查及供需格局分析预测报告
- 水电天然气的费用——小数应用
- 劳动合同书英文版范文
- 估计不规则图形的面积(新人教版)教学设计
- 污水排放记录(共3页)
- 第五章成本论解析
- 舞台照明灯项目申请报告(模板范文)
- 如何激发小学生学习英语的兴趣(课堂PPT)
- DSM-5量表(ADHD)
- BAS建筑设备监控系统施工调试方法及程序
- 物业合资公司组建方案
- 《埃博霉素生物合成》PPT课件.ppt
评论
0/150
提交评论