matlab 温度计,Matlab与Excel相结合实现标准铂电阻温度计检定数据处理

1990年国际温标(ITS-90)规定在-1893442~660323温区内的温度值由在一组规定的定义固定点分度的标准铂电阻温度计确定,使用规定定义的参考函数和偏差函数内插计算定义固定点之间的温度值[1]。由于参考函数高度非线性,偏差函数的系数需要根据定义固定点的检定数据反求,且各温区的参考函数和偏差函数有区别,出具分度表时,还需迭代计算,因此,标准铂电阻温度计检定数据处理非常复杂。常规的计量检定数据处理通常采用Excel电子表格,具有操作界面友好、数据输入直观的优点,能够满足大部分计量数据处理的需要[2~5]。但是Excel实现高阶方程计算过程复杂,且不易进行迭代数值计算,更不能直接求解偏差函数系数方程组,因此很难实现标准铂电阻温度计计量检定数据的有效处理。而Matlab在复杂数据处理方面具有明显优势,且Matlab提供了ExcelLink插件,通过该插件,可以在Excel操作环境中调用Matlab的强大计算功能,实现复杂的计量检定数据计算。本文整合Matlab和Excel的优点,以简洁高效地解决标准铂电阻温度计检定数据处理问题。1检定数据处理问题以0~660323温区为例,讨论标准铂电阻温度计检定数据处理问题。0~660323温区的参考函数定义为:Wr(t)=C0+9i=1Cit-481481i(1)其逆函数为9t=D0+i=1DiWr(t)-264164i(2)偏差函数为W(t)=W(t)-Wr(t)(3)式(2)和式(3)中:t为温度值,Wr(t)为参考函数,定义固定点上的Wr(t)函数值和C0、D0、Ci、Di均在ITS-90中给出,i=1,2,…,9。W(t)为被检铂电阻的电阻比函数,是温度点t的电阻值Rt与水三相点的电阻值Rtp之比,即W(t)=RtRtp(4)0~660323温区的偏差函数为:W7(t)=a7[W(t)-1]+b7[W(t)-1]2+c7[W(t)-1]3(5)式中:系数a7、b7、c7由锡、锌、铝凝固点测得的W(t)联立求出。0~419527温区内,偏差函数为:W8(t)=a8[W(t)-1]+b8[W(t)-1]2(6)系数a8、b8由锡、锌凝固点测得的W(t)联立方程求出。0~231928温区内,偏差函数为:W9(t)=a9[W(t)-1]+b9[W(t)-1]2(7)系数a9、b9由温度计分别在锡、铟凝固点测得的INFORMATIONTECHNIQUE&CONTROL信息化技术与控制工业计量2011年第21卷第3期W(t)联立方程求出。0~1565985温区内,偏差函数为:W10(t)=a10[W(t)-1](8)系数a10由温度计在铟凝固点测得的W(t)求出。0~297646温区内,偏差函数为:W11(t)=a11[W(t)-1](9)系数a11由温度计在镓熔点测得的W(t)求出。根据JJG1602007要求,检定证书上需根据不同温区给出相应的检定数据,如Rtp、WAl、WZn、WSn、WIn、WGa、WHg、WAr,在水三相点测得的自热效应,以及偏差函数系数ai、bi、ci。对首次检定的温度计除给出以上数据外,还需给出稳定性试验前后的Wt及Rtp。此外,检定证书要求各分区的定义固定点的电阻比数据位数要求保留小数点后4位以上,系数ai、bi、ci的数据位数要求保留小数点后6位以上。当需要出具分度表时,还需进行迭代数值计算,数据处理工作量将大幅增加。2检定数据处理方法21ExcelLink的使用方法ExcelLink插件提供了11个核心函数,对于计量检定数据处理而言,表1中的三个函数基本能实现所有数据处理功