本发明公开一种基于多元线性回归模型的皮带秤主要误差因素分析的方法，首先搭建皮带秤实验平台，根据实际情况改变主要误差因素张力、温度、等效流量的大小，记录各个测试点的张力、温度、等效流量、标定值及料斗称，求出标定值与料斗称的相对误差；然后计算皮带秤主要误差因素之间的相关系数，判断彼此之间的相关性；再根据相关性确定各因素之间的关系，设定多元线性模型并求解；最后以残差平方和、决定系数和MS残差方差作为检验指标，评价所建立的模型的回归效果；运用拟合模型预测各个测试点的测试值，并与实际值进行对比计算出预测误差，确定拟合结果的准确性。为定量分析各个主要误差影响因素对皮带秤精度的影响程度提供理论基础。
1.基于多元线性回归模型的皮带秤主要误差因素分析的方法，其特征在于，包括如下步骤：步骤1、搭建皮带秤实验平台，根据实际情况改变主要误差因素张力、温度、等效流量的大小，记录各个测试点的张力、温度、等效流量、标定值及料斗称，求出标定值与料斗称的相对误差，即皮带秤精度值；
步骤2、计算皮带秤主要误差因素之间的相关系数，判断彼此之间的相关性，具体的选用皮尔森相关系数统计主要误差因素之间的相关系数r，用来描述两个变量之间的线性相关程度，计算公式如下：式中，n为样本个数，X、Y为两变量的观测值， 为两变量的均值，SX、SY为两变量的标准差；
步骤3、当各主要误差因素成弱相关或不相关，则设定多元线性回归模型并求解，求解多元线性回归模型的具体方法为：步骤3.1、设随机变量y与m个自变量存在线性关系，表示为：式中，β0,β1,…,βm为回归系数，ξ是随机误差，σ是标准差，x1，x2…xm表示自变量；
y～(β0+β1x1+β2x2+…+βmxm,σ2)设有n组样本观测数据，即：
y＝(y1,y2,…,yn)T, β＝(β0,β1,…,βm)T,ξ＝(ξ1,ξ2,…,ξn)T式中，xi1 xi2 … xim表示第n个样本的第m个自变量；
步骤3.2、设b0,b1,…,bm分别为β0,β1,…,βm的最小二乘估计值，则y的观测值为：yk＝b0+b1xk1+b2xk2+…+bmxkm+ek,k＝1,2,…,nyk的估计值为：
为了使实际值与估计值的拟合程度最好，应使下式值达到最小：即转换为求解下面的方程组：
求解出方程组即得b0,b1,…,bm，从而求得回归模型：步骤4、以残差平方和、决定系数和MS残差方差作为检验指标，评价所建立的模型的回归效果；
步骤5、运用拟合模型预测各个测试点的测试值，并与实际值进行对比计算出预测误差，确定拟合结果的准确性。
2.根据权利要求1所述的基于多元线性回归模型的皮带秤主要误差因素分析的方法，其特征在于，所述步骤2中相关系数r的绝对值越接近1，相关度越强，规定：相关系数在0.8-
1.0范围，表示极强相关；0.6-0.8，表示强相关；在0.4-0.6范围，表示中等程度相关；在0.2-
0.4范围，表示弱相关；在0-0.2范围，表示极弱相关或无相关。
3.根据权利要求1所述的基于多元线性回归模型的皮带秤主要误差因素分析的方法，其特征在于，所述步骤4中MS残差方差值、残差平方和SS残、决定系数R2的计算公式如下：式中，Yi表示数据点，表示回归直线上相应的位置值；
式中，决定系数R2的值接近于1，说明回归方程的效果越好，表示样本数据能较好地拟合选用的线性回归模型；SS总为总平方和，SS回为回归平方和，df残是残差自由度。
基于多元线性回归模型的皮带秤主要误差因素分析的方法
技术领域
[0001] 本发明属于制造装备检测、控制、诊断与维护技术领域，特别涉及一种基于多元线性回归模型的皮带秤主要误差因素分析的方法。
背景技术
[0002] 在动态称重中发挥着重要作用的电子皮带秤可以对散状物料进行自动连续累计计量，因这些独有的特点使其在工业生产以及港口贸易中被广泛使用。国内电子皮带秤制造与生产从20世纪70年代开始进入了一个新的发展时期，国内许多厂家通过学习与交流，开发了全悬浮式皮带秤，使得皮带秤的计量精度可以达到了0.25％。
[0003] 对于电子皮带秤而言，“皮带效应”带来的计量误差是影响计量精度的一个重要方面。“皮带效应”受输送带的张力、称重托辊的非准直度、输送带刚度、物料的流量、环境温度等因素的影响。目前对于皮带秤主要误差影响因素的分析主要是对单一的影响因素进行分析，缺少综合的进行分析。
[0004] 在实际工作中影响精度的因素是多样的，通过误差分析可以得知，张力与温度对精度的影响较大，流量是引起“皮带效应”的寄生载荷的重要影响因素，研究张力、温度和流量对精度的混合影响，可以更贴近实际的了解精度变化规律。温度、张力、等效流量等参数均没有严密的推导公式，但是均对精度有一定的影响。
发明内容
[0005] 本发明的目的是提供一种基于多元线性回归模型的皮带秤主要误差因素分析的方法，为研究各个主要误差影响因素对皮带秤精度的影响程度，以及为企业进行精度补偿提供了一定的依据。
[0006] 实现本发明目的的技术解决方案为：基于多元线性回归分析的皮带秤主要影响因素关系的方法主要步骤如下：
[0007] 步骤1、搭建皮带秤实验平台，根据实际情况改变主要误差因素张力、温度、等效流量的大小，记录各个测试点的张力、温度、等效流量、标定值及料斗称，求出标定值与料斗称的相对误差，即皮带秤精度值；
[0008] 步骤2、计算皮带秤主要误差因素之间的相关系数，判断彼此之间的相关性；
[0009] 步骤3、当各主要误差因素成弱相关或不相关，则设定多元线性模型并求解；
[0010] 步骤4、以残差平方和、决定系数和MS残差方差作为检验指标，评价所建立的模型的回归效果；
[0011] 步骤5、运用拟合模型预测各个测试点的测试值，并与实际值进行对比计算出预测误差，确定拟合结果的准确性。
[0012] 与现有发明相比，本发明具有如下优势：1)可以将复杂的皮带秤工况数据转化为简单的数学变量关系，再结合计算机技术，使求参解算精确而且迅速。2)并且使用多元线性回归分析，可以分析各个因素对结果影响的权重，针对影响较大的因素进行分析研究，减小误差，补偿精度。
附图说明
[0013] 图1为本发明的一种基于多元线性回归模型的皮带秤主要误差因素分析的方法流程图。
具体实施方式
[0014] 结合附图对本发明的技术路线做出进一步说明，此处所应用的实例适用于电子皮带秤的误差影响因素关系分析。
[0015] 参照图1，本发明的一种基于多元线性回归模型的皮带秤主要误差因素分析的方法包括以下步骤：
[0016] 步骤1、搭建皮带秤实验平台，根据实际情况改变主要误差因素张力、温度、等效流量的大小，记录各个测试点的张力、温度、等效流量、标定值及料斗称，求出标定值与料斗称的相对误差，即皮带秤精度值，实验数据如表1所示：
[0017] 表1实验数据
[0018]
[0019] 步骤2、计算各个皮带秤主要误差因素之间的相关系数，判断彼此之间的相关性能。
[0020] 所述步骤2中各个皮带秤主要误差因素之间的相关系数选用皮尔森相关系数进行统计。相关系数用r表示，其中n为样本量，X、Y和 分别为两个变量的观测值和均值。r用来描述两个变量之间的线性相关程度，r的绝对值越大相关性越强。计算公式如下：
[0021]
[0022] 式中，n为样本个数，X、Y为两变量的观测值， 为两变量的均值，SX、SY为两变量的标准差；通常规定：相关系数在0.8-1.0范围，表示极强相关；在0.6-0.8范围，表示强相关；在0.4-0.6范围，表示中等程度相关；在0.2-0.4范围，表示弱相关；在0-0.2范围，表示极弱相关或无相关，计算出的张力、温度、等效流量的相关性整理如表2所示：
[0023] 表2张力、温度、等效流量的相关性
[0024]
[0025] 步骤3、根据上述相关性的判断确定各因素之间的关系，由表2可知各因素之间呈现弱相关，彼此互的影响较小；设定多元线性模型并求解，并运用MATLAB编程进行拟合计算。
[0026] 所述步骤3中模型的确定以及计算，主要包括以下步骤：
[0027] 步骤3.1，设随机变量y与m个自变量存在线性关系：
[0028]
[0029] 式中，β0,β1,…,βm为回归系数，ξ是随机误差，σ是标准差，于是线性关系的另一种表达方式为：
[0030] y～(β0+β1x1+β2x2+…+βmxm,σ2)
[0031] 设有n组样本观测数据，即：
[0032]
[0033] 则上述关系式表示为：
[0034] y＝Xβ+ξ
[0035] 步骤3.2，设b0,b1,…,bm分别为β0,β1,…,βm的最小二乘估计值，于是y的观测值为：
[0036] yk＝b0+b1xk1+b2xk2+…+bmxkm+ek,k＝1,2,…,n
[0037] yk的估计值为：
[0038]
[0039] 其中，ek为误差ξk的估计值，称为残差。为了使得实际值与估计值的拟合程度最好，应使得下式值达到最小：
[0040]
[0041] 即可以转换为求解下面的方程组：
[0042]
[0043] 化简整理后得到正规方程组如下
[0044]
[0045] 记
[0046]
[0047]
[0048] 则正规方程的矩阵形式为：(X′X)b＝X′y，求解出方程组可得b0,b1,…,bm，从而求得回归模型为：
[0049]
[0050] 通过MATLAB编程，对实验数据进行拟合计算，拟合的结果为：
[0051] δ＝-2.18*10-4x1-4.473*10-2x2-1.472*10-4x3+1.6537
[0052] 其中，x1表示张力值，x2表示温度值，x3表示等效流量值，δ表示皮带秤精度值。
[0053] 步骤4、以残差平方和、决定系数和MS残差方差作为检验指标，评价所建立的模型的回归效果，用Spass数据处理软件处理计算获得各个评价指标的结果。
[0054] 所述步骤4中残差平方和SS残、决定系数R2的计算公式如下：
[0055]
[0056] 式中，Yi表示数据点，表示回归直线上相应的位置值；
[0057]
[0058] 式中，R2的值接近于1，说明回归方程的效果越好，表示样本数据能较好地拟合了选用的线性回归模型。用Spass数据处理软件处理计算获得各个评价指标的结果，结果见表3所示。
[0059] 表3各个评价指标的结果
[0060]
[0061]
[0062] 由计算结果可知：决定系数R2>0.9并且SS残差与MS均很小，表明拟合结果较好。
[0063] 步骤5、运用拟合模型预测各个测试点的测试值，并与实际值进行对比计算出预测误差，确定拟合结果的准确性。运用拟合结果公式预测出各个测试点的测试值，并与实际值进行对比计算出预测误差，整理如表4所示。
[0064] 表4各个测试点的测试值与实际值
[0065]
[0066] 通过以上表述，可以看出本发明有以下优点：
[0067] (1)可以将复杂的皮带秤工况数据转化为简单的数学变量关系，再结合计算机技术，使求参解算精确而且迅速。
[0068] (2)并且使用多元线性回归分析，可以分析各个因素对结果影响的权重，针对影响较大的因素进行分析研究，减小误差，补偿精度。
[0069] 应当理解，虽然本说明书按照实施方式加以描述，但这种叙述方式仅为清楚可见，对于该技术领域的技术人员来说，在不脱离本发明原理的前提下，做出的改进也应视为本发明的保护范围。
