发明名称 运输安全记录仪数据分析方法
摘要 本发明涉及一种运输安全记录仪数据分析方法,其方法的步骤主要包括采集数据导入、采集数据的时域预处理,剔除数据野点,去除均值、低通滤波滤掉高频信号、计算加速度均方根、计算冲击持续时间、得出功率谱密度、得出冲击持续时间速度变化量~加速度峰值统计点散图及数据结果应用。本发明方法对采集的时域信号首先进行时域预处理,低通滤波处理,为采集数据的分析打下了基础,本发明方法全面分析运输过程中货物的动态状况,及时判断货物的受损情况,为提高物流管理水平的提高打下了基础。
申请公布号 CN102890750B 申请公布日期 2015.06.10
申请号 CN201210333408.1 申请日期 2012.09.11
申请人 中国包装科研测试中心 发明人 韩雪山;王刚;张子汉;许萍萍;李宇杰;何渊井;陈志强;陈振强;张卫红
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 天津盛理知识产权代理有限公司 12209 代理人 韩奎勇
主权项 一种运输安全记录仪数据分析方法,该方法运行的硬件系统包括安全记录仪及分析主机,安全记录仪通过USB端口将采集到的加速度数据传送给分析主机,其特征在于:分析主机对采集数据的分析方法步骤如下:⑴将安全记录仪数据线插入运行分析主机USB端口,将安全记录仪采集的数据导入分析主机;⑵采集数据的时域预处理,剔除数据野点,零均值处理;其中,所述采集数据的时域预处理,剔除数据野点,零均值处理的具体步骤是:①根据带限比较法,将安全记录仪加速度采集的数据依时间顺序设置采样值时间序列{x<sub>i</sub>},i=1,2,3,...,N,计算采样值时间序列的方均与均方差值:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mn>2</mn></msup><mo>=</mo><mover><msubsup><mi>x</mi><mi>i</mi><mn>2</mn></msubsup><mo>&OverBar;</mo></mover><mo>-</mo><msup><mover><msub><mi>x</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000660762130000011.GIF" wi="297" he="84" /></maths>②数据野点判断:如果,x<sub>i</sub>‑kS<sub>i</sub><x<sub>i+1</sub><x<sub>i</sub>+kS<sub>i</sub>则x<sub>i+1</sub>是正常点,否则,x<sub>i+1</sub>为野点,对于判定为野点的数据将其从采集数据中剔除,系数k值在3至5之间取值;③零均值处理:均值即采集信号的直流干扰项,根据信号的均值估计,在进行统计量计算之前,将均值去除,设{x<sub>i</sub>},i=1,2,3,...,N,为加速度的采样值时间序列,记采样值均值为<img file="FDA0000660762130000012.GIF" wi="73" he="76" />有:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>x</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000660762130000013.GIF" wi="250" he="143" /></maths>零均值处理后的新信号序列为{u<sub>i</sub>},i=1,2,3,...,N,有<img file="FDA0000660762130000014.GIF" wi="243" he="83" />且新信号序列均值为零,即<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>u</mi><mo>&OverBar;</mo></mover><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000660762130000015.GIF" wi="177" he="81" /></maths>⑶低通滤波:采用巴特沃思低通滤波器滤除对货物损害极小的高频信号;⑷计算加速度均方根;⑸得出功率谱密度:运用周期图法得出振动加速度信号的功率谱密度;其中,所述运用周期图法得出振动加速度信号的功率谱密度为:①加窗处理:为抑制功率泄漏,采用汉宁窗对采集的加速度时域信号进行采样截断;②傅立叶变换:将加窗后的加速度信号进行傅立叶变换,得到幅值谱;③得出功率谱密度:将傅立叶变换后幅值谱中的每个幅值数据平方后除以频率变换量,既得出最终的功率谱密度结果;⑹计算冲击持续时间:取加速度信号冲击响应峰值的第一个5%‑20%信号所对应的时间为冲击计时的起点,取加速度信号冲击响应峰值的第二个5%‑20%信号所对应的时间为冲击计时的终点,计算出冲击持续时间;⑺得出冲击持续时间速度变化量~加速度峰值统计点散图:以冲击持续时间内速度变化量为横轴、冲击峰值即每个加速度信号的最大值为纵轴统计出点散图;⑻数据结果应用:包括加速度均方根分析、功率谱密度应用及冲击持续时间速度变化量~加速度峰值统计点散图的应用;其中,所述加速度均方根分析、功率谱密度应用及冲击持续时间速度变化量~加速度峰值统计点散图的应用为:①加速度均方根分析:通过加速度均方根数值的大小判断货物的损坏程度;②功率谱密度应用:将功率谱密度值输入实验室振动台的程序控制计算机,实现等效路试的实验室随机振动试验,达到准确模拟测试的目的;③冲击持续时间速度变化量~加速度峰值统计点散图的应用:对比冲击持续时间速度变化量~加速度峰值统计点散图与产品的破损边界曲线图,如果点散图中的统计点大量落在产品的破损边界曲线的破损区域内,那么产品就会受到损害,如果点散图中的统计点几乎落在产品的破损边界曲线的破损区域外,则产品状态良好。
地址 300457 天津市滨海新区黄海路海川街2号