滤波程序全局库

软件版本：TIA15.1

功能：可实现对采集数据的滤波，根据现场不同情况，设置滤波参数，达到理想曲线。

移植性：功能块编写语言均为SCL,方便移植到其它品牌plc。

内容：此库内包含三中滤波程序，一阶滞后滤波法、中位值滤波法、堆栈平均滤波法。

滤波效果如下图：

①一阶滞后滤波法封装： ↓

代码：

IF #"a (0-1)" >= 1.0 THEN

    #"a (0-1)" := 0.5;

END_IF;

#ms_time := INT_TO_TIME(#time_ms);

#IEC_Timer_0_Instance(IN := NOT #p,

                      PT := #ms_time,

                      Q => #p);

IF #p THEN

    #new_value := (1 - #"a (0-1)") * #IN + #"a (0-1)" * #last_value;

    #last_value := #new_value;

    #OUT := #new_value;

END_IF;

②堆栈平均滤波法封装  ↓

代码：

#IEC_Timer_0_Instance(IN:=NOT #TIMER_INOUT,

                      PT:=#time_采集时间,

                      Q=>#TIMER_INOUT);

#S_数组[#"N_采集数量(3-999)"] := #In_输入值;

IF #TIMER_INOUT THEN

    FOR #S_j := 1 TO  #"N_采集数量(3-999)" DO

        #S_数组[#S_j-1] := #S_数组[#S_j];

    END_FOR;

END_IF;

#T_max := #S_数组[0];

#T_min := #S_数组[0];

FOR #S_j := 0 TO #"N_采集数量(3-999)"-1 DO

    IF #S_数组[#S_j] > #T_max THEN

        #T_max := #S_数组[#S_j];

    END_IF;

    IF #S_数组[#S_j] < #T_min THEN

        #T_min := #S_数组[#S_j];

    END_IF;

    #T_sum := #T_sum + #S_数组[#S_j];

END_FOR;

IF #"MAX/MIN_EN" THEN

#Out_输出结果:=(#T_sum-#T_max- #T_min )/ (#"N_采集数量(3-999)"-2);

ELSE

    #Out_输出结果 := #T_sum / #"N_采集数量(3-999)";

END_IF;

#Out_max := #T_max;

#Out_min := #T_min;

③中位值滤波法封装  ↓

代码：

//n 须为3<n<999的奇数

IF #N < 3 THEN

    #N := 3;

END_IF;

IF #N > 999 THEN

    #N := 999;

END_IF;

#MOD :=#N MOD 2;

IF #MOD = 0 THEN

    #N := #N + 1;

END_IF;

//转换成时间格式

#time_temp := INT_TO_TIME(#time_ms);

//生成脉冲，对n个数堆栈

#IEC_Timer_0_Instance(IN:=NOT #p,

                      PT:=#time_temp,

                      Q=>#p );

IF #p THEN

    FOR #x := 0 TO #N - 2 DO

        #value[#x] := #value[#x + 1];

    END_FOR;

   #value[#N-1] := #IN;

END_IF;

 // 将value数组中的数复制到temp_value数组中

MOVE_BLK(IN:=#value[0],

             COUNT:=1000,

             OUT=>#temp_value[0]);

    //对temp_value数组中数从小到大排序

    FOR #j := 0 TO #N - 3 DO

        FOR #i := #j TO #N - 2 DO

            IF #temp_value[#j] > #temp_value[#i + 1] THEN

                #value_temp := #temp_value[#i + 1];

                #temp_value[#i + 1] := #temp_value[#j];

                #temp_value[#j] := #value_temp;

            END_IF;

        END_FOR;

    END_FOR;

    //在排好序的数组中,取中间数做为输出

     #OUT := #temp_value[(#N - 1) / 2];

三种滤波方法全局库.rar

编辑不易给个赞呗

点赞是必须的，辛苦了。

我一直以来都喜欢一阶滤波，简单方便。

 论坛有你更精彩！

谢谢分享，我先去学习下基础知识

能分享知识就值得点赞！！

谢谢分享，点赞！

下载不了啊

能下载呀，点赞！

点赞，这个是真的干货

感谢分享~~

太专业了，有点看不懂了，，，

，，赞一个

感谢分享！

感谢分享！

感觉使用1500的基本都是高手。

有没有一种滤波方法可以把输入信号中断近似n赫兹成分的波动去除？

重点，只对n赫兹附近的波动做去除。去除深度可调，附近半径可调。其他频率不去除。

感谢分享，必须赞

为大神分享干货点赞!

非常不错图文并茂，三种滤波方法对比效果明显。

,做成全局库确实不错，可用性与可移植性都不错

谢谢大神分享，必须点赞。

绝对是干货

一阶滤波公式写错了吧，应该是 #new_value := (1 - #"a (0-1)") * #last_value + #"a (0-1)" * #IN;

你说的这种叫陷波器，以前模拟电路中经常见到的，原理上就是低通滤波+高通滤波。用电阻、电容、电感拼凑就能够做到。

其实，现在用电阻、电容、运放组成0衰减陷波器的更好。

说到电阻电容，想到的肯定是积分、微分了，所以某一个频率陷波器就是根据频率周期判断2个周期采集的2个值之间相减的绝对值必须小于(低通)或者大于(高通)某个范围值。

嗯，如果学过傅立叶的应该都知道这个，FFT运算就是太费时了，不过有个查表法可以提高效率...

好嘞，以上纯粹我瞎扯淡的。

必须给赞，感谢楼主的分享。

有时滤波对于信号处理是必需的算法，谢谢分享！

楼主，谢谢你的分享。真心感谢！

感谢分享！！

。下载下来仔细看一下。

牛 逼格拉斯！！！

简单的信号不需要滤波;复杂的信号又过滤不了;滤波权重的算法比滤波算法本身还复杂.

这个用编程计算的话， 首先，缓存一段数据， 然后，用福利也变化求出对应频率的分量。

用缓存减去分量。

有，西门子1500自带，带通滤波

西门子库自带，带通滤波

hahahahhahahahahaaaa

谢谢分享，给楼主点赞。

谢谢分享！

我感觉也写错了,你不说我都不知道,,百度了一下:一阶滤波算法是比较常用的滤波算法，它的滤波结果=a*本次采样值+（1-a）*上次滤波结果，其中，a为0~1之间的数。一阶滤波相当于是将新的采样值与上次的滤波结果计算一个加权平均值。a的取值决定了算法的灵敏度，a越大，新采集的值占的权重越大，算法越灵敏，但平顺性差；相反，a越小，新采集的值占的权重越小，灵敏度差，但平顺性好。

感谢分享。

真是好东西，辛苦了

感谢楼主分享

还是喜欢看这样的贴子，而不是一些假大空的开天辟地的

有图，滤波后效果一目了然

仔细看了一下堆栈平均滤波法的代码，发现了一个小Bug，不知道对不对。

n个数据没采集完时，也就是数组里数据量小与n时，算平均值时还是除以n，会导致输出结果不对。

等到采集的数据量大于等于n时就没问题了。

V16怎么打开这个程序啊之前我记得可以打开的找不到了

，学习下，感谢

很好的分享。

一阶滞后滤波  公式错了  楼主改下吧。

他们说公式有错？