技术领域
[0001] 本
发明涉及一种称重滤波方法，尤其涉及一种稳定可靠的动态称重滤波方法，并涉及采用了该动态称重滤波方法的动态称重滤波系统。
背景技术
[0002] 动态
电子秤是流
水线上用的称重工具，现在使用已经越来越广泛，其工作原理是将需要称重的物体不断地由传送带经过动态电子秤，进而使得动态电子秤能够不停地检测重量，
现有技术普遍使用的是300Hz的
采样频率，这种动态电子秤虽然能够满足最基本的流水称重工作，但是由于运动中的物体会振动，那么，必然会给称重数据造成很大的噪声进而使得这种动态的称重结果满足不了稳定可靠的要求，存在一定的误差。
[0003] 现有的动态电子秤等仪器虽然也存在采用模拟滤波或数字滤波的方法，模拟滤波采用贝塞尔(Bessel)
滤波器，通过设计模拟滤波
电路实现；模拟滤波电路的主要器件包括高
精度运算放大器和精密
电阻；通常采用多阶贝塞尔滤波器实现稳定的滤波效果，对应的滤波电路需要较多的高精度
运算放大器和精密电阻，不仅难以实现，而且成本高昂。
[0004] 随着
采样频率、采样精度的提高以及数字处理器性能的提高，数字滤波方式逐渐得到推广，现有技术中是将贝塞尔滤波器数字化转为
数字滤波器，可实现与贝塞尔模
块滤波器相近的效果。但是，这种数字化过程因离散采样和计算位数有限产生的舍入误差将降低滤波器的性能。此外，不同流水线上的物品重量及称重采样变化规律并不相同，因此，这种传统的、单一的滤波器是难以满足实际需求的。
发明内容
[0005] 本发明所要解决的技术问题是需要提供一种稳定可靠、速度快且精度高，能够尽量降低因为物体振动而带来的称重误差，可以平滑稳定过渡的动态称重滤波方法，并提供采用了该动态称重滤波方法的动态称重滤波系统。
[0006] 对此，本发明提供一种动态称重滤波方法，包括以下步骤：
[0007] 步骤S1，对传送带上的物体重量进行连续等周期的采样；
[0008] 步骤S2，对物体重量的采样结果进行过滤，得到滤波后的重量值；
[0009] 步骤S3，根据滤波后的重量值计算重量变化斜率；
[0010] 步骤S4，通过所述重量变化斜率判断所述物体的称量状态；
[0011] 步骤S5，根据称量状态确定是否采用
低通滤波器进行滤波，并输出最终的重量结果。
[0012] 本发明的进一步改进在于，所述步骤S1中，对所述物体重量进行连续等周期的采样过程中，采样周期时间为Ts，采样结果为ms,k，k表示第k个采样周期，k为自然数。
[0013] 本发明的进一步改进在于，所述步骤S2中，通过贝塞尔数字滤波器对所述采样结果进行过滤。
[0014] 本发明的进一步改进在于，所述步骤S2中，通过公式mb,k＝fB(mb,k-1,…,mb,k-M-1,ms,k,…,ms,k-M-1)对所述采样结果进行过滤，进而得到滤波后的重量值mb,k，其中，fB表示贝塞尔数字滤波器的数学函数，M表示贝塞尔数字滤波器的阶数，mb,k-1和mb,k-M-1分别表示第k-1个周期和第k-M-1个周期的贝塞尔数字滤波后的重量值；若k-x<1(x＝k-1,…,k-M-1)，则令mb,k-x＝0，ms,k-x＝0。
[0015] 本发明的进一步改进在于，所述步骤S3中，通过公式 计算重量变化斜率Sk，其中，p表示斜率间隔点。
[0016] 本发明的进一步改进在于，所述步骤S4中，设定重量变化斜率极大
阈值为Smax、重量变化斜率平衡阈值为Sbal和重量变化斜率最小阈值为Smin，然后将所述重量变化斜率Sk分别与所述重量变化斜率极大阈值Smax、重量变化斜率平衡阈值Sbal和重量变化斜率最小阈值Smin对比，进而判断所述物体的称量状态。
[0017] 本发明的进一步改进在于，所述步骤S4中，若所述重量变化斜率Sk大于所述重量变化斜率极大阈值Smax，则判定所述物体的称量状态为上升状态；若当前称量状态为上升状态且所述重量变化斜率Sk小于所述重量变化斜率平衡阈值Sbal，则判定所述物体的称量状态为稳定状态；若当前称量状态为稳定状态且所述重量变化斜率Sk小于所述重量变化斜率最小阈值Smin，则判定所述物体的称量状态为下降状态。
[0018] 本发明的进一步改进在于，所述步骤S5中，若称量状态为稳定状态，则采用所述低通滤波器进行滤波；若称量状态为初始状态、上升状态或下降状态中的任意一种，则不采用所述低通滤波器进行滤波，所述初始状态为最初的称量状态。
[0019] 本发明的进一步改进在于，所述步骤S5中，若采用所述低通滤波器进行滤波，则通过公式mf,k＝(1-α)·mf,k-1+α·mb,k计算低通滤波后的重量值mf,k，其中，α表示低通滤波因子，f表示低通滤波；若称量状态为稳定状态，则输出最终的重量结果mO,k＝mf,k；若称量状态为初始状态、上升状态或下降状态中的任意一种，则输出最终的重量结果mO,k＝mb,k。
[0020] 本发明还提供一种动态称重滤波系统，所述动态称重滤波系统采用了如上所述的动态称重滤波方法。
[0021] 与现有技术相比，本发明的有益效果在于：通过对传送带上的物体重量进行连续等周期的采样，然后对物体重量的采样结果进行过滤，并通过滤波后的重量值计算其重量变化斜率，进而能够根据重量变化斜率来判断所述物体的称量状态，以此决定是否采用低通滤波器进行滤波，因此，本发明符合物体在传送带上运动的特点，能够有效降低因为物体在传送带上运动而带来振动噪声和称重误差；本发明的滤波过程速度快且精度高，滤波平滑，过渡稳定，使得动态称重的过程更加稳定可靠，并使得动态称重的准确率和工作效率都得到了很大的提升。
附图说明
[0022] 图1是本发明一种
实施例的工作流程示意图；
[0023] 图2是本发明一种实施例中2阶贝塞尔数字滤波器的幅频特性仿真示意图；
[0024] 图3是本发明一种实施例中4阶贝塞尔数字滤波器的幅频特性仿真示意图；
[0025] 图4是本发明一种实施例中6阶贝塞尔数字滤波器的幅频特性仿真示意图；
[0026] 图5是本发明一种实施例中只使用贝塞尔数字滤波器以及使用了贝塞尔数字滤波器和低通滤波器的双重滤波之间的对比仿真效果示意图；
[0027] 图6是采用了阶数为4阶、截止频率为2Hz、采样频率为300Hz的贝塞尔数字滤波前和贝塞尔数字滤波后的对比仿真效果示意图；
[0028] 图7是采用了阶数为4阶、截止频率为4Hz、采样频率为300Hz的贝塞尔数字滤波前和贝塞尔数字滤波后的对比仿真效果示意图；
[0029] 图8是采用了阶数为4阶、截止频率为6Hz、采样频率为300Hz的贝塞尔数字滤波前和贝塞尔数字滤波后的对比仿真效果示意图；
[0030] 图9是采用了阶数为4阶、截止频率为8Hz、采样频率为300Hz的贝塞尔数字滤波前和贝塞尔数字滤波后的对比仿真效果示意图；
[0031] 图10是采用了阶数为2阶、截止频率为8Hz、采样频率为300Hz的贝塞尔数字滤波前和贝塞尔数字滤波后的对比仿真效果示意图；
[0032] 图11是采用了阶数为6阶、截止频率为8Hz、采样频率为300Hz的贝塞尔数字滤波前和贝塞尔数字滤波后的对比仿真效果示意图。
具体实施方式
[0033] 下面结合附图，对本发明的较优的实施例作进一步的详细说明。
[0034] 如图1所示，本例提供一种动态称重滤波方法，包括以下步骤：
[0035] 步骤S1，对传送带上的物体重量进行连续等周期的采样；
[0036] 步骤S2，对物体重量的采样结果进行过滤，得到滤波后的重量值；
[0037] 步骤S3，根据滤波后的重量值计算重量变化斜率；
[0038] 步骤S4，通过所述重量变化斜率判断所述物体的称量状态；
[0039] 步骤S5，根据称量状态确定是否采用低通滤波器进行滤波，并输出最终的重量结果。
[0040] 本例所述步骤S1中，对所述物体重量进行连续等周期的采样过程中，采样周期时间为Ts，采样结果为ms,k，k表示第k个采样周期，k为自然数，即k为大于或等于1的整数。
[0041] 本例所述步骤S2中，通过贝塞尔数字滤波器对所述采样结果进行过滤。所述步骤S2中，通过公式mb,k＝fB(mb,k-1,…,mb,k-M-1,ms,k,…,ms,k-M-1)对所述采样结果进行过滤，进而得到滤波后的重量值mb,k，其中，fB表示贝塞尔数字滤波器的数学函数，M表示贝塞尔数字滤波器的阶数，阶数M的优选取值范围为2～8的整数，mb,k-1和mb,k-M-1分别表示第k-1个周期和第k-M-1个周期的贝塞尔数字滤波后的重量值；若k-x<1(x＝k-1,…,k-M-1)，则令mb,k-x＝0，ms,k-x＝0。
[0042] 贝塞尔数字滤波器具有最大平坦的群延迟(线性
相位响应)的线性
过滤器。具有单位群时延的n阶贝塞尔滤波器的传递函数为 式中 为n阶贝塞尔数字滤波器的多项式。
[0043]
[0044]
[0045]
[0046]
[0047] 阶数越高或截止频率越低，贝塞尔数字滤波器的响应速度越慢，但滤波效果越好。
[0048] 利用数学辅助
软件Matlab可以求解不同阶数以及不同截止频率的贝塞尔模拟低通滤波器的传递函数，然后进行Z变换，转为
状态方程[0049] 获得状态方程系数ak(k＝1,2,…,n)和bk(k＝1,2,…,n-1)后，即可以在
单片机中设置贝塞尔数字滤波器，进行
迭代计算，完成贝塞尔数字滤波。
[0050] 本例所述贝塞尔数字滤波器的幅频特性如图2至图4所示，设定截止频率为8Hz(截止频率对幅频特性没有影响)，观察2、4、6阶贝塞尔数字滤波器的伯德(Bode)图，即通过仿真示意图比较不同阶数的贝塞尔数字滤波器的幅频特性。
[0051] 由图2至图4可以看出，阶数为2阶的贝塞尔数字滤波器的分频斜率约为9.8dB/oct，而阶数为4阶的贝塞尔数字滤波器其分频斜率增加至18.8dB/oct。根据原有设备的指标，4阶(或3阶)的贝塞尔数字滤波器能满足要求。
[0052] 本例所述步骤S3中，通过公式 计算重量变化斜率Sk，其中，p表示斜率间隔点，斜率间隔点p的优选取值范围为1～6的整数。
[0053] 本例所述步骤S4中，设定重量变化斜率极大阈值为Smax、重量变化斜率平衡阈值为Sbal和重量变化斜率最小阈值为Smin，然后将所述重量变化斜率Sk分别与所述重量变化斜率极大阈值Smax、重量变化斜率平衡阈值Sbal和重量变化斜率最小阈值Smin对比，进而判断所述物体的称量状态。
[0054] 本例所述步骤S4中，若所述重量变化斜率Sk大于所述重量变化斜率极大阈值Smax，则判定所述物体的称量状态为上升状态；若当前称量状态为上升状态且所述重量变化斜率Sk小于所述重量变化斜率平衡阈值Sbal，则判定所述物体的称量状态为稳定状态；若当前称量状态为稳定状态且所述重量变化斜率Sk小于所述重量变化斜率最小阈值Smin，则判定所述物体的称量状态为下降状态。
[0055] 本例所述步骤S5中，若称量状态为稳定状态，则采用所述低通滤波器进行滤波；若称量状态为初始状态、上升状态或下降状态中的任意一种，则不采用所述低通滤波器进行滤波，所述初始状态为最初的称量状态。
[0056] 本例所述步骤S5中，若采用所述低通滤波器进行滤波，通过公式mf,k＝(1-α)·mf,k-1+α·mb,k计算低通滤波后的重量值mf,k，其中，α表示低通滤波因子，低通滤波因子α的优选取值范围为0～0.5的实数，f表示低通滤波；若称量状态为稳定状态，则输出最终的重量结果mO,k＝mf,k；若称量状态为初始状态、上升状态或下降状态中的任意一种，则输出最终的重量结果mO,k＝mb,k。
[0057] 本例还提供一种动态称重滤波系统，所述动态称重滤波系统采用了如上所述的动态称重滤波方法。
[0058] 值得一提的是，本例所述低通滤波器为改进型的低通滤波器。
[0059] 本例引入改进型的低通滤波器(MLP)的原因是仅仅只有贝塞尔数字滤波器(Bessel滤波器)在平台段的滤波器是不够理想，如图6至图11所示，即使将截止频率设置为2Hz，在平台段仍存在2.5Hz的频率振荡。采用Butterworth
带阻滤波器或贝塞尔数字滤波器二次滤波均无法解决低频率振荡的问题。至于振荡的原因，很可能是物体在等间距传送带的转动轴(轮)上运动造成的。由图6至图11可知，截止频率越高，或贝塞尔数字滤波器的阶数越低，响应特性越好，但是平滑程度越差。
[0060] 因此，本例在步骤S4和步骤S5中引入改进型的低通滤波器，其设计出发点为判断重量检测数据的斜率变化以确定是否需要进行低通滤波。因为物体在传送带上运动时，重量检测数值先逐渐增大，然后保持稳定，最后逐渐减小。所以，在重量检测数值逐渐增大过程中，重量变化的斜率先增大后减小至0甚至为负值，当斜率减小至某个值时，即若称量状态为稳定状态时，则启动所述低通滤波器；在重量检测数值逐渐减小过程中，重量变化的斜率由接近0开始减小，当斜率减小至某个值时，即当称量状态为下降状态时，停止低通滤波器；在初始状态和上升状态，同样也不需要启动低通滤波器。
[0061] 本例引入的改进型的所述低通滤波器的滤波效果如图5所示；图5中，双重滤波指的是经过了步骤S2的贝塞尔数字滤波器之后，还经过了步骤S4和步骤S5的低通滤波器，即本例所述的动态称重滤波方法其实就是双重滤波。在实验过程中，对采样频率为300Hz的600个点进行处理，结果如图5所示，图5的双重滤波指的是先用贝塞尔数字滤波器(n＝4，fs＝2)滤波，然后在步骤S4和S5中再根据实际情况用改进型的低通滤波器进行滤波。由图5可以看出，经过本例所述动态称重滤波方法的双重滤波后，平台段数值变化范围小于5；而仅用贝塞尔数字滤波器(n＝4，fs＝2)滤波的结果数值变化范围将近25，如图6至图11所示。
[0062] 综上所述，本例通过对传送带上的物体重量进行连续等周期的采样，然后对物体重量的采样结果进行过滤，并通过滤波后的重量值计算其重量变化斜率，进而能够根据重量变化斜率来判断所述物体的称量状态，以此决定是否采用低通滤波器进行滤波，因此，本例符合物体在传送带上运动的特点，能够有效降低因为物体在传送带上运动而带来振动噪声和称重误差；本例的滤波过程速度快且精度高，滤波平滑，过渡稳定，使得动态称重的过程更加稳定可靠，并使得动态称重的准确率和工作效率都得到了很大的提升。
[0063] 以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。
