发明名称 水果高速分选系统动态称重信号的智能处理方法
摘要 本发明公开了一种水果高速分选系统动态称重信号的智能处理方法,基于改进AR模型,实现水果重量的高精度获取与预测。包括以下步骤:实时获取分级速度,采样存储一组标准果数据;滤波消除信号采样值中的低频噪声;操作站系统读取初始数据,用改进的AR模型建模并计算模型阶数和参数;将不同分选速度下训练得到的模型参数和阶数存储于文件中;控制站系统根据分选速度利用训练好的模型参数预测水果重量信息;操作站系统可实现对预测结果的精度分析。本发明在每秒20个果的分级速度下,称重精度可达(1.5±0.5)%以内,改变分级速度后可重新获取最优的模型阶数和参数,对硬件要求低,结构简单,扩展性好,且可方便移植用于其它领域高速称重信号的处理。
申请公布号 CN102176117B 申请公布日期 2012.12.05
申请号 CN201110022875.8 申请日期 2011.01.20
申请人 浙江大学 发明人 何慧梅;黄平捷;张光新;侯迪波;刘喆;蔡文;田径;周泽魁
分类号 G05B13/04(2006.01)I;B07C5/28(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 张法高
主权项 1.一种水果高速分选系统动态称重信号的智能处理方法,其特征在于它的步骤如下:1)水果输送装置静止时,采样获得一组标准果的稳定值,然后传送装置开始转动,根据分级速度,确定采样周期,采样获得每个标准果经过称重传感器时的一组数据;2)采用一阶滞后滤波,消除采样信号中的低频噪声,获得预处理后的信号;3)操作站系统读取预处理数据,分析数据的特点,根据实时性要求,得到改进的AR模型;4)确定改进AR模型的参数和阶数,将不同分级速度下训练得到的模型参数和阶数存储于文件中;5)控制站系统利用训练好的模型参数预测水果重量信息;6)对预测结果进行精度评价;所述步骤1)为:需要n个重量不同的标准果,在静态条件下,对每一个标准果多次采样后的值取平均,作为果子最终的静态采样值,分别记为:b<sub>1</sub>,b<sub>2</sub>,…,b<sub>n</sub>;在高速情况下,根据分级速度,确定采样周期:设称重传感器的承重脚长度为L,分级速度为V,信号从开始上升到开始下降的时间间隔t<sub>AC</sub>=L/V,采样只能在t<sub>AC</sub>时间长度内完成,若每个果子的采样次数为m,则采样周期t<sub>s</sub><L/(V*m);最后按设定的采样周期采样,存储每个标准果的采样数据y<sub>1</sub>(1),y<sub>1</sub>(2),…,y<sub>1</sub>(m);y<sub>2</sub>(1),y<sub>2</sub>(2),…,y<sub>2</sub>(m);……;y<sub>n</sub>(1),y<sub>n</sub>(2),…,y<sub>n</sub>(m);所述步骤2)为:利用一阶滞后滤波消除采样数据中的低频噪声,定义y(k)为第k次采样值,x(k-1)为第k-1次的滤波输出值,α为滤波系数,其值小于1,代入一阶滞后滤波公式x(k)=a*y(k)+(1-a)*y(k-1),得到第k次采样数据对应的滤波输出x(k);将每组采样数据滤波处理后,得到一组数据x<sub>1</sub>(1),x<sub>1</sub>(2),…,x<sub>1</sub>(m);x<sub>2</sub>(1),x<sub>2</sub>(2),…,x<sub>2</sub>(m);……;x<sub>n</sub>(1),x<sub>n</sub>(2),…,x<sub>n</sub>(m),作为模型训练的数据;所述步骤3)为:称重传感器的输出信号近似为平稳随机的时间序列,一个时间序列的AR模型描述为:序列的当前值x(n)被表示为过去值的线性函数加上一个误差项e(n),用如下方程式(1)表示:x(n)=a<sub>1</sub>x(n-1)+a<sub>2</sub>x(n-2)+...+a<sub>p</sub>x(n-p)+e(n)    (1)方程式(1)中a<sub>1</sub>,a<sub>2</sub>,…,a<sub>p</sub>为模型系数;x(n-1),x(n-2),…,x(n-p)为信号p个过去值,构成p阶AR模型;现假设已知p个初始采样值x(1),x(2),…,x(p),信号N个时刻后达到稳定,要预测信号稳定值x(p+N),利用方程式(1)递推得到如下方程式(2)为:x(p+N)=a<sub>1</sub>x(N)+a<sub>2</sub>x(N+1)+...+a<sub>p</sub>x(N+p-1)=A<sub>1</sub>x(1)+A<sub>1</sub>x(2)+...+A<sub>p</sub>x(p)        (2)方程式(2)中系数A<sub>1</sub>,A<sub>2</sub>,…,A<sub>p</sub>为常数,信号稳定值x(p+N)由初始采样值x(1),x(2),…,x(p)线性表示;对于不同重量的果子,经过同一个称重传感器,信号稳定值与初始采样值的线性关系一致;所述步骤4)为:利用操作站系统计算改进AR模型阶数和参数,首先操作站系统读取模型训练数据x<sub>1</sub>(1),x<sub>1</sub>(2),…,x<sub>1</sub>(m);x<sub>2</sub>(1),x<sub>2</sub>(2),…,x<sub>2</sub>(m);……;x<sub>n</sub>(1),x<sub>n</sub>(2),…,x<sub>n</sub>(m),以及静态稳定值b<sub>1</sub>,b<sub>2</sub>,…,b<sub>n</sub>,存储在矩阵变量中,记为:<img file="FDA00001629685600021.GIF" wi="1162" he="287" />其中m代表每个果子的采样次数;n代表用于模型参数预测的标准果数量,A矩阵每一行代表一组模型训练的初始数据,A矩阵第i行对应的信号稳定值为矩阵B的第i行;先假设模型的阶数为p,定义另外一组变量,分别为:<img file="FDA00001629685600022.GIF" wi="810" he="293" /><img file="FDA00001629685600023.GIF" wi="913" he="299" /><maths num="0001"><![CDATA[<math><mrow><msub><mi>X</mi><mi>p</mi></msub><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>p</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>Y</mi><mi>p</mi></msub><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mi>p</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>E</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>e</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>e</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,A<sub>p1</sub>用于确定模型参数的训练数据,B<sub>p1</sub>为训练数据A<sub>p1</sub>对应的稳定值,A<sub>p2</sub>用于模型阶数确定时估计模型误差,B<sub>p2</sub>为A<sub>p2</sub>数据对应的稳定值,X<sub>p</sub>为阶数为P时的模型参数,Y<sub>p</sub>用于存储A<sub>p2</sub>计算得到的预测数据,e<sub>i</sub>用于存储阶次为i时预测值与实际稳定值的标准差,作为模型误差;将方程式(2)得到的改进AR模型用于模型参数计算如下:A<sub>p1</sub>*X<sub>p</sub>=B<sub>p1</sub>    (3)计算线性系数X<sub>p</sub>,然后用该系数对A<sub>p2</sub>中数据进行预测,令Y<sub>p</sub>=A<sub>p2</sub>*X<sub>p</sub>     (4)得到一组预测值Y<sub>p</sub>,再比较和Y<sub>p</sub>和B<sub>p2</sub>的差异,计算p阶模型时的标准差,计算公式为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>e</mi><mi>p</mi></msub><mo>=</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>q</mi></munderover><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>b</mi><mrow><mi>p</mi><mo>+</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><mi>q</mi></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>计算得到一组模型阶数对应的模型误差,取e<sub>i</sub>最小时对应的阶数i,作为模型的阶数,同时计算阶数为i时的模型参数X<sub>i</sub>;最后将模型的阶数和参数,作为当前采样速度下的最优系数,传给控制站系统,用于水果重量信号采样稳定值的预测;在系统初始化时,重复上述步骤进行多次训练,获得不同分级速度下对应的模型参数和阶数,用文件存储;当根据现场条件需要改变分级速度时,操作站系统选择直接将文件中相应分级速度下的模型参数和阶数传给控制站系统,或重新计算模型的参数和阶数,传给控制站系统;若现场条件改变,重新计算的值与文件中存储的值不一致,可选择是否更新文件;所述步骤5)为:预测称重信号稳定值,首先获取分级速度,确定该分级速度对应的模型参数和阶数已经存储在控制站系统中;当水果经过称重传感器时,由外部同步信号上升沿触发控制站系统的中断,采样次数由模型阶数确定;当采样获得的一组数据经过滤波后,控制站系统用训练好的AR模型进行预测,得到信号稳定值,标定后转化为重量信息;所述步骤6)为:随机取K个果子,确定采样周期,采样获得K组高速动态数据,用训练好的模型进行稳定值预测,得到K个预测数据,标定后转化为K个预测重量,计算得到预测重量与实际重量的偏差,相对误差,预测精度。
地址 310027 浙江省杭州市西湖区浙大路38号