基于Python的最小二乘法在电子天平中的应用.pdf

175lceApl安民军：基于Python的最小乘法在电子天平中的应用1业内思考基于Python的最小二乘法在电子天平中的应用安民军（甘肃省计量研究院，甘肃兰州7 30 0 50）【摘要】本文通过分析在不同电子天平上应用最小二乘法的算法来修正电子天平的误差，使其能够满足计量检定要求，通过读取不同电子天平的检定数据后在Python上对其进行各种数据分析和计算，利用最小二乘法进行曲线拟合得出相对应的系数，然后代人公式计算出其修正后的误差。通过对不同电子天平的对比来发现问题来解决电子天平应用中的实际问题【关键词】电子天平；最小二乘法;Python【D O 1 编码】1 0.39 6 9/j.issn.1674-4977.2023.06.056olication of Least Squares Method Based on Python in Electronic BalanANMinjun(Gansu Institute of Metrology,Lanzhou 730050,China)Abstract:This article analyzes the application of the least squares method algorithm on different electronic scales to correct the error ofthe electronic balance,so that it can meet the requirements of metrological verification.By reading the calibration data of differentelectronic balances and conducting various data analysis and calculations on Python,the least squares method is used to fit the curve toobtain the corresponding coefficients,and then the corrected error is calculated by incorporating the formula.By comparing differentelectronic scales to identify problems and solve practical problems in the application of electronic scales.Key words:electronic balance;least squares method;python电子天平作为非自动衡器，特别是作为机械天平的替代物。一般采用应变式传感器、电容式传感器、电磁平衡式传感器。应变式传感器，结构简单、造价低，但精度有限；电容式传感器称量速度快，性价比较高，但无法达到很高的精度；采用电磁平衡传感器的电子天平，其特点是称量准确可靠、显示快速清晰并具有自动检测系统、简便的自动校准装置以及超载保护等装置。电子天平被称作电子和机械相结合的产物。随着现代工业的发展，称量系统的应用越来越广泛，特别是电子天平的应用，电子天平作为现代称量系统的核心部件，起到了不可或缺的作用。作为电子和机械相结合的产物，电子天平在使用一段时间后，由于受机械变形或者电子元器件老化等各种原因，导致电子天平各种性能下降，一般通过两点标定来修正其系统误差，但是有些电子天平的传感器其线性比较差，不能很好地解决此类问题，本文通过利用最小二乘法的算法，通过多点采集可以解决传感器线性不太理想的此类问题，分析最小二乘法最适合于在哪种情况下的应用。1最小二乘法最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合，最小二乘法是一种广泛应用的算法，在最优控制等各种控制系统和预测中都有应用。本文通过其可用于曲线拟合来计算修正系数。2具体过程本文以分别用厂家为赛多利斯科学仪器（北京）有限公【基金项目】甘肃省市场监督管理局科技计划资助项目：电子天平检定管理系统研究与开发(SSCJG-JL-B202202。176deferror(p,X,y):99.9995,199.9995,209.99931)Yi-np.array0，0.1 0 0 0，1 0.0 0 0 1，2 0.0 0 0 2,49.9 9 9 9Xi00,200,220Enp.arra#采集点BRAND&STANDARDIZATION品牌与标准化2023年第6 期司、梅特勒-托利多仪器（上海）有限公司、RADWAG,型号为BSA224S、M S2 0 4T S/0 2、A S2 2 0.R 2 的三台2 2 0 g的万分之一电子天平为例，分别采集0 g、0.1 g、1 0 g、2 0 g、50 g、1 0 0 g、2 0 0 g220g八个点的检定数据，数据如表1。表1三台电子天平的采集数据检定点/g天平一/g天平二/g天平三/g00.00000.00000.00000.10.10000.10000.10001010.000110.000210.00032020.000220.000520.00055049.999950.000950.001010099.9995100.0021100.0040200199.9995200.0044200.0050220219.9993220.0045220.0055不合格不合格备注合格（线性可以）（线性不行）由于python语言具有开放性，可以灵活应用，很多源代码都是公开的，适合于各个学习阶层，利用网络资源上的开源代码就可以解决。利用其开源代码，调用numpy库，就可以利用python对其进行数据分析及计算 2。其中的matplotlib库可以将数据可视化，便于观察发现问题。scipy是一个高级的科学计算库，可以很方便对数据进行送代运算。通过调用这几个库，特别是numpy和matplotlib库，实现了数据分析及其可视化，便于发现问题，使其更加直观。其程序如下。2.1程序设计import numpy as npimport matplotlib.pyplot as pltimport scipy as spfrom scipy.optimize import leastsq#天平数据#需要拟合的函数func(）指定函数的形状deffunc(p,x):k,b=preturn k*x+b#定义偏差函数，x，y 为数组中对应Xi,Yi的值return func(p,x)-y#设置k,b的初始值，可以任意设定，经过实验，发现po的值会影响cost的值：Para1p0-1,20#把error函数中除了po以外的参数打包到args中，leastsq（）为最小二乘法函数Para=leastsq(error,p0,args=(Xi,Yi)#读取结果k,b=Para0print(k=,k,b=,b)#画样本点plt.figure(figsize=(12,12)plt.scatter(Xi,Yi,color=red,label=Sample datal,linewidth=2)#画拟合直线x=np.linspace(0,220,220)y=k*x+b#绘制拟合曲线plt.plot(x,y,类 color-blue,label-Fitting Curve,linewidth=2)plt.legend()#绘制图例plt.xlabel(Height:g,fontproperties=simHei,fontsize=12)plt.ylabel(Weight:g,fontproperties=simHei,fontsize=12)plt.show()我们通过调用scipy库中的最小二乘法算法来分别得到k和b，k 和b的值分别在表2 中，如下表。表2 k和b值电子天平kb电子天平一0.99999660403972256.723942665942673e5电子天平二1.0000211375776342-1.0582585648677423es电子天平三1.0000257213092980.000108080242970487462.2数据分析对于电子天平来说，作为一个称重系统，出厂最理想的建模输出函数应该为y=x，但是由于在使用过程中由于其传感器的变形、电子元器件的老化以及使用中电压和环境等各方面因素的影响,其输出函数已经不是理想状态下的y=x,本文我们通过一阶函数y=kx+b对其进行拟合，分别对3个天平的输人和输出进行一阶函数y=kx+b的曲线拟合，如图1 所示。177（编辑：侯睿琪）安民军：基于Python的最小乘法在电子天平中的应用月/业内思考FengCrmSample detalFitngCurveSaneledeta2Sonipledotas2020015020131001001005000150200150503003502图1三台电子天平的拟合曲线图然后我们通过对已采集到的八个点输人到一阶函数y=kx+b进行验证，并且同时在将修正值加上后得到其修正后的误差，其修正值为测量值减去标准值。如表3。表3电子天平的拟合值和修正后的误差电子天平一电子天平二电子天平三采集点/g电子天平一电子天平二电子天平三（修正后误差）（修正后误差）（修正后误差）00.00010.00030.00010.00010.00030.00010.10.10010.10000.10010.00010.00000.00011010.000010.000210.0004-0.00010.00000.00012020.000020.000420.0006-0.0002-0.00010.00015049.999950.001050.00140.00000.00010.000410099.9997100.0021100.00270.00020.0000-0.0013200199.9994200.0042200.0053-0.0001-0.00020.0003220219.9993220.0046220.00580.00000.00010.00033结束语通过对表3的电子天平修正后误差可以看出来，电子天平修正后电子天平三在1 0 0 g这个点不合格，但是对于线性可以的电子天平具有较好的补偿作用，所以最小二乘法对于线性好的电子天平具有较好的补偿修正作用。电子天平一的检定点是检定合格的，其修正以后误差非常小；电子天平二的检定点是不合格，但是其传感器线性可以，通过修正以后误差也是非常小；电子天平三的检定点也不合格，其传感器线性在小量程不好，经过修正后其1 0 0 g这个点不合格，通过对比分析，发现传感器在局部非线性时，不能很好地解决此类问题，这就需要控制中的鲁棒控制、神经网络、专家系统等各种先进控制算法。品【参考文献】1曹连江.电子信息测量及其误差分析校正的研究 M.长春：东北师范大学出版社，2 0 1 7-0 9.2 王立友.基于Python的线性回归最小二乘算法的设计与实现.电脑知识与技术,2 0 2 1,1 7(35)：1 2 3-1 2 4,1 33.3 丁克良，欧吉坤，陈义.整体最小二乘法及其在测量数据处理中的应用 C/中国测绘学会，中国测绘学会第九次全国会员代表大会暨学会成立50 周年纪念大会论文集.2 0 0 9:40 6-41 2.4 姜伟，付佳媛.最小二乘法及其应用.中国传媒大学学报（自然科学版），2 0 2 0,2 7（5）：7 2-7 8.5 李蓉，张一帆，杨薇薇.基于最小二乘法和卡尔曼滤波算法的平面电子罗盘误差补偿方法研究.宿州学院学报,2 0 1 9,34(1 1):7 7-8 1.【作者简介】安民军，男，1 9 8 1 年出生，高级工程师，硕士，研究方向为计量检定。