CN109145398B - 一种基于称重的货品数量计算方法 - Google Patents
一种基于称重的货品数量计算方法 Download PDFInfo
- Publication number
- CN109145398B CN109145398B CN201810853293.6A CN201810853293A CN109145398B CN 109145398 B CN109145398 B CN 109145398B CN 201810853293 A CN201810853293 A CN 201810853293A CN 109145398 B CN109145398 B CN 109145398B
- Authority
- CN
- China
- Prior art keywords
- weight
- goods
- calculating
- weighing
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005303 weighing Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000006243 chemical reaction Methods 0.000 claims abstract description 32
- RCKMWOKWVGPNJF-UHFFFAOYSA-N diethylcarbamazine Chemical compound CCN(CC)C(=O)N1CCN(C)CC1 RCKMWOKWVGPNJF-UHFFFAOYSA-N 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 241000282326 Felis catus Species 0.000 claims description 7
- 230000002019 anti-mutation Effects 0.000 abstract description 15
- 238000012937 correction Methods 0.000 abstract description 15
- 238000004364 calculation method Methods 0.000 abstract description 11
- 238000001914 filtration Methods 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000035772 mutation Effects 0.000 abstract description 6
- 238000007726 management method Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 description 20
- 238000005070 sampling Methods 0.000 description 18
- 230000008859 change Effects 0.000 description 17
- 238000005259 measurement Methods 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 235000008446 instant noodles Nutrition 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002790 anti-mutagenic effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- QVFWZNCVPCJQOP-UHFFFAOYSA-N chloralodol Chemical compound CC(O)(C)CC(C)OC(O)C(Cl)(Cl)Cl QVFWZNCVPCJQOP-UHFFFAOYSA-N 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010926 purge Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000013585 weight reducing agent Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G23/00—Auxiliary devices for weighing apparatus
- G01G23/01—Testing or calibrating of weighing apparatus
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G23/00—Auxiliary devices for weighing apparatus
- G01G23/14—Devices for determining tare weight or for cancelling out the tare by zeroising, e.g. mechanically operated
Abstract
本发明公开一种基于称重的货品数量计算方法,包括以下步骤:步骤一、重量转换;将称重读数准确地转换为对应的重量值W;步骤二、重量缩放;即对称重单位的转换;步骤三、零点纠偏;对称重读数进行校准;步骤四、去皮;去除重量值中非货品物体的重量;步骤五、抗突变滤波;对重量值进行抗突变滤波;步骤六、数量转换;将经过了上述步骤一至五校准和纠偏后的重量值转换为货品大体数量;步骤七、数量校准;对上述货品大体数量进行校准。本发明计算精准,著提升了基于称重的货位管理技术在实际工作场景中的精度和稳定性。
Description
技术领域
本发明涉及一种货品数量计算方法,尤其涉及一种基于称重的货品数量计算方法。
背景技术
随着物联网(IoT)技术和智能系统的逐步普及,为每个货位加装称重传感器,并通过重量来实时统计该货位中的货品数量的方法,正成为近年来较流行的发展方向。此类方案通常由称重传感器、数模转换(ADC)模块,以及工控机等部分组成。其中称重传感器通常安装在货位托盘下方,将重量转换为相应的电压或电流值;数模转换模块则实时测量称重传感器的电压变化,并将其转换为数字信号;工控机从数模转换模块读取最新的称重裸读数(ADC裸读数),将其换算为重量单位,并进一步完成数量统计。
上述方案具备成本低、可靠性高、环境耐受性好、维护简单等优点。但与此同时,该方案也存在如下不足:
1.重量转换的精度问题:由ADC裸读数转换为重量时,通常仅采用两点线性标定算法,转换精度不足。
2.重量缩放问题:重量转换后的单位可能不便于人员实际操作(如:重量转换标定时使用“克”为单位,但实际操作是可能使用“吨”更为方便)。
3.零点偏移问题:称重传感器可能由于长期负重而产生永久性或半永久性机械形变(蠕变),或由于结露、温度变化、落尘、沙土等环境因素导致其零点测量发生偏移,致使测量读数失准。现有的手动校准技术需要人工操作和测量,浪费大量人工成本,并使得业务中断。同时还有校准间隔大,难以保持称重读数长期稳定和准确,以及人工操作容易失误引起更大误差等问题。
4.去皮问题:系统需要能够感应并自动去除货位上用于盛放货品的蓝框等容器的重量。
5.突变问题:现代称重传感器由于制作材质和工艺的原因,其满量程电压通常也只有几个毫伏(mV)的范围,非常微弱。因此极易受到外部环境的电磁干扰而产生较高的抖动和误差;另一方面,货位本身也可能被移动,而移动过程中产生的惯性和离心力等因素也会导致称重读数瞬间发生变化。这些突变因素均可导致称重读数以及货品数量计量的不准确。
例如:在一个典型的现代仓库(一个仓库内可同时容纳多个货架,每个货架内又可包含多个货位)场景中,充斥着WIFI、手机信号(3G/4G/5G)、Lora、Zigbee等射频干扰,以及各类机电设备和交流输电线路产生的工频干扰等电磁干扰。同时,现代仓储中普遍使用的AGV车在托举或放下货架时,会产生惯性(电梯效应),并且其在背负货架运动时,也会由于颠簸,转向等原因产生惯性和离心力上的变化。这些都会对相应货位上称重传感器的读数产生影响。
6.公差问题:由于一个货位中可同时容纳多件货品,因此其中每件货品的公差将可能产生累积,导致最终的货品数量计算错误。
例如:考虑一种货品A的重量为100g,其公差范围为正负2%。即:出厂合格品的单件重量范围在98g到102g之间。假如一个货位同时堆放了200件此类A货品,则计入公差后,其合法的重量范围在19600g到20400g之间。意即:即使我们克服了之前的5点问题,获得了精确且稳定的称重读数,以此最终计算出的商品数量仍然可能有高达8件(正负4件)的偏差。
显然,上述问题严重阻碍了基于称重的货位管理技术在实际应用场景中的效果。
发明内容
本发明的目的是提供一种基于称重的货品数量计算方法,显著提升了基于称重的货位管理技术在实际工作场景中的精度和稳定性。
为了实现上述目的,本发明的技术方案是:
一种基于称重的货品数量计算方法,包括以下步骤:
步骤一、重量转换;将称重读数准确地转换为对应的重量值W;
步骤二、重量缩放;即对称重单位的转换;
步骤三、零点纠偏;对称重读数进行校准;
步骤四、去皮;去除重量值中非货品物体的重量;
步骤五、抗突变滤波;对重量值进行抗突变滤波;
步骤六、数量转换;将经过了上述步骤一至五校准和纠偏后的重量值转换为货品大体数量;
步骤七、数量校准;对上述货品大体数量进行校准。
进一步的,所述步骤一中的转化公式为:W=f({S},ADC),其中W为重量值,f表示要使用的线性或非线性插值算法,{S}表示事先经标定的采样集合,ADC表示的称重数据。
进一步的,所述步骤二的计算公式为:SW=WS*W(WU);其中SW为缩放后重量,WS为事先指定的可变系数,WU为单位。
进一步的,所述步骤三中零点纠偏为零点自动校准,当系统检测到当前货位已被清空,但称重读数却并未归零时,在至少有预设数值的连续采样读数落入称重读数阈值ZPT的范围内后,将连续采样数ZPCS个的采样按照加权的形式求取平均数,并将该数值作为新的零点偏移补偿量保存到零点偏移补偿量ZPO中。
进一步的,所述步骤三中零点纠偏为零点自动追随校准,在货位处于稳定状态时,且至少有预设数值的连续采样读数的变化量落入了[ZPTMIN,ZPTMAX]范围内后,将连续采样数(ZPCS)个的采样数学平均数或加权平均值累加到零点偏移补偿量(ZPO)中,其中ZPTMAX为最大跟随步进值,ZPTMIN为最小跟随步进值。
进一步的,所述步骤三中零点纠偏为零点自动追随校准,在货位处于不稳定状态时,自动记录并调整货品出货、进货等业务引起的正常重量变化,(出货导致对应货位正常减重,进货导致对应货位正常增重)。在此基础上,若检测到至少预设数值的连续采样读数,其正常变化之外的差量落入了[ZPTMIN,ZPTMAX]范围内后,将连续采样数ZPCS个的采样数学平均数或加权平均值累加到零点偏移补偿量ZPO中,其中ZPTMAX为最大跟随步进值,ZPTMIN为最小跟随步进值。
进一步的,所述步骤四旨在去除重量值中篮筐容器等非货品物体的重量。去皮操作通过在当前重量值上,减去一个事先指定的皮重参数BW来实现。减去BW参数后的值即为当前货位上的纯货品重量值。
进一步的,所述步骤五从采样中排除由电磁干扰、惯性、离心力等现象带来的称重读数突变问题。抗突变滤波过程主要由:抗突变最大范围AMMAX、抗突变最大范围最小范围AMMIN、抗突变比例范围AMR以及抗突变采样数AMS等参数控制。
进一步的,所述步骤六中的数量转换公式为:货品大体数量P=经过步骤五后的重量值/货品的标准毛重GSW
进一步的,所述步骤七中,通过对满仓容量CAP、货品标准毛重GSW、公差范围APU、最大累计误差CAMAX、最小累计误差CAMIN以及误差容限ATF的限定,计算累计误差容忍度CAT的范围;进而使用累计误差容忍度CAT对计算出的货品大体数量进一步校准。
本发明相对于现有技术具有以下优点:
本发明显著提升了基于称重的货位管理技术在实际工作场景中的精度和稳定性。
附图说明
图1是本发明流程示意图。
图2是本发明重量转换的线性插值算法示例。
图3是本发明重量转换的非线性插值算法示例。
图4是本发明数量转换中标准毛重GSW在重量分布中的位置示意图。
图5是本发明数量校准中优选实施例“倒S型”(inverse sigmoid)平滑调整曲线的示意图。
具体实施方式
以下结合附图进一步说明本发明的实施例。
实施例:请参见图1所示,
一种基于称重的货品数量计算方法,包括以下步骤:
步骤一、重量转换;将称重读数准确地转换为对应的重量值W;
步骤二、重量缩放;即对称重单位的转换;
步骤三、零点纠偏;对称重读数进行校准;
步骤四、去皮;去除重量值中非货品物体的重量;
步骤五、抗突变滤波;对重量值进行抗突变滤波;
步骤六、数量转换;将经过了上述步骤一至五校准和纠偏后的重量值转换为货品大体数量;
步骤七、数量校准;对上述货品大体数量进行校准;
其中,所述步骤六中的数量转换公式为:货品大体数量P=经过步骤一至五后的重量值/货品的标准毛重GSW。
优选的,由于传统的两点标定线性转换无法真实反映称重传感器的物理特性,导致重量转换结果准确度不足。为此,步骤一中的转化公式为:W=f({S},ADC),其中W为重量值,f表示要使用的线性或非线性插值算法,{S}表示事先经标定的采样集合,ADC表示的称重数据(裸读数)。请参见图2和图3所示,该转换方法支持任意多个称重数据采样S(S1、S2、S3等),其中每个采样均包含ADC裸读数X及与之对应的实际重量值Y(S1=[X1,Y1]、S2=[X2,Y2]、S3=[X3,Y3]……),并结合线性插值(分段函数)或非线性插值(平滑曲线)算法,可将ADC裸读数准确地转换为对应重量值。
优选的,步骤二中重量缩放本质上是对称重单位的转换,可通过将当前称重值乘以一个事先指定的可变系数WS来实现的计算公式为:SW=WS*W(WU);其中SW为缩放后重量,WS为事先指定的可变系数,WU为单位(如:吨、公斤、市斤、克等等)。
优选的,步骤三中零点纠偏为零点自动校准,零点自动校准指系统检测到当前货位已被清空,但经过重量转换和重量缩放等步骤后的称重读数却并未归零时,自动重新校准称重读数的零点补偿偏移值,以帮助称重系统回零的过程。
具体的零点自动校准主要由两个参数进行控制:称重读数阈值ZPT和连续采样数ZPCS。其中称重读数阈值ZPT用来保证当货位已清空时的称重读数最大范围(绝对值,ZPT也可被分开指定为ZPT+和ZPT-:即正值最大范围和负值最大范围);而连续采样数ZPCS则规定了至少有多少个连续采样读数完全落入ZPT参数所规定的范围内后,才执行零点自动校准过程(可以为1,表示不需要执行连续采样)。
当系统检测到至少有预设数值(如连续采样数ZPCS)的连续采样读数落入称重读数阈值ZPT的范围内(说明当前货位已被清空),但称重读数却并未归零时,将连续采样数ZPCS个的采样按照加权的形式求取平均数,并将该数值作为新的零点偏移补偿量保存到零点偏移补偿量ZPO中,并在未来的重量计算过程中,应用该偏移值,以对传感器的零点漂移(蠕变)进行补偿。
例如:ZPT=30、ZPCS=20表示至少有连续20个采样的重量值落入±30的范围内时,才开始一次零点自动校准过程。具体的自动校准过程可以为:将最近20次的采样按照加权(例如:越新的采样权重越高;或简单的数学平均数,即:所有采样的权重均为1;或以第一个采样为准,即:第一个采样的权重为1,后续采样的权重均为0,等等)的形式求取平均数,并将该数值作为新的零点偏移补偿量保存到ZPO中。未来的每次新采样在完成重量转换和缩放等步骤后,均会追加ZPO值以补偿零点漂移。
显然,ZPT值应该远小于当前货位上篮筐等容器或单件货品(无容器情况下)的重量,否则可能导致误判(将篮筐或货品的重量错误地判定为零点漂移)。推荐的ZPT值应小于篮筐或单件货品重量的50%。较大的ZPCS值则可保护系统不会受到由电磁干扰、惯性、离心力等情况引起的突变干扰。
优选的,步骤三中零点纠偏为零点自动追随校准,与零点自动校准不同,零点自动追随功能无需清空货位即可被自动激活。零点自动追随功能可以在货位处于稳定状态(例如:货品数量未发生变化)以及非稳定状态(动态,例如:持续进行进货、出货等操作,货品数量频繁发生变化等情况)的前提下,实时、动态地对当前货位上,由蠕变、环境变化等原因引起的零点漂移现象进行实时校准。
零点自动跟随过程由最大跟随步进值ZPTMAX、最小跟随步进值ZPTMIN、以及连续采样数ZPTCS等参数控制。
其中最大跟随步进值ZPTMAX决定了每次发生零点自动追随过程时,可接受的最大重量变化值(Delta,绝对值,变化值是指当前重量值与上次进入稳态或上次完成零点自动追随后的重量值之差,ZPTMAX也可被分开指定为ZPTMAX+和ZPTMAX-:即正值最大范围和负值最大范围)。
而最小跟随步进值ZPTMIN则决定了要激活一次零点自动追随过程需要的最低重量变化值范围(绝对值,“变化值”的定义见上文;可以为0,表示不限制最小波动范围,ZPTMIN也可被分开指定为ZPTMIN+和ZPTMIN-:即正值最小范围和负值最小范围)。
最后,连续采样数ZPTCS规定了至少有多少个连续采样读数的变化(差量)值(Delta)的绝对值完全落入了[ZPTMIN..ZPTMAX](在不使用绝对值时,则是:[ZPTMIN+..ZPTMAX+]或[ZPTMAX-..ZPTMIN-])参数所规定的范围内后,才执行零点自动追踪过程(可以为1,表示不需要执行连续采样)。
具体的校准方式为:在货位处于稳定状态时,且至少有预设数值(如连续采样数ZPCS)的连续采样读数落入了[ZPTMIN,ZPTMAX]范围内后,将连续采样数ZPCS个的采样数学平均数或加权平均值累加到零点偏移补偿量ZPO中,并在未来的重量计算过程中,应用该偏移值,以对传感器的零点漂移(蠕变)进行补偿。
例如:假设ZPTMIN=1.5、ZPTMAX=5、ZPTCS=100、ZPO=10,此时若在上次货品数量发生变化以后,连续100次采样的变化值(Delta)范围在[1.5..5]以内,则将这100个变化值的加权平均值(假设为2.5)累加到ZPO(新的ZPO=12.5)参数上。未来的每次新采样在完成重量转换和缩放等步骤后,均会追加该ZPO值以补偿零点漂移。
反之,在货位处于非稳定状态时,则首先要对状态的变化进行补偿。方法是:将当前状态与上一状态的差值逐一累加到当前已保存的每一个连续采样中,或将其累加到当前已保存的基准值中。
例如:假设开始零点自动追随过程,并且已存在10次有效的连续采样时,当前货位的称重读数(经过步骤一到步骤五清洗)为200.5g,货品数量为10。而从第11次采样开始,其称重读数变更为240.75g,货品数量为12。则此时状态发生了变化(货品数量由10变为了12,同时重量增加了40.25g)。为了能够动态地继续跟踪零点偏移值,此时应对前10个采样的数值逐一增加40.25,即: 以补偿其前10个采样中的非偏移性称重差异。或者,若算法仅保留一个基准重量值(而不是为每个采样保留完整的称重读数),则应对该基准重量值追加40.25。
由于蠕变或环境变化导致的零点漂移通常是以微小步进量逐步发生,因此ZPTMAX值通常可设置为一个相对当前传感器有效量程来说非常小的值,这也有助于防止系统发生误判。将ZPTMIN设置为一个非零数值既可以防止系统过于频繁地启动零点自动跟随过程,也可以保护本过程免受信号等干扰,或其它随机波动等因素的影响。较大的ZPTCS值则可保护系统不会受到由电磁干扰、惯性、离心力等情况引起的突变干扰。
优选的,所述步骤四旨在去除重量值中篮筐容器等非货品物体的重量。去皮操作通过在当前重量值上,减去一个事先指定的皮重参数BW来实现。减去BW参数后的值即为当前货位上的纯货品重量值。
优选的,所述步骤五从采样中排除由电磁干扰、惯性、离心力等现象带来的称重读数突变问题。
抗突变滤波过程主要由:抗突变最大范围AMMAX、抗突变最大范围最小范围AMMIN、抗突变比例范围AMR以及抗突变采样数AMS等参数控制。
抗突变比例范围AMR是一组连续到达的稳定采样之间的最大波动比例范围。同属于一个采样组内,不同采样之间的变化波动比例在AMR指定范围内时,才被认为该组采样是稳定的,没有发生突变。
抗突变最大范围AMMAX决定了AMR的绝对上限值,若当前重量值与AMR的积大于此值,则以此值为准(绝对值,AMMAX也可被分开指定为AMMAX+和AMMAX-:即正值最大范围和负值最大范围;可以为0,表示最大范围不限)。
抗突变最小范围AMMIN决定了AMR的绝对下限值,若当前重量值与AMR的积小于此值,则以此值为准(绝对值,AMMIN也可被分开指定为AMMIN+和AMMIN-:即正值最小范围和负值最小范围;可以为0,表示最小范围不限)。
抗突变采样数AMS规定了至少有多少个连续采样读数的绝对值完全落入了受AMMAX和AMMIN约束的AMR比例范围内后,才认为没有发生突变(可以为1,表示不需要执行连续采样)。
具体滤波方式为:将AMS次采样的数学平均数或加权平均值作为该组采样的最终重量值输出至下一步。
例如:AMR=0.1、AMMAX=30、AMMIN=1.5、AMS=5表示每5个连续采样为一组,若这5个采样的波动范围不超过±10%(同时,对于该组采样的重量值来说,其与0.1的积若小于1.5则强制设置为1.5、若大于30则强制设定为30),则认为该组采样稳定,该组采样的最终重量值为上述5个重量值的加权平均数。
显然,AMR限制的波动比例范围保证了数据的相对稳定。AMMAX则在货位重载(货品数量多,重量值读数大)的情况下帮助限制了可允许的读数最大波动范围,防止按比例限制后绝对值仍然过大(过于迟钝)的现象发生;相反AMMIN在货位轻载(货品数量少,重量值读数小)的情况下帮助限制了可允许的读数最小波动范围,防止按比例限制后绝对值过小(过于敏感)的现象发生。较大的AMR值则可保护系统不会受到由电磁干扰、惯性、离心力等情况引起的突变干扰。
优选的,步骤六数量转换过程具体为:将经过了上述步骤一至五步骤校准和纠偏后的重量值转换为大体数量。其转换过程为:
当前货位货品大体数量P=当前货位重量值/当前货位上货品的标准毛重GSW。
请参见图4所示,其中货品的标准毛重GSW是指,通过实际称量测试或通过可信渠道获得的,货品毛重标准值。它应该在计入公差影响后,能够满足正态分布条件下的中位数。意即:若随机从N件该型号产品中抽取一件,则该件货品:1.重量高于GSW的概率与重量低于GSW的概率上应当基本相等(即:均为50%);2.同时,其重量高于GSW的幅度在概率上也应与其低于GSW的幅度相当。
例如:若一件产品的GSW为100g,其标准公差为2%。从1000件该型号产品中随机抽取一件,则:1.该件产品重量大于100g或小于100g的概率均为50%;2.并且其重量为99g的概率应与其重量为101g的概率相当。
GSW的近似值可以通过从一批或多批相同型号的货品中,随机抽取M件进行称重,并求其平均值来得出,即:
GSW≈M件随机货品的总重/M
显而易见,GSW是保证数量转换准确的关键指标之一。它越接近前文中描述的正态分布中位数,就在“一个货位内,同时盛放多件相同型号货品”的情况下,消除其累积误差,使得数量计算更为精确。
可以想象,若某品牌方便面,其包装标称保重为100g,标准公差3%。但由于厂商追逐利益,实际测量后其所有产品重量均分布在97g到98g之间。此时若将该货品的GSW值设置为厂方标称的100g,则在该货位盛放100包方便面时,实际测得的结果可能只有97包。
因此GSW通常需要通过自行实际测量或从其它可信的权威渠道获取。
虽然在步骤六中引入了GSW后,已能做到比较准确地完成数量计算。但仍无法避免如±1件之类的微小偏差。为了进一步增强数量统计的准确度,在执行数量转换后,仍需进行进一步的数量精细校准过程。数量校准过程较为复杂,大体上来讲,该过程的主要动作是对当前货位中所有货品的可接受累计公差(重量值)进行合理的控制,并据此对上一步产生的数量转换结果进行必要的修正。
优选的,步骤七中,数量校准过程中涉及到的主要参数包括:
满仓容量CAP:当前货位的满仓容量(货品件数,为0表示无件数限制)。
公差范围APU:每件货品相对于其中位数的最大可接受公差范围(可以为百分比或重量绝对值)。
最大累计误差CAMAX:当前货位所有货品的最大可接受累计误差上限(可以为相对于GSW的百分比,或重量绝对值,为0表示无最大限制)。
最小累计误差CAMIN:当前货位所有货品的最小可接受累计误差下限(可以为相对于GSW的百分比,或重量绝对值,为0表示无最小限制)。
误差容限ATF:当前货位对误差的容忍度。从几何意义上说,此参数将与当前货位满仓容量一起决定逐步限制累计误差的曲线或直线陡峭程度。此值越大,对误差的宽容度越高(曲线越平滑、或直线斜率的绝对值越小,为1表示无调整)。
以及在步骤6中得到的大体数量P和标准毛重GSW。
通过对上述几个方面施加合理限制,即可计算出合理的累计误差容忍度CAT范围。进而使用CAT对实际货品数量执行进一步校准,显著提升系统准确度。
举例来说,此过程可通过综合多中参数和数据,对当前货位的可接受累积误差给出如图5所示的“倒S型”(inverse sigmoid)平滑调整曲线,基于“倒S型平滑曲线”的CAT计算方法:
1.计算“倒S型曲线”平滑下降因子ISF:ISF=1/(CAP*ATF)。
2.计算“倒S型曲线”缩放补偿值ISC:ISC=1/(1/(1+exp(ISF)))。
3.将由步骤六中得到的货品大体数量值P带入,计算其误差容忍中间系数ISITF:ISITF=(1/(1+exp(P*ISF)))*ISC。
以及其在最差情况下,可能造成的累积误差最大单位(件数)比例ISWCT:ISWCT=P*APU。
4.计算其受限后的最大累计误差单位比例ISCWCT:ISCWCT=ISITF*ISWCT。
5.计算累计公差容忍度CAT:
a)若CAMIN不为零,且ISCWCT<CAMIN(此时CAMIN应转换为相对于GSW的百分比),则:CAT=CAMIN/ISWCT。此举确保累计公差范围满足CAMIN指定的最小范围约束。
b)若CAMAX不为零,且ISCWCT>CAMAX(此时CAMAX应转换为相对于GSW的百分比),则:CAT=CAMAX/ISWCT。此举确保累计公差范围满足CAMAX指定的最大范围约束。
c)否则说明ISCWCT介于[CAMIN..CAMAX]之间,此时CAT受上述“倒S型曲线”控制:CAT=ISITF。
6.计算实际允许的货位总公差(百分比)ISRPA:ISRPA=CAT*APU。
7.计算实际误差容忍范围(重量值)ISAAR:ISAAR=经过了上述步骤一至五校准和纠偏后的重量值*ISRPA。
8.计算步骤六中得到的货品大体数量值P(P为整数)后的称重值余数ISWCO:ISWCO=经过了上述步骤一至五校准和纠偏后的重量值–P*GSW。
9.若ISWCO+ISAAR≥GSW,则执行补偿:P=P+1,为货品数量值加1。
注:本文中用到的数学算式与C语言中的数学表达式语法基本相同,即:“+”表示加法、“-”表示减法、“*”表示乘法、“/”表示除法、exp表示“以自然数e为底的幂运算”,例如:exp(10)表示e的10次方、exp(N)表示e的N次方。
应当指出的是,上例中仅列出了对包括P、GSW、CAP、APU、CAMAX、CAMIN、ATF在内的各项校准方法的一种合理应用。实际上,我们可以有任意多中方式对上述限制方法进行合理应用。例如:将上述算式中的自然数e替换为圆周率π或任意其它实数、将上述算式中的幂运算替换为对数运算、将上述非线性(曲线)计算替换为线性(直线)计算等等。
因此上例中提到的所有算式以及“倒S型曲线”算法均仅为举例,并不对本发明的范围有任何限制作用。
可以看出,CAP和ATF从几何意义上来说,主要用来控制容忍度曲线的平滑程度(或容忍度直线的斜率绝对值)。而CAMIN和CAMAX则用于对曲线进行适当缩放,保证货位在轻载和重载条件下也能保持合理的CAT值。而GSW和APU则决定了最终误差范围的实际值。
综上,本发明计算精准,著提升了基于称重的货位管理技术在实际工作场景中的精度和稳定性。
Claims (5)
1.一种基于称重的货品数量计算方法,其特征在于,包括数量转换方法:货品大体数量P=重量值W/货品的标准毛重GSW,P为整数;以及数量校准方法:对所述货品大体数量P进行校准;通过对满仓容量CAP、货品的标准毛重GSW、公差范围APU、最大累计误差CAMAX、最小累计误差CAMIN以及误差容限ATF的限定,计算累计误差容忍度CAT的范围;进而使用累计误差容忍度CAT对计算出的货品大体数量进一步校准,具体步骤如下:
(1).计算“倒S型曲线”平滑下降因子ISF,ISF=1/(CAP*ATF);
(2).计算“倒S型曲线”缩放补偿值ISC,ISC=1/(1/(1+exp(ISF)));
(3).将P带入,计算误差容忍中间系数ISITF,ISITF=(1/(1+exp(P*ISF)))*ISC,以及累积误差最大单位比例ISWCT,ISWCT=P*APU;
(4).计算受限后的最大累计误差单位比例ISCWCT,ISCWCT=ISITF*ISWCT;
(5).计算累计公差容忍度CAT:若最小累计误差CAMIN不为零,且ISCWCT<CAMIN,则CAT=CAMIN/ISWCT;若最大累计误差CAMAX不为零,且ISCWCT>CAMAX,则CAT=CAMAX/ISWCT;否则说明ISCWCT介于CAMIN与CAMAX之间,此时CAT=ISITF;
(6).计算实际允许的货位总公差ISRPA,ISRPA=CAT*APU;
(7).计算实际误差容忍范围重量值ISAAR,ISAAR=W*ISRPA;
(8).计算称重值余数ISWCO,ISWCO=W–P*GSW;
(9).若ISWCO+ISAAR≥GSW,则执行补偿:P=P+1,为货品数量值加1。
2.根据权利要求1所述的一种基于称重的货品数量计算方法,其特征在于:还包括重量缩放方法,对称重单位进行转换;转换公式为:SW=WS*W(WU);其中SW为缩放后重量,WS为事先指定的可变系数,WU为单位。
3.根据权利要求1所述的一种基于称重的货品数量计算方法,其特征在于:还包括去皮方法,去除重量值中非货品物体的重量,将称得的重量减去预设的皮重参数BW。
4.根据权利要求1所述的一种基于称重的货品数量计算方法,其特征在于:将步骤(2)算式中的自然数e替换为任意其它实数、将步骤(2)算式中的幂运算替换为对数运算或线性运算。
5.根据权利要求1所述的一种基于称重的货品数量计算方法,其特征在于:将步骤(3)算式中的自然数e替换为任意其它实数、将步骤(3)算式中的幂运算替换为对数运算或线性运算。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810853293.6A CN109145398B (zh) | 2018-07-30 | 2018-07-30 | 一种基于称重的货品数量计算方法 |
PCT/CN2019/097853 WO2020024872A1 (zh) | 2018-07-30 | 2019-07-26 | 一种基于称重的货品数量计算方法 |
US17/263,400 US20210156727A1 (en) | 2018-07-30 | 2019-07-26 | Method for calculating the quantity of goods based on weighing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810853293.6A CN109145398B (zh) | 2018-07-30 | 2018-07-30 | 一种基于称重的货品数量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109145398A CN109145398A (zh) | 2019-01-04 |
CN109145398B true CN109145398B (zh) | 2024-01-05 |
Family
ID=64798448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810853293.6A Active CN109145398B (zh) | 2018-07-30 | 2018-07-30 | 一种基于称重的货品数量计算方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20210156727A1 (zh) |
CN (1) | CN109145398B (zh) |
WO (1) | WO2020024872A1 (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063325B (zh) * | 2018-07-30 | 2024-02-23 | 白杨 | 一种基于称重的货品数量计算方法 |
CN109145398B (zh) * | 2018-07-30 | 2024-01-05 | 白杨 | 一种基于称重的货品数量计算方法 |
CN110059998A (zh) * | 2019-03-11 | 2019-07-26 | 北京辰森世纪科技股份有限公司 | 自动盘点货物的方法及装置、系统、存储介质、电子装置 |
CN109883529A (zh) * | 2019-04-04 | 2019-06-14 | 杭州点微科技有限公司 | 一种称重存取计数的抗抖方法 |
CN110095175A (zh) * | 2019-05-08 | 2019-08-06 | 北京汇通天下物联科技有限公司 | 载重设备参数校准方法、装置及系统 |
CN113554802B (zh) * | 2021-04-21 | 2023-06-02 | 浙江星星冷链集成股份有限公司 | 一种称重装置及无人售货系统 |
CN116412892B (zh) * | 2023-03-17 | 2024-04-16 | 江苏云涌电子科技股份有限公司 | 一种智能称重货架零点漂移补偿系统及方法 |
Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1025781B (de) * | 1955-06-09 | 1958-03-06 | Pneumatic Scale Corp | Kontinuierlich arbeitende Verpackungsmaschine fuer Schuettgut mit Mengen- und Gewichtskorrektureinrichtung |
US4157738A (en) * | 1977-05-15 | 1979-06-12 | Shinko Denshi Company Limited | Method for counting the number of articles using a weighing machine |
US4316516A (en) * | 1979-12-20 | 1982-02-23 | Mettler Instrumente Ag | Weighing balance with automatic zero correction |
US4412298A (en) * | 1979-09-20 | 1983-10-25 | Pitney Bowes Inc. | Method for tracking creep and drift in a digital scale under full load |
US4645019A (en) * | 1980-07-01 | 1987-02-24 | Kabushiki Kaisha Ishida Koki Seisakusho | Method of weighting and counting |
EP0312272A2 (en) * | 1987-10-09 | 1989-04-19 | Ishida Scales Mfg. Co. Ltd. | Combinatorial counting method and weighing apparatus therefor |
US5193629A (en) * | 1988-06-07 | 1993-03-16 | Percell Group Limited | Method and apparatus for determining the number and/or value of weighed items |
JP2010008214A (ja) * | 2008-06-26 | 2010-01-14 | Teraoka Seiko Co Ltd | 秤装置 |
JP2010025687A (ja) * | 2008-07-17 | 2010-02-04 | Teraoka Seiko Co Ltd | 秤装置 |
FR2938644A1 (fr) * | 2008-11-14 | 2010-05-21 | Balea | Procede et dispositif de controle de comptage d'articles par pesage |
CN202119518U (zh) * | 2011-05-26 | 2012-01-18 | 深圳市合智创盈电子有限公司 | 一种货品实时管理系统及电子计算装置 |
CN102359813A (zh) * | 2011-08-15 | 2012-02-22 | 孙冉 | 皮带秤称重传感器校准方法 |
EP2469248A1 (fr) * | 2010-12-27 | 2012-06-27 | Balea | Procédé et dispositif de comptage d'articles |
CN102997983A (zh) * | 2012-10-12 | 2013-03-27 | 孙新利 | 基于电子天平的硅片自动计片装置 |
CN103148920A (zh) * | 2012-08-17 | 2013-06-12 | 苏州艾隆科技股份有限公司 | 盘点机 |
CN104570839A (zh) * | 2014-12-03 | 2015-04-29 | 杭州信雅达科技有限公司 | 托盘设备、装配线系统及补料方法 |
CN204807179U (zh) * | 2015-06-27 | 2015-11-25 | 东莞市华兰海电子有限公司 | 可显示货物数量和重量的装置 |
CN107402059A (zh) * | 2017-07-28 | 2017-11-28 | 杭州高达软件系统股份有限公司 | 一种行车货品称重方法 |
CN108268930A (zh) * | 2018-01-31 | 2018-07-10 | 华北理工大学 | 棒材称重计数器以及称重计数方法 |
Family Cites Families (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1816676A1 (de) * | 1967-12-23 | 1969-07-31 | Zaklady Mekh Precyzyjnej Przed | Einrichtung zur Ablesung der Waagebalkenstellung sowie zum Auflegen und Abnehmen der Gewichte in einer selbsttaetigen Neigungswaage |
CA991661A (en) * | 1973-12-07 | 1976-06-22 | Joseph T. Sniezek | Endless conveyor belt load measurement system and method of automatically calibrating same |
DE3176893D1 (en) * | 1980-05-20 | 1988-11-03 | Percell Group Ltd | Method of measuring by weighing |
US4393950A (en) * | 1981-06-26 | 1983-07-19 | Triangle Package Machinery Company | Tare weighing apparatus and method therefor |
JPS5850429A (ja) * | 1981-09-21 | 1983-03-24 | Ishida Scales Mfg Co Ltd | 計数方法 |
US4406338A (en) * | 1982-02-05 | 1983-09-27 | Dunbar Glenn G | Weight cart |
US4535854A (en) * | 1982-11-16 | 1985-08-20 | Pennsylvania Scale Company | Calibration method and apparatus for an electronic weight indicator |
JPS59200923A (ja) * | 1983-04-28 | 1984-11-14 | Ishida Scales Mfg Co Ltd | スパン調整方法 |
JPS6011119A (ja) * | 1983-06-30 | 1985-01-21 | Shimadzu Corp | 電子天びん |
DE3340421A1 (de) * | 1983-11-09 | 1985-05-15 | Sartorius GmbH, 3400 Göttingen | Elektrische waage |
US4637481A (en) * | 1983-11-17 | 1987-01-20 | Shoemaker Roger E | Self-propelled scale testing vehicle |
US4691290A (en) * | 1984-08-06 | 1987-09-01 | Reliance Electric Company | Creep-compensated weighing apparatus |
JPS61178617A (ja) * | 1985-02-01 | 1986-08-11 | Yamato Scale Co Ltd | 秤のスパン自動調整方法及びその装置 |
US4667153A (en) * | 1985-05-17 | 1987-05-19 | Doyle James H | Automatic calibrator |
DE3530170C2 (de) * | 1985-08-23 | 1994-03-10 | Sartorius Gmbh | Zählwaage |
US4727947A (en) * | 1985-09-03 | 1988-03-01 | Ishida Scales Manufacturing Company, Ltd. | Span adjustment system and method |
US4660662A (en) * | 1985-09-11 | 1987-04-28 | National Controls, Inc. | Digital electronic scale with stabilized display |
US4630696A (en) * | 1985-10-25 | 1986-12-23 | Azionaria Costruzioni Macchine Automatiche | Apparatus and method for automatic system calibration to provide enhanced resolution in computerized weighing systems |
CH670508A5 (zh) * | 1986-05-23 | 1989-06-15 | Mettler Instrumente Ag | |
CH670507A5 (zh) * | 1986-05-23 | 1989-06-15 | Mettler Instrumente Ag | |
US4751661A (en) * | 1986-05-27 | 1988-06-14 | Ncr Corporation | Automatic zero balancing of a weighing apparatus |
US4715457A (en) * | 1986-05-27 | 1987-12-29 | Ncr Corporation | Method for fast weighing |
US4909338A (en) * | 1989-06-12 | 1990-03-20 | Ncr Corporation | Method and apparatus for scale calibration and weighing |
JPH0652190B2 (ja) * | 1989-07-24 | 1994-07-06 | 株式会社島津製作所 | 電子天びん |
US6794586B1 (en) * | 2002-06-27 | 2004-09-21 | Ncr Corporation | System and method for controlling weighing operations at a scale of a self-checkout terminal |
PL1736745T3 (pl) * | 2005-06-21 | 2012-12-31 | Mettler Toledo Gmbh | Sposób adaptacyjnej korekcji zjawisk dryftu w urządzeniu do pomiaru siły oraz urządzenie do pomiaru siły dla realizacji sposobu |
US7441436B2 (en) * | 2005-09-27 | 2008-10-28 | Mettler-Toledo, Inc. | Method for weighing apparatus calibration management |
JP4940026B2 (ja) * | 2007-06-15 | 2012-05-30 | 株式会社タニタ | 生体測定装置 |
GB2519295A (en) * | 2013-10-15 | 2015-04-22 | Cashmaster Internat Ltd | Improvements in or relating to counting by weight counting systems |
CN106979812A (zh) * | 2016-01-15 | 2017-07-25 | 广西大学 | 一种高速公路车辆动态称重系统 |
CN106208050B (zh) * | 2016-08-17 | 2018-12-18 | 华北电力大学 | 一种基于pmu的电网支路静态参数检测辨识方法 |
CN109145398B (zh) * | 2018-07-30 | 2024-01-05 | 白杨 | 一种基于称重的货品数量计算方法 |
CN109063325B (zh) * | 2018-07-30 | 2024-02-23 | 白杨 | 一种基于称重的货品数量计算方法 |
-
2018
- 2018-07-30 CN CN201810853293.6A patent/CN109145398B/zh active Active
-
2019
- 2019-07-26 US US17/263,400 patent/US20210156727A1/en active Pending
- 2019-07-26 WO PCT/CN2019/097853 patent/WO2020024872A1/zh active Application Filing
Patent Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1025781B (de) * | 1955-06-09 | 1958-03-06 | Pneumatic Scale Corp | Kontinuierlich arbeitende Verpackungsmaschine fuer Schuettgut mit Mengen- und Gewichtskorrektureinrichtung |
US4157738A (en) * | 1977-05-15 | 1979-06-12 | Shinko Denshi Company Limited | Method for counting the number of articles using a weighing machine |
US4412298A (en) * | 1979-09-20 | 1983-10-25 | Pitney Bowes Inc. | Method for tracking creep and drift in a digital scale under full load |
US4316516A (en) * | 1979-12-20 | 1982-02-23 | Mettler Instrumente Ag | Weighing balance with automatic zero correction |
US4645019A (en) * | 1980-07-01 | 1987-02-24 | Kabushiki Kaisha Ishida Koki Seisakusho | Method of weighting and counting |
EP0312272A2 (en) * | 1987-10-09 | 1989-04-19 | Ishida Scales Mfg. Co. Ltd. | Combinatorial counting method and weighing apparatus therefor |
US5193629A (en) * | 1988-06-07 | 1993-03-16 | Percell Group Limited | Method and apparatus for determining the number and/or value of weighed items |
JP2010008214A (ja) * | 2008-06-26 | 2010-01-14 | Teraoka Seiko Co Ltd | 秤装置 |
JP2010025687A (ja) * | 2008-07-17 | 2010-02-04 | Teraoka Seiko Co Ltd | 秤装置 |
FR2938644A1 (fr) * | 2008-11-14 | 2010-05-21 | Balea | Procede et dispositif de controle de comptage d'articles par pesage |
EP2469248A1 (fr) * | 2010-12-27 | 2012-06-27 | Balea | Procédé et dispositif de comptage d'articles |
CN202119518U (zh) * | 2011-05-26 | 2012-01-18 | 深圳市合智创盈电子有限公司 | 一种货品实时管理系统及电子计算装置 |
CN102359813A (zh) * | 2011-08-15 | 2012-02-22 | 孙冉 | 皮带秤称重传感器校准方法 |
CN103148920A (zh) * | 2012-08-17 | 2013-06-12 | 苏州艾隆科技股份有限公司 | 盘点机 |
CN102997983A (zh) * | 2012-10-12 | 2013-03-27 | 孙新利 | 基于电子天平的硅片自动计片装置 |
CN104570839A (zh) * | 2014-12-03 | 2015-04-29 | 杭州信雅达科技有限公司 | 托盘设备、装配线系统及补料方法 |
CN204807179U (zh) * | 2015-06-27 | 2015-11-25 | 东莞市华兰海电子有限公司 | 可显示货物数量和重量的装置 |
CN107402059A (zh) * | 2017-07-28 | 2017-11-28 | 杭州高达软件系统股份有限公司 | 一种行车货品称重方法 |
CN108268930A (zh) * | 2018-01-31 | 2018-07-10 | 华北理工大学 | 棒材称重计数器以及称重计数方法 |
Non-Patent Citations (1)
Title |
---|
《基因单片机的智能型称重计数器的研究》;秦乐园;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20150615;C030-7 * |
Also Published As
Publication number | Publication date |
---|---|
US20210156727A1 (en) | 2021-05-27 |
WO2020024872A1 (zh) | 2020-02-06 |
CN109145398A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109063325B (zh) | 一种基于称重的货品数量计算方法 | |
CN109145398B (zh) | 一种基于称重的货品数量计算方法 | |
CN203241127U (zh) | 温度自动补偿电子秤及校准系统 | |
US8063320B2 (en) | Counting scale with a moving average start calculating section, for determining whether or not to calculate a moving average | |
EP2336735B1 (en) | Measuring device | |
US5300736A (en) | Adaptive timing in-motion checkweigher | |
FI70485B (fi) | Maetningsfoerfarande foer impedanser saerskilt smao kapacitanser vid vilket man anvaender en eller flera referenser | |
CN1884982B (zh) | 测量传感器的输出信号的处理方法以及测力装置 | |
CN110849459B (zh) | 一种称重传感器的蠕变修正方法 | |
CN112304412B (zh) | 称重装置的滞后补偿方法 | |
CN113701859A (zh) | 称重计数方法、装置、系统、计算机设备及可读存储介质 | |
KR101239354B1 (ko) | 글로브박스 내의 기압변화 보정기능을 가지는 전자저울 | |
US4703815A (en) | Span adjusting device for weigher | |
CN114910146B (zh) | 一种猪场料塔称重模拟量传感器失效后自动测算重量方法 | |
CN109870227A (zh) | 一种货品存取称重计数的自动去皮方法 | |
CN109696897B (zh) | 一种多点顶升机构载荷检测方法及其检测装置 | |
CN112082633A (zh) | 一种在人体晃动情况下稳定体重值的测重方法 | |
CN113945262A (zh) | 偏载误差修正方法和系统 | |
JPS6314885B2 (zh) | ||
US4417632A (en) | Automatic weighing method and device | |
US4792003A (en) | Span adjusting device for weigher | |
CN109489792A (zh) | 一种原料配料精度在线监控方法 | |
CN217358734U (zh) | 一种用于电子秤全周倾斜称重误差补偿的装置 | |
CN113790791A (zh) | 一种重力售货柜的重力传感器标定方法 | |
CN117782274A (zh) | 一种应变片电子天平及其修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |